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