diff options
| author | knarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-10-05 11:07:13 +0000 |
|---|---|---|
| committer | knarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-10-05 11:07:13 +0000 |
| commit | 748569b8a25ce2d5d3b774f18df0345a3fe68768 (patch) | |
| tree | d432e9d6a51ae8a8021f85f40ce2ec85053ee764 /src/TwoPunctures.h | |
| parent | 4b903632cbc6158f80bb0b064ff7fcf45539ac39 (diff) | |
- add source terms provided by aliased function
- small fix suggested by Marcus concerning the initialisation of the dvs
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@19 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
Diffstat (limited to 'src/TwoPunctures.h')
| -rw-r--r-- | src/TwoPunctures.h | 18 |
1 files changed, 11 insertions, 7 deletions
diff --git a/src/TwoPunctures.h b/src/TwoPunctures.h index 9915181..19763e7 100644 --- a/src/TwoPunctures.h +++ b/src/TwoPunctures.h @@ -30,7 +30,8 @@ int Index (int ivar, int i, int j, int k, int nvar, int n1, int n2, int n3); void allocate_derivs (derivs * v, int n); void free_derivs (derivs * v, int n); void Derivatives_AB3 (int nvar, int n1, int n2, int n3, derivs v); -void F_of_v (int nvar, int n1, int n2, int n3, derivs v, +void F_of_v (CCTK_POINTER_TO_CONST cctkGH, + int nvar, int n1, int n2, int n3, derivs v, double *F, derivs u); void J_times_dv (int nvar, int n1, int n2, int n3, derivs dv, double *Jdv, derivs u); @@ -61,7 +62,8 @@ void rx3_To_xyz (int nvar, double x, double r, double phi, double *y, // Routines in "Equations.c" double BY_KKofxyz (double x, double y, double z); void BY_Aijofxyz (double x, double y, double z, double Aij[3][3]); -void NonLinEquations (double A, double B, double X, double R, +void NonLinEquations (CCTK_REAL rho_adm, + double A, double B, double X, double R, double x, double r, double phi, double y, double z, derivs U, double *values); void LinEquations (double A, double B, double X, double R, @@ -77,12 +79,14 @@ void LineRelax_be (double *dv, int i, int k, int nvar, int n1, int n2, int n3, double *rhs, int *ncols, int **cols, double **JFD); void relax (double *dv, int nvar, int n1, int n2, int n3, double *rhs, int *ncols, int **cols, double **JFD); -void TestRelax (int nvar, int n1, int n2, int n3, derivs v, - double *dv); -int bicgstab (int nvar, int n1, int n2, int n3, derivs v, +void TestRelax (CCTK_POINTER_TO_CONST cctkGH, + int nvar, int n1, int n2, int n3, derivs v, double *dv); +int bicgstab (CCTK_POINTER_TO_CONST cctkGH, + int nvar, int n1, int n2, int n3, derivs v, derivs dv, int output, int itmax, double tol, double *normres); -void Newton (int nvar, int n1, int n2, int n3, derivs v, - double tol, int itmax); +void Newton (CCTK_POINTER_TO_CONST cctkGH, + int nvar, int n1, int n2, int n3, derivs v, + double tol, int itmax); /* |