diff options
Diffstat (limited to 'src/sphericalharmonic.cc')
| -rw-r--r-- | src/sphericalharmonic.cc | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/sphericalharmonic.cc b/src/sphericalharmonic.cc index 116613c..949d261 100644 --- a/src/sphericalharmonic.cc +++ b/src/sphericalharmonic.cc @@ -88,8 +88,16 @@ extern "C" void Multipole_SetHarmonic(CCTK_ARGUMENTS) CCTK_REAL theta = acos(z[index]/r[index]); CCTK_REAL phi = atan2(y[index],x[index]); + CCTK_REAL re = 0; + CCTK_REAL im = 0; + Multipole_SphericalHarmonic(test_sw,test_l,test_m,theta,phi, - &harmonic_re[index], &harmonic_im[index]); + &re, &im); + + CCTK_REAL fac = test_mode_proportional_to_r ? r[index] : 1.0; + + harmonic_re[index] = re * fac; + harmonic_im[index] = im * fac; } } } |