Chombo + EB + MF  3.2
AdvectPhysicsF_F.H
Go to the documentation of this file.
1 #ifndef _ADVECTPHYSICSF_F_H_
2 #define _ADVECTPHYSICSF_F_H_
3 
4 #include "FORT_PROTO.H"
5 #include "CH_Timer.H"
6 #include "REAL.H"
7 
8 extern "C"
9 {
10 
11 #ifndef GUARDAVEFACESCALTOFACEVECT
12 #define GUARDAVEFACESCALTOFACEVECT
13 // Prototype for Fortran procedure avefacescaltofacevect ...
14 //
15 void FORTRAN_NAME( AVEFACESCALTOFACEVECT ,avefacescaltofacevect )(
16  CHFp_FRA(facevect)
17  ,CHFp_CONST_FRA1(facescal)
18  ,CHFp_CONST_INT(facedir)
19  ,CHFp_CONST_INT(vectdir)
20  ,CHFp_BOX(dcalcface)
21  ,CHFp_BOX(offbox) );
22 
23 #define FORT_AVEFACESCALTOFACEVECT FORTRAN_NAME( inlineAVEFACESCALTOFACEVECT, inlineAVEFACESCALTOFACEVECT)
24 #define FORTNT_AVEFACESCALTOFACEVECT FORTRAN_NAME( AVEFACESCALTOFACEVECT, avefacescaltofacevect)
25 
26 inline void FORTRAN_NAME(inlineAVEFACESCALTOFACEVECT, inlineAVEFACESCALTOFACEVECT)(
27  CHFp_FRA(facevect)
28  ,CHFp_CONST_FRA1(facescal)
29  ,CHFp_CONST_INT(facedir)
30  ,CHFp_CONST_INT(vectdir)
31  ,CHFp_BOX(dcalcface)
32  ,CHFp_BOX(offbox) )
33 {
34  CH_TIMELEAF("FORT_AVEFACESCALTOFACEVECT");
35  FORTRAN_NAME( AVEFACESCALTOFACEVECT ,avefacescaltofacevect )(
36  CHFt_FRA(facevect)
37  ,CHFt_CONST_FRA1(facescal)
38  ,CHFt_CONST_INT(facedir)
39  ,CHFt_CONST_INT(vectdir)
40  ,CHFt_BOX(dcalcface)
41  ,CHFt_BOX(offbox) );
42 }
43 #endif // GUARDAVEFACESCALTOFACEVECT
44 
45 #ifndef GUARDAVESCALTOFACE
46 #define GUARDAVESCALTOFACE
47 // Prototype for Fortran procedure avescaltoface ...
48 //
49 void FORTRAN_NAME( AVESCALTOFACE ,avescaltoface )(
50  CHFp_FRA1(facescal)
51  ,CHFp_CONST_FRA1(cellscal)
52  ,CHFp_CONST_INT(idir)
53  ,CHFp_BOX(dcalcface) );
54 
55 #define FORT_AVESCALTOFACE FORTRAN_NAME( inlineAVESCALTOFACE, inlineAVESCALTOFACE)
56 #define FORTNT_AVESCALTOFACE FORTRAN_NAME( AVESCALTOFACE, avescaltoface)
57 
58 inline void FORTRAN_NAME(inlineAVESCALTOFACE, inlineAVESCALTOFACE)(
59  CHFp_FRA1(facescal)
60  ,CHFp_CONST_FRA1(cellscal)
61  ,CHFp_CONST_INT(idir)
62  ,CHFp_BOX(dcalcface) )
63 {
64  CH_TIMELEAF("FORT_AVESCALTOFACE");
65  FORTRAN_NAME( AVESCALTOFACE ,avescaltoface )(
66  CHFt_FRA1(facescal)
67  ,CHFt_CONST_FRA1(cellscal)
68  ,CHFt_CONST_INT(idir)
69  ,CHFt_BOX(dcalcface) );
70 }
71 #endif // GUARDAVESCALTOFACE
72 
73 #ifndef GUARDAVECELLTOFACE
74 #define GUARDAVECELLTOFACE
75 // Prototype for Fortran procedure avecelltoface ...
76 //
77 void FORTRAN_NAME( AVECELLTOFACE ,avecelltoface )(
78  CHFp_FRA1(facevel)
79  ,CHFp_CONST_FRA(cellvel)
80  ,CHFp_CONST_INT(idir)
81  ,CHFp_BOX(dcalcface) );
82 
83 #define FORT_AVECELLTOFACE FORTRAN_NAME( inlineAVECELLTOFACE, inlineAVECELLTOFACE)
84 #define FORTNT_AVECELLTOFACE FORTRAN_NAME( AVECELLTOFACE, avecelltoface)
85 
86 inline void FORTRAN_NAME(inlineAVECELLTOFACE, inlineAVECELLTOFACE)(
87  CHFp_FRA1(facevel)
88  ,CHFp_CONST_FRA(cellvel)
89  ,CHFp_CONST_INT(idir)
90  ,CHFp_BOX(dcalcface) )
91 {
92  CH_TIMELEAF("FORT_AVECELLTOFACE");
93  FORTRAN_NAME( AVECELLTOFACE ,avecelltoface )(
94  CHFt_FRA1(facevel)
95  ,CHFt_CONST_FRA(cellvel)
96  ,CHFt_CONST_INT(idir)
97  ,CHFt_BOX(dcalcface) );
98 }
99 #endif // GUARDAVECELLTOFACE
100 
101 #ifndef GUARDAVEFACETOCELL
102 #define GUARDAVEFACETOCELL
103 // Prototype for Fortran procedure avefacetocell ...
104 //
105 void FORTRAN_NAME( AVEFACETOCELL ,avefacetocell )(
106  CHFp_FRA(cellvel)
107  ,CHFp_CONST_FRA1(facevel)
108  ,CHFp_CONST_INT(idir)
109  ,CHFp_BOX(dcalccell) );
110 
111 #define FORT_AVEFACETOCELL FORTRAN_NAME( inlineAVEFACETOCELL, inlineAVEFACETOCELL)
112 #define FORTNT_AVEFACETOCELL FORTRAN_NAME( AVEFACETOCELL, avefacetocell)
113 
114 inline void FORTRAN_NAME(inlineAVEFACETOCELL, inlineAVEFACETOCELL)(
115  CHFp_FRA(cellvel)
116  ,CHFp_CONST_FRA1(facevel)
117  ,CHFp_CONST_INT(idir)
118  ,CHFp_BOX(dcalccell) )
119 {
120  CH_TIMELEAF("FORT_AVEFACETOCELL");
121  FORTRAN_NAME( AVEFACETOCELL ,avefacetocell )(
122  CHFt_FRA(cellvel)
123  ,CHFt_CONST_FRA1(facevel)
124  ,CHFt_CONST_INT(idir)
125  ,CHFt_BOX(dcalccell) );
126 }
127 #endif // GUARDAVEFACETOCELL
128 
129 #ifndef GUARDMAGNITUDEF
130 #define GUARDMAGNITUDEF
131 // Prototype for Fortran procedure MAGNITUDEF ...
132 //
133 void FORTRAN_NAME( MAGNITUDEF ,magnitudef )(
134  CHFp_FRA1(magdata)
135  ,CHFp_CONST_FRA(data)
136  ,CHFp_BOX(box) );
137 
138 #define FORT_MAGNITUDEF FORTRAN_NAME( inlineMAGNITUDEF, inlineMAGNITUDEF)
139 #define FORTNT_MAGNITUDEF FORTRAN_NAME( MAGNITUDEF, magnitudef)
140 
141 inline void FORTRAN_NAME(inlineMAGNITUDEF, inlineMAGNITUDEF)(
142  CHFp_FRA1(magdata)
143  ,CHFp_CONST_FRA(data)
145 {
146  CH_TIMELEAF("FORT_MAGNITUDEF");
147  FORTRAN_NAME( MAGNITUDEF ,magnitudef )(
148  CHFt_FRA1(magdata)
149  ,CHFt_CONST_FRA(data)
150  ,CHFt_BOX(box) );
151 }
152 #endif // GUARDMAGNITUDEF
153 
154 #ifndef GUARDGETRELGRADF
155 #define GUARDGETRELGRADF
156 // Prototype for Fortran procedure GETRELGRADF ...
157 //
158 void FORTRAN_NAME( GETRELGRADF ,getrelgradf )(
159  CHFp_FRA1(du)
160  ,CHFp_CONST_FRA1(u)
161  ,CHFp_CONST_INT(idir)
162  ,CHFp_BOX(loBox)
163  ,CHFp_CONST_INT(hasLo)
164  ,CHFp_BOX(hiBox)
165  ,CHFp_CONST_INT(hasHi)
166  ,CHFp_BOX(centerBox) );
167 
168 #define FORT_GETRELGRADF FORTRAN_NAME( inlineGETRELGRADF, inlineGETRELGRADF)
169 #define FORTNT_GETRELGRADF FORTRAN_NAME( GETRELGRADF, getrelgradf)
170 
171 inline void FORTRAN_NAME(inlineGETRELGRADF, inlineGETRELGRADF)(
172  CHFp_FRA1(du)
173  ,CHFp_CONST_FRA1(u)
174  ,CHFp_CONST_INT(idir)
175  ,CHFp_BOX(loBox)
176  ,CHFp_CONST_INT(hasLo)
177  ,CHFp_BOX(hiBox)
178  ,CHFp_CONST_INT(hasHi)
179  ,CHFp_BOX(centerBox) )
180 {
181  CH_TIMELEAF("FORT_GETRELGRADF");
182  FORTRAN_NAME( GETRELGRADF ,getrelgradf )(
183  CHFt_FRA1(du)
184  ,CHFt_CONST_FRA1(u)
185  ,CHFt_CONST_INT(idir)
186  ,CHFt_BOX(loBox)
187  ,CHFt_CONST_INT(hasLo)
188  ,CHFt_BOX(hiBox)
189  ,CHFt_CONST_INT(hasHi)
190  ,CHFt_BOX(centerBox) );
191 }
192 #endif // GUARDGETRELGRADF
193 
194 #ifndef GUARDPOSTNORMALSOURCE
195 #define GUARDPOSTNORMALSOURCE
196 // Prototype for Fortran procedure POSTNORMALSOURCE ...
197 //
198 void FORTRAN_NAME( POSTNORMALSOURCE ,postnormalsource )(
199  CHFp_FRA(dWminus)
200  ,CHFp_FRA(dWplus)
201  ,CHFp_CONST_FRA(W)
202  ,CHFp_CONST_FRA1(advVel)
203  ,CHFp_CONST_REAL(dt)
204  ,CHFp_CONST_REAL(dx)
205  ,CHFp_CONST_INT(idir)
206  ,CHFp_BOX(box) );
207 
208 #define FORT_POSTNORMALSOURCE FORTRAN_NAME( inlinePOSTNORMALSOURCE, inlinePOSTNORMALSOURCE)
209 #define FORTNT_POSTNORMALSOURCE FORTRAN_NAME( POSTNORMALSOURCE, postnormalsource)
210 
211 inline void FORTRAN_NAME(inlinePOSTNORMALSOURCE, inlinePOSTNORMALSOURCE)(
212  CHFp_FRA(dWminus)
213  ,CHFp_FRA(dWplus)
214  ,CHFp_CONST_FRA(W)
215  ,CHFp_CONST_FRA1(advVel)
216  ,CHFp_CONST_REAL(dt)
217  ,CHFp_CONST_REAL(dx)
218  ,CHFp_CONST_INT(idir)
219  ,CHFp_BOX(box) )
220 {
221  CH_TIMELEAF("FORT_POSTNORMALSOURCE");
222  FORTRAN_NAME( POSTNORMALSOURCE ,postnormalsource )(
223  CHFt_FRA(dWminus)
224  ,CHFt_FRA(dWplus)
225  ,CHFt_CONST_FRA(W)
226  ,CHFt_CONST_FRA1(advVel)
227  ,CHFt_CONST_REAL(dt)
228  ,CHFt_CONST_REAL(dx)
229  ,CHFt_CONST_INT(idir)
230  ,CHFt_BOX(box) );
231 }
232 #endif // GUARDPOSTNORMALSOURCE
233 
234 #ifndef GUARDRIEMANNF
235 #define GUARDRIEMANNF
236 // Prototype for Fortran procedure RIEMANNF ...
237 //
238 void FORTRAN_NAME( RIEMANNF ,riemannf )(
239  CHFp_FRA(Wgdnv)
240  ,CHFp_CONST_FRA(WLeft)
241  ,CHFp_CONST_FRA(WRight)
242  ,CHFp_CONST_FRA1(advVel)
243  ,CHFp_CONST_INT(idir)
244  ,CHFp_BOX(box) );
245 
246 #define FORT_RIEMANNF FORTRAN_NAME( inlineRIEMANNF, inlineRIEMANNF)
247 #define FORTNT_RIEMANNF FORTRAN_NAME( RIEMANNF, riemannf)
248 
249 inline void FORTRAN_NAME(inlineRIEMANNF, inlineRIEMANNF)(
250  CHFp_FRA(Wgdnv)
251  ,CHFp_CONST_FRA(WLeft)
252  ,CHFp_CONST_FRA(WRight)
253  ,CHFp_CONST_FRA1(advVel)
254  ,CHFp_CONST_INT(idir)
255  ,CHFp_BOX(box) )
256 {
257  CH_TIMELEAF("FORT_RIEMANNF");
258  FORTRAN_NAME( RIEMANNF ,riemannf )(
259  CHFt_FRA(Wgdnv)
260  ,CHFt_CONST_FRA(WLeft)
261  ,CHFt_CONST_FRA(WRight)
262  ,CHFt_CONST_FRA1(advVel)
263  ,CHFt_CONST_INT(idir)
264  ,CHFt_BOX(box) );
265 }
266 #endif // GUARDRIEMANNF
267 
268 #ifndef GUARDQUASILINEARUPDATE
269 #define GUARDQUASILINEARUPDATE
270 // Prototype for Fortran procedure QUASILINEARUPDATE ...
271 //
272 void FORTRAN_NAME( QUASILINEARUPDATE ,quasilinearupdate )(
273  CHFp_FRA(AdWdx)
274  ,CHFp_CONST_FRA(WHalf)
275  ,CHFp_CONST_FRA1(cellVel)
277  ,CHFp_CONST_INT(idir)
278  ,CHFp_BOX(box) );
279 
280 #define FORT_QUASILINEARUPDATE FORTRAN_NAME( inlineQUASILINEARUPDATE, inlineQUASILINEARUPDATE)
281 #define FORTNT_QUASILINEARUPDATE FORTRAN_NAME( QUASILINEARUPDATE, quasilinearupdate)
282 
283 inline void FORTRAN_NAME(inlineQUASILINEARUPDATE, inlineQUASILINEARUPDATE)(
284  CHFp_FRA(AdWdx)
285  ,CHFp_CONST_FRA(WHalf)
286  ,CHFp_CONST_FRA1(cellVel)
288  ,CHFp_CONST_INT(idir)
289  ,CHFp_BOX(box) )
290 {
291  CH_TIMELEAF("FORT_QUASILINEARUPDATE");
292  FORTRAN_NAME( QUASILINEARUPDATE ,quasilinearupdate )(
293  CHFt_FRA(AdWdx)
294  ,CHFt_CONST_FRA(WHalf)
295  ,CHFt_CONST_FRA1(cellVel)
297  ,CHFt_CONST_INT(idir)
298  ,CHFt_BOX(box) );
299 }
300 #endif // GUARDQUASILINEARUPDATE
301 
302 #ifndef GUARDAVECELLVECTOFACEVEC
303 #define GUARDAVECELLVECTOFACEVEC
304 // Prototype for Fortran procedure avecellvectofacevec ...
305 //
306 void FORTRAN_NAME( AVECELLVECTOFACEVEC ,avecellvectofacevec )(
307  CHFp_FRA(facevec)
308  ,CHFp_CONST_FRA(cellvec)
309  ,CHFp_CONST_INT(facedir)
310  ,CHFp_BOX(dcalcface) );
311 
312 #define FORT_AVECELLVECTOFACEVEC FORTRAN_NAME( inlineAVECELLVECTOFACEVEC, inlineAVECELLVECTOFACEVEC)
313 #define FORTNT_AVECELLVECTOFACEVEC FORTRAN_NAME( AVECELLVECTOFACEVEC, avecellvectofacevec)
314 
315 inline void FORTRAN_NAME(inlineAVECELLVECTOFACEVEC, inlineAVECELLVECTOFACEVEC)(
316  CHFp_FRA(facevec)
317  ,CHFp_CONST_FRA(cellvec)
318  ,CHFp_CONST_INT(facedir)
319  ,CHFp_BOX(dcalcface) )
320 {
321  CH_TIMELEAF("FORT_AVECELLVECTOFACEVEC");
322  FORTRAN_NAME( AVECELLVECTOFACEVEC ,avecellvectofacevec )(
323  CHFt_FRA(facevec)
324  ,CHFt_CONST_FRA(cellvec)
325  ,CHFt_CONST_INT(facedir)
326  ,CHFt_BOX(dcalcface) );
327 }
328 #endif // GUARDAVECELLVECTOFACEVEC
329 
330 }
331 
332 #endif
#define CHFt_CONST_INT(i)
Definition: FORT_PROTO.H:989
void CHFp_CONST_FRA1(facescal)
#define CHFp_FRA1(a)
Definition: FORT_PROTO.H:913
void CHFp_CONST_FRA(cellvel)
void CHFp_BOX(dcalcface)
Definition: AdvectPhysicsF_F.H:62
#define CHFt_FRA(a)
Definition: FORT_PROTO.H:1053
IndexTM< T, N > scale(const IndexTM< T, N > &a_p, T a_s)
Definition: IndexTMI.H:382
void CHFp_CONST_REAL(dt)
Definition: EBPatchGodunovF_F.H:704
Definition: EBInterface.H:45
#define CHFt_CONST_FRA(a)
Definition: FORT_PROTO.H:1095
void CHFp_CONST_INT(facedir)
void FORTRAN_NAME(AVEFACESCALTOFACEVECT, avefacescaltofacevect)(CHFp_FRA(facevect)
void CHFp_FRA(dWplus)
#define CHFt_FRA1(a)
Definition: FORT_PROTO.H:1136
#define CHFt_CONST_FRA1(a)
Definition: FORT_PROTO.H:1188
#define CHFt_BOX(b)
Definition: FORT_PROTO.H:1026
#define CH_TIMELEAF(name)
Definition: CH_Timer.H:100
#define CHFt_CONST_REAL(r)
Definition: FORT_PROTO.H:991