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
|
/* default order */
#ifndef p_N
#define p_N 6
#endif
#define NODETOL 1e-4
#ifdef NDG2d
#define p_Nfp (p_N+1)
#define p_Np ((p_N+1)*(p_N+2)/2)
#define p_Nfields 3
#define p_Nfaces 3
#endif
#ifdef NDG3d
#define p_Nfp ((p_N+1)*(p_N+2)/2)
#define p_Np ((p_N+1)*(p_N+2)*(p_N+3)/6)
#define p_Nfields 6
#define p_Nfaces 4
#endif
#define BSIZE (16*((p_Np+15)/16))
//#define BSIZE p_Np
#define max(a,b) ( (a>b)?a:b )
#define min(a,b) ( (a<b)?a:b )
/* list of vertices on each edge */
typedef struct foo {
int nprocs; /* number of processes */
int procid; /* number of this process */
int Nverts; /* number of vertices per element */
int Nfaces; /* number of faces per element */
int Nedges; /* number of faces per element (3d only) */
int Nv; /* number of mesh vertices */
int K; /* number of mesh elements */
int **EToV; /* element to global vertex list */
int **EToG; /* element to global face number list */
int **EToS; /* element to global edge number list */
int **EToE; /* element to neighbor element (elements numbered by their proc) */
int **EToF; /* element to neighbor face (element local number 0,1,2) */
int **EToP; /* element to neighbor proc (element prc number 0,1,.., nprocs-1) */
int **localEToV; /* element to process local vertex list */
int localNunique; /* number of unique nodes on this process */
int *bcflag; /* vector. entry n is 1 if vertex n is on a boundary */
double **GX; /* x-coordinates of element vertices */
double **GY; /* y-coordinates of element vertices */
double **GZ; /* z-coordinates of element vertices (3d) */
/* high order node info */
int **Fmask; /* face node numbers in element volume data */
double *r, *s, *t; /* (r,s) coordinates of reference nodes */
double **Dr, **Ds, **Dt; /* local nodal derivative matrices */
double **LIFT; /* local lift matrix */
double **x; /* x-coordinates of element nodes */
double **y; /* y-coordinates of element nodes */
double **z; /* z-coordinates of element nodes (3d) */
/* DG STUFF (EXPERIMENTAL) */
int *vmapM; /* volume id of -ve trace of face node */
int *vmapP; /* volume id of +ve trace of face node */
int *Npar; /* # of nodes to send recv to each proc */
int **parK; /* element of parallel nodes */
int **parF; /* face of parallel nodes */
/* MPI stuff */
int parNtotalout; /* total number of nodes to send */
int *parmapOUT; /* list of nodes to send out */
int *c_parmapOUT; /* device list */
float *f_outQ, *f_inQ; /* MPI data buffers */
/* float version of operators */
float *f_Dr, *f_Ds, *f_Dt;
float *f_LIFT;
/* float geometric info */
float *vgeo; /* geometric factors */
float *surfinfo;
/* float field storage (CPU) */
float *f_Q, *f_rhsQ, *f_resQ;
/* time stepping constants */
double *rk4a, *rk4b, *rk4c;
}Mesh;
|
