-
Notifications
You must be signed in to change notification settings - Fork 0
/
riemann_toro.h
48 lines (33 loc) · 907 Bytes
/
riemann_toro.h
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
#ifndef RIEMANN_TORO_BRANT
#define RIEMANN_TORO_BRANT
//#define LOUD_RIEMANN
struct RiemannStructure
{
double g1,g2,g3,g4,g5,g6,g7,g8;
double rho_L;
double u_L;
double p_L;
double c_sound_L;
double rho_R;
double u_R;
double p_R;
double c_sound_R;
double gamma;
double mpa;
};
class RiemannProblem
{
public:
RiemannStructure r;
int riemann(double s, double stl[5], double str[5], double sta[5]);
double guess_p_riemann(void);
void star_pu_riemann(double *pget, double *uget);
void prefun_riemann(double *f, double *fd, double p, double dk, double pk, double ck);
void sample_riemann(double s, double pm, double um, double *rho_s, double *u_s, double *p_s);
void Test_Problem_Toro_1(void);
void Test_Problem_Toro_2(void);
void Test_Problem_Toro_3(void);
void Test_Problem_Toro_4(void);
void Test_Problem_Toro_5(void);
};
#endif //RIEMANN_TORO_BRANT