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
|
/* Dfinition d'une structure metrique pour sauver les coefs de la metrique */
struct point
{
double rho;
double gamma;
double alpha;
double omega;
double drhodr;
double dgammadr;
double dalphadr;
double domegadr;
double drhodtheta;
double dgammadtheta;
double dalphadtheta;
double domegadtheta;
};
typedef struct point point;
int max( int i, int j)
{
if ( i<=j ) return j;
else return i;
}
int min( int i, int j)
{
if ( i>=j ) return j;
else return i;
}
double sign( double x )
{
if ( x >= 0 ) return 1. ;
else return -1.;
}
|
