1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
|
/* TwoPunctures: File "TwoPunctures.h"*/
#define StencilSize 19
#define N_PlaneRelax 1
#define NRELAX 200
#define Step_Relax 1
typedef struct DERIVS
{
CCTK_REAL *d0, *d1, *d2, *d3, *d11, *d12, *d13, *d22, *d23, *d33;
} derivs;
/*
Files of "TwoPunctures":
TwoPunctures.c
FuncAndJacobian.c
CoordTransf.c
Equations.c
Newton.c
utilities.c (see utilities.h)
**************************
*/
/* Routines in "TwoPunctures.c"*/
CCTK_REAL TestSolution (CCTK_REAL A, CCTK_REAL B, CCTK_REAL X, CCTK_REAL R, CCTK_REAL phi);
void TestVector_w (CCTK_REAL *par, int nvar, int n1, int n2, int n3, CCTK_REAL *w);
/* Routines in "FuncAndJacobian.c"*/
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 (CCTK_POINTER_TO_CONST cctkGH,
int nvar, int n1, int n2, int n3, derivs v,
CCTK_REAL *F, derivs u);
void J_times_dv (int nvar, int n1, int n2, int n3, derivs dv,
CCTK_REAL *Jdv, derivs u);
void JFD_times_dv (int i, int j, int k, int nvar, int n1,
int n2, int n3, derivs dv, derivs u, CCTK_REAL *values);
void SetMatrix_JFD (int nvar, int n1, int n2, int n3,
derivs u, int *ncols, int **cols, CCTK_REAL **Matrix);
CCTK_REAL PunctEvalAtArbitPosition (CCTK_REAL *v, int ivar, CCTK_REAL A, CCTK_REAL B, CCTK_REAL phi,
int nvar, int n1, int n2, int n3);
void calculate_derivs (int i, int j, int k, int ivar, int nvar, int n1,
int n2, int n3, derivs v, derivs vv);
CCTK_REAL interpol (CCTK_REAL a, CCTK_REAL b, CCTK_REAL c, derivs v);
CCTK_REAL PunctTaylorExpandAtArbitPosition (int ivar, int nvar, int n1,
int n2, int n3, derivs v, CCTK_REAL x,
CCTK_REAL y, CCTK_REAL z);
CCTK_REAL PunctIntPolAtArbitPosition (int ivar, int nvar, int n1,
int n2, int n3, derivs v, CCTK_REAL x,
CCTK_REAL y, CCTK_REAL z);
void SpecCoef(int n1, int n2, int n3, int ivar, CCTK_REAL *v, CCTK_REAL *cf);
CCTK_REAL PunctEvalAtArbitPositionFast (CCTK_REAL *v, int ivar, CCTK_REAL A, CCTK_REAL B, CCTK_REAL phi,
int nvar, int n1, int n2, int n3);
CCTK_REAL PunctIntPolAtArbitPositionFast (int ivar, int nvar, int n1,
int n2, int n3, derivs v, CCTK_REAL x,
CCTK_REAL y, CCTK_REAL z);
/* Routines in "CoordTransf.c"*/
void AB_To_XR (int nvar, CCTK_REAL A, CCTK_REAL B, CCTK_REAL *X,
CCTK_REAL *R, derivs U);
void C_To_c (int nvar, CCTK_REAL X, CCTK_REAL R, CCTK_REAL *x,
CCTK_REAL *r, derivs U);
void rx3_To_xyz (int nvar, CCTK_REAL x, CCTK_REAL r, CCTK_REAL phi, CCTK_REAL *y,
CCTK_REAL *z, derivs U);
/* Routines in "Equations.c"*/
CCTK_REAL BY_KKofxyz (CCTK_REAL x, CCTK_REAL y, CCTK_REAL z);
void BY_Aijofxyz (CCTK_REAL x, CCTK_REAL y, CCTK_REAL z, CCTK_REAL Aij[3][3]);
void NonLinEquations (CCTK_REAL rho_adm,
CCTK_REAL A, CCTK_REAL B, CCTK_REAL X, CCTK_REAL R,
CCTK_REAL x, CCTK_REAL r, CCTK_REAL phi,
CCTK_REAL y, CCTK_REAL z, derivs U, CCTK_REAL *values);
void LinEquations (CCTK_REAL A, CCTK_REAL B, CCTK_REAL X, CCTK_REAL R,
CCTK_REAL x, CCTK_REAL r, CCTK_REAL phi,
CCTK_REAL y, CCTK_REAL z, derivs dU, derivs U, CCTK_REAL *values);
/* Routines in "Newton.c"*/
void TestRelax (CCTK_POINTER_TO_CONST cctkGH,
int nvar, int n1, int n2, int n3, derivs v, CCTK_REAL *dv);
void Newton (CCTK_POINTER_TO_CONST cctkGH,
int nvar, int n1, int n2, int n3, derivs v,
CCTK_REAL tol, int itmax);
/*
27: -1.325691774825335e-03
37: -1.325691778944117e-03
47: -1.325691778942711e-03
17: -1.510625972641537e-03
21: -1.511443006977708e-03
27: -1.511440785153687e-03
37: -1.511440809549005e-03
39: -1.511440809597588e-03
*/
|
