program test_mkl
   
      
      IMPLICIT NONE
      INTEGER          N
      PARAMETER        ( N = 2 )
      INTEGER          LDA
      PARAMETER        ( LDA = N )
      INTEGER          LWMAX
      PARAMETER        ( LWMAX = 1000 )

      INTEGER          INFO, LWORK, LIWORK

      
      INTEGER          IWORK( LWMAX )
      DOUBLE PRECISION A( LDA, N ), W( N ), WORK( LWMAX )
      

      DATA A/1,3,3,4/
      
      WRITE(*,*)'DSYEVD Example Program Results'
      
      LWORK = -1
      LIWORK = -1
     
      CALL DSYEVD ( 'V', 'U', N, A, LDA, W, WORK, LWORK,IWORK, LIWORK, INFO )
      
         
     ! Print eigenvalues.
      PRINT *, 'Eigenvalues', W

     ! Print eigenvectors.

      PRINT *, 'Eigenvectors (stored columnwise)', A
      
      STOP
   
    end

2楼: Originally posted by wlscj at 2014-11-06 19:17:22
调用了MKL里面的DSYEVD函数?

DSYEVD 是个啥?

3楼: Originally posted by 徐谊 at 2014-11-06 19:23:50
嗯，是的。求本征值的
?syevd
Computes all eigenvalues and (optionally) all
eigenvectors of a real symmetric matrix using divide
and conquer algorithm....