chombo-discharge
CD_AmrMesh.H
Go to the documentation of this file.
1 /* chombo-discharge
2  * Copyright © 2021 SINTEF Energy Research.
3  * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
4  */
5 
14 #ifndef CD_AmrMesh_H
15 #define CD_AmrMesh_H
16 
17 // Chombo includes
18 #include <DisjointBoxLayout.H>
19 #include <ProblemDomain.H>
20 
21 // Our includes
22 #include <CD_EBAMRData.H>
23 #include <CD_EBCoarAve.H>
24 #include <CD_ComputationalGeometry.H>
25 #include <CD_MultiFluidIndexSpace.H>
26 #include <CD_IrregAmrStencil.H>
27 #include <CD_LoadBalancing.H>
28 #include <CD_MFLevelGrid.H>
29 #include <CD_EBReflux.H>
30 #include <CD_DomainFluxIFFAB.H>
31 #include <CD_DomainFluxIFFABFactory.H>
32 #include <CD_ParticleContainer.H>
33 #include <CD_NonConservativeDivergenceStencil.H>
34 #include <CD_Realm.H>
35 #include <CD_CopyStrategy.H>
36 #include <CD_LoadBalancing.H>
37 #include <CD_NamespaceHeader.H>
38 
43 class AmrMesh
44 {
45 public:
49  AmrMesh();
50 
55  AmrMesh(const AmrMesh& a_other) = delete;
56 
61  AmrMesh(const AmrMesh&& a_other) = delete;
62 
66  virtual ~AmrMesh();
67 
72  AmrMesh&
73  operator=(const AmrMesh& a_other) = delete;
74 
79  AmrMesh&
80  operator=(const AmrMesh&& a_other) = delete;
81 
92  template <typename T>
93  void
94  copyData(EBAMRData<T>& a_dst,
95  const EBAMRData<T>& a_src,
96  const CopyStrategy& a_toRegion = CopyStrategy::Valid,
97  const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
98 
110  template <typename T>
111  void
112  copyData(EBAMRData<T>& a_dst,
113  const EBAMRData<T>& a_src,
114  const Interval& a_dstComps,
115  const Interval& a_srcComps,
116  const CopyStrategy& a_toRegion = CopyStrategy::Valid,
117  const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
118 
132  template <typename T>
133  void
134  copyData(LevelData<T>& a_dst,
135  const LevelData<T>& a_src,
136  const int a_level,
137  const std::string a_toRealm,
138  const std::string a_fromRealm,
139  const CopyStrategy& a_toRegion = CopyStrategy::Valid,
140  const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
141 
157  template <typename T>
158  void
159  copyData(LevelData<T>& a_dst,
160  const LevelData<T>& a_src,
161  const int a_level,
162  const std::string a_toRealm,
163  const std::string a_fromRealm,
164  const Interval& a_dstComps,
165  const Interval& a_srcComps,
166  const CopyStrategy& a_toRegion = CopyStrategy::Valid,
167  const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
168 
173  template <typename T>
174  void
175  deallocate(Vector<T*>& a_data) const;
176 
181  template <typename T>
182  void
183  deallocate(Vector<RefCountedPtr<T>>& a_data) const;
184 
189  template <typename T>
190  void
191  deallocate(EBAMRData<T>& a_data) const;
192 
199  template <typename T>
200  void
201  alias(Vector<T*>& a_alias, const Vector<RefCountedPtr<T>>& a_data) const;
202 
209  template <typename T, typename S>
210  void
211  alias(Vector<T*>& a_alias, const EBAMRData<S>& a_data) const;
212 
218  template <typename T>
219  void
220  allocate(Vector<RefCountedPtr<ParticleData<T>>>& a_particles, const std::string a_realm) const;
221 
227  template <typename T>
228  void
229  allocate(ParticleContainer<T>& a_container, const std::string a_realm) const;
230 
236  template <typename T>
237  void
238  allocatePointer(Vector<RefCountedPtr<T>>& a_data) const;
239 
246  template <typename T>
247  void
248  allocatePointer(Vector<RefCountedPtr<T>>& a_data, const int a_finestLevel) const;
249 
256  template <typename T>
257  void
258  allocatePointer(EBAMRData<T>& a_data, const std::string a_realm) const;
259 
267  template <typename T>
268  void
269  allocatePointer(EBAMRData<T>& a_data, const std::string a_realm, const int a_finestLevel) const;
270 
277  EBAMRCellData
278  slice(EBAMRCellData& a_original, const Interval a_variables) const noexcept;
279 
286  const EBAMRCellData
287  slice(const EBAMRCellData& a_original, const Interval a_variables) const noexcept;
288 
292  void
293  parseOptions();
294 
299  void
300  parseRuntimeOptions();
301 
309  void
310  registerOperator(const std::string a_operator, const std::string a_realm, const phase::which_phase a_phase);
311 
320  void
321  registerMask(const std::string a_mask, const int a_buffer, const std::string a_realm);
322 
326  void
327  buildDomains();
328 
334  void
335  setFinestLevel(const int a_finestLevel);
336 
340  void
341  setCoarsestGrid(const IntVect& a_nCells);
342 
348  bool
349  queryRealm(const std::string a_realm) const;
350 
355  bool
356  getEbCf() const;
357 
363  void
364  registerRealm(const std::string a_realm);
365 
373  void
374  alias(EBAMRCellData& a_data,
375  const phase::which_phase a_phase,
376  const MFAMRCellData& a_mfdata,
377  const int a_finestLevel) const;
378 
386  void
387  alias(EBAMRFluxData& a_data,
388  const phase::which_phase a_phase,
389  const MFAMRFluxData& a_mfdata,
390  const int a_finestLevel) const;
391 
398  void
399  alias(EBAMRCellData& a_data, const phase::which_phase a_phase, const MFAMRCellData& a_mfdata) const;
400 
407  void
408  alias(EBAMRFluxData& a_data, const phase::which_phase a_phase, const MFAMRFluxData& a_mfdata) const;
409 
416  void
417  alias(EBAMRIVData& a_data, const phase::which_phase a_phase, const MFAMRIVData& a_mfdata) const;
418 
425  EBAMRCellData
426  alias(const phase::which_phase a_phase, const MFAMRCellData& a_mfdata) const;
427 
434  EBAMRFluxData
435  alias(const phase::which_phase a_phase, const MFAMRFluxData& a_mfdata) const;
436 
443  EBAMRIVData
444  alias(const phase::which_phase a_phase, const MFAMRIVData& a_mfdata) const;
445 
455  void
456  computeGradient(EBAMRCellData& a_gradient,
457  const EBAMRCellData& a_phi,
458  const std::string a_realm,
459  const phase::which_phase a_phase) const;
460 
470  void
471  computeGradient(EBAMRFluxData& a_gradient,
472  const EBAMRCellData& a_phi,
473  const std::string a_realm,
474  const phase::which_phase a_phase) const;
475 
484  void
485  computeGradient(MFAMRCellData& a_gradient, const MFAMRCellData& a_phi, const std::string a_realm) const;
486 
495  void
496  computeGradient(MFAMRFluxData& a_gradient, const MFAMRCellData& a_phi, const std::string a_realm) const;
497 
507  void
508  allocate(EBAMRCellData& a_data,
509  const std::string a_realm,
510  const phase::which_phase a_phase,
511  const int a_nComp,
512  const int a_nGhost = -1) const;
513 
524  void
525  allocate(LevelData<EBCellFAB>& a_data,
526  const std::string a_realm,
527  const phase::which_phase a_phase,
528  const int a_level,
529  const int a_nComp,
530  const int a_nGhost = -1) const;
531 
541  void
542  allocate(EBAMRFluxData& a_data,
543  const std::string a_realm,
544  const phase::which_phase a_phase,
545  const int a_nComp,
546  const int a_nGhost = -1) const;
547 
557  void
558  allocate(EBAMRIVData& a_data,
559  const std::string a_realm,
560  const phase::which_phase a_phase,
561  const int a_nComp,
562  const int a_nGhost = -1) const;
563 
572  void
573  allocate(EBAMRIFData& a_data,
574  const std::string a_realm,
575  const phase::which_phase a_phase,
576  const int a_nComp,
577  const int a_nGhost = -1) const;
578 
587  void
588  allocate(EBAMRBool& a_data, const std::string a_realm, const int a_nComp, const int a_nGhost = -1) const;
589 
598  void
599  allocate(MFAMRCellData& a_data, const std::string a_realm, const int a_nComp, const int a_nGhost = -1) const;
600 
609  void
610  allocate(MFAMRFluxData& a_data, const std::string a_realm, const int a_nComp, const int a_ghost = -1) const;
611 
620  void
621  allocate(MFAMRIVData& a_data, const std::string a_realm, const int a_nComp, const int a_ghost = 0) const;
622 
630  void
631  reallocate(EBAMRCellData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
632 
640  void
641  reallocate(EBAMRFluxData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
642 
650  void
651  reallocate(EBAMRIVData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
652 
660  void
661  reallocate(EBAMRIFData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
662 
669  void
670  reallocate(EBAMRBool& a_data, const int a_lmin) const;
671 
678  void
679  reallocate(MFAMRCellData& a_data, const int a_lmin) const;
680 
687  void
688  reallocate(MFAMRFluxData& a_data, const int a_lmin) const;
689 
696  void
697  reallocate(MFAMRIVData& a_data, const int a_lmin) const;
698 
705  void
706  average(MFAMRCellData& a_data, const std::string a_realm, const Average& a_average) const;
707 
713  void
714  arithmeticAverage(MFAMRCellData& a_data, const std::string a_realm) const;
715 
721  void
722  harmonicAverage(MFAMRCellData& a_data, const std::string a_realm) const;
723 
729  void
730  conservativeAverage(MFAMRCellData& a_data, const std::string a_realm) const;
731 
738  void
739  average(MFAMRFluxData& a_data, const std::string a_realm, const Average& a_average) const;
740 
747  void
748  arithmeticAverage(MFAMRFluxData& a_data, const std::string a_realm) const;
749 
756  void
757  harmonicAverage(MFAMRFluxData& a_data, const std::string a_realm) const;
758 
764  void
765  conservativeAverage(MFAMRFluxData& a_data, const std::string a_realm) const;
766 
774  void
775  average(EBAMRCellData& a_data,
776  const std::string a_realm,
777  const phase::which_phase a_phase,
778  const Average& a_average) const;
779 
788  void
789  average(EBAMRCellData& a_data,
790  const std::string a_realm,
791  const phase::which_phase a_phase,
792  const Average& a_average,
793  const Interval& a_variables) const;
794 
801  void
802  arithmeticAverage(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
803 
811  void
812  arithmeticAverage(EBAMRCellData& a_data,
813  const std::string a_realm,
814  const phase::which_phase a_phase,
815  const Interval& a_variables) const;
816 
823  void
824  harmonicAverage(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
825 
833  void
834  harmonicAverage(EBAMRCellData& a_data,
835  const std::string a_realm,
836  const phase::which_phase a_phase,
837  const Interval& a_variables) const;
838 
845  void
846  conservativeAverage(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
847 
855  void
856  conservativeAverage(EBAMRCellData& a_data,
857  const std::string a_realm,
858  const phase::which_phase a_phase,
859  const Interval& a_variables) const;
860 
868  void
869  average(EBAMRFluxData& a_data,
870  const std::string a_realm,
871  const phase::which_phase a_phase,
872  const Average& a_average) const;
873 
881  void
882  arithmeticAverage(EBAMRFluxData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
883 
891  void
892  harmonicAverage(EBAMRFluxData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
893 
900  void
901  conservativeAverage(EBAMRFluxData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
902 
910  void
911  average(EBAMRIVData& a_data,
912  const std::string a_realm,
913  const phase::which_phase a_phase,
914  const Average& a_average) const;
915 
922  void
923  arithmeticAverage(EBAMRIVData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
924 
931  void
932  harmonicAverage(EBAMRIVData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
933 
940  void
941  conservativeAverage(EBAMRIVData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
942 
964  template <class P, const Real& (P::*particleScalarField)() const>
965  void
966  depositParticles(EBAMRCellData& a_meshData,
967  const std::string& a_realm,
968  const phase::which_phase& a_phase,
969  const ParticleContainer<P>& a_particles,
970  const DepositionType a_depositionType,
971  const CoarseFineDeposition a_coarseFineDeposition,
972  const bool a_forceIrregNGP = false);
973 
985  template <class P, Real (P::*particleScalarField)() const>
986  void
987  depositParticles(EBAMRCellData& a_meshData,
988  const std::string& a_realm,
989  const phase::which_phase& a_phase,
990  const ParticleContainer<P>& a_particles,
991  const DepositionType a_depositionType,
992  const CoarseFineDeposition a_coarseFineDeposition,
993  const bool a_forceIrregNGP = false);
994 
1016  template <class P, const RealVect& (P::*particleVectorField)() const>
1017  void
1018  depositParticles(EBAMRCellData& a_meshData,
1019  const std::string& a_realm,
1020  const phase::which_phase& a_phase,
1021  const ParticleContainer<P>& a_particles,
1022  const DepositionType a_depositionType,
1023  const CoarseFineDeposition a_coarseFineDeposition,
1024  const bool a_forceIrregNGP = false);
1025 
1037  template <class P, RealVect (P::*particleVectorField)() const>
1038  void
1039  depositParticles(EBAMRCellData& a_meshData,
1040  const std::string& a_realm,
1041  const phase::which_phase& a_phase,
1042  const ParticleContainer<P>& a_particles,
1043  const DepositionType a_depositionType,
1044  const CoarseFineDeposition a_coarseFineDeposition,
1045  const bool a_forceIrregNGP = false);
1046 
1058  template <class P, const Real& (P::*particleScalarField)() const>
1059  void
1060  depositParticles(EBAMRIVData& a_meshData,
1061  const std::string& a_realm,
1062  const phase::which_phase& a_phase,
1063  const ParticleContainer<P>& a_particles) const noexcept;
1064 
1076  template <class P, Real (P::*particleScalarField)() const>
1077  void
1078  depositParticles(EBAMRIVData& a_meshData,
1079  const std::string& a_realm,
1080  const phase::which_phase& a_phase,
1081  const ParticleContainer<P>& a_particles) const noexcept;
1082 
1101  template <class P, Real& (P::*particleScalarField)()>
1102  void
1103  interpolateParticles(ParticleContainer<P>& a_particles,
1104  const std::string& a_realm,
1105  const phase::which_phase& a_phase,
1106  const EBAMRCellData& a_meshScalarField,
1107  const DepositionType a_interpType,
1108  const bool a_forceIrregNGP = false) const;
1109 
1128  template <class P, RealVect& (P::*particleVectorField)()>
1129  void
1130  interpolateParticles(ParticleContainer<P>& a_particles,
1131  const std::string& a_realm,
1132  const phase::which_phase& a_phase,
1133  const EBAMRCellData& a_meshVectorField,
1134  const DepositionType a_interpType,
1135  const bool a_forceIrregNGP = false) const;
1136 
1147  template <class P>
1148  void
1149  removeCoveredParticlesIF(ParticleContainer<P>& a_particles,
1150  const phase::which_phase& a_phase,
1151  const Real a_tolerance = 0.0) const;
1152 
1163  template <class P>
1164  void
1165  removeCoveredParticlesDiscrete(ParticleContainer<P>& a_particles,
1166  const phase::which_phase& a_phase,
1167  const Real a_tolerance = 0.0) const;
1168 
1177  template <class P>
1178  void
1179  removeCoveredParticlesVoxels(ParticleContainer<P>& a_particles, const phase::which_phase& a_phase) const;
1180 
1192  template <class P>
1193  void
1194  transferCoveredParticlesIF(ParticleContainer<P>& a_particlesFrom,
1195  ParticleContainer<P>& a_particlesTo,
1196  const phase::which_phase& a_phase,
1197  const Real a_tolerance = 0.0) const;
1198 
1210  template <class P>
1211  void
1212  transferCoveredParticlesDiscrete(ParticleContainer<P>& a_particlesFrom,
1213  ParticleContainer<P>& a_particlesTo,
1214  const phase::which_phase& a_phase,
1215  const Real a_tolerance = 0.0) const;
1216 
1227  template <class P>
1228  void
1229  transferCoveredParticlesVoxels(ParticleContainer<P>& a_particlesFrom,
1230  ParticleContainer<P>& a_particlesTo,
1231  const phase::which_phase& a_phase) const;
1232 
1240  template <class P>
1241  void
1242  transferIrregularParticles(
1243  ParticleContainer<P>& a_dstParticles,
1244  ParticleContainer<P>& a_srcParticles,
1245  const phase::which_phase a_phase,
1246  const std::function<void(P&)> a_transferModifier = [](P&) -> void {
1247  return;
1248  }) const noexcept;
1249 
1269  template <class P>
1270  void
1271  intersectParticlesRaycastIF(
1272  ParticleContainer<P>& a_activeParticles,
1273  ParticleContainer<P>& a_ebParticles,
1274  ParticleContainer<P>& a_domainParticles,
1275  const phase::which_phase a_phase,
1276  const Real a_tolerance,
1277  const bool a_deleteParticles,
1278  const std::function<void(P&)> a_nonDeletionModifier = [](P&) -> void {
1279  return;
1280  }) const noexcept;
1281 
1297  template <class P>
1298  void
1299  intersectParticlesBisectIF(
1300  ParticleContainer<P>& a_activeParticles,
1301  ParticleContainer<P>& a_ebParticles,
1302  ParticleContainer<P>& a_domainParticles,
1303  const phase::which_phase a_phase,
1304  const Real a_bisectionStep,
1305  const bool a_deleteParticles,
1306  const std::function<void(P&)> a_nonDeletionModifier = [](P&) -> void {
1307  return;
1308  }) const noexcept;
1309 
1316  void
1317  interpGhost(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
1318 
1327  void
1328  interpGhost(LevelData<EBCellFAB>& a_fineData,
1329  const LevelData<EBCellFAB>& a_coarData,
1330  const int a_level,
1331  const std::string a_realm,
1332  const phase::which_phase a_phase) const;
1333 
1339  void
1340  interpGhost(MFAMRCellData& a_data, const std::string a_realm) const;
1341 
1348  void
1349  interpGhostPwl(MFAMRCellData& a_data, const std::string a_realm) const;
1350 
1357  void
1358  interpGhostPwl(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
1359 
1368  void
1369  interpGhostMG(MFAMRCellData& a_data, const std::string a_realm) const;
1370 
1379  void
1380  interpGhostMG(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
1381 
1393  void
1394  interpToNewGrids(MFAMRCellData& a_newData,
1395  const MFAMRCellData& a_oldData,
1396  const int a_lmin,
1397  const int a_oldFinestLevel,
1398  const int a_newFinestLevel,
1399  const EBCoarseToFineInterp::Type a_type);
1400 
1413  void
1414  interpToNewGrids(EBAMRCellData& a_newData,
1415  const EBAMRCellData& a_oldData,
1416  const phase::which_phase a_phase,
1417  const int a_lmin,
1418  const int a_oldFinestLevel,
1419  const int a_newFinestLevel,
1420  const EBCoarseToFineInterp::Type a_type);
1421 
1434  void
1435  interpToNewGrids(EBAMRIVData& a_newData,
1436  const EBAMRIVData& a_oldData,
1437  const phase::which_phase a_phase,
1438  const int a_lmin,
1439  const int a_oldFinestLevel,
1440  const int a_newFinestLevel,
1441  const EBCoarseToFineInterp::Type a_type);
1442 
1449  template <class P>
1450  void
1451  remapToNewGrids(ParticleContainer<P>& a_particles, const int a_lmin, const int a_newFinestLevel) const noexcept;
1452 
1460  void
1461  interpToCentroids(EBAMRCellData& a_data, const std::string a_realm, const phase::which_phase a_phase) const;
1462 
1470  void
1471  interpToCentroids(LevelData<EBCellFAB>& a_data,
1472  const std::string a_realm,
1473  const phase::which_phase a_phase,
1474  const int a_level) const;
1475 
1480  void
1481  setMultifluidIndexSpace(const RefCountedPtr<MultiFluidIndexSpace>& a_multiFluidIndexSpace);
1482 
1488  void
1489  setBaseImplicitFunction(const phase::which_phase a_phase, const RefCountedPtr<BaseIF>& a_baseIF);
1490 
1494  void
1495  preRegrid();
1496 
1503  void
1504  regridAmr(const Vector<IntVectSet>& a_tags, const int a_lmin, const int a_hardcap = -1);
1505 
1513  void
1514  regridRealm(const std::string a_realm,
1515  const Vector<Vector<int>>& a_procs,
1516  const Vector<Vector<Box>>& a_boxes,
1517  const int a_lmin);
1518 
1522  void
1523  postRegrid();
1524 
1531  void
1532  setGrids(const Vector<Vector<Box>>& a_boxes, const std::map<std::string, Vector<Vector<long int>>>& a_realmsAndLoads);
1533 
1538  void
1539  regridOperators(const int a_lmin);
1540 
1546  void
1547  regridOperators(const std::string a_realm, const int a_lmin);
1548 
1552  void
1553  sanityCheck() const;
1554 
1558  int
1559  getFinestLevel() const;
1560 
1564  int
1565  getMaxAmrDepth() const;
1566 
1570  int
1571  getMaxSimulationDepth() const;
1572 
1576  int
1577  getBlockingFactor() const;
1578 
1582  int
1583  getMaxBoxSize() const;
1584 
1588  int
1589  getBrBuffer() const;
1590 
1594  int
1595  getMaxEbisBoxSize() const;
1596 
1600  int
1601  getNumberOfGhostCells() const;
1602 
1606  int
1607  getNumberOfEbGhostCells() const;
1608 
1612  int
1613  getRedistributionRadius() const;
1614 
1618  Real
1619  getFinestDx() const;
1620 
1624  RealVect
1625  getProbLo() const;
1626 
1630  RealVect
1631  getProbHi() const;
1632 
1636  ProblemDomain
1637  getFinestDomain() const;
1638 
1642  const Vector<Real>&
1643  getDx() const;
1644 
1648  const Vector<int>&
1649  getRefinementRatios() const;
1650 
1657  int
1658  getRefinementRatio(const int a_level1, const int a_level2) const;
1659 
1664  const RefCountedPtr<BaseIF>&
1665  getBaseImplicitFunction(const phase::which_phase a_phase) const;
1666 
1671  const Vector<DisjointBoxLayout>&
1672  getProxyGrids() const;
1673 
1677  const Vector<DisjointBoxLayout>&
1678  getGrids(const std::string a_realm) const;
1679 
1683  const Vector<ProblemDomain>&
1684  getDomains() const;
1685 
1692  const AMRMask&
1693  getMask(const std::string a_mask, const int a_buffer, const std::string a_realm) const;
1694 
1699  const AMRMask&
1700  getValidCells(const std::string a_realm) const;
1701 
1706  const Vector<RefCountedPtr<LevelTiles>>&
1707  getLevelTiles(const std::string a_realm) const;
1708 
1714  const Vector<RefCountedPtr<EBLevelGrid>>&
1715  getEBLevelGrid(const std::string a_realm, const phase::which_phase a_phase) const;
1716 
1722  const Vector<RefCountedPtr<EBLevelGrid>>&
1723  getEBLevelGridCoFi(const std::string a_realm, const phase::which_phase a_phase) const;
1724 
1730  const Vector<EBISLayout>&
1731  getEBISLayout(const std::string a_realm, const phase::which_phase a_phase) const;
1732 
1737  const Vector<RefCountedPtr<MFLevelGrid>>&
1738  getMFLevelGrid(const std::string a_realm) const;
1739 
1744  Vector<RefCountedPtr<LayoutData<VoFIterator>>>&
1745  getVofIterator(const std::string a_realm, const phase::which_phase a_phase) const;
1746 
1752  const EBAMRFAB&
1753  getLevelset(const std::string a_realm, const phase::which_phase a_phase) const;
1754 
1760  EBAMRParticleMesh&
1761  getParticleMesh(const std::string a_realm, const phase::which_phase a_phase) const;
1762 
1768  EBAMRSurfaceDeposition&
1769  getSurfaceDeposition(const std::string a_realm, const phase::which_phase a_phase) const;
1770 
1774  Vector<RefCountedPtr<EBCoarAve>>&
1775  getCoarseAverage(const std::string a_realm, const phase::which_phase a_phase) const;
1776 
1782  Vector<RefCountedPtr<EBMultigridInterpolator>>&
1783  getMultigridInterpolator(const std::string a_realm, const phase::which_phase a_phase) const;
1784 
1790  Vector<RefCountedPtr<EBCoarseToFineInterp>>&
1791  getFineInterp(const std::string a_realm, const phase::which_phase a_phase) const;
1792 
1798  Vector<RefCountedPtr<EBReflux>>&
1799  getFluxRegister(const std::string a_realm, const phase::which_phase a_phase) const;
1800 
1806  Vector<RefCountedPtr<EBFluxRedistribution>>&
1807  getRedistributionOp(const std::string a_realm, const phase::which_phase a_phase) const;
1808 
1814  const IrregAmrStencil<CentroidInterpolationStencil>&
1815  getCentroidInterpolationStencils(const std::string a_realm, const phase::which_phase a_phase) const;
1816 
1822  const IrregAmrStencil<EbCentroidInterpolationStencil>&
1823  getEbCentroidInterpolationStencils(const std::string a_realm, const phase::which_phase a_phase) const;
1824 
1830  const IrregAmrStencil<NonConservativeDivergenceStencil>&
1831  getNonConservativeDivergenceStencils(const std::string a_realm, const phase::which_phase a_phase) const;
1832 
1837  std::vector<std::string>
1838  getRealms() const;
1839 
1844  BoxSorting
1845  getBoxSorting() const;
1846 
1851  RefCountedPtr<EBIndexSpace>&
1852  getEBIndexSpace(const phase::which_phase a_phase)
1853  {
1854  return m_multifluidIndexSpace->getEBIndexSpace(a_phase);
1855  }
1856 
1857 protected:
1861  enum class GridGenerationMethod
1862  {
1863  BergerRigoutsous,
1864  Tiled,
1865  };
1866 
1870  mutable std::map<std::string, RefCountedPtr<Realm>> m_realms;
1871 
1875  std::map<phase::which_phase, RefCountedPtr<BaseIF>> m_baseif;
1876 
1880  std::map<std::string, Vector<DisjointBoxLayout>> m_oldGrids;
1881 
1886  std::map<std::string, Vector<Copier>> m_oldToNewCellCopiers;
1887 
1892  std::map<std::string, Vector<Copier>> m_oldToNewEBCopiers;
1893 
1898  std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validToValidRealmCopiers;
1899 
1904  std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validToValidGhostRealmCopiers;
1905 
1910  std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validGhostToValidRealmCopiers;
1911 
1916  std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validGhostToValidGhostRealmCopiers;
1917 
1921  GridGenerationMethod m_gridGenerationMethod;
1922 
1926  BoxSorting m_boxSort;
1927 
1931  RefCountedPtr<MultiFluidIndexSpace> m_multifluidIndexSpace;
1932 
1936  IrregStencil::StencilType m_stencilType;
1937 
1941  IrregStencil::StencilType m_centroidStencilType;
1942 
1946  IrregStencil::StencilType m_ebCentroidStencilType;
1947 
1951  IntVect m_numCells;
1952 
1956  Real m_fillRatioBR;
1957 
1961  RealVect m_probLo;
1962 
1966  RealVect m_probHi;
1967 
1971  int m_refRatio;
1972 
1976  int m_verbosity;
1977 
1981  int m_finestLevel;
1982 
1986  int m_oldFinestLevel;
1987 
1991  int m_maxAmrDepth;
1992 
1996  int m_maxSimulationDepth;
1997 
2001  int m_maxBoxSize;
2002 
2006  int m_maxEbisBoxSize;
2007 
2011  int m_bufferSizeBR;
2012 
2016  int m_blockingFactor;
2017 
2021  int m_numEbGhostsCells;
2022 
2026  int m_numGhostCells;
2027 
2031  int m_numLsfGhostCells;
2032 
2036  int m_multigridInterpOrder;
2037 
2041  int m_multigridInterpRadius;
2042 
2046  int m_multigridInterpWeight;
2047 
2051  int m_redistributionRadius;
2052 
2056  int m_centroidStencilOrder;
2057 
2061  int m_centroidStencilRadius;
2062 
2066  int m_ebCentroidStencilOrder;
2067 
2071  int m_ebCentroidStencilRadius;
2072 
2076  int m_irregStenOrder;
2077 
2081  int m_irregStenRadius;
2082 
2086  bool m_hasGrids;
2087 
2091  bool m_hasRegridCopiers;
2092 
2096  Vector<DisjointBoxLayout> m_grids;
2097 
2101  Vector<ProblemDomain> m_domains;
2102 
2106  Vector<int> m_refinementRatios;
2107 
2111  Vector<Real> m_dx;
2112 
2116  void
2117  defineRealms();
2118 
2126  void
2127  buildGrids(const Vector<IntVectSet>& a_tags, const int a_lmin, const int a_hardcap = -1);
2128 
2132  void
2133  buildCopiers();
2134 
2138  void
2139  parseProbLoHiCorners();
2140 
2144  void
2145  parseGridGeneration();
2146 
2150  void
2151  parseVerbosity();
2152 
2156  void
2157  parseCoarsestLevelNumCells();
2158 
2162  void
2163  parseMaxAmrDepth();
2164 
2168  void
2169  parseMaxSimulationDepth();
2170 
2174  void
2175  parseRefinementRatios();
2176 
2180  void
2181  parseBrBufferSize();
2182 
2186  void
2187  parseBrFillRatio();
2188 
2192  void
2193  parseMaxBoxSize();
2194 
2199  void
2200  parseMaxEbisBoxSize();
2201 
2205  void
2206  parseBlockingFactor();
2207 
2212  void
2213  parseEbGhostCells();
2214 
2218  void
2219  parseNumGhostCells();
2220 
2224  void
2225  parseMultigridInterpolator();
2226 
2230  void
2231  parseRedistributionRadius();
2232 
2236  void
2237  parseCentroidStencils();
2238 
2242  void
2243  parseEbCentroidStencils();
2244 
2255  void
2256  computeGradient(LevelData<EBCellFAB>& a_gradient,
2257  const LevelData<EBCellFAB>& a_phi,
2258  const std::string a_realm,
2259  const phase::which_phase a_phase,
2260  const int a_lvl) const;
2261 };
2262 
2263 #include <CD_NamespaceFooter.H>
2264 
2265 #include <CD_AmrMeshImplem.H>
2266 
2267 #endif
CD_AmrMeshImplem.H
Implementation of CD_AmrMesh.H.
Average
Average
Various averaging methods.
Definition: CD_Average.H:24
BoxSorting
BoxSorting
Enum for sorting boxes.
Definition: CD_BoxSorting.H:21
CoarseFineDeposition
CoarseFineDeposition
Coarse-fine deposition types (see CD_EBAMRParticleMesh for how these are handled).
Definition: CD_CoarseFineDeposition.H:26
CD_ComputationalGeometry.H
Declaration of base class for defining geometries.
CD_CopyStrategy.H
Simple enum for distinguishing copying strategies.
CopyStrategy
CopyStrategy
Enum for distinguishing how we copy data Valid => valid region ValidGhost => valid+ghost region.
Definition: CD_CopyStrategy.H:23
DepositionType
DepositionType
Deposition types.
Definition: CD_DepositionType.H:23
CD_DomainFluxIFFABFactory.H
Declaration of factory class for DomainFluxIFFAB.
CD_DomainFluxIFFAB.H
Declaration of a BaseIFFAB wrapper that holds domain fluxes.
CD_EBAMRData.H
Class for holding data across EBAMR hierarchies.
CD_EBCoarAve.H
Declaration of conservative coarsening utility.
CD_EBReflux.H
Declaration of a class which can reflux over the coarse-fine interface.
CD_IrregAmrStencil.H
Class for holding stencils on irregular over an entire AMR hierarchy.
CD_LoadBalancing.H
Declaration of a static class for various load balancing operations.
CD_MFLevelGrid.H
Declaration of a wrapper for wrapping multifluid EBLevelGrids.
CD_MultiFluidIndexSpace.H
Multi-fluid index space.
CD_NonConservativeDivergenceStencil.H
Implementation of IrregStencil that can perform the nonconservative divergence averaging.
CD_ParticleContainer.H
Declaration of a class for holding particles on an AMR hierarchy.
CD_Realm.H
Declaration of the Realm class.
AMRMask
Vector< RefCountedPtr< LevelData< BaseFab< bool > >> > AMRMask
Alias for cutting down on the typic of booleans defined over AMR grids.
Definition: CD_Realm.H:36
AmrMesh
Class for handling spatial operations.
Definition: CD_AmrMesh.H:44
AmrMesh::removeCoveredParticlesVoxels
void removeCoveredParticlesVoxels(ParticleContainer< P > &a_particles, const phase::which_phase &a_phase) const
Function which removes particles from the domain if they fall inside the EB.
Definition: CD_AmrMeshImplem.H:695
AmrMesh::computeGradient
void computeGradient(EBAMRCellData &a_gradient, const EBAMRCellData &a_phi, const std::string a_realm, const phase::which_phase a_phase) const
Compute cell-centered gradient over an AMR hierarchy.
Definition: CD_AmrMesh.cpp:1337
AmrMesh::transferCoveredParticlesDiscrete
void transferCoveredParticlesDiscrete(ParticleContainer< P > &a_particlesFrom, ParticleContainer< P > &a_particlesTo, const phase::which_phase &a_phase, const Real a_tolerance=0.0) const
Function which transferse particles from one particle container to another if they fall inside the EB...
Definition: CD_AmrMeshImplem.H:828
AmrMesh::m_refRatio
int m_refRatio
Refinement ratio.
Definition: CD_AmrMesh.H:1971
AmrMesh::parseRuntimeOptions
void parseRuntimeOptions()
Parse runtime options.
Definition: CD_AmrMesh.cpp:949
AmrMesh::m_oldToNewCellCopiers
std::map< std::string, Vector< Copier > > m_oldToNewCellCopiers
Storage for copiers from the old grids to the new ones.
Definition: CD_AmrMesh.H:1886
AmrMesh::getRealms
std::vector< std::string > getRealms() const
Get the name of all Realms.
Definition: CD_AmrMesh.cpp:3471
AmrMesh::getProbHi
RealVect getProbHi() const
Get upper-left corner of computational domain.
Definition: CD_AmrMesh.cpp:2755
AmrMesh::parseBrBufferSize
void parseBrBufferSize()
Parse buffer size for Berger-Rigoutsous grid algorithm.
Definition: CD_AmrMesh.cpp:2499
AmrMesh::getFinestDx
Real getFinestDx() const
Get resolution on the finest grid level.
Definition: CD_AmrMesh.cpp:2906
AmrMesh::getValidCells
const AMRMask & getValidCells(const std::string a_realm) const
Get a map of all valid cells on a specified realm.
Definition: CD_AmrMesh.cpp:3038
AmrMesh::getBaseImplicitFunction
const RefCountedPtr< BaseIF > & getBaseImplicitFunction(const phase::which_phase a_phase) const
Get implicit function for a specific phase.
Definition: CD_AmrMesh.cpp:2939
AmrMesh::AmrMesh
AmrMesh(const AmrMesh &&a_other)=delete
Disallowed move constructor.
AmrMesh::harmonicAverage
void harmonicAverage(MFAMRCellData &a_data, const std::string a_realm) const
Harmonic coarsening of multifluid data.
Definition: CD_AmrMesh.cpp:1523
AmrMesh::m_fillRatioBR
Real m_fillRatioBR
Fill ratio.
Definition: CD_AmrMesh.H:1956
AmrMesh::getFinestDomain
ProblemDomain getFinestDomain() const
Get finest domain.
Definition: CD_AmrMesh.cpp:2895
AmrMesh::m_oldFinestLevel
int m_oldFinestLevel
Finest level before a regrid.
Definition: CD_AmrMesh.H:1986
AmrMesh::copyData
void copyData(EBAMRData< T > &a_dst, const EBAMRData< T > &a_src, const CopyStrategy &a_toRegion=CopyStrategy::Valid, const CopyStrategy &a_fromRegion=CopyStrategy::Valid) const noexcept
Method for copying from a source container to a destination container. User supplies information abou...
Definition: CD_AmrMeshImplem.H:22
AmrMesh::getProxyGrids
const Vector< DisjointBoxLayout > & getProxyGrids() const
Get the "proxy" grids in AmrMesh.
Definition: CD_AmrMesh.cpp:2961
AmrMesh::transferCoveredParticlesIF
void transferCoveredParticlesIF(ParticleContainer< P > &a_particlesFrom, ParticleContainer< P > &a_particlesTo, const phase::which_phase &a_phase, const Real a_tolerance=0.0) const
Function which transferse particles from one particle container to another if they fall inside the EB...
Definition: CD_AmrMeshImplem.H:754
AmrMesh::m_numCells
IntVect m_numCells
Coarsest box where we compute.
Definition: CD_AmrMesh.H:1951
AmrMesh::postRegrid
void postRegrid()
Run post-regrid operations.
Definition: CD_AmrMesh.cpp:1042
AmrMesh::parseMaxAmrDepth
void parseMaxAmrDepth()
Parse the maximum permitted AMR depth.
Definition: CD_AmrMesh.cpp:2335
AmrMesh::getMultigridInterpolator
Vector< RefCountedPtr< EBMultigridInterpolator > > & getMultigridInterpolator(const std::string a_realm, const phase::which_phase a_phase) const
Get multigrid interpolation utility.
Definition: CD_AmrMesh.cpp:3188
AmrMesh::m_probHi
RealVect m_probHi
Domain simulation corner.
Definition: CD_AmrMesh.H:1966
AmrMesh::m_maxBoxSize
int m_maxBoxSize
Max box size.
Definition: CD_AmrMesh.H:2001
AmrMesh::registerOperator
void registerOperator(const std::string a_operator, const std::string a_realm, const phase::which_phase a_phase)
Register an operator over a realm and a phase.
Definition: CD_AmrMesh.cpp:3351
AmrMesh::m_centroidStencilOrder
int m_centroidStencilOrder
Order for interpolation stencils to centroids.
Definition: CD_AmrMesh.H:2056
AmrMesh::m_grids
Vector< DisjointBoxLayout > m_grids
Grids.
Definition: CD_AmrMesh.H:2096
AmrMesh::m_multifluidIndexSpace
RefCountedPtr< MultiFluidIndexSpace > m_multifluidIndexSpace
MultiFluidIndexSpace.
Definition: CD_AmrMesh.H:1931
AmrMesh::regridRealm
void regridRealm(const std::string a_realm, const Vector< Vector< int >> &a_procs, const Vector< Vector< Box >> &a_boxes, const int a_lmin)
Regrid a realm. This generates the grids for the realm, but does not do the operators on the realm.
Definition: CD_AmrMesh.cpp:3414
AmrMesh::m_probLo
RealVect m_probLo
Domain simulation corner.
Definition: CD_AmrMesh.H:1961
AmrMesh::remapToNewGrids
void remapToNewGrids(ParticleContainer< P > &a_particles, const int a_lmin, const int a_newFinestLevel) const noexcept
Regrid particle to new grids.
Definition: CD_AmrMeshImplem.H:362
AmrMesh::m_bufferSizeBR
int m_bufferSizeBR
Set buffer size.
Definition: CD_AmrMesh.H:2011
AmrMesh::setCoarsestGrid
void setCoarsestGrid(const IntVect &a_nCells)
Set the coarsest grid cells.
Definition: CD_AmrMesh.cpp:2294
AmrMesh::parseMaxSimulationDepth
void parseMaxSimulationDepth()
Set maximum simulation depth.
Definition: CD_AmrMesh.cpp:2354
AmrMesh::getMaxSimulationDepth
int getMaxSimulationDepth() const
Get maximum permitted simulation depth.
Definition: CD_AmrMesh.cpp:2788
AmrMesh::interpGhostPwl
void interpGhostPwl(MFAMRCellData &a_data, const std::string a_realm) const
Interpolate ghost cells over a realm and phase. This uses piecewise linear interpolation (with limite...
Definition: CD_AmrMesh.cpp:1953
AmrMesh::m_maxSimulationDepth
int m_maxSimulationDepth
Maximum allowed depth for simulation.
Definition: CD_AmrMesh.H:1996
AmrMesh::setFinestLevel
void setFinestLevel(const int a_finestLevel)
Set the finest level.
Definition: CD_AmrMesh.cpp:2411
AmrMesh::getEBIndexSpace
RefCountedPtr< EBIndexSpace > & getEBIndexSpace(const phase::which_phase a_phase)
Get EBIndexSpace corresponding to a particular phase.
Definition: CD_AmrMesh.H:1852
AmrMesh::transferIrregularParticles
void transferIrregularParticles(ParticleContainer< P > &a_dstParticles, ParticleContainer< P > &a_srcParticles, const phase::which_phase a_phase, const std::function< void(P &)> a_transferModifier=[](P &) -> void { return;}) const noexcept
Transfer particles that are on the wrong side of the EB to a different container.
Definition: CD_AmrMeshImplem.H:987
AmrMesh::queryRealm
bool queryRealm(const std::string a_realm) const
Query if a realm exists.
Definition: CD_AmrMesh.cpp:3310
AmrMesh::~AmrMesh
virtual ~AmrMesh()
Destructor.
Definition: CD_AmrMesh.cpp:46
AmrMesh::parseEbCentroidStencils
void parseEbCentroidStencils()
Parse EB interpolation (or extrapolation) stencils.
Definition: CD_AmrMesh.cpp:2682
AmrMesh::removeCoveredParticlesIF
void removeCoveredParticlesIF(ParticleContainer< P > &a_particles, const phase::which_phase &a_phase, const Real a_tolerance=0.0) const
Function which removes particles from the domain if they fall inside the EB.
Definition: CD_AmrMeshImplem.H:539
AmrMesh::m_multigridInterpRadius
int m_multigridInterpRadius
Multigrid interpolation radius.
Definition: CD_AmrMesh.H:2041
AmrMesh::arithmeticAverage
void arithmeticAverage(MFAMRCellData &a_data, const std::string a_realm) const
Arithmetic coarsening of multifluid data.
Definition: CD_AmrMesh.cpp:1512
AmrMesh::buildDomains
void buildDomains()
Build domains.
Definition: CD_AmrMesh.cpp:985
AmrMesh::m_ebCentroidStencilRadius
int m_ebCentroidStencilRadius
Radius for EB interpolation stencil.
Definition: CD_AmrMesh.H:2071
AmrMesh::m_gridGenerationMethod
GridGenerationMethod m_gridGenerationMethod
Grid generation method.
Definition: CD_AmrMesh.H:1921
AmrMesh::interpGhostMG
void interpGhostMG(MFAMRCellData &a_data, const std::string a_realm) const
Interpolate ghost cells over a realm and phase.
Definition: CD_AmrMesh.cpp:2016
AmrMesh::setGrids
void setGrids(const Vector< Vector< Box >> &a_boxes, const std::map< std::string, Vector< Vector< long int >>> &a_realmsAndLoads)
Set grids from boxes and computational loads.
Definition: CD_AmrMesh.cpp:2424
AmrMesh::parseRefinementRatios
void parseRefinementRatios()
Parse refinement ratios.
Definition: CD_AmrMesh.cpp:2373
AmrMesh::m_refinementRatios
Vector< int > m_refinementRatios
AMR resolutions.
Definition: CD_AmrMesh.H:2106
AmrMesh::depositParticles
void depositParticles(EBAMRCellData &a_meshData, const std::string &a_realm, const phase::which_phase &a_phase, const ParticleContainer< P > &a_particles, const DepositionType a_depositionType, const CoarseFineDeposition a_coarseFineDeposition, const bool a_forceIrregNGP=false)
Deposit scalar particle quantities on the mesh.
Definition: CD_AmrMeshImplem.H:383
AmrMesh::m_verbosity
int m_verbosity
Verbosity.
Definition: CD_AmrMesh.H:1976
AmrMesh::conservativeAverage
void conservativeAverage(MFAMRCellData &a_data, const std::string a_realm) const
Conservative coarsening of multifluid data.
Definition: CD_AmrMesh.cpp:1534
AmrMesh::m_hasRegridCopiers
bool m_hasRegridCopiers
Has regrid copiers or not.
Definition: CD_AmrMesh.H:2091
AmrMesh::getFineInterp
Vector< RefCountedPtr< EBCoarseToFineInterp > > & getFineInterp(const std::string a_realm, const phase::which_phase a_phase) const
Get interpolator.
Definition: CD_AmrMesh.cpp:3205
AmrMesh::parseMaxBoxSize
void parseMaxBoxSize()
Parse the maximum permitted box size.
Definition: CD_AmrMesh.cpp:2461
AmrMesh::interpToCentroids
void interpToCentroids(EBAMRCellData &a_data, const std::string a_realm, const phase::which_phase a_phase) const
Interpolate data to centroids on realm and phase.
Definition: CD_AmrMesh.cpp:2255
AmrMesh::m_multigridInterpOrder
int m_multigridInterpOrder
Multigrid interpolation order.
Definition: CD_AmrMesh.H:2036
AmrMesh::m_centroidStencilRadius
int m_centroidStencilRadius
Radius for centroid interpolation stencil.
Definition: CD_AmrMesh.H:2061
AmrMesh::m_validToValidRealmCopiers
std::map< std::pair< std::string, std::string >, Vector< Copier > > m_validToValidRealmCopiers
Map for copying between various Realms. First index is the "from" realm and second index is the "to" ...
Definition: CD_AmrMesh.H:1898
AmrMesh::m_validGhostToValidRealmCopiers
std::map< std::pair< std::string, std::string >, Vector< Copier > > m_validGhostToValidRealmCopiers
Map for copying between various Realms. First index is the "from" realm and second index is the "to" ...
Definition: CD_AmrMesh.H:1910
AmrMesh::intersectParticlesRaycastIF
void intersectParticlesRaycastIF(ParticleContainer< P > &a_activeParticles, ParticleContainer< P > &a_ebParticles, ParticleContainer< P > &a_domainParticles, const phase::which_phase a_phase, const Real a_tolerance, const bool a_deleteParticles, const std::function< void(P &)> a_nonDeletionModifier=[](P &) -> void { return;}) const noexcept
Particle intersection algorithm based on ray-casting.
Definition: CD_AmrMeshImplem.H:1071
AmrMesh::parseNumGhostCells
void parseNumGhostCells()
Parse the number of ghost cells.
Definition: CD_AmrMesh.cpp:2589
AmrMesh::operator=
AmrMesh & operator=(const AmrMesh &&a_other)=delete
Disallowed move assignment.
AmrMesh::sanityCheck
void sanityCheck() const
Do a sanity check to make sure everything is set up correctly.
Definition: CD_AmrMesh.cpp:2717
AmrMesh::m_numGhostCells
int m_numGhostCells
Number of ghost cells.
Definition: CD_AmrMesh.H:2026
AmrMesh::buildGrids
void buildGrids(const Vector< IntVectSet > &a_tags, const int a_lmin, const int a_hardcap=-1)
Build new internal AMR grids.
Definition: CD_AmrMesh.cpp:1131
AmrMesh::parseVerbosity
void parseVerbosity()
Parse the verbosity for AmrMesh.
Definition: CD_AmrMesh.cpp:2305
AmrMesh::getEBLevelGrid
const Vector< RefCountedPtr< EBLevelGrid > > & getEBLevelGrid(const std::string a_realm, const phase::which_phase a_phase) const
Get the EBLevelGrid for a Realm and phase.
Definition: CD_AmrMesh.cpp:3070
AmrMesh::GridGenerationMethod
GridGenerationMethod
Enum for having understandable notation for grid generation.
Definition: CD_AmrMesh.H:1862
AmrMesh::m_validGhostToValidGhostRealmCopiers
std::map< std::pair< std::string, std::string >, Vector< Copier > > m_validGhostToValidGhostRealmCopiers
Map for copying between various Realms. First index is the "from" realm and second index is the "to" ...
Definition: CD_AmrMesh.H:1916
AmrMesh::registerRealm
void registerRealm(const std::string a_realm)
Register a new realm.
Definition: CD_AmrMesh.cpp:3338
AmrMesh::removeCoveredParticlesDiscrete
void removeCoveredParticlesDiscrete(ParticleContainer< P > &a_particles, const phase::which_phase &a_phase, const Real a_tolerance=0.0) const
Function which removes particles from the domain if they fall inside the EB.
Definition: CD_AmrMeshImplem.H:605
AmrMesh::getEbCf
bool getEbCf() const
Check if mesh has an EBCF.
Definition: CD_AmrMesh.cpp:3327
AmrMesh::m_dx
Vector< Real > m_dx
Level resolutions.
Definition: CD_AmrMesh.H:2111
AmrMesh::getGrids
const Vector< DisjointBoxLayout > & getGrids(const std::string a_realm) const
Get the grids.
Definition: CD_AmrMesh.cpp:2972
AmrMesh::intersectParticlesBisectIF
void intersectParticlesBisectIF(ParticleContainer< P > &a_activeParticles, ParticleContainer< P > &a_ebParticles, ParticleContainer< P > &a_domainParticles, const phase::which_phase a_phase, const Real a_bisectionStep, const bool a_deleteParticles, const std::function< void(P &)> a_nonDeletionModifier=[](P &) -> void { return;}) const noexcept
Particle intersection algorithm based on bisection.
Definition: CD_AmrMeshImplem.H:1243
AmrMesh::m_numEbGhostsCells
int m_numEbGhostsCells
Number of ghost cells to use for eb stuff.
Definition: CD_AmrMesh.H:2021
AmrMesh::parseOptions
void parseOptions()
Parse options. Called during the constructor.
Definition: CD_AmrMesh.cpp:915
AmrMesh::buildCopiers
void buildCopiers()
Build copiers for copying between realms.
Definition: CD_AmrMesh.cpp:1272
AmrMesh::m_stencilType
IrregStencil::StencilType m_stencilType
Stencil type.
Definition: CD_AmrMesh.H:1936
AmrMesh::getProbLo
RealVect getProbLo() const
Get lower-left corner of computational domain.
Definition: CD_AmrMesh.cpp:2744
AmrMesh::m_boxSort
BoxSorting m_boxSort
Box sorting.
Definition: CD_AmrMesh.H:1926
AmrMesh::setMultifluidIndexSpace
void setMultifluidIndexSpace(const RefCountedPtr< MultiFluidIndexSpace > &a_multiFluidIndexSpace)
Sets multifluid index space.
Definition: CD_AmrMesh.cpp:891
AmrMesh::m_maxAmrDepth
int m_maxAmrDepth
Maximum amr depth.
Definition: CD_AmrMesh.H:1991
AmrMesh::defineRealms
void defineRealms()
Define Realms.
Definition: CD_AmrMesh.cpp:3384
AmrMesh::m_finestLevel
int m_finestLevel
Finest level.
Definition: CD_AmrMesh.H:1981
AmrMesh::reallocate
void reallocate(EBAMRCellData &a_data, const phase::which_phase a_phase, const int a_lmin) const
Reallocate data.
Definition: CD_AmrMesh.cpp:582
AmrMesh::m_oldToNewEBCopiers
std::map< std::string, Vector< Copier > > m_oldToNewEBCopiers
Storage for copiers from the old grids to the new ones.
Definition: CD_AmrMesh.H:1892
AmrMesh::operator=
AmrMesh & operator=(const AmrMesh &a_other)=delete
Disallowed copy assignment.
AmrMesh::getEBISLayout
const Vector< EBISLayout > & getEBISLayout(const std::string a_realm, const phase::which_phase a_phase) const
Get EBISLayouts for a Realm and phase.
Definition: CD_AmrMesh.cpp:2988
AmrMesh::getFluxRegister
Vector< RefCountedPtr< EBReflux > > & getFluxRegister(const std::string a_realm, const phase::which_phase a_phase) const
Get flux register.
Definition: CD_AmrMesh.cpp:3222
AmrMesh::getRedistributionOp
Vector< RefCountedPtr< EBFluxRedistribution > > & getRedistributionOp(const std::string a_realm, const phase::which_phase a_phase) const
Get the redistribution operators.
Definition: CD_AmrMesh.cpp:3239
AmrMesh::m_baseif
std::map< phase::which_phase, RefCountedPtr< BaseIF > > m_baseif
Implicit functions.
Definition: CD_AmrMesh.H:1875
AmrMesh::getNumberOfGhostCells
int getNumberOfGhostCells() const
Get the default number of ghost cells.
Definition: CD_AmrMesh.cpp:2799
AmrMesh::interpGhost
void interpGhost(EBAMRCellData &a_data, const std::string a_realm, const phase::which_phase a_phase) const
Interpolate ghost vectors over a realm, using the default ghost cell interpolation method.
Definition: CD_AmrMesh.cpp:1863
AmrMesh::getMFLevelGrid
const Vector< RefCountedPtr< MFLevelGrid > > & getMFLevelGrid(const std::string a_realm) const
Get EBISLayouts for a Realm.
Definition: CD_AmrMesh.cpp:3104
AmrMesh::regridAmr
void regridAmr(const Vector< IntVectSet > &a_tags, const int a_lmin, const int a_hardcap=-1)
Regrid AMR. This versions generates the grids and Realms, but not the operator.
Definition: CD_AmrMesh.cpp:1079
AmrMesh::parseCoarsestLevelNumCells
void parseCoarsestLevelNumCells()
Parse the coarsest domain grid.
Definition: CD_AmrMesh.cpp:2318
AmrMesh::getEbCentroidInterpolationStencils
const IrregAmrStencil< EbCentroidInterpolationStencil > & getEbCentroidInterpolationStencils(const std::string a_realm, const phase::which_phase a_phase) const
Get interpolation stencils for going from cell center to cell boundary centroid.
Definition: CD_AmrMesh.cpp:3274
AmrMesh::parseRedistributionRadius
void parseRedistributionRadius()
Parse the default redistribution radius.
Definition: CD_AmrMesh.cpp:2631
AmrMesh::interpolateParticles
void interpolateParticles(ParticleContainer< P > &a_particles, const std::string &a_realm, const phase::which_phase &a_phase, const EBAMRCellData &a_meshScalarField, const DepositionType a_interpType, const bool a_forceIrregNGP=false) const
Interpolate a scalar field onto the particle position.
Definition: CD_AmrMeshImplem.H:495
AmrMesh::m_ebCentroidStencilOrder
int m_ebCentroidStencilOrder
Order for interpolation stencilst to EB.
Definition: CD_AmrMesh.H:2066
AmrMesh::getNonConservativeDivergenceStencils
const IrregAmrStencil< NonConservativeDivergenceStencil > & getNonConservativeDivergenceStencils(const std::string a_realm, const phase::which_phase a_phase) const
Get the nonconservative divergence stencils.
Definition: CD_AmrMesh.cpp:3292
AmrMesh::m_validToValidGhostRealmCopiers
std::map< std::pair< std::string, std::string >, Vector< Copier > > m_validToValidGhostRealmCopiers
Map for copying between various Realms. First index is the "from" realm and second index is the "to" ...
Definition: CD_AmrMesh.H:1904
AmrMesh::getDx
const Vector< Real > & getDx() const
Get spatial resolutions.
Definition: CD_AmrMesh.cpp:2917
AmrMesh::getBrBuffer
int getBrBuffer() const
Return buffer for B-R mesh refinement algorithm.
Definition: CD_AmrMesh.cpp:2854
AmrMesh::getMaxBoxSize
int getMaxBoxSize() const
Get maximum permitted box size.
Definition: CD_AmrMesh.cpp:2843
AmrMesh::m_realms
std::map< std::string, RefCountedPtr< Realm > > m_realms
These are all the Realms.
Definition: CD_AmrMesh.H:1870
AmrMesh::getMaxEbisBoxSize
int getMaxEbisBoxSize() const
Get maximum box size for EBIS generation.
Definition: CD_AmrMesh.cpp:2865
AmrMesh::AmrMesh
AmrMesh()
Default constructor.
Definition: CD_AmrMesh.cpp:31
AmrMesh::parseEbGhostCells
void parseEbGhostCells()
Parse number of ghost cells for eb stuff.
Definition: CD_AmrMesh.cpp:2573
AmrMesh::parseMultigridInterpolator
void parseMultigridInterpolator()
Parse settings for the multigrid interpolator.
Definition: CD_AmrMesh.cpp:2609
AmrMesh::getCoarseAverage
Vector< RefCountedPtr< EBCoarAve > > & getCoarseAverage(const std::string a_realm, const phase::which_phase a_phase) const
Get the coarsening utility.
Definition: CD_AmrMesh.cpp:3171
AmrMesh::slice
EBAMRCellData slice(EBAMRCellData &a_original, const Interval a_variables) const noexcept
Slice cell-centered data in order to fetch a subset of components.
Definition: CD_AmrMesh.cpp:50
AmrMesh::m_blockingFactor
int m_blockingFactor
Blocking factor.
Definition: CD_AmrMesh.H:2016
AmrMesh::getRefinementRatio
int getRefinementRatio(const int a_level1, const int a_level2) const
Get refinement factor between two levels.
Definition: CD_AmrMesh.cpp:2876
AmrMesh::getCentroidInterpolationStencils
const IrregAmrStencil< CentroidInterpolationStencil > & getCentroidInterpolationStencils(const std::string a_realm, const phase::which_phase a_phase) const
Get interpolation stencils for going from cell center to cell centroid.
Definition: CD_AmrMesh.cpp:3256
AmrMesh::parseBlockingFactor
void parseBlockingFactor()
Parse the B-R blocking factor. For tiled mesh refinement this parses the tile size.
Definition: CD_AmrMesh.cpp:2557
AmrMesh::interpolateParticles
void interpolateParticles(ParticleContainer< P > &a_particles, const std::string &a_realm, const phase::which_phase &a_phase, const EBAMRCellData &a_meshVectorField, const DepositionType a_interpType, const bool a_forceIrregNGP=false) const
Interpolate a vector field onto the particle position.
AmrMesh::m_irregStenOrder
int m_irregStenOrder
Stencil order for irregular stencils.
Definition: CD_AmrMesh.H:2076
AmrMesh::m_irregStenRadius
int m_irregStenRadius
radius for irregular stencils
Definition: CD_AmrMesh.H:2081
AmrMesh::getBlockingFactor
int getBlockingFactor() const
Get blocking factor (i.e. the smallest possible box).
Definition: CD_AmrMesh.cpp:2832
AmrMesh::allocatePointer
void allocatePointer(Vector< RefCountedPtr< T >> &a_data) const
Allocate pointer but not any memory blocks.
Definition: CD_AmrMeshImplem.H:303
AmrMesh::m_multigridInterpWeight
int m_multigridInterpWeight
Multigrid interpolation weights.
Definition: CD_AmrMesh.H:2046
AmrMesh::m_oldGrids
std::map< std::string, Vector< DisjointBoxLayout > > m_oldGrids
Old grids.
Definition: CD_AmrMesh.H:1880
AmrMesh::getMask
const AMRMask & getMask(const std::string a_mask, const int a_buffer, const std::string a_realm) const
Get a registered mask.
Definition: CD_AmrMesh.cpp:3022
AmrMesh::depositParticles
void depositParticles(EBAMRIVData &a_meshData, const std::string &a_realm, const phase::which_phase &a_phase, const ParticleContainer< P > &a_particles) const noexcept
Deposit scalar particle quantities on the mesh.
AmrMesh::average
void average(MFAMRCellData &a_data, const std::string a_realm, const Average &a_average) const
Average multifluid data over a specified realm.
Definition: CD_AmrMesh.cpp:1472
AmrMesh::getSurfaceDeposition
EBAMRSurfaceDeposition & getSurfaceDeposition(const std::string a_realm, const phase::which_phase a_phase) const
Get EBAMRSurfaceDeposition surface deposition operator.
Definition: CD_AmrMesh.cpp:3154
AmrMesh::preRegrid
void preRegrid()
Run pre-regrid operations.
Definition: CD_AmrMesh.cpp:1010
AmrMesh::registerMask
void registerMask(const std::string a_mask, const int a_buffer, const std::string a_realm)
Register a boolean mask over a realm.
Definition: CD_AmrMesh.cpp:3368
AmrMesh::AmrMesh
AmrMesh(const AmrMesh &a_other)=delete
Disallowed copy constructor.
AmrMesh::parseBrFillRatio
void parseBrFillRatio()
Parse the Berger-Rigoutsos fill ratio.
Definition: CD_AmrMesh.cpp:2395
AmrMesh::m_domains
Vector< ProblemDomain > m_domains
Problem domains.
Definition: CD_AmrMesh.H:2101
AmrMesh::parseProbLoHiCorners
void parseProbLoHiCorners()
Parse the low/high corners of the computational domain.
Definition: CD_AmrMesh.cpp:967
AmrMesh::parseGridGeneration
void parseGridGeneration()
Parse the grid generation algorithm.
Definition: CD_AmrMesh.cpp:2518
AmrMesh::allocate
void allocate(Vector< RefCountedPtr< ParticleData< T >>> &a_particles, const std::string a_realm) const
Template class for generic allocation of particle data.
Definition: CD_AmrMeshImplem.H:237
AmrMesh::interpToNewGrids
void interpToNewGrids(MFAMRCellData &a_newData, const MFAMRCellData &a_oldData, const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel, const EBCoarseToFineInterp::Type a_type)
Interpolate data to new grids.
Definition: CD_AmrMesh.cpp:2080
AmrMesh::transferCoveredParticlesVoxels
void transferCoveredParticlesVoxels(ParticleContainer< P > &a_particlesFrom, ParticleContainer< P > &a_particlesTo, const phase::which_phase &a_phase) const
Function which transferse particles from one particle container to another if they fall inside the EB...
Definition: CD_AmrMeshImplem.H:922
AmrMesh::parseMaxEbisBoxSize
void parseMaxEbisBoxSize()
Parse the maximum permitted box size.
Definition: CD_AmrMesh.cpp:2480
AmrMesh::getParticleMesh
EBAMRParticleMesh & getParticleMesh(const std::string a_realm, const phase::which_phase a_phase) const
Get EBAMRParticleMesh operator.
Definition: CD_AmrMesh.cpp:3137
AmrMesh::alias
void alias(Vector< T * > &a_alias, const Vector< RefCountedPtr< T >> &a_data) const
Turn smart-pointer data structure into regular-pointer data structure.
Definition: CD_AmrMeshImplem.H:209
AmrMesh::m_hasGrids
bool m_hasGrids
Has grids or not.
Definition: CD_AmrMesh.H:2086
AmrMesh::getVofIterator
Vector< RefCountedPtr< LayoutData< VoFIterator > > > & getVofIterator(const std::string a_realm, const phase::which_phase a_phase) const
Get vof iterators for a Realm and phase. This has the capability of iterating through cut-cells.
Definition: CD_AmrMesh.cpp:3005
AmrMesh::m_ebCentroidStencilType
IrregStencil::StencilType m_ebCentroidStencilType
Stencil type for EB interpolation.
Definition: CD_AmrMesh.H:1946
AmrMesh::deallocate
void deallocate(Vector< T * > &a_data) const
Deallocate data.
Definition: CD_AmrMeshImplem.H:161
AmrMesh::getMaxAmrDepth
int getMaxAmrDepth() const
Get maximum permitted amr depth.
Definition: CD_AmrMesh.cpp:2777
AmrMesh::getRefinementRatios
const Vector< int > & getRefinementRatios() const
Get refinement ratios.
Definition: CD_AmrMesh.cpp:2928
AmrMesh::regridOperators
void regridOperators(const int a_lmin)
Regrid AMR operators. This is done for all realms.
Definition: CD_AmrMesh.cpp:1100
AmrMesh::getRedistributionRadius
int getRedistributionRadius() const
Get default redistribution radius.
Definition: CD_AmrMesh.cpp:2821
AmrMesh::getLevelTiles
const Vector< RefCountedPtr< LevelTiles > > & getLevelTiles(const std::string a_realm) const
Get the tiled space representation.
Definition: CD_AmrMesh.cpp:3054
AmrMesh::getBoxSorting
BoxSorting getBoxSorting() const
Get box sorting method.
Definition: CD_AmrMesh.cpp:3488
AmrMesh::parseCentroidStencils
void parseCentroidStencils()
Parse centroid interpolation stencils.
Definition: CD_AmrMesh.cpp:2647
AmrMesh::m_maxEbisBoxSize
int m_maxEbisBoxSize
Maximum box size for EBIS generation.
Definition: CD_AmrMesh.H:2006
AmrMesh::getNumberOfEbGhostCells
int getNumberOfEbGhostCells() const
Get number of ghost cells used for EB grid generation.
Definition: CD_AmrMesh.cpp:2810
AmrMesh::getLevelset
const EBAMRFAB & getLevelset(const std::string a_realm, const phase::which_phase a_phase) const
Get levelset function, allocated over a grid for a Realm and phase.
Definition: CD_AmrMesh.cpp:3121
AmrMesh::m_centroidStencilType
IrregStencil::StencilType m_centroidStencilType
Stencil type for centroid interpolation.
Definition: CD_AmrMesh.H:1941
AmrMesh::setBaseImplicitFunction
void setBaseImplicitFunction(const phase::which_phase a_phase, const RefCountedPtr< BaseIF > &a_baseIF)
Set implicit function for a specific phase. Need e.g. for level-sets.
Definition: CD_AmrMesh.cpp:904
AmrMesh::getFinestLevel
int getFinestLevel() const
Get finest grid level.
Definition: CD_AmrMesh.cpp:2766
AmrMesh::getEBLevelGridCoFi
const Vector< RefCountedPtr< EBLevelGrid > > & getEBLevelGridCoFi(const std::string a_realm, const phase::which_phase a_phase) const
Get the coarsened fine EBLevelGrid for a Realm and phase.
Definition: CD_AmrMesh.cpp:3087
AmrMesh::m_redistributionRadius
int m_redistributionRadius
Redistribution radius.
Definition: CD_AmrMesh.H:2051
AmrMesh::getDomains
const Vector< ProblemDomain > & getDomains() const
Get domains.
Definition: CD_AmrMesh.cpp:2950
AmrMesh::m_numLsfGhostCells
int m_numLsfGhostCells
Number of ghost cells to use when writing level-set to grid.
Definition: CD_AmrMesh.H:2031
EBAMRData
Default class for holding LevelData<T> data across an EBAMR realm.
Definition: CD_EBAMRData.H:40
EBAMRParticleMesh
Class for handling particle-mesh operations with AMR.
Definition: CD_EBAMRParticleMesh.H:47
EBAMRSurfaceDeposition
class for handling surface deposition of particles with EB and AMR.
Definition: CD_EBAMRSurfaceDeposition.H:29
EBCoarseToFineInterp::Type
Type
Type of interpolation methods supported. PWC = Piecewise constant, ignoring the embedded boundary....
Definition: CD_EBCoarseToFineInterp.H:42
IrregAmrStencil
Class for holding VoFStencils on irregular cells over an entire AMR hierarchy. The template parameter...
Definition: CD_IrregAmrStencil.H:29
IrregStencil::StencilType
StencilType
Enum for identifying stencil – only meant for enhancing code visibility.
Definition: CD_IrregStencil.H:44
ParticleContainer
Templated class for holding particles on an AMR hierarchy with particle remapping.
Definition: CD_ParticleContainer.H:50
Generated by doxygen 1.9.1