User: hinder Date: 2010/06/14 01:12 AM
Modified: /trunk/ param.ccl /trunk/src/ sphericalharmonic.cc
Log: Add testing code for making grid function harmonic coefficient proportional to radius
File Changes:
Directory: /trunk/src/ ======================
File [modified]: sphericalharmonic.cc Delta lines: +9 -1 =================================================================== --- trunk/src/sphericalharmonic.cc 2010-06-13 23:13:33 UTC (rev 57) +++ trunk/src/sphericalharmonic.cc 2010-06-14 06:12:17 UTC (rev 58) @@ -88,8 +88,16 @@ 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; } } }
Directory: /trunk/ ==================
File [modified]: param.ccl Delta lines: +5 -0 =================================================================== --- trunk/param.ccl 2010-06-13 23:13:33 UTC (rev 57) +++ trunk/param.ccl 2010-06-14 06:12:17 UTC (rev 58) @@ -94,6 +94,11 @@ * :: "Any integer" } -2
+CCTK_BOOLEAN test_mode_proportional_to_r "whether the test spherical harmonic coefficient is proportional to the radial coordinate" +{ +} "no" + + # This parameter is deprecated and is no longer used KEYWORD mode_type "Which type of mode extraction do we have" STEERABLE=always {