1

サイズの異なる行列の固有値と固有ベクトルを計算しようとしています。私は非常に単純な Fortran90 コードを使用しており、Ubuntu で動作する私のマシンで利用可能な Intel MKL パッケージに含まれる適切な Lapack ライブラリにリンクしてコンパイルしています。メッセージの末尾に「matrix_diag_01.f90」というコードが添付されています。「random」モジュールには、Numerical Recipes の「ran」乱数ジェネレーターが含まれています。コードは次を使用して適切にコンパイルされます

ifort -I $(MKLPATH) -o matrix_diag_01 matrix_diag_01.f90 
      random.f90 $(MKLPATH)libmkl_lapack95.a -Wl,--start-group  
      $(MKLPATH)libmkl_intel_lp64.a $(MKLPATH)libmkl_lapack.a 
      $(MKLPATH)libmkl_intel_thread.a $(MKLPATH)libmkl_core.a
       -Wl,--end-group -lguide -lpthread

小さい行列が与えられた場合、実行可能ファイルはうまく機能します。ただし、サイズが 3000x3000 の行列の場合、奇妙な動作が発生します。最初にこのエラーが発生します

MKL ERROR : Parameter 8 was incorrect on entry to SSYEVD

ただし、SSYEVD の呼び出しには 3 つのパラメーターしかありません。次に、固有ベクトルを返しますが、固有値は返しません。より大きなメモリを搭載した別のマシンでコンパイルして確認しましたが、結果は同じでした。

誰でも助けてもらえますか？

ありがとう！

PROGRAM matrix_diag_01

USE random

IMPLICIT NONE

INTERFACE 
   SUBROUTINE diag(mat,n)
      INTEGER n
      REAL,DIMENSION(n,n) :: mat
   END SUBROUTINE
END INTERFACE 

INTEGER n,i,j,iseed
REAL, DIMENSION(:), ALLOCATABLE  :: w
REAL, DIMENSION(:,:), ALLOCATABLE :: mat

write (*,*) ' Please enter size of matrix'
read (*,*) n
write (*,*) ' Please type seed'
read (*,*) iseed

allocate (mat(n,n))

do i = 1,n
   do j = 1,n
      mat(i,j) = ran(iseed)
   end do
end do

call diag(mat,n)

stop
END
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SUBROUTINE diag(mat,n)

  USE mkl95_lapack
  USE mkl95_precision

  IMPLICIT NONE

  CHARACTER(len=1) :: jobz = 'V'

  INTEGER n,i

  REAL,DIMENSION(n,n) :: mat
  REAL,DIMENSION(:,:),ALLOCATABLE :: matt,a
  REAL,DIMENSION(:),ALLOCATABLE :: w

  allocate (matt(n,n),a(n,n),w(n))

  matt = mat*transpose(mat)
  a = sqrt(matt)
  open (unit=7,file="matrix.dat",status="unknown")
  do i = 1,n
     write (7,100) a(i,:)
  end do
  close (unit=7)

  call syevd(a,w,jobz)

  open (unit=8,file="eig_val.dat",status="unknown")
  do i = 1,n
     write (8,100) w(i)
  end do
  close (unit=8)

  open (unit=9,file="eig_vec.dat",status="unknown")
  do i = 1,n
     write (9,100) a(i,:)
  end do
  close (unit=9)

  return
 100 format(5000f16.5)
end
4

1 に答える 1

2

呼び出している関数のソースを見ると、それ自体が次の呼び出しを発行します: http://software.intel.com/sites/products/documentation/hpc/mkl/lapack/mkl_lapack_examples/ssyevd_ex.f.htm

 LWORK = MIN( LWMAX, INT( WORK( 1 ) ) )

 CALL SSYEVD( 'Vectors', 'Upper', N, A, LDA, W, WORK, LWORK,
         $             IWORK, LIWORK, INFO )

おそらく、パラメータ LWORK が間違っています。

于 2011-07-06T07:58:41.293 に答える