为什么CVF6.6的IMSL不提供MPI
不能用,比如
USE linear_operators
USE mpi_setup_int
   integer,parameter :: n=3, k=16
   real, dimension(n,n,k) :: A,U,V,R,S(n,k), units_of_error(k)
   MP_NPROCS=MP_SETUP()       ! Set up MPI.
   A=rand(A); S=SVD(A, U=U, V=V)
   R = U .x. diag(S) .xt. V; 
   units_of_error = norm(A-R)/S(1,1:k)/epsilon(A)
   MP_NPROCS=MP_SETUP('Final') ! Shut down MPI.
   print*,units_of_error
   end
这样的代码无法运行。
去掉  MP_NPROCS=MP_SETUP()       ! Set up MPI.
    MP_NPROCS=MP_SETUP('Final') ! Shut down MPI.
就能运行
请高手指点!