g***i
发帖数: 90 | 1 multiple right hand sides CG method
算法或程序, thanks | g***i
发帖数: 90 | 2 finally got some algorithm and programmed in fortran
using lapack. but it turns out that the block CG doesn't
work well
here is the program if someone is interested:
subroutine bcg(A,B,X,n,r,eps,kmax)
implicit none
integer n,r,k,kmax,i,j,INFO,LWORK
real(8) A(n,n),B(n,r), X(n,r), R0(n,r), D0(r,r), D1(r,r), &
Al(r,r), Be(r,r), P(n,r), &
Q(n,r), RR(r,r), IPIV(r), NR(n,r), &
eps,errm,err(r),err0(r)
real(8), dimen
【在 g***i 的大作中提到】
: multiple right hand sides CG method
: 算法或程序, thanks
|
|
