PLUTO
four_vel.c
Go to the documentation of this file.
1 #include "pluto.h"
2 
3 /* ********************************************************************* */
4 void ConvertTo4vel (double **v, int beg, int end)
5 /*
6  *
7  *
8  *********************************************************************** */
9 {
10  int i;
11  double lor;
12 
13  for (i = beg; i <= end; i++){
14  lor = EXPAND( v[i][VX1]*v[i][VX1],
15  + v[i][VX2]*v[i][VX2],
16  + v[i][VX3]*v[i][VX3]);
17  lor = 1.0/sqrt(1.0 - lor);
18 
19  EXPAND(v[i][VX1] *= lor; ,
20  v[i][VX2] *= lor; ,
21  v[i][VX3] *= lor;)
22  }
23 }
24 /* ********************************************************************* */
25 void ConvertTo3vel (double **v, int beg, int end)
26 /*
27  *
28  *
29  *********************************************************************** */
30 {
31  int i;
32  double lor;
33 
34  for (i = beg; i <= end; i++){
35  lor = EXPAND( v[i][VX1]*v[i][VX1],
36  + v[i][VX2]*v[i][VX2],
37  + v[i][VX3]*v[i][VX3]);
38  lor = sqrt(1.0 + lor);
39 
40  EXPAND(v[i][VX1] /= lor; ,
41  v[i][VX2] /= lor; ,
42  v[i][VX3] /= lor;)
43  }
44 }
VX2
#define VX2
Definition: mod_defs.h:29
VX1
#define VX1
Definition: mod_defs.h:28
pluto.h
PLUTO main header file.
i
int i
Definition: analysis.c:2
ConvertTo4vel
void ConvertTo4vel(double **v, int beg, int end)
Definition: four_vel.c:4
VX3
#define VX3
Definition: mod_defs.h:30
ConvertTo3vel
void ConvertTo3vel(double **v, int beg, int end)
Definition: four_vel.c:25