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