aboutsummaryrefslogtreecommitdiff
path: root/src/Equations.c
diff options
context:
space:
mode:
authorknarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf>2005-02-26 14:05:39 +0000
committerknarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf>2005-02-26 14:05:39 +0000
commit0b28e2b6563dcc87e3e8ae0dd325f7b9dce76493 (patch)
tree6b6a655aaa9d8170b62f0f73ed0c130bfb1c0f01 /src/Equations.c
parent4e6de8925c991fa2ba8cdb41cb4c3ed2044330e0 (diff)
transform // style comments into /* */ style comments because some compilers
do not understand C99 yet, unfortunately :( no other change done git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@40 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
Diffstat (limited to 'src/Equations.c')
-rw-r--r--src/Equations.c48
1 files changed, 24 insertions, 24 deletions
diff --git a/src/Equations.c b/src/Equations.c
index ed228e7..1c9c5f4 100644
--- a/src/Equations.c
+++ b/src/Equations.c
@@ -1,4 +1,4 @@
-// TwoPunctures: File "Equations.c"
+/* TwoPunctures: File "Equations.c"*/
#include <stdio.h>
#include <stdlib.h>
@@ -9,16 +9,16 @@
#include "TP_utilities.h"
#include "TwoPunctures.h"
-// U.d0[ivar] = U[ivar]; (ivar = 0..nvar-1)
-// U.d1[ivar] = U[ivar]_x;
-// U.d2[ivar] = U[ivar]_y;
-// U.d3[ivar] = U[ivar]_z;
-// U.d11[ivar] = U[ivar]_xx;
-// U.d12[ivar] = U[ivar]_xy;
-// U.d13[ivar] = U[ivar]_xz;
-// U.d22[ivar] = U[ivar]_yy;
-// U.d23[ivar] = U[ivar]_yz;
-// U.d33[ivar] = U[ivar]_zz;
+/* U.d0[ivar] = U[ivar]; (ivar = 0..nvar-1) */
+/* U.d1[ivar] = U[ivar]_x; */
+/* U.d2[ivar] = U[ivar]_y; */
+/* U.d3[ivar] = U[ivar]_z; */
+/* U.d11[ivar] = U[ivar]_xx; */
+/* U.d12[ivar] = U[ivar]_xy; */
+/* U.d13[ivar] = U[ivar]_xz;*/
+/* U.d22[ivar] = U[ivar]_yy;*/
+/* U.d23[ivar] = U[ivar]_yz;*/
+/* U.d33[ivar] = U[ivar]_zz;*/
double
BY_KKofxyz (double x, double y, double z)
@@ -42,7 +42,7 @@ BY_KKofxyz (double x, double y, double z)
n_plus[2] = z / r_plus;
n_minus[2] = z / r_minus;
- // dot product: np_Pp = (n_+).(P_+); nm_Pm = (n_-).(P_-)
+ /* dot product: np_Pp = (n_+).(P_+); nm_Pm = (n_-).(P_-) */
np_Pp = 0;
nm_Pm = 0;
for (i = 0; i < 3; i++)
@@ -50,7 +50,7 @@ BY_KKofxyz (double x, double y, double z)
np_Pp += n_plus[i] * par_P_plus[i];
nm_Pm += n_minus[i] * par_P_minus[i];
}
- // cross product: np_Sp[i] = [(n_+) x (S_+)]_i; nm_Sm[i] = [(n_-) x (S_-)]_i
+ /* cross product: np_Sp[i] = [(n_+) x (S_+)]_i; nm_Sm[i] = [(n_-) x (S_-)]_i*/
np_Sp[0] = n_plus[1] * par_S_plus[2] - n_plus[2] * par_S_plus[1];
np_Sp[1] = n_plus[2] * par_S_plus[0] - n_plus[0] * par_S_plus[2];
np_Sp[2] = n_plus[0] * par_S_plus[1] - n_plus[1] * par_S_plus[0];
@@ -61,7 +61,7 @@ BY_KKofxyz (double x, double y, double z)
for (i = 0; i < 3; i++)
{
for (j = 0; j < 3; j++)
- { // Bowen-York-Curvature :
+ { /* Bowen-York-Curvature :*/
Aij =
+ 1.5 * (par_P_plus[i] * n_plus[j] + par_P_plus[j] * n_plus[i]
+ np_Pp * n_plus[i] * n_plus[j]) / r2_plus
@@ -100,7 +100,7 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3])
n_plus[2] = z / r_plus;
n_minus[2] = z / r_minus;
- // dot product: np_Pp = (n_+).(P_+); nm_Pm = (n_-).(P_-)
+ /* dot product: np_Pp = (n_+).(P_+); nm_Pm = (n_-).(P_-) */
np_Pp = 0;
nm_Pm = 0;
for (i = 0; i < 3; i++)
@@ -108,7 +108,7 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3])
np_Pp += n_plus[i] * par_P_plus[i];
nm_Pm += n_minus[i] * par_P_minus[i];
}
- // cross product: np_Sp[i] = [(n_+) x (S_+)]_i; nm_Sm[i] = [(n_-) x (S_-)]_i
+ /* cross product: np_Sp[i] = [(n_+) x (S_+)]_i; nm_Sm[i] = [(n_-) x (S_-)]_i*/
np_Sp[0] = n_plus[1] * par_S_plus[2] - n_plus[2] * par_S_plus[1];
np_Sp[1] = n_plus[2] * par_S_plus[0] - n_plus[0] * par_S_plus[2];
np_Sp[2] = n_plus[0] * par_S_plus[1] - n_plus[1] * par_S_plus[0];
@@ -118,7 +118,7 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3])
for (i = 0; i < 3; i++)
{
for (j = 0; j < 3; j++)
- { // Bowen-York-Curvature :
+ { /* Bowen-York-Curvature :*/
Aij[i][j] =
+ 1.5 * (par_P_plus[i] * n_plus[j] + par_P_plus[j] * n_plus[i]
+ np_Pp * n_plus[i] * n_plus[j]) / r2_plus
@@ -132,9 +132,9 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3])
}
}
-//---------------------------------------------------------------
-//******* Nonlinear Equations **********
-//---------------------------------------------------------------
+/*-----------------------------------------------------------*/
+/******** Nonlinear Equations ***********/
+/*-----------------------------------------------------------*/
void
NonLinEquations (CCTK_REAL rho_adm,
double A, double B, double X, double R,
@@ -160,9 +160,9 @@ NonLinEquations (CCTK_REAL rho_adm,
}
-//---------------------------------------------------------------
-//******* Linear Equations **********
-//---------------------------------------------------------------
+/*-----------------------------------------------------------*/
+/******** Linear Equations ***********/
+/*-----------------------------------------------------------*/
void
LinEquations (double A, double B, double X, double R,
double x, double r, double phi,
@@ -184,4 +184,4 @@ LinEquations (double A, double B, double X, double R,
- 0.875 * BY_KKofxyz (x, y, z) / psi8 * dU.d0[0];
}
-//---------------------------------------------------------------
+/*-----------------------------------------------------------*/