[Commits] [svn:einsteintoolkit] Multipole/trunk/ (Rev. 58)

[ian.hinder at aei.mpg.de]

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
 {