我对C编程很陌生,并为Jacobi符号的计算编写了以下算法。虽然有一些更快的版本,这个算法可用,我只是寻找一些反馈,我的编码风格。+*****************************************************************************
* Notes: This function evaluates the Jacobi
我正在用Jacobi方法求解泊松方程(在2d轴对称圆柱坐标系中)。L2范数从第一次迭代时的~1E3 (我有一个非常糟糕的猜测)下降到~0.2,非常慢。然后,L2范数开始在多次迭代中增加。这是我的Jacobi方法算法。然后把它封装在一个while循环中。subroutine Jacobi(PoissonRHS, V, resid)real, dimension(0:,0:) :: PoissonRHS, V
REAL residcompare to old V values to check for conver