[Wien] mmatrix

[Yongbin Lee]

Dear WIEN-user

 I must be comfused somthing.
There are 4 inline functions in mmatrix.f for SRC_optic.

 F0(L)=(2*L+1)*(2*L+3)
 F1(L,M)= -SQRT((L+M+1)*(L+M+2)/F0(L))
 F3(L,M)=  SQRT((L-M+1)*(L-M+2)/F0(L))
 F5(L,M)=  SQRT((L-M+1)*(L+M+1)/F0(L))

Instead of these functions, I got 

 F0(L)= 2*(2*L+1)*(2*L+3)
 F1(L,M)= ((-1)**(L+M))*SQRT((L+M+1)*(L+M+2)/F0(L))
 F3(L,M)= ((-1)**(L+M))*SQRT((L-M+1)*(L-M+2)/F0(L))
 F5(L,M)= ((-1)**(L+M))*SQRT(2*(L-M+1)*(L+M+1)/F0(L))

from  EDMONDS's gradient formula.
But I could not see the way to go the fuctions in the mmatrix.f.
Any suggestion ?

 Thank you

  Yongbin