blob: f4715e26e17e426a8bda53d6c2b2d2beec7fea7c (
plain)
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
|
/*@@
@file util.c
@date 23 Oct 2001
@author Jonathan Thornburg <jthorn@aei.mpg.de>
@desc
Utility routines for generalized interpolation.
This file contains various utility routines for the interpolator.
@enddesc
@version $Id$
@@*/
#include <math.h>
#include <limits.h>
#include <stdlib.h>
#include <stdio.h>
#ifndef LOCALINTERP_STANDALONE_BUILD
#include "cctk.h"
#endif
#include "InterpLocalUniform.h"
/* the rcs ID and its dummy function to use it */
static const char *rcsid = "$Header$";
#ifndef LOCALINTERP_STANDALONE_BUILD
CCTK_FILEVERSION(CactusBase_LocalInterp_src_GeneralizedPolynomialUniform_util_c)
#endif
/******************************************************************************/
/*@@
@routine LocalInterp_zero_int_array
@date 23 Oct 2001
@author Jonathan Thornburg <jthorn@aei.mpg.de>
@desc This function sets a CCTK_INT array to all zeros.
@enddesc
@var N
@vdesc The number of elements in the array.
@vtype CCTK_INT N
@endvar
@var array
@vdesc (A pointer to) the array
@vtype CCTK_INT array[]
@vio out
@endvar
@@*/
void LocalInterp_zero_int_array(int N, CCTK_INT array[])
{
int i;
for (i = 0 ; i < N ; ++i)
{
array[i] = 0;
}
}
/******************************************************************************/
/*@@
@routine LocalInterp_decode_N_parts
@date 22 Jan 2002
@author Jonathan Thornburg <jthorn@aei.mpg.de>
@desc This function decodes a CCTK_VARIABLE_* variable type code
(cctk_Constants.h) to determine whether the type is real
or complex.
@enddesc
@var type_code
@vdesc The type code to be decoded
@vtype int
@endvar
@returntype int
@returndesc This function returns
1 if the data type represents a real number of some sort
(includes integers)
2 if the data type represents a complex number
0 if the data type doesn't represent a number,
eg strings and pointers
-1 if the data type is invalid
@endreturndesc
@@*/
int LocalInterp_decode_N_parts(int type_code)
{
switch (type_code)
{
case CCTK_VARIABLE_VOID: return 0;
case CCTK_VARIABLE_BYTE: return 1;
case CCTK_VARIABLE_INT: return 1;
case CCTK_VARIABLE_INT2: return 1;
case CCTK_VARIABLE_INT4: return 1;
case CCTK_VARIABLE_INT8: return 1;
case CCTK_VARIABLE_REAL: return 1;
case CCTK_VARIABLE_REAL4: return 1;
case CCTK_VARIABLE_REAL8: return 1;
case CCTK_VARIABLE_REAL16: return 1;
case CCTK_VARIABLE_COMPLEX: return 2;
case CCTK_VARIABLE_COMPLEX8: return 2;
case CCTK_VARIABLE_COMPLEX16: return 2;
case CCTK_VARIABLE_COMPLEX32: return 2;
case CCTK_VARIABLE_STRING: return 0;
case CCTK_VARIABLE_POINTER: return 0;
case CCTK_VARIABLE_FPOINTER: return 0;
default: return -1;
}
}
|
