Athena++/Atmosphere
Planetary Atmosphere Simulator
coordinates.hpp
Go to the documentation of this file.
1 #ifndef COORDINATES_COORDINATES_HPP_
2 #define COORDINATES_COORDINATES_HPP_
3 //========================================================================================
4 // Athena++ astrophysical MHD code
5 // Copyright(C) 2014 James M. Stone <jmstone@princeton.edu> and other code contributors
6 // Licensed under the 3-clause BSD License, see LICENSE file for details
7 //========================================================================================
9 // \brief defines abstract base and derived classes for coordinates. These classes
10 // provide data and functions to compute/store coordinate positions and spacing, as well
11 // as geometrical factors (areas, volumes, coordinate source terms) for various
12 // coordinate systems.
13 
14 // C headers
15 
16 // C++ headers
17 #include <iostream>
18 
19 // Athena++ headers
20 #include "../athena.hpp"
21 #include "../athena_arrays.hpp"
22 #include "../hydro/srcterms/hydro_srcterms.hpp"
23 #include "../mesh/mesh.hpp"
24 
25 // forward declarations
26 class MeshBlock;
27 class ParameterInput;
28 
29 //----------------------------------------------------------------------------------------
31 // \brief abstract base class for all coordinate derived classes
32 
33 class Coordinates {
34  public:
35  friend class HydroSourceTerms;
36  Coordinates(MeshBlock *pmb, ParameterInput *pin, bool flag = false);
37  virtual ~Coordinates() = default;
38 
39  // data
40  MeshBlock *pmy_block; // ptr to MeshBlock containing this Coordinates
41  AthenaArray<Real> dx1f, dx2f, dx3f, x1f, x2f, x3f; // face spacing and positions
42  AthenaArray<Real> dx1v, dx2v, dx3v, x1v, x2v, x3v; // volume spacing and positions
43  AthenaArray<Real> x1s2, x1s3, x2s1, x2s3, x3s1, x3s2; // area averaged positions for AMR
44  // geometry coefficients (only used in SphericalPolar, Cylindrical, Cartesian)
45  AthenaArray<Real> h2f, dh2fd1, h31f, h32f, dh31fd1, dh32fd2;
46  AthenaArray<Real> h2v, dh2vd1, h31v, h32v, dh31vd1, dh32vd2;
47 
48  // functions...
49  // ...to compute length of edges
50  virtual void Edge1Length(const int k, const int j, const int il, const int iu,
51  AthenaArray<Real> &len);
52  virtual void Edge2Length(const int k, const int j, const int il, const int iu,
53  AthenaArray<Real> &len);
54  virtual void Edge3Length(const int k, const int j, const int il, const int iu,
55  AthenaArray<Real> &len);
56  virtual Real GetEdge1Length(const int k, const int j, const int i);
57  virtual Real GetEdge2Length(const int k, const int j, const int i);
58  virtual Real GetEdge3Length(const int k, const int j, const int i);
59  // ...to compute length connecting cell centers (for non-ideal MHD)
60  virtual void VolCenter1Length(const int k, const int j, const int il, const int iu,
61  AthenaArray<Real> &len);
62  virtual void VolCenter2Length(const int k, const int j, const int il, const int iu,
63  AthenaArray<Real> &len);
64  virtual void VolCenter3Length(const int k, const int j, const int il, const int iu,
65  AthenaArray<Real> &len);
66  // ...to compute physical width at cell center
67  virtual void CenterWidth1(const int k, const int j, const int il, const int iu,
68  AthenaArray<Real> &dx1);
69  virtual void CenterWidth2(const int k, const int j, const int il, const int iu,
70  AthenaArray<Real> &dx2);
71  virtual void CenterWidth3(const int k, const int j, const int il, const int iu,
72  AthenaArray<Real> &dx3);
73 
74  // ...to compute area of faces
75  virtual void Face1Area(const int k, const int j, const int il, const int iu,
76  AthenaArray<Real> &area);
77  virtual void Face2Area(const int k, const int j, const int il, const int iu,
78  AthenaArray<Real> &area);
79  virtual void Face3Area(const int k, const int j, const int il, const int iu,
80  AthenaArray<Real> &area);
81  virtual Real GetFace1Area(const int k, const int j, const int i);
82  virtual Real GetFace2Area(const int k, const int j, const int i);
83  virtual Real GetFace3Area(const int k, const int j, const int i);
84  // ...to compute area of faces joined by cell centers (for non-ideal MHD)
85  virtual void VolCenterFace1Area(const int k, const int j, const int il, const int iu,
86  AthenaArray<Real> &area);
87  virtual void VolCenterFace2Area(const int k, const int j, const int il, const int iu,
88  AthenaArray<Real> &area);
89  virtual void VolCenterFace3Area(const int k, const int j, const int il, const int iu,
90  AthenaArray<Real> &area);
91 
92  // ...to compute Laplacian of quantities in the coord system and orthogonal subspaces
93  virtual void Laplacian(
94  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
95  const int il, const int iu, const int jl, const int ju, const int kl, const int ku,
96  const int nl, const int nu);
97  virtual void LaplacianX1(
98  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
99  const int n, const int k, const int j, const int il, const int iu);
100  virtual void LaplacianX1All(
101  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
102  const int nl, const int nu, const int kl, const int ku,
103  const int jl, const int ju, const int il, const int iu);
104  virtual void LaplacianX2(
105  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
106  const int n, const int k, const int j, const int il, const int iu);
107  virtual void LaplacianX2All(
108  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
109  const int nl, const int nu, const int kl, const int ku,
110  const int jl, const int ju, const int il, const int iu);
111  virtual void LaplacianX3(
112  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
113  const int n, const int k, const int j, const int il, const int iu);
114  virtual void LaplacianX3All(
115  const AthenaArray<Real> &s, AthenaArray<Real> &delta_s,
116  const int nl, const int nu, const int kl, const int ku,
117  const int jl, const int ju, const int il, const int iu);
118 
119  // ...to compute volume of cells
120  virtual void CellVolume(const int k, const int j, const int il, const int iu,
121  AthenaArray<Real> &vol);
122  virtual Real GetCellVolume(const int k, const int j, const int i);
123 
124  // ...to compute geometrical source terms
125  virtual void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
126  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
127  AthenaArray<Real> &u);
128 
129  // ...to determine if index is a pole
130  bool IsPole(int j);
131 
132 
133  // In GR, functions...
134  // ...to return private variables
135  Real GetMass() const {return bh_mass_;}
136  Real GetSpin() const {return bh_spin_;}
137 
138  // ...to compute metric
139  void Metric(
140  Real x1, Real x2, Real x3, ParameterInput *pin, AthenaArray<Real> &g,
141  AthenaArray<Real> &g_inv, AthenaArray<Real> &dg_dx1, AthenaArray<Real> &dg_dx2,
142  AthenaArray<Real> &dg_dx3);
143  virtual void CellMetric(const int k, const int j, const int il, const int iu,
144  AthenaArray<Real> &g, AthenaArray<Real> &gi) {}
145  virtual void Face1Metric(const int k, const int j, const int il, const int iu,
146  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) {}
147  virtual void Face2Metric(const int k, const int j, const int il, const int iu,
148  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) {}
149  virtual void Face3Metric(const int k, const int j, const int il, const int iu,
150  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) {}
151 
152  // ...to transform primitives to locally flat space
153  virtual void PrimToLocal1(
154  const int k, const int j, const int il, const int iu,
155  const AthenaArray<Real> &b1_vals, AthenaArray<Real> &prim_left,
156  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) {}
157  virtual void PrimToLocal2(
158  const int k, const int j, const int il, const int iu,
159  const AthenaArray<Real> &b2_vals, AthenaArray<Real> &prim_left,
160  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) {}
161  virtual void PrimToLocal3(
162  const int k, const int j, const int il, const int iu,
163  const AthenaArray<Real> &b3_vals, AthenaArray<Real> &prim_left,
164  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) {}
165 
166  // ...to transform fluxes in locally flat space to global frame
167  virtual void FluxToGlobal1(
168  const int k, const int j, const int il, const int iu,
169  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
170  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) {}
171  virtual void FluxToGlobal2(
172  const int k, const int j, const int il, const int iu,
173  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
174  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) {}
175  virtual void FluxToGlobal3(
176  const int k, const int j, const int il, const int iu,
177  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
178  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) {}
179 
180  // ...to raise (lower) covariant (contravariant) components of a vector
181  virtual void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j,
182  int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3) {}
183  virtual void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i,
184  Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) {}
185 
186  protected:
187  bool coarse_flag; // true if this coordinate object is parent (coarse) mesh in AMR
188  Mesh *pm;
189  int il, iu, jl, ju, kl, ku, ng; // limits of indices of arrays (normal or coarse)
190  int nc1, nc2, nc3; // # cells in each dir of arrays (normal or coarse)
191  // Scratch arrays for coordinate factors
192  // Format: coord_<type>[<direction>]_<index>[<count>]_
193  // type: vol[ume], area, etc.
194  // direction: 1/2/3 depending on which face, edge, etc. is in play
195  // index: i/j/k indicating which coordinates index array
196  // count: 1/2/... in case multiple arrays are needed for different terms
197  AthenaArray<Real> coord_vol_i_, coord_vol_i1_, coord_vol_i2_;
198  AthenaArray<Real> coord_vol_j_, coord_vol_j1_, coord_vol_j2_;
199  AthenaArray<Real> coord_vol_k1_;
200  AthenaArray<Real> coord_vol_kji_;
201  AthenaArray<Real> coord_area1_i_, coord_area1_i1_;
202  AthenaArray<Real> coord_area1_j_, coord_area1_j1_, coord_area1_j2_;
203  AthenaArray<Real> coord_area1_k1_;
204  AthenaArray<Real> coord_area1_kji_;
205  AthenaArray<Real> coord_area2_i_, coord_area2_i1_, coord_area2_i2_;
206  AthenaArray<Real> coord_area2_j_, coord_area2_j1_, coord_area2_j2_;
207  AthenaArray<Real> coord_area2_k1_;
208  AthenaArray<Real> coord_area2_kji_;
209  AthenaArray<Real> coord_area3_i_, coord_area3_i1_, coord_area3_i2_;
210  AthenaArray<Real> coord_area3_j1_, coord_area3_j2_;
211  AthenaArray<Real> coord_area3_kji_;
212  AthenaArray<Real> coord_area1vc_i_,coord_area1vc_j_; //nonidealmhd additions
213  AthenaArray<Real> coord_area2vc_i_,coord_area2vc_j_; //nonidealmhd additions
214  AthenaArray<Real> coord_area3vc_i_; //nonidealmhd addition
215  AthenaArray<Real> coord_len1_i1_, coord_len1_i2_;
216  AthenaArray<Real> coord_len1_j1_, coord_len1_j2_;
217  AthenaArray<Real> coord_len1_kji_;
218  AthenaArray<Real> coord_len2_i1_;
219  AthenaArray<Real> coord_len2_j1_, coord_len2_j2_;
220  AthenaArray<Real> coord_len2_kji_;
221  AthenaArray<Real> coord_len3_i1_;
222  AthenaArray<Real> coord_len3_j1_, coord_len3_j2_;
223  AthenaArray<Real> coord_len3_k1_;
224  AthenaArray<Real> coord_len3_kji_;
225  AthenaArray<Real> coord_width1_i1_;
226  AthenaArray<Real> coord_width1_kji_;
227  AthenaArray<Real> coord_width2_i1_;
228  AthenaArray<Real> coord_width2_j1_;
229  AthenaArray<Real> coord_width2_kji_;
230  AthenaArray<Real> coord_width3_j1_, coord_width3_j2_, coord_width3_j3_;
231  AthenaArray<Real> coord_width3_k1_;
232  AthenaArray<Real> coord_width3_ji1_;
233  AthenaArray<Real> coord_width3_kji_;
234  AthenaArray<Real> coord_src_j1_, coord_src_j2_;
235  AthenaArray<Real> coord_src_kji_;
236  AthenaArray<Real> coord_src1_i_;
237  AthenaArray<Real> coord_src1_j_;
238  AthenaArray<Real> coord_src2_i_;
239  AthenaArray<Real> coord_src2_j_;
240  AthenaArray<Real> coord_src3_j_;
241 
242  // Scratch arrays for physical source terms
243  AthenaArray<Real> phy_src1_i_, phy_src2_i_;
244 
245  // GR-specific scratch arrays
246  AthenaArray<Real> metric_cell_i1_, metric_cell_i2_;
247  AthenaArray<Real> metric_cell_j1_, metric_cell_j2_;
248  AthenaArray<Real> metric_cell_kji_;
249  AthenaArray<Real> metric_face1_i1_, metric_face1_i2_;
250  AthenaArray<Real> metric_face1_j1_, metric_face1_j2_;
251  AthenaArray<Real> metric_face1_kji_;
252  AthenaArray<Real> metric_face2_i1_, metric_face2_i2_;
253  AthenaArray<Real> metric_face2_j1_, metric_face2_j2_;
254  AthenaArray<Real> metric_face2_kji_;
255  AthenaArray<Real> metric_face3_i1_, metric_face3_i2_;
256  AthenaArray<Real> metric_face3_j1_, metric_face3_j2_;
257  AthenaArray<Real> metric_face3_kji_;
258  AthenaArray<Real> trans_face1_i1_, trans_face1_i2_;
259  AthenaArray<Real> trans_face1_j1_;
260  AthenaArray<Real> trans_face1_ji1_, trans_face1_ji2_, trans_face1_ji3_,
261  trans_face1_ji4_, trans_face1_ji5_, trans_face1_ji6_, trans_face1_ji7_;
262  AthenaArray<Real> trans_face1_kji_;
263  AthenaArray<Real> trans_face2_i1_, trans_face2_i2_;
264  AthenaArray<Real> trans_face2_j1_;
265  AthenaArray<Real> trans_face2_ji1_, trans_face2_ji2_, trans_face2_ji3_,
266  trans_face2_ji4_, trans_face2_ji5_, trans_face2_ji6_;
267  AthenaArray<Real> trans_face2_kji_;
268  AthenaArray<Real> trans_face3_i1_, trans_face3_i2_;
269  AthenaArray<Real> trans_face3_j1_;
270  AthenaArray<Real> trans_face3_ji1_, trans_face3_ji2_, trans_face3_ji3_,
271  trans_face3_ji4_, trans_face3_ji5_, trans_face3_ji6_;
272  AthenaArray<Real> trans_face3_kji_;
273  AthenaArray<Real> g_, gi_;
274 
275  // GR-specific variables
276  Real bh_mass_;
277  Real bh_spin_;
278 };
279 
280 //----------------------------------------------------------------------------------------
282 // \brief derived class for Cartesian coordinates. None of the virtual funcs
283 // in the Coordinates abstract base class need to be overridden.
284 
285 class Cartesian : public Coordinates {
286  friend class HydroSourceTerms;
287 
288  public:
289  Cartesian(MeshBlock *pmb, ParameterInput *pin, bool flag);
290 };
291 
292 //----------------------------------------------------------------------------------------
294 // \brief derived class for Cylindrical coordinates. Some of the length, area,
295 // and volume functions in the Coordinates abstract base class are overridden
296 
297 class Cylindrical : public Coordinates {
298  friend class HydroSourceTerms;
299 
300  public:
301  Cylindrical(MeshBlock *pmb, ParameterInput *pin, bool flag);
302 
303  // functions...
304  // ...to compute length of edges
305  void Edge2Length(const int k, const int j, const int il, const int iu,
306  AthenaArray<Real> &len) final;
307  Real GetEdge2Length(const int k, const int j, const int i) final;
308  // ...to compute length connecting cell centers (for non-ideal MHD)
309  void VolCenter2Length(const int k, const int j, const int il, const int iu,
310  AthenaArray<Real> &len) final;
311  // ...to compute physical width at cell center
312  void CenterWidth2(const int k, const int j, const int il, const int iu,
313  AthenaArray<Real> &dx2) final;
314 
315  // ...to compute area of faces
316  void Face1Area(const int k, const int j, const int il, const int iu,
317  AthenaArray<Real> &area) final;
318  void Face3Area(const int k, const int j, const int il, const int iu,
319  AthenaArray<Real> &area) final;
320  Real GetFace1Area(const int k, const int j, const int i) final;
321  Real GetFace3Area(const int k, const int j, const int i) final;
322  // ...to compute area of faces joined by cell centers (for non-ideal MHD)
323  void VolCenterFace1Area(const int k, const int j, const int il, const int iu,
324  AthenaArray<Real> &area) final;
325  void VolCenterFace3Area(const int k, const int j, const int il, const int iu,
326  AthenaArray<Real> &area) final;
327  // ...to compute volumes of cells
328  void CellVolume(const int k, const int j, const int il, const int iu,
329  AthenaArray<Real> &vol) final;
330  Real GetCellVolume(const int k, const int j, const int i) final;
331 
332  // ...to compute geometrical source terms
333  void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
334  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
335  AthenaArray<Real> &u) final;
336 };
337 
338 //----------------------------------------------------------------------------------------
340 // \brief derived class for spherical polar coordinates. Many of the length, area,
341 // and volume functions in the Coordinates abstract base class are overridden.
342 
343 class SphericalPolar : public Coordinates {
344  friend class HydroSourceTerms;
345 
346  public:
347  SphericalPolar(MeshBlock *pmb, ParameterInput *pin, bool flag);
348 
349  // functions...
350  // ...to compute length of edges
351  void Edge2Length(const int k, const int j, const int il, const int iu,
352  AthenaArray<Real> &len) final;
353  void Edge3Length(const int k, const int j, const int il, const int iu,
354  AthenaArray<Real> &len) final;
355  Real GetEdge2Length(const int k, const int j, const int i) final;
356  Real GetEdge3Length(const int k, const int j, const int i) final;
357  // ...to compute length connecting cell centers (for non-ideal MHD)
358  void VolCenter2Length(const int k, const int j, const int il, const int iu,
359  AthenaArray<Real> &len) final;
360  void VolCenter3Length(const int k, const int j, const int il, const int iu,
361  AthenaArray<Real> &len) final;
362  // ...to compute physical width at cell center
363  void CenterWidth2(const int k, const int j, const int il, const int iu,
364  AthenaArray<Real> &dx2) final;
365  void CenterWidth3(const int k, const int j, const int il, const int iu,
366  AthenaArray<Real> &dx3) final;
367 
368  // ...to compute area of faces
369  void Face1Area(const int k, const int j, const int il, const int iu,
370  AthenaArray<Real> &area) final;
371  void Face2Area(const int k, const int j, const int il, const int iu,
372  AthenaArray<Real> &area) final;
373  void Face3Area(const int k, const int j, const int il, const int iu,
374  AthenaArray<Real> &area) final;
375  Real GetFace1Area(const int k, const int j, const int i) final;
376  Real GetFace2Area(const int k, const int j, const int i) final;
377  Real GetFace3Area(const int k, const int j, const int i) final;
378  // ...to compute area of faces joined by cell centers (for non-ideal MHD)
379  void VolCenterFace1Area(const int k, const int j, const int il, const int iu,
380  AthenaArray<Real> &area) final;
381  void VolCenterFace2Area(const int k, const int j, const int il, const int iu,
382  AthenaArray<Real> &area) final;
383  void VolCenterFace3Area(const int k, const int j, const int il, const int iu,
384  AthenaArray<Real> &area) final;
385  // ...to compute volumes of cells
386  void CellVolume(const int k, const int j, const int il, const int iu,
387  AthenaArray<Real> &vol) final;
388  Real GetCellVolume(const int k, const int j, const int i) final;
389 
390  // ...to compute geometrical source terms
391  void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
392  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
393  AthenaArray<Real> &u) final;
394 };
395 
396 //----------------------------------------------------------------------------------------
398 // \brief derived class for Minkowski (flat) spacetime and Cartesian coordinates in GR.
399 // None of the length, area, and volume functions in the abstract base class need to be
400 // overridden, but all the metric and transforms functions are.
401 
402 class Minkowski : public Coordinates {
403  friend class HydroSourceTerms;
404 
405  public:
406  Minkowski(MeshBlock *pmb, ParameterInput *pin, bool flag);
407 
408  // In GR, functions...
409  // ...to compute metric
410  void CellMetric(const int k, const int j, const int il, const int iu,
411  AthenaArray<Real> &g, AthenaArray<Real> &gi) final;
412  void Face1Metric(const int k, const int j, const int il, const int iu,
413  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
414  void Face2Metric(const int k, const int j, const int il, const int iu,
415  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
416  void Face3Metric(const int k, const int j, const int il, const int iu,
417  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
418 
419  // ...to transform primitives to locally flat space
420  void PrimToLocal1(const int k, const int j, const int il, const int iu,
421  const AthenaArray<Real> &b1_vals, AthenaArray<Real> &prim_left,
422  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
423  void PrimToLocal2(const int k, const int j, const int il, const int iu,
424  const AthenaArray<Real> &b2_vals, AthenaArray<Real> &prim_left,
425  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
426  void PrimToLocal3(const int k, const int j, const int il, const int iu,
427  const AthenaArray<Real> &b3_vals, AthenaArray<Real> &prim_left,
428  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
429 
430  // ...to transform fluxes in locally flat space to global frame
431  void FluxToGlobal1(
432  const int k, const int j, const int il, const int iu,
433  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
434  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
435  void FluxToGlobal2(
436  const int k, const int j, const int il, const int iu,
437  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
438  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
439  void FluxToGlobal3(
440  const int k, const int j, const int il, const int iu,
441  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
442  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
443 
444  // for raising (lowering) covariant (contravariant) components of a vector
445  void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i,
446  Real *pa0, Real *pa1, Real *pa2, Real *pa3) final;
447  void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i,
448  Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final;
449 };
450 
451 //----------------------------------------------------------------------------------------
453 // \brief derived class for Schwarzschild spacetime and spherical polar coordinates in GR
454 // Nearly every function in the abstract base class need to be overridden.
455 
456 class Schwarzschild : public Coordinates {
457  friend class HydroSourceTerms;
458 
459  public:
460  Schwarzschild(MeshBlock *pmb, ParameterInput *pin, bool flag);
461 
462  // functions...
463  // ...to compute length of edges
464  void Edge1Length(const int k, const int j, const int il, const int iu,
465  AthenaArray<Real> &len) final;
466  void Edge2Length(const int k, const int j, const int il, const int iu,
467  AthenaArray<Real> &len) final;
468  void Edge3Length(const int k, const int j, const int il, const int iu,
469  AthenaArray<Real> &len) final;
470  Real GetEdge1Length(const int k, const int j, const int i) final;
471  Real GetEdge2Length(const int k, const int j, const int i) final;
472  Real GetEdge3Length(const int k, const int j, const int i) final;
473 
474  // ...to compute physical width at cell center
475  void CenterWidth1(const int k, const int j, const int il, const int iu,
476  AthenaArray<Real> &dx1) final;
477  void CenterWidth2(const int k, const int j, const int il, const int iu,
478  AthenaArray<Real> &dx2) final;
479  void CenterWidth3(const int k, const int j, const int il, const int iu,
480  AthenaArray<Real> &dx3) final;
481 
482  // ...to compute area of faces
483  void Face1Area(const int k, const int j, const int il, const int iu,
484  AthenaArray<Real> &area) final;
485  void Face2Area(const int k, const int j, const int il, const int iu,
486  AthenaArray<Real> &area) final;
487  void Face3Area(const int k, const int j, const int il, const int iu,
488  AthenaArray<Real> &area) final;
489  Real GetFace1Area(const int k, const int j, const int i) final;
490  Real GetFace2Area(const int k, const int j, const int i) final;
491  Real GetFace3Area(const int k, const int j, const int i) final;
492 
493  // ...to compute volumes of cells
494  void CellVolume(const int k, const int j, const int il, const int iu,
495  AthenaArray<Real> &vol) final;
496  Real GetCellVolume(const int k, const int j, const int i) final;
497 
498  // ...to compute geometrical source terms
499  void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
500  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
501  AthenaArray<Real> &u) final;
502 
503  // In GR, functions...
504  // ...to compute metric
505  void CellMetric(const int k, const int j, const int il, const int iu,
506  AthenaArray<Real> &g, AthenaArray<Real> &gi) final;
507  void Face1Metric(const int k, const int j, const int il, const int iu,
508  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
509  void Face2Metric(const int k, const int j, const int il, const int iu,
510  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
511  void Face3Metric(const int k, const int j, const int il, const int iu,
512  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
513 
514  // ...to transform primitives to locally flat space
515  void PrimToLocal1(const int k, const int j, const int il, const int iu,
516  const AthenaArray<Real> &b1_vals, AthenaArray<Real> &prim_left,
517  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
518  void PrimToLocal2(const int k, const int j, const int il, const int iu,
519  const AthenaArray<Real> &b2_vals, AthenaArray<Real> &prim_left,
520  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
521  void PrimToLocal3(const int k, const int j, const int il, const int iu,
522  const AthenaArray<Real> &b3_vals, AthenaArray<Real> &prim_left,
523  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
524 
525  // ...to transform fluxes in locally flat space to global frame
526  void FluxToGlobal1(
527  const int k, const int j, const int il, const int iu,
528  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
529  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
530  void FluxToGlobal2(
531  const int k, const int j, const int il, const int iu,
532  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
533  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
534  void FluxToGlobal3(
535  const int k, const int j, const int il, const int iu,
536  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
537  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
538 
539  // for raising (lowering) covariant (contravariant) components of a vector
540  void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i,
541  Real *pa0, Real *pa1, Real *pa2, Real *pa3) final;
542  void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i,
543  Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final;
544 };
545 
546 //----------------------------------------------------------------------------------------
548 // \brief derived class for Kerr spacetime and Kerr-Schild coordinates in GR.
549 // Nearly every function in the abstract base class need to be overridden.
550 
551 class KerrSchild : public Coordinates {
552  friend class HydroSourceTerms;
553 
554  public:
555  KerrSchild(MeshBlock *pmb, ParameterInput *pin, bool flag);
556 
557  // functions...
558  // ...to compute length of edges
559  void Edge1Length(const int k, const int j, const int il, const int iu,
560  AthenaArray<Real> &len) final;
561  void Edge2Length(const int k, const int j, const int il, const int iu,
562  AthenaArray<Real> &len) final;
563  void Edge3Length(const int k, const int j, const int il, const int iu,
564  AthenaArray<Real> &len) final;
565  Real GetEdge1Length(const int k, const int j, const int i) final;
566  Real GetEdge2Length(const int k, const int j, const int i) final;
567  Real GetEdge3Length(const int k, const int j, const int i) final;
568 
569  // ...to compute physical width at cell center
570  void CenterWidth1(const int k, const int j, const int il, const int iu,
571  AthenaArray<Real> &dx1) final;
572  void CenterWidth2(const int k, const int j, const int il, const int iu,
573  AthenaArray<Real> &dx2) final;
574  void CenterWidth3(const int k, const int j, const int il, const int iu,
575  AthenaArray<Real> &dx3) final;
576 
577  // ...to compute area of faces
578  void Face1Area(const int k, const int j, const int il, const int iu,
579  AthenaArray<Real> &area) final;
580  void Face2Area(const int k, const int j, const int il, const int iu,
581  AthenaArray<Real> &area) final;
582  void Face3Area(const int k, const int j, const int il, const int iu,
583  AthenaArray<Real> &area) final;
584  Real GetFace1Area(const int k, const int j, const int i) final;
585  Real GetFace2Area(const int k, const int j, const int i) final;
586  Real GetFace3Area(const int k, const int j, const int i) final;
587 
588  // ...to compute volumes of cells
589  void CellVolume(const int k, const int j, const int il, const int iu,
590  AthenaArray<Real> &vol) final;
591  Real GetCellVolume(const int k, const int j, const int i) final;
592 
593  // ...to compute geometrical source terms
594  void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
595  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
596  AthenaArray<Real> &u) final;
597 
598  // In GR, functions...
599  // ...to compute metric
600  void CellMetric(const int k, const int j, const int il, const int iu,
601  AthenaArray<Real> &g, AthenaArray<Real> &gi) final;
602  void Face1Metric(const int k, const int j, const int il, const int iu,
603  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
604  void Face2Metric(const int k, const int j, const int il, const int iu,
605  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
606  void Face3Metric(const int k, const int j, const int il, const int iu,
607  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
608 
609  // ...to transform primitives to locally flat space
610  void PrimToLocal1(const int k, const int j, const int il, const int iu,
611  const AthenaArray<Real> &b1_vals, AthenaArray<Real> &prim_left,
612  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
613  void PrimToLocal2(const int k, const int j, const int il, const int iu,
614  const AthenaArray<Real> &b2_vals, AthenaArray<Real> &prim_left,
615  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
616  void PrimToLocal3(const int k, const int j, const int il, const int iu,
617  const AthenaArray<Real> &b3_vals, AthenaArray<Real> &prim_left,
618  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
619 
620  // ...to transform fluxes in locally flat space to global frame
621  void FluxToGlobal1(
622  const int k, const int j, const int il, const int iu,
623  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
624  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
625  void FluxToGlobal2(
626  const int k, const int j, const int il, const int iu,
627  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
628  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
629  void FluxToGlobal3(
630  const int k, const int j, const int il, const int iu,
631  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
632  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
633 
634  // for raising (lowering) covariant (contravariant) components of a vector
635  void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i,
636  Real *pa0, Real *pa1, Real *pa2, Real *pa3) final;
637  void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i,
638  Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final;
639 };
640 
641 //----------------------------------------------------------------------------------------
643 // \brief derived class for arbitrary (stationary) user-defined coordinates in GR.
644 // Nearly every function in the abstract base class need to be overridden.
645 
646 class GRUser : public Coordinates {
647  friend class HydroSourceTerms;
648 
649  public:
650  GRUser(MeshBlock *pmb, ParameterInput *pin, bool flag);
651 
652  // functions...
653  // ...to compute length of edges
654  void Edge1Length(const int k, const int j, const int il, const int iu,
655  AthenaArray<Real> &len) final;
656  void Edge2Length(const int k, const int j, const int il, const int iu,
657  AthenaArray<Real> &len) final;
658  void Edge3Length(const int k, const int j, const int il, const int iu,
659  AthenaArray<Real> &len) final;
660  Real GetEdge1Length(const int k, const int j, const int i) final;
661  Real GetEdge2Length(const int k, const int j, const int i) final;
662  Real GetEdge3Length(const int k, const int j, const int i) final;
663 
664  // ...to compute physical width at cell center
665  void CenterWidth1(const int k, const int j, const int il, const int iu,
666  AthenaArray<Real> &dx1) final;
667  void CenterWidth2(const int k, const int j, const int il, const int iu,
668  AthenaArray<Real> &dx2) final;
669  void CenterWidth3(const int k, const int j, const int il, const int iu,
670  AthenaArray<Real> &dx3) final;
671 
672  // ...to compute area of faces
673  void Face1Area(const int k, const int j, const int il, const int iu,
674  AthenaArray<Real> &area) final;
675  void Face2Area(const int k, const int j, const int il, const int iu,
676  AthenaArray<Real> &area) final;
677  void Face3Area(const int k, const int j, const int il, const int iu,
678  AthenaArray<Real> &area) final;
679  Real GetFace1Area(const int k, const int j, const int i) final;
680  Real GetFace2Area(const int k, const int j, const int i) final;
681  Real GetFace3Area(const int k, const int j, const int i) final;
682 
683  // ...to compute volumes of cells
684  void CellVolume(const int k, const int j, const int il, const int iu,
685  AthenaArray<Real> &vol) final;
686  Real GetCellVolume(const int k, const int j, const int i) final;
687 
688  // ...to compute geometrical source terms
689  void AddCoordTermsDivergence(const Real dt, const AthenaArray<Real> *flux,
690  const AthenaArray<Real> &prim, const AthenaArray<Real> &bcc,
691  AthenaArray<Real> &u) final;
692 
693  // In GR, functions...
694  // ...to compute metric
695  void CellMetric(const int k, const int j, const int il, const int iu,
696  AthenaArray<Real> &g, AthenaArray<Real> &gi) final;
697  void Face1Metric(const int k, const int j, const int il, const int iu,
698  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
699  void Face2Metric(const int k, const int j, const int il, const int iu,
700  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
701  void Face3Metric(const int k, const int j, const int il, const int iu,
702  AthenaArray<Real> &g, AthenaArray<Real> &g_inv) final;
703 
704  // ...to transform primitives to locally flat space
705  void PrimToLocal1(const int k, const int j, const int il, const int iu,
706  const AthenaArray<Real> &b1_vals, AthenaArray<Real> &prim_left,
707  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
708  void PrimToLocal2(const int k, const int j, const int il, const int iu,
709  const AthenaArray<Real> &b2_vals, AthenaArray<Real> &prim_left,
710  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
711  void PrimToLocal3(const int k, const int j, const int il, const int iu,
712  const AthenaArray<Real> &b3_vals, AthenaArray<Real> &prim_left,
713  AthenaArray<Real> &prim_right, AthenaArray<Real> &bx) final;
714 
715  // ...to transform fluxes in locally flat space to global frame
716  void FluxToGlobal1(
717  const int k, const int j, const int il, const int iu,
718  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
719  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
720  void FluxToGlobal2(
721  const int k, const int j, const int il, const int iu,
722  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
723  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
724  void FluxToGlobal3(
725  const int k, const int j, const int il, const int iu,
726  const AthenaArray<Real> &cons, const AthenaArray<Real> &bbx,
727  AthenaArray<Real> &flux, AthenaArray<Real> &ey, AthenaArray<Real> &ez) final;
728 
729  // ...for raising (lowering) covariant (contravariant) components of a vector
730  void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i,
731  Real *pa0, Real *pa1, Real *pa2, Real *pa3) final;
732  void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i,
733  Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final;
734 };
735 
736 #endif // COORDINATES_COORDINATES_HPP_
Real
double Real
Definition: athena.hpp:29
Cartesian::Cartesian
Cartesian(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: cartesian.cpp:24
Coordinates::trans_face3_i1_
AthenaArray< Real > trans_face3_i1_
Definition: coordinates.hpp:268
Coordinates::coord_area2_k1_
AthenaArray< Real > coord_area2_k1_
Definition: coordinates.hpp:207
Coordinates::x1s2
AthenaArray< Real > x1s2
Definition: coordinates.hpp:43
Coordinates::coord_len2_i1_
AthenaArray< Real > coord_len2_i1_
Definition: coordinates.hpp:218
Coordinates::coord_area2vc_i_
AthenaArray< Real > coord_area2vc_i_
Definition: coordinates.hpp:213
Coordinates::metric_face1_i2_
AthenaArray< Real > metric_face1_i2_
Definition: coordinates.hpp:249
Coordinates::coord_width3_j2_
AthenaArray< Real > coord_width3_j2_
Definition: coordinates.hpp:230
Coordinates::trans_face3_ji4_
AthenaArray< Real > trans_face3_ji4_
Definition: coordinates.hpp:271
Coordinates::metric_face1_j1_
AthenaArray< Real > metric_face1_j1_
Definition: coordinates.hpp:250
Coordinates::GetSpin
Real GetSpin() const
Definition: coordinates.hpp:136
Coordinates::LowerVectorCell
virtual void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i, Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3)
Definition: coordinates.hpp:183
Coordinates::LaplacianX1
virtual void LaplacianX1(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int n, const int k, const int j, const int il, const int iu)
Definition: coordinates.cpp:577
Coordinates::trans_face2_i2_
AthenaArray< Real > trans_face2_i2_
Definition: coordinates.hpp:263
Coordinates::trans_face3_ji3_
AthenaArray< Real > trans_face3_ji3_
Definition: coordinates.hpp:270
Coordinates::coord_area1_j2_
AthenaArray< Real > coord_area1_j2_
Definition: coordinates.hpp:202
Coordinates::dx2v
AthenaArray< Real > dx2v
Definition: coordinates.hpp:42
Coordinates::x2f
AthenaArray< Real > x2f
Definition: coordinates.hpp:41
Coordinates::phy_src2_i_
AthenaArray< Real > phy_src2_i_
Definition: coordinates.hpp:243
Coordinates::coord_len3_j2_
AthenaArray< Real > coord_len3_j2_
Definition: coordinates.hpp:222
Coordinates::trans_face3_ji1_
AthenaArray< Real > trans_face3_ji1_
Definition: coordinates.hpp:270
Coordinates::h31f
AthenaArray< Real > h31f
Definition: coordinates.hpp:45
Coordinates::dx1f
AthenaArray< Real > dx1f
Definition: coordinates.hpp:41
Coordinates::metric_face3_j2_
AthenaArray< Real > metric_face3_j2_
Definition: coordinates.hpp:256
Coordinates::x3s2
AthenaArray< Real > x3s2
Definition: coordinates.hpp:43
Coordinates::trans_face1_kji_
AthenaArray< Real > trans_face1_kji_
Definition: coordinates.hpp:262
Coordinates::metric_face2_i1_
AthenaArray< Real > metric_face2_i1_
Definition: coordinates.hpp:252
Coordinates::Face1Area
virtual void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:436
Coordinates::dx2f
AthenaArray< Real > dx2f
Definition: coordinates.hpp:41
Coordinates::trans_face3_kji_
AthenaArray< Real > trans_face3_kji_
Definition: coordinates.hpp:272
Coordinates::GetEdge3Length
virtual Real GetEdge3Length(const int k, const int j, const int i)
Definition: coordinates.cpp:371
Coordinates::Edge2Length
virtual void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:339
Coordinates::trans_face1_ji6_
AthenaArray< Real > trans_face1_ji6_
Definition: coordinates.hpp:261
Coordinates::coord_area2_j1_
AthenaArray< Real > coord_area2_j1_
Definition: coordinates.hpp:206
Coordinates::dx3v
AthenaArray< Real > dx3v
Definition: coordinates.hpp:42
Coordinates::coord_area2_i_
AthenaArray< Real > coord_area2_i_
Definition: coordinates.hpp:205
Coordinates::trans_face3_j1_
AthenaArray< Real > trans_face3_j1_
Definition: coordinates.hpp:269
Coordinates::coord_area1vc_j_
AthenaArray< Real > coord_area1vc_j_
Definition: coordinates.hpp:212
Coordinates::coord_len1_i2_
AthenaArray< Real > coord_len1_i2_
Definition: coordinates.hpp:215
Coordinates::trans_face2_j1_
AthenaArray< Real > trans_face2_j1_
Definition: coordinates.hpp:264
Coordinates::dx3f
AthenaArray< Real > dx3f
Definition: coordinates.hpp:41
Coordinates::CenterWidth2
virtual void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2)
Definition: coordinates.cpp:415
Coordinates::metric_face2_j1_
AthenaArray< Real > metric_face2_j1_
Definition: coordinates.hpp:253
Coordinates::coord_width2_i1_
AthenaArray< Real > coord_width2_i1_
Definition: coordinates.hpp:227
Coordinates::FluxToGlobal1
virtual void FluxToGlobal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez)
Definition: coordinates.hpp:167
Coordinates::coord_area1_kji_
AthenaArray< Real > coord_area1_kji_
Definition: coordinates.hpp:204
Coordinates::trans_face1_i1_
AthenaArray< Real > trans_face1_i1_
Definition: coordinates.hpp:258
Coordinates::metric_face3_i1_
AthenaArray< Real > metric_face3_i1_
Definition: coordinates.hpp:255
Coordinates::CellMetric
virtual void CellMetric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &gi)
Definition: coordinates.hpp:143
Coordinates::trans_face2_ji6_
AthenaArray< Real > trans_face2_ji6_
Definition: coordinates.hpp:266
Coordinates::coord_src2_j_
AthenaArray< Real > coord_src2_j_
Definition: coordinates.hpp:239
Coordinates::coord_width2_kji_
AthenaArray< Real > coord_width2_kji_
Definition: coordinates.hpp:229
Coordinates::coord_width2_j1_
AthenaArray< Real > coord_width2_j1_
Definition: coordinates.hpp:228
Coordinates::coord_src_j1_
AthenaArray< Real > coord_src_j1_
Definition: coordinates.hpp:234
Coordinates::trans_face2_ji5_
AthenaArray< Real > trans_face2_ji5_
Definition: coordinates.hpp:266
Coordinates::coord_len2_kji_
AthenaArray< Real > coord_len2_kji_
Definition: coordinates.hpp:220
Coordinates::Coordinates
Coordinates(MeshBlock *pmb, ParameterInput *pin, bool flag=false)
Definition: coordinates.cpp:26
Coordinates::trans_face1_ji2_
AthenaArray< Real > trans_face1_ji2_
Definition: coordinates.hpp:260
Coordinates::LaplacianX3All
virtual void LaplacianX3All(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int nl, const int nu, const int kl, const int ku, const int jl, const int ju, const int il, const int iu)
Definition: coordinates.cpp:701
Coordinates::GetEdge1Length
virtual Real GetEdge1Length(const int k, const int j, const int i)
Definition: coordinates.cpp:363
Coordinates::metric_face2_i2_
AthenaArray< Real > metric_face2_i2_
Definition: coordinates.hpp:252
Coordinates::~Coordinates
virtual ~Coordinates()=default
Coordinates::gi_
AthenaArray< Real > gi_
Definition: coordinates.hpp:273
Coordinates::PrimToLocal1
virtual void PrimToLocal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b1_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx)
Definition: coordinates.hpp:153
Coordinates::h32v
AthenaArray< Real > h32v
Definition: coordinates.hpp:46
Coordinates::coord_len3_i1_
AthenaArray< Real > coord_len3_i1_
Definition: coordinates.hpp:221
Coordinates::coord_vol_j2_
AthenaArray< Real > coord_vol_j2_
Definition: coordinates.hpp:198
Coordinates::coord_area1_j1_
AthenaArray< Real > coord_area1_j1_
Definition: coordinates.hpp:202
Coordinates::coord_src3_j_
AthenaArray< Real > coord_src3_j_
Definition: coordinates.hpp:240
Coordinates::coord_vol_i2_
AthenaArray< Real > coord_vol_i2_
Definition: coordinates.hpp:197
Coordinates::trans_face1_j1_
AthenaArray< Real > trans_face1_j1_
Definition: coordinates.hpp:259
Coordinates::h2v
AthenaArray< Real > h2v
Definition: coordinates.hpp:46
Coordinates::coord_len2_j1_
AthenaArray< Real > coord_len2_j1_
Definition: coordinates.hpp:219
Coordinates::VolCenter3Length
virtual void VolCenter3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:394
Coordinates::coord_area3_i_
AthenaArray< Real > coord_area3_i_
Definition: coordinates.hpp:209
Coordinates::GetEdge2Length
virtual Real GetEdge2Length(const int k, const int j, const int i)
Definition: coordinates.cpp:367
Coordinates::trans_face1_i2_
AthenaArray< Real > trans_face1_i2_
Definition: coordinates.hpp:258
Coordinates::coord_len3_j1_
AthenaArray< Real > coord_len3_j1_
Definition: coordinates.hpp:222
Coordinates::trans_face1_ji3_
AthenaArray< Real > trans_face1_ji3_
Definition: coordinates.hpp:260
Coordinates::metric_face2_j2_
AthenaArray< Real > metric_face2_j2_
Definition: coordinates.hpp:253
Coordinates::Edge1Length
virtual void Edge1Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:328
Coordinates::h2f
AthenaArray< Real > h2f
Definition: coordinates.hpp:45
Coordinates::coord_len2_j2_
AthenaArray< Real > coord_len2_j2_
Definition: coordinates.hpp:219
Coordinates::coord_area1vc_i_
AthenaArray< Real > coord_area1vc_i_
Definition: coordinates.hpp:212
Coordinates::coord_vol_j1_
AthenaArray< Real > coord_vol_j1_
Definition: coordinates.hpp:198
Coordinates::GetMass
Real GetMass() const
Definition: coordinates.hpp:135
Coordinates::x1f
AthenaArray< Real > x1f
Definition: coordinates.hpp:41
Coordinates::coord_src_j2_
AthenaArray< Real > coord_src_j2_
Definition: coordinates.hpp:234
Coordinates::coord_area3_i2_
AthenaArray< Real > coord_area3_i2_
Definition: coordinates.hpp:209
Coordinates::coord_len1_j2_
AthenaArray< Real > coord_len1_j2_
Definition: coordinates.hpp:216
Coordinates::metric_face1_i1_
AthenaArray< Real > metric_face1_i1_
Definition: coordinates.hpp:249
Coordinates::coord_len3_k1_
AthenaArray< Real > coord_len3_k1_
Definition: coordinates.hpp:223
Coordinates::coord_vol_i_
AthenaArray< Real > coord_vol_i_
Definition: coordinates.hpp:197
Coordinates::coord_len1_i1_
AthenaArray< Real > coord_len1_i1_
Definition: coordinates.hpp:215
Coordinates::coord_area2_i2_
AthenaArray< Real > coord_area2_i2_
Definition: coordinates.hpp:205
Coordinates::dh31vd1
AthenaArray< Real > dh31vd1
Definition: coordinates.hpp:46
Coordinates::x3f
AthenaArray< Real > x3f
Definition: coordinates.hpp:41
Coordinates::coord_area3_i1_
AthenaArray< Real > coord_area3_i1_
Definition: coordinates.hpp:209
Coordinates::LaplacianX2All
virtual void LaplacianX2All(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int nl, const int nu, const int kl, const int ku, const int jl, const int ju, const int il, const int iu)
Definition: coordinates.cpp:650
Coordinates::g_
AthenaArray< Real > g_
Definition: coordinates.hpp:273
Coordinates::coord_len1_j1_
AthenaArray< Real > coord_len1_j1_
Definition: coordinates.hpp:216
Coordinates::coord_area3_kji_
AthenaArray< Real > coord_area3_kji_
Definition: coordinates.hpp:211
Coordinates::pmy_block
MeshBlock * pmy_block
Definition: coordinates.hpp:40
Coordinates::dh32fd2
AthenaArray< Real > dh32fd2
Definition: coordinates.hpp:45
Coordinates::x1v
AthenaArray< Real > x1v
Definition: coordinates.hpp:42
Coordinates::VolCenter1Length
virtual void VolCenter1Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:378
Coordinates::coord_area2_j2_
AthenaArray< Real > coord_area2_j2_
Definition: coordinates.hpp:206
Coordinates::VolCenterFace1Area
virtual void VolCenterFace1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:488
Coordinates::Face1Metric
virtual void Face1Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv)
Definition: coordinates.hpp:145
Coordinates::GetFace1Area
virtual Real GetFace1Area(const int k, const int j, const int i)
Definition: coordinates.cpp:472
Coordinates::coord_area2_i1_
AthenaArray< Real > coord_area2_i1_
Definition: coordinates.hpp:205
Coordinates::coord_len3_kji_
AthenaArray< Real > coord_len3_kji_
Definition: coordinates.hpp:224
Coordinates::x2v
AthenaArray< Real > x2v
Definition: coordinates.hpp:42
Coordinates::Face2Metric
virtual void Face2Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv)
Definition: coordinates.hpp:147
Coordinates::VolCenterFace2Area
virtual void VolCenterFace2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:498
Coordinates::phy_src1_i_
AthenaArray< Real > phy_src1_i_
Definition: coordinates.hpp:243
Coordinates::dx1v
AthenaArray< Real > dx1v
Definition: coordinates.hpp:42
Coordinates::coord_area2vc_j_
AthenaArray< Real > coord_area2vc_j_
Definition: coordinates.hpp:213
Coordinates::coord_src1_j_
AthenaArray< Real > coord_src1_j_
Definition: coordinates.hpp:237
Coordinates::metric_cell_i2_
AthenaArray< Real > metric_cell_i2_
Definition: coordinates.hpp:246
Coordinates::trans_face3_ji2_
AthenaArray< Real > trans_face3_ji2_
Definition: coordinates.hpp:270
Coordinates::coord_area3_j1_
AthenaArray< Real > coord_area3_j1_
Definition: coordinates.hpp:210
Coordinates::trans_face1_ji5_
AthenaArray< Real > trans_face1_ji5_
Definition: coordinates.hpp:261
Coordinates::GetFace3Area
virtual Real GetFace3Area(const int k, const int j, const int i)
Definition: coordinates.cpp:480
Coordinates::IsPole
bool IsPole(int j)
Definition: coordinates.cpp:749
Coordinates::coord_width1_i1_
AthenaArray< Real > coord_width1_i1_
Definition: coordinates.hpp:225
Coordinates::x1s3
AthenaArray< Real > x1s3
Definition: coordinates.hpp:43
Coordinates::trans_face2_i1_
AthenaArray< Real > trans_face2_i1_
Definition: coordinates.hpp:263
Coordinates::CenterWidth3
virtual void CenterWidth3(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx3)
Definition: coordinates.cpp:424
Coordinates::trans_face3_ji5_
AthenaArray< Real > trans_face3_ji5_
Definition: coordinates.hpp:271
Coordinates::x3v
AthenaArray< Real > x3v
Definition: coordinates.hpp:42
Coordinates::h32f
AthenaArray< Real > h32f
Definition: coordinates.hpp:45
Coordinates::coord_src_kji_
AthenaArray< Real > coord_src_kji_
Definition: coordinates.hpp:235
Coordinates::coord_width1_kji_
AthenaArray< Real > coord_width1_kji_
Definition: coordinates.hpp:226
Coordinates::trans_face1_ji1_
AthenaArray< Real > trans_face1_ji1_
Definition: coordinates.hpp:260
Coordinates::coord_area1_j_
AthenaArray< Real > coord_area1_j_
Definition: coordinates.hpp:202
Coordinates::coord_vol_k1_
AthenaArray< Real > coord_vol_k1_
Definition: coordinates.hpp:199
Coordinates::coord_vol_kji_
AthenaArray< Real > coord_vol_kji_
Definition: coordinates.hpp:200
Coordinates::trans_face1_ji7_
AthenaArray< Real > trans_face1_ji7_
Definition: coordinates.hpp:261
Coordinates::trans_face1_ji4_
AthenaArray< Real > trans_face1_ji4_
Definition: coordinates.hpp:261
Coordinates::VolCenterFace3Area
virtual void VolCenterFace3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:508
Coordinates::FluxToGlobal2
virtual void FluxToGlobal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez)
Definition: coordinates.hpp:171
Coordinates::x3s1
AthenaArray< Real > x3s1
Definition: coordinates.hpp:43
Coordinates::RaiseVectorCell
virtual void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3)
Definition: coordinates.hpp:181
Coordinates::metric_cell_j2_
AthenaArray< Real > metric_cell_j2_
Definition: coordinates.hpp:247
Coordinates::coord_vol_i1_
AthenaArray< Real > coord_vol_i1_
Definition: coordinates.hpp:197
Coordinates::dh31fd1
AthenaArray< Real > dh31fd1
Definition: coordinates.hpp:45
Coordinates::coord_width3_kji_
AthenaArray< Real > coord_width3_kji_
Definition: coordinates.hpp:233
Coordinates::coord_vol_j_
AthenaArray< Real > coord_vol_j_
Definition: coordinates.hpp:198
Coordinates::coord_area1_i_
AthenaArray< Real > coord_area1_i_
Definition: coordinates.hpp:201
Coordinates::x2s3
AthenaArray< Real > x2s3
Definition: coordinates.hpp:43
Coordinates::coord_area1_k1_
AthenaArray< Real > coord_area1_k1_
Definition: coordinates.hpp:203
Coordinates::dh2vd1
AthenaArray< Real > dh2vd1
Definition: coordinates.hpp:46
Coordinates::coord_area2_j_
AthenaArray< Real > coord_area2_j_
Definition: coordinates.hpp:206
Coordinates::trans_face2_ji3_
AthenaArray< Real > trans_face2_ji3_
Definition: coordinates.hpp:265
Coordinates::metric_face3_j1_
AthenaArray< Real > metric_face3_j1_
Definition: coordinates.hpp:256
Coordinates::coord_width3_j3_
AthenaArray< Real > coord_width3_j3_
Definition: coordinates.hpp:230
Coordinates::metric_face2_kji_
AthenaArray< Real > metric_face2_kji_
Definition: coordinates.hpp:254
Coordinates::metric_face3_kji_
AthenaArray< Real > metric_face3_kji_
Definition: coordinates.hpp:257
Coordinates::coord_area2_kji_
AthenaArray< Real > coord_area2_kji_
Definition: coordinates.hpp:208
Coordinates::Face2Area
virtual void Face2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:447
Coordinates::trans_face2_ji1_
AthenaArray< Real > trans_face2_ji1_
Definition: coordinates.hpp:265
Coordinates::PrimToLocal2
virtual void PrimToLocal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b2_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx)
Definition: coordinates.hpp:157
Coordinates::LaplacianX3
virtual void LaplacianX3(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int n, const int k, const int j, const int il, const int iu)
Definition: coordinates.cpp:683
Coordinates::h31v
AthenaArray< Real > h31v
Definition: coordinates.hpp:46
Coordinates::coord_width3_j1_
AthenaArray< Real > coord_width3_j1_
Definition: coordinates.hpp:230
Coordinates::Edge3Length
virtual void Edge3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:350
Coordinates::trans_face3_ji6_
AthenaArray< Real > trans_face3_ji6_
Definition: coordinates.hpp:271
Coordinates::coord_src1_i_
AthenaArray< Real > coord_src1_i_
Definition: coordinates.hpp:236
Coordinates::trans_face2_ji4_
AthenaArray< Real > trans_face2_ji4_
Definition: coordinates.hpp:266
Coordinates::FluxToGlobal3
virtual void FluxToGlobal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez)
Definition: coordinates.hpp:175
Coordinates::CellVolume
virtual void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol)
Definition: coordinates.cpp:521
Coordinates::coord_area1_i1_
AthenaArray< Real > coord_area1_i1_
Definition: coordinates.hpp:201
Coordinates::LaplacianX2
virtual void LaplacianX2(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int n, const int k, const int j, const int il, const int iu)
Definition: coordinates.cpp:633
Coordinates::Laplacian
virtual void Laplacian(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int il, const int iu, const int jl, const int ju, const int kl, const int ku, const int nl, const int nu)
Definition: coordinates.cpp:543
Coordinates::dh2fd1
AthenaArray< Real > dh2fd1
Definition: coordinates.hpp:45
Coordinates::dh32vd2
AthenaArray< Real > dh32vd2
Definition: coordinates.hpp:46
Coordinates::trans_face3_i2_
AthenaArray< Real > trans_face3_i2_
Definition: coordinates.hpp:268
Coordinates::metric_face1_j2_
AthenaArray< Real > metric_face1_j2_
Definition: coordinates.hpp:250
Coordinates::coord_width3_k1_
AthenaArray< Real > coord_width3_k1_
Definition: coordinates.hpp:231
Coordinates::coord_area3_j2_
AthenaArray< Real > coord_area3_j2_
Definition: coordinates.hpp:210
Coordinates::CenterWidth1
virtual void CenterWidth1(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx1)
Definition: coordinates.cpp:406
Coordinates::metric_face1_kji_
AthenaArray< Real > metric_face1_kji_
Definition: coordinates.hpp:251
Coordinates::metric_face3_i2_
AthenaArray< Real > metric_face3_i2_
Definition: coordinates.hpp:255
Coordinates::coord_len1_kji_
AthenaArray< Real > coord_len1_kji_
Definition: coordinates.hpp:217
Coordinates::GetFace2Area
virtual Real GetFace2Area(const int k, const int j, const int i)
Definition: coordinates.cpp:476
Coordinates::trans_face2_ji2_
AthenaArray< Real > trans_face2_ji2_
Definition: coordinates.hpp:265
Coordinates::Metric
void Metric(Real x1, Real x2, Real x3, ParameterInput *pin, AthenaArray< Real > &g, AthenaArray< Real > &g_inv, AthenaArray< Real > &dg_dx1, AthenaArray< Real > &dg_dx2, AthenaArray< Real > &dg_dx3)
Definition: coordinates.cpp:772
Coordinates::coord_src2_i_
AthenaArray< Real > coord_src2_i_
Definition: coordinates.hpp:238
Coordinates::Face3Metric
virtual void Face3Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv)
Definition: coordinates.hpp:149
Coordinates::coord_width3_ji1_
AthenaArray< Real > coord_width3_ji1_
Definition: coordinates.hpp:232
Coordinates::VolCenter2Length
virtual void VolCenter2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len)
Definition: coordinates.cpp:386
Coordinates::metric_cell_kji_
AthenaArray< Real > metric_cell_kji_
Definition: coordinates.hpp:248
Coordinates::trans_face2_kji_
AthenaArray< Real > trans_face2_kji_
Definition: coordinates.hpp:267
Coordinates::x2s1
AthenaArray< Real > x2s1
Definition: coordinates.hpp:43
Coordinates::GetCellVolume
virtual Real GetCellVolume(const int k, const int j, const int i)
Definition: coordinates.cpp:535
Coordinates::PrimToLocal3
virtual void PrimToLocal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b3_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx)
Definition: coordinates.hpp:161
Coordinates::LaplacianX1All
virtual void LaplacianX1All(const AthenaArray< Real > &s, AthenaArray< Real > &delta_s, const int nl, const int nu, const int kl, const int ku, const int jl, const int ju, const int il, const int iu)
Definition: coordinates.cpp:597
Coordinates::metric_cell_i1_
AthenaArray< Real > metric_cell_i1_
Definition: coordinates.hpp:246
Coordinates::metric_cell_j1_
AthenaArray< Real > metric_cell_j1_
Definition: coordinates.hpp:247
Coordinates::Face3Area
virtual void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area)
Definition: coordinates.cpp:458
Coordinates::AddCoordTermsDivergence
virtual void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u)
Definition: coordinates.cpp:736
Coordinates::coord_area3vc_i_
AthenaArray< Real > coord_area3vc_i_
Definition: coordinates.hpp:214
Cylindrical::AddCoordTermsDivergence
void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u) final
Definition: cylindrical.cpp:286
Cylindrical::Cylindrical
Cylindrical(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: cylindrical.cpp:28
Cylindrical::GetCellVolume
Real GetCellVolume(const int k, const int j, const int i) final
Definition: cylindrical.cpp:279
Cylindrical::CenterWidth2
void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2) final
Definition: cylindrical.cpp:195
Cylindrical::GetFace3Area
Real GetFace3Area(const int k, const int j, const int i) final
Definition: cylindrical.cpp:235
Cylindrical::GetFace1Area
Real GetFace1Area(const int k, const int j, const int i) final
Definition: cylindrical.cpp:230
Cylindrical::Edge2Length
void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: cylindrical.cpp:162
Cylindrical::VolCenter2Length
void VolCenter2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: cylindrical.cpp:182
Cylindrical::VolCenterFace3Area
void VolCenterFace3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: cylindrical.cpp:255
Cylindrical::Face1Area
void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: cylindrical.cpp:207
Cylindrical::Face3Area
void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: cylindrical.cpp:217
Cylindrical::CellVolume
void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol) final
Definition: cylindrical.cpp:266
Cylindrical::VolCenterFace1Area
void VolCenterFace1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: cylindrical.cpp:245
Cylindrical::GetEdge2Length
Real GetEdge2Length(const int k, const int j, const int i) final
Definition: cylindrical.cpp:174
GRUser::Face1Metric
void Face1Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: gr_user.cpp:737
GRUser::GetFace2Area
Real GetFace2Area(const int k, const int j, const int i) final
Definition: gr_user.cpp:536
GRUser::Edge1Length
void Edge1Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: gr_user.cpp:405
GRUser::CenterWidth3
void CenterWidth3(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx3) final
Definition: gr_user.cpp:476
GRUser::FluxToGlobal1
void FluxToGlobal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: gr_user.cpp:1304
GRUser::Edge3Length
void Edge3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: gr_user.cpp:425
GRUser::GetEdge1Length
Real GetEdge1Length(const int k, const int j, const int i) final
Definition: gr_user.cpp:438
GRUser::GRUser
GRUser(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: gr_user.cpp:42
GRUser::Face2Area
void Face2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: gr_user.cpp:504
GRUser::PrimToLocal3
void PrimToLocal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b3_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: gr_user.cpp:1136
GRUser::CellVolume
void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol) final
Definition: gr_user.cpp:554
GRUser::FluxToGlobal3
void FluxToGlobal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: gr_user.cpp:1528
GRUser::Face1Area
void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: gr_user.cpp:494
GRUser::GetFace1Area
Real GetFace1Area(const int k, const int j, const int i) final
Definition: gr_user.cpp:531
GRUser::GetEdge3Length
Real GetEdge3Length(const int k, const int j, const int i) final
Definition: gr_user.cpp:448
GRUser::CellMetric
void CellMetric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &gi) final
Definition: gr_user.cpp:725
GRUser::CenterWidth1
void CenterWidth1(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx1) final
Definition: gr_user.cpp:456
GRUser::Edge2Length
void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: gr_user.cpp:415
GRUser::Face2Metric
void Face2Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: gr_user.cpp:749
GRUser::GetFace3Area
Real GetFace3Area(const int k, const int j, const int i) final
Definition: gr_user.cpp:541
GRUser::PrimToLocal1
void PrimToLocal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b1_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: gr_user.cpp:794
GRUser::GetEdge2Length
Real GetEdge2Length(const int k, const int j, const int i) final
Definition: gr_user.cpp:443
GRUser::CenterWidth2
void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2) final
Definition: gr_user.cpp:466
GRUser::LowerVectorCell
void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i, Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final
Definition: gr_user.cpp:1666
GRUser::PrimToLocal2
void PrimToLocal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b2_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: gr_user.cpp:965
GRUser::FluxToGlobal2
void FluxToGlobal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: gr_user.cpp:1416
GRUser::GetCellVolume
Real GetCellVolume(const int k, const int j, const int i) final
Definition: gr_user.cpp:571
GRUser::Face3Area
void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: gr_user.cpp:514
GRUser::RaiseVectorCell
void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3) final
Definition: gr_user.cpp:1630
GRUser::Face3Metric
void Face3Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: gr_user.cpp:761
GRUser::AddCoordTermsDivergence
void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u) final
Definition: gr_user.cpp:586
KerrSchild::CenterWidth1
void CenterWidth1(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx1) final
Definition: kerr-schild.cpp:536
KerrSchild::Face1Area
void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: kerr-schild.cpp:576
KerrSchild::GetFace3Area
Real GetFace3Area(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:633
KerrSchild::PrimToLocal2
void PrimToLocal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b2_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: kerr-schild.cpp:1171
KerrSchild::FluxToGlobal2
void FluxToGlobal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: kerr-schild.cpp:1554
KerrSchild::Face2Metric
void Face2Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: kerr-schild.cpp:922
KerrSchild::LowerVectorCell
void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i, Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final
Definition: kerr-schild.cpp:1791
KerrSchild::CenterWidth2
void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2) final
Definition: kerr-schild.cpp:547
KerrSchild::GetEdge1Length
Real GetEdge1Length(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:515
KerrSchild::GetEdge3Length
Real GetEdge3Length(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:528
KerrSchild::Face3Area
void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: kerr-schild.cpp:599
KerrSchild::RaiseVectorCell
void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3) final
Definition: kerr-schild.cpp:1744
KerrSchild::AddCoordTermsDivergence
void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u) final
Definition: kerr-schild.cpp:689
KerrSchild::PrimToLocal3
void PrimToLocal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b3_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: kerr-schild.cpp:1313
KerrSchild::FluxToGlobal3
void FluxToGlobal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: kerr-schild.cpp:1654
KerrSchild::CellMetric
void CellMetric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &gi) final
Definition: kerr-schild.cpp:828
KerrSchild::GetCellVolume
Real GetCellVolume(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:674
KerrSchild::Edge3Length
void Edge3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: kerr-schild.cpp:501
KerrSchild::PrimToLocal1
void PrimToLocal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b1_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: kerr-schild.cpp:1028
KerrSchild::Face3Metric
void Face3Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: kerr-schild.cpp:969
KerrSchild::Edge2Length
void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: kerr-schild.cpp:490
KerrSchild::KerrSchild
KerrSchild(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: kerr-schild.cpp:48
KerrSchild::GetFace1Area
Real GetFace1Area(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:619
KerrSchild::Edge1Length
void Edge1Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: kerr-schild.cpp:478
KerrSchild::Face1Metric
void Face1Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: kerr-schild.cpp:875
KerrSchild::Face2Area
void Face2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: kerr-schild.cpp:587
KerrSchild::FluxToGlobal1
void FluxToGlobal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: kerr-schild.cpp:1452
KerrSchild::CellVolume
void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol) final
Definition: kerr-schild.cpp:653
KerrSchild::CenterWidth3
void CenterWidth3(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx3) final
Definition: kerr-schild.cpp:557
KerrSchild::GetFace2Area
Real GetFace2Area(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:626
KerrSchild::GetEdge2Length
Real GetEdge2Length(const int k, const int j, const int i) final
Definition: kerr-schild.cpp:522
Mesh
Definition: mesh.hpp:203
Minkowski::Face1Metric
void Face1Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: minkowski.cpp:127
Minkowski::Face2Metric
void Face2Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: minkowski.cpp:143
Minkowski::PrimToLocal2
void PrimToLocal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b2_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: minkowski.cpp:203
Minkowski::CellMetric
void CellMetric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &gi) final
Definition: minkowski.cpp:102
Minkowski::Minkowski
Minkowski(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: minkowski.cpp:33
Minkowski::PrimToLocal1
void PrimToLocal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b1_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: minkowski.cpp:190
Minkowski::FluxToGlobal1
void FluxToGlobal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: minkowski.cpp:243
Minkowski::LowerVectorCell
void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i, Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final
Definition: minkowski.cpp:307
Minkowski::PrimToLocal3
void PrimToLocal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b3_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: minkowski.cpp:216
Minkowski::FluxToGlobal3
void FluxToGlobal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: minkowski.cpp:269
Minkowski::Face3Metric
void Face3Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: minkowski.cpp:159
Minkowski::FluxToGlobal2
void FluxToGlobal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: minkowski.cpp:256
Minkowski::RaiseVectorCell
void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3) final
Definition: minkowski.cpp:290
Schwarzschild::Face1Area
void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: schwarzschild.cpp:370
Schwarzschild::GetFace3Area
Real GetFace3Area(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:414
Schwarzschild::CenterWidth3
void CenterWidth3(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx3) final
Definition: schwarzschild.cpp:352
Schwarzschild::Face3Metric
void Face3Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: schwarzschild.cpp:709
Schwarzschild::AddCoordTermsDivergence
void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u) final
Definition: schwarzschild.cpp:461
Schwarzschild::GetFace2Area
Real GetFace2Area(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:409
Schwarzschild::Edge1Length
void Edge1Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: schwarzschild.cpp:281
Schwarzschild::CellMetric
void CellMetric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &gi) final
Definition: schwarzschild.cpp:574
Schwarzschild::CellVolume
void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol) final
Definition: schwarzschild.cpp:429
Schwarzschild::FluxToGlobal1
void FluxToGlobal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: schwarzschild.cpp:1168
Schwarzschild::GetFace1Area
Real GetFace1Area(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:404
Schwarzschild::Schwarzschild
Schwarzschild(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: schwarzschild.cpp:38
Schwarzschild::Face3Area
void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: schwarzschild.cpp:388
Schwarzschild::FluxToGlobal3
void FluxToGlobal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: schwarzschild.cpp:1334
Schwarzschild::CenterWidth1
void CenterWidth1(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx1) final
Definition: schwarzschild.cpp:332
Schwarzschild::GetEdge2Length
Real GetEdge2Length(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:319
Schwarzschild::PrimToLocal1
void PrimToLocal1(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b1_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: schwarzschild.cpp:766
Schwarzschild::PrimToLocal2
void PrimToLocal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b2_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: schwarzschild.cpp:901
Schwarzschild::LowerVectorCell
void LowerVectorCell(Real a0, Real a1, Real a2, Real a3, int k, int j, int i, Real *pa_0, Real *pa_1, Real *pa_2, Real *pa_3) final
Definition: schwarzschild.cpp:1437
Schwarzschild::Edge3Length
void Edge3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: schwarzschild.cpp:301
Schwarzschild::FluxToGlobal2
void FluxToGlobal2(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &cons, const AthenaArray< Real > &bbx, AthenaArray< Real > &flux, AthenaArray< Real > &ey, AthenaArray< Real > &ez) final
Definition: schwarzschild.cpp:1251
Schwarzschild::GetEdge1Length
Real GetEdge1Length(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:314
Schwarzschild::RaiseVectorCell
void RaiseVectorCell(Real a_0, Real a_1, Real a_2, Real a_3, int k, int j, int i, Real *pa0, Real *pa1, Real *pa2, Real *pa3) final
Definition: schwarzschild.cpp:1407
Schwarzschild::Face2Metric
void Face2Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: schwarzschild.cpp:664
Schwarzschild::Face1Metric
void Face1Metric(const int k, const int j, const int il, const int iu, AthenaArray< Real > &g, AthenaArray< Real > &g_inv) final
Definition: schwarzschild.cpp:619
Schwarzschild::Face2Area
void Face2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: schwarzschild.cpp:379
Schwarzschild::CenterWidth2
void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2) final
Definition: schwarzschild.cpp:342
Schwarzschild::GetCellVolume
Real GetCellVolume(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:447
Schwarzschild::Edge2Length
void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: schwarzschild.cpp:291
Schwarzschild::PrimToLocal3
void PrimToLocal3(const int k, const int j, const int il, const int iu, const AthenaArray< Real > &b3_vals, AthenaArray< Real > &prim_left, AthenaArray< Real > &prim_right, AthenaArray< Real > &bx) final
Definition: schwarzschild.cpp:1036
Schwarzschild::GetEdge3Length
Real GetEdge3Length(const int k, const int j, const int i) final
Definition: schwarzschild.cpp:324
SphericalPolar::Face1Area
void Face1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:341
SphericalPolar::GetCellVolume
Real GetCellVolume(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:434
SphericalPolar::CenterWidth3
void CenterWidth3(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx3) final
Definition: spherical_polar.cpp:329
SphericalPolar::VolCenter2Length
void VolCenter2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: spherical_polar.cpp:296
SphericalPolar::GetFace3Area
Real GetFace3Area(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:382
SphericalPolar::VolCenterFace3Area
void VolCenterFace3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:408
SphericalPolar::CellVolume
void CellVolume(const int k, const int j, const int il, const int iu, AthenaArray< Real > &vol) final
Definition: spherical_polar.cpp:421
SphericalPolar::SphericalPolar
SphericalPolar(MeshBlock *pmb, ParameterInput *pin, bool flag)
Definition: spherical_polar.cpp:33
SphericalPolar::AddCoordTermsDivergence
void AddCoordTermsDivergence(const Real dt, const AthenaArray< Real > *flux, const AthenaArray< Real > &prim, const AthenaArray< Real > &bcc, AthenaArray< Real > &u) final
Definition: spherical_polar.cpp:441
SphericalPolar::CenterWidth2
void CenterWidth2(const int k, const int j, const int il, const int iu, AthenaArray< Real > &dx2) final
Definition: spherical_polar.cpp:320
SphericalPolar::GetEdge3Length
Real GetEdge3Length(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:287
SphericalPolar::VolCenter3Length
void VolCenter3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: spherical_polar.cpp:307
SphericalPolar::GetFace1Area
Real GetFace1Area(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:374
SphericalPolar::Edge3Length
void Edge3Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: spherical_polar.cpp:270
SphericalPolar::Face3Area
void Face3Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:361
SphericalPolar::Face2Area
void Face2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:351
SphericalPolar::GetFace2Area
Real GetFace2Area(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:378
SphericalPolar::VolCenterFace1Area
void VolCenterFace1Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:390
SphericalPolar::Edge2Length
void Edge2Length(const int k, const int j, const int il, const int iu, AthenaArray< Real > &len) final
Definition: spherical_polar.cpp:258
SphericalPolar::GetEdge2Length
Real GetEdge2Length(const int k, const int j, const int i) final
Definition: spherical_polar.cpp:283
SphericalPolar::VolCenterFace2Area
void VolCenterFace2Area(const int k, const int j, const int il, const int iu, AthenaArray< Real > &area) final
Definition: spherical_polar.cpp:399
Weno3Coeff::a2
double a2[2]
Definition: interp_weno3.hpp:8
Weno3Coeff::a1
double a1[2]
Definition: interp_weno3.hpp:6
Weno5Coeff::a3
double a3[3]
Definition: interp_weno5.hpp:18