|
地板

楼主 |
发表于 2018-1-15 10:18:07
|
只看该作者
错误行应该是第20行(x=work(a,c,d,y,nmax))
[Fortran] 纯文本查看 复制代码 01 | subroutine absorption ( p 0 , p , m , lamda , kmax , mmax , nmax ) |
05 | integer :: kmax , mmax , nmax |
06 | real , dimension ( kmax , mmax , nmax ) :: p , p 0 |
07 | real , dimension ( nmax ) :: x , y |
08 | real , dimension ( nmax ) :: a , c , d |
23 | end subroutine absorption |
另外work函数的代码是:
[Fortran] 纯文本查看 复制代码 01 | function work ( a , c , d , b , ranka ) |
05 | real , dimension ( ranka ) :: a , c , d , work , y , b , r , delta , beta |
15 | r ( i ) = d ( i ) - beta ( i ) * delta ( i -1 ) |
17 | y ( i ) = ( b ( i ) - beta ( i ) * y ( i -1 ) ) / r ( i ) |
23 | work ( i ) = y ( i ) - delta ( i ) * work ( i +1 ) |
|
|