PLUTO
four_vel.c File Reference
#include "pluto.h"
Include dependency graph for four_vel.c:

Go to the source code of this file.

Functions

void ConvertTo4vel (double **v, int beg, int end)
 
void ConvertTo3vel (double **v, int beg, int end)
 

Function Documentation

void ConvertTo3vel ( double **  v,
int  beg,
int  end 
)

Definition at line 25 of file four_vel.c.

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 }
int end
Global end index for the local array.
Definition: structs.h:116
#define VX2
Definition: mod_defs.h:29
#define VX1
Definition: mod_defs.h:28
int beg
Global start index for the local array.
Definition: structs.h:115
int i
Definition: analysis.c:2
#define VX3
Definition: mod_defs.h:30

Here is the caller graph for this function:

void ConvertTo4vel ( double **  v,
int  beg,
int  end 
)

Definition at line 4 of file four_vel.c.

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 }
int end
Global end index for the local array.
Definition: structs.h:116
#define VX2
Definition: mod_defs.h:29
#define VX1
Definition: mod_defs.h:28
int beg
Global start index for the local array.
Definition: structs.h:115
int i
Definition: analysis.c:2
#define VX3
Definition: mod_defs.h:30

Here is the caller graph for this function: