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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
|
/*@@
@file Coord.c
@date Indeterminate 11-12th April 1999
@author Gabrielle Allen
@desc
Routines to deal with cooordinates and coordinate registration
@enddesc
@@*/
/*#define DEBUG_COORD*/
#include <stdlib.h>
#include <stdio.h>
#include "cctk.h"
#include "StoreHandledData.h"
#include "WarnLevel.h"
#include "Coord.h"
#include "ErrorCodes.h"
#include "Groups.h"
#include "Misc.h"
#include "FortranString.h"
static cHandledData *coordinates = NULL;
static int num_coords = 0;
/*@@
@routine RegisterCoord_ByIndex
@date 11-12th April 1999
@author Gabrielle Allen
@desc
Register a GF as a coordinate with a name, and index
and a direction
@enddesc
@calls Util_GetHandle, Util_NewHandle, CCTK_Warn
@var name
@vdesc Name coordinate is registered as
@vtype const char *
@vio in
@vcomment The name must not be the same as the variable name
@endvar
@var index
@vdesc The index of the variable providing the coordinate
@vtype int
@vio in
@vcomment
@endvar
@var dir
@vdesc An index giving the coordinate direction
@vtype int
@vio in
@vcomment This should be an optional argument
@endvar
@returntype int
@returndesc
-1 = Duplicate coordinate name, all duplicates
will be ignored
ERROR_MEMORY = Failure with malloc
>= 0 Handle returned by coordinate registration
@endreturndesc
@@*/
int CCTK_RegisterCoordI(const char *name, int index, int dir)
{
int handle;
struct Coordprops *new_coord;
/* Check that the method hasn't already been registered */
handle = Util_GetHandle(coordinates, name, NULL);
if(handle < 0)
{
/* New extension. */
new_coord = (struct Coordprops *)malloc(sizeof(struct Coordprops));
if(new_coord)
{
/* Get a handle for it. */
handle = Util_NewHandle(&coordinates, name, new_coord);
/* Initialise the coordinate properties structure */
new_coord->name = (char *)name;
new_coord->index = index;
new_coord->direction = dir;
new_coord->origin = 0;
/* Remember how many methods there are */
num_coords++;
#ifdef DEBUG_COORD
printf(" In CCTK_RegisterCoordI\n");
printf(" -----------------------------\n");
printf(" handle %d, name %s,\n",handle,name);
printf(" index %d, direction %d\n",index,dir);
#endif
}
else
{
/* Memory failure. */
handle = ERROR_MEMORY;
}
}
else
{
/* Method already exists. */
char *msg;
msg = (char *)malloc(200*sizeof(char)+sizeof(name));
sprintf(msg,"Coordinate with name -%s- already registered",name);
CCTK_WARN(1,msg);
if (msg) free(msg);
handle = -1;
}
return handle;
}
/*@@
@routine RegisterCoord
@date 11-12th April 1999
@author Gabrielle Allen
@desc
Register a GF as a coordinate with a name and
a direction.
@enddesc
@calls CCTK_RegisterCoordI, CCTK_VarIndex
@var name
@vdesc Name coordinate is registered as
@vtype const char *
@vio in
@vcomment The name must not be the same as the variable name
@endvar
@var gfname
@vdesc The name of the variable providing the coordinate
@vtype const char *
@vio in
@vcomment
@endvar
@var dir
@vdesc An index giving the coordinate direction
@vtype int
@vio in
@vcomment This should be an optional argument
@endvar
@returntype int
@returndesc
-1 = Duplicate coordinate name, all duplicates
will be ignored
ERROR_MEMORY = Failure with malloc
>= 0 Handle returned by coordinate registration
@endreturndesc
@@*/
int CCTK_RegisterCoord(const char *coordname,
const char *gfname,
int dir)
{
int retval;
int index;
index = CCTK_VarIndex(gfname);
if (index >= 0)
{
retval = CCTK_RegisterCoordI(coordname,index,dir);
}
else
{
/* FIXME Should this register this error or a failed to register coord error */
retval = ERROR_INDEXNOTINRANGE;
}
return retval;
}
int CCTK_RegisterCoordRange(cGH *GH, CCTK_REAL min, CCTK_REAL max,
const char *coordname)
{
int retval;
return retval;
}
int CCTK_CoordIndex(const char *name)
{
int handle;
struct Coordprops *coord;
for (handle = 0;;handle++)
{
coord = (struct Coordprops *)Util_GetHandledData(coordinates, handle);
if (coord)
{
if (CCTK_Equals(name,(const char *)coord->name))
return coord->index;
}
else
{
char *msg;
msg = (char *)malloc( 100*sizeof(char)+sizeof(name) );
sprintf(msg,"Could not find registered coordinate %s",name);
CCTK_WARN(2,msg);
if (msg) free(msg);
return ERROR_COORDNOTFOUND;
}
}
}
void FMODIFIER FORTRAN_NAME(CCTK_CoordIndex)(int *handle, ONE_FORTSTRING_ARG)
{
ONE_FORTSTRING_CREATE(name)
*handle = CCTK_CoordIndex (name);
free(name);
}
CCTK_REAL CCTK_CoordOrigin(const char *name)
{
int handle;
struct Coordprops *coord;
for (handle = 0;;handle++)
{
coord = (struct Coordprops *)Util_GetHandledData(coordinates, handle);
if (coord)
{
if (CCTK_Equals(name,(const char *)coord->name))
return coord->origin;
}
else
{
char *msg;
msg = (char *)malloc( 100*sizeof(char)+sizeof(name) );
sprintf(msg,"Could not find registered coordinate %s",name);
CCTK_WARN(2,msg);
if (msg) free(msg);
return ERROR_COORDNOTFOUND;
}
}
}
|