Changeset 1621 in Sophya
- Timestamp:
- Aug 6, 2001, 7:32:33 PM (24 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/Samba/alm.cc
r865 r1621 64 64 int_4 nlmax=Lmax(); 65 65 66 TVector<T> pow(nlmax+1); 67 pow.SetTemp(true); 66 TVector<T> powsp(nlmax+1); 68 67 69 68 for (int l=0; l<=nlmax;l++) 70 69 { 71 pow (l)=( (*this)(l,0) ).real()*( (*this)(l,0) ).real();70 powsp(l)=( (*this)(l,0) ).real()*( (*this)(l,0) ).real(); 72 71 for (int m=1; m<=l; m++) 73 72 { 74 pow (l)+=2.*( (*this)(l,m).real()*(*this)(l,m).real()+73 powsp(l)+=2.*( (*this)(l,m).real()*(*this)(l,m).real()+ 75 74 (*this)(l,m).imag()*(*this)(l,m).imag() ); 76 75 } 77 pow (l)/=(2.*l+1.);76 powsp(l)/=(2.*l+1.); 78 77 } 79 return pow ;78 return powsp; 80 79 } 81 80
Note:
See TracChangeset
for help on using the changeset viewer.