Quantcast
Channel: Intel® Software - Intel® Visual Fortran Compiler for Windows*
Viewing all articles
Browse latest Browse all 5691

Problem with MKL dpotrf function

$
0
0

Can anyone see why I am getting all zeroes back in [A] when calling the MKL dpotrf function below? Info comes back as 0 which indicates success, but A is just full of zeroes. What am I missing?

I have set it up in Visual Studio 2010 as a 32-bit application as follows:
Fortran => Libraries => Use Intel Math Kernel Library = Sequential
Linker => General => Additional Library Directories = C:\Program Files (x86)\Intel\Composer XE 2013 SP1.139\mkl\lib\ia32
Linker => Input => Additional Dependencies = mkl_blas95.lib mkl_lapack95.lib mkl_intel_c.lib mkl_core.lib mkl_sequential.lib

      SUBROUTINE MatrixTest

      IMPLICIT NONE

      INCLUDE 'mkl.fi'

      INTEGER    n
      PARAMETER (n=4)

      REAL*8     A(n,n)
      INTEGER    info,i,j

c     Fill the upper triangle
      A(1,1)=11.0d0
      A(1,2)=-5.0d0
      A(1,3)=3.0d0
      A(1,4)=2.0d0
      A(2,2)=13.0d0
      A(2,3)=4.0d0
      A(2,4)=-6.0d0
      A(3,3)=16.0d0
      A(3,4)=7.0d0
      A(4,4)=19.0d0

c     Make it symmetrical
      do i=2,n
         do j=1,i-1
            A(i,j)=A(j,i)
         enddo
      enddo

      call dpotrf ('L',n,A,n,info)

      end


Viewing all articles
Browse latest Browse all 5691

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>