diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/AHFinder_fun.F | 19 |
1 files changed, 12 insertions, 7 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F index 0af446c..2790e56 100644 --- a/src/AHFinder_fun.F +++ b/src/AHFinder_fun.F @@ -287,8 +287,8 @@ ! the recursion relations requires 2 starting values. do m=1,lmax-2 - - legen_base = -legen_base*aux2*lm_factor3(m) ! L=M +! L=M + legen_base = -legen_base*aux2*lm_factor3(m) aux1 = aux1 + legen_base*cc_cosa(m,m) @@ -297,7 +297,8 @@ end if local_legen_old_old = legen_base - local_legen_old = cost*lm_factor3(m+1) ! L=M+1 +! L=M+1 + local_legen_old = cost*lm_factor3(m+1) . *legen_base aux1 = aux1 + local_legen_old*cc_cosa(m+1,m) @@ -306,7 +307,8 @@ aux1 = aux1 + local_legen_old*cs_sina(m+1,m) end if - do l=m+2,lmax ! L=M+2,lmax +! L=M+2,lmax + do l=m+2,lmax local_legen = cost*lm_factor1(l,m)*local_legen_old . - lm_factor2(l,m)*local_legen_old_old @@ -326,7 +328,8 @@ if (lmax.gt.1) then - m=lmax-1 ! L=lmax-1 +! L=lmax-1 + m=lmax-1 legen_base = -legen_base*aux2*lm_factor3(m) aux1 = aux1 + legen_base*cc_cosa(m,m) @@ -334,7 +337,8 @@ aux1 = aux1 + legen_base*cs_sina(m,m) end if - local_legen = cost*lm_factor3(lmax)* ! L=lmax +! L=lmax + local_legen = cost*lm_factor3(lmax)* . legen_base aux1 = aux1 + local_legen*cc_cosa(lmax,m) @@ -346,7 +350,8 @@ ! And finally do M=lmax. - m=lmax ! L=lmax +! L=lmax + m=lmax legen_base = -legen_base*aux2*lm_factor3(m) aux1 = aux1 + legen_base*cc_cosa(m,m) |