主题:为什么CVF6.6的IMSL不提供MPI
为什么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.
就能运行
请高手指点!
不能用,比如
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.
就能运行
请高手指点!