chombo-discharge
Loading...
Searching...
No Matches
CD_AmrMesh.H
Go to the documentation of this file.
1/*
2 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
3 *
4 * SPDX-License-Identifier: GPL-3.0-or-later
5 */
6
15#ifndef CD_AMRMESH_H
16#define CD_AMRMESH_H
17
18// Chombo includes
19#include <DisjointBoxLayout.H>
20#include <ProblemDomain.H>
21
22// Our includes
23#include <CD_EBAMRData.H>
24#include <CD_EBCoarAve.H>
25#include <CD_ComputationalGeometry.H>
26#include <CD_MultiFluidIndexSpace.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_Realm.H>
34#include <CD_CopyStrategy.H>
35#include <CD_LoadBalancing.H>
36#include <CD_NamespaceHeader.H>
37
42class AmrMesh
43{
44public:
48 AmrMesh();
49
54 AmrMesh(const AmrMesh& a_other) = delete;
55
60 AmrMesh(const AmrMesh&& a_other) = delete;
61
65 virtual ~AmrMesh();
66
71 AmrMesh&
72 operator=(const AmrMesh& a_other) = delete;
73
78 AmrMesh&
79 operator=(const AmrMesh&& a_other) = delete;
80
91 template <typename T>
92 void
93 copyData(EBAMRData<T>& a_dst,
94 const EBAMRData<T>& a_src,
95 const CopyStrategy& a_toRegion = CopyStrategy::Valid,
96 const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
97
109 template <typename T>
110 void
111 copyData(EBAMRData<T>& a_dst,
112 const EBAMRData<T>& a_src,
113 const Interval& a_dstComps,
114 const Interval& a_srcComps,
115 const CopyStrategy& a_toRegion = CopyStrategy::Valid,
116 const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
117
131 template <typename T>
132 void
133 copyData(LevelData<T>& a_dst,
134 const LevelData<T>& a_src,
135 const int a_level,
136 const std::string a_toRealm,
137 const std::string a_fromRealm,
138 const CopyStrategy& a_toRegion = CopyStrategy::Valid,
139 const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
140
156 template <typename T>
157 void
158 copyData(LevelData<T>& a_dst,
159 const LevelData<T>& a_src,
160 const int a_level,
161 const std::string a_toRealm,
162 const std::string a_fromRealm,
163 const Interval& a_dstComps,
164 const Interval& a_srcComps,
165 const CopyStrategy& a_toRegion = CopyStrategy::Valid,
166 const CopyStrategy& a_fromRegion = CopyStrategy::Valid) const noexcept;
167
172 template <typename T>
173 void
174 deallocate(Vector<T*>& a_data) const;
175
180 template <typename T>
181 void
182 deallocate(Vector<RefCountedPtr<T>>& a_data) const;
183
188 template <typename T>
189 void
190 deallocate(EBAMRData<T>& a_data) const;
191
198 template <typename T>
199 void
200 alias(Vector<T*>& a_alias, const Vector<RefCountedPtr<T>>& a_data) const;
201
208 template <typename T, typename S>
209 void
210 alias(Vector<T*>& a_alias, const EBAMRData<S>& a_data) const;
211
217 template <typename T>
218 void
219 allocate(Vector<RefCountedPtr<ParticleData<T>>>& a_particles, const std::string& a_realm) const;
220
226 template <typename T>
227 void
228 allocate(ParticleContainer<T>& a_container, const std::string& a_realm) const;
229
235 template <typename T>
236 void
237 allocatePointer(Vector<RefCountedPtr<T>>& a_data) const;
238
245 template <typename T>
246 void
247 allocatePointer(Vector<RefCountedPtr<T>>& a_data, const int a_finestLevel) const;
248
255 template <typename T>
256 void
257 allocatePointer(EBAMRData<T>& a_data, const std::string& a_realm) const;
258
266 template <typename T>
267 void
268 allocatePointer(EBAMRData<T>& a_data, const std::string& a_realm, const int a_finestLevel) const;
269
276 EBAMRCellData
277 slice(EBAMRCellData& a_original, const Interval a_variables) const noexcept;
278
285 EBAMRCellData
286 slice(const EBAMRCellData& a_original, const Interval a_variables) const noexcept;
287
291 void
292 parseOptions();
293
298 void
299 parseRuntimeOptions();
300
308 void
309 registerOperator(const std::string& a_operator, const std::string& a_realm, const phase::which_phase a_phase);
310
319 void
320 registerMask(const std::string& a_mask, const int a_buffer, const std::string& a_realm);
321
325 void
326 buildDomains();
327
333 void
334 setFinestLevel(const int a_finestLevel);
335
340 void
341 setCoarsestGrid(const IntVect& a_nCells);
342
348 bool
349 queryRealm(const std::string& a_realm) const;
350
356 bool
357 getEbCf() const;
358
364 void
365 registerRealm(const std::string& a_realm);
366
374 void
375 alias(EBAMRCellData& a_data,
376 const phase::which_phase a_phase,
377 const MFAMRCellData& a_mfdata,
378 const int a_finestLevel) const;
379
387 void
388 alias(EBAMRFluxData& a_data,
389 const phase::which_phase a_phase,
390 const MFAMRFluxData& a_mfdata,
391 const int a_finestLevel) const;
392
399 void
400 alias(EBAMRCellData& a_data, const phase::which_phase a_phase, const MFAMRCellData& a_mfdata) const;
401
408 void
409 alias(EBAMRFluxData& a_data, const phase::which_phase a_phase, const MFAMRFluxData& a_mfdata) const;
410
417 void
418 alias(EBAMRIVData& a_data, const phase::which_phase a_phase, const MFAMRIVData& a_mfdata) const;
419
426 EBAMRCellData
427 alias(const phase::which_phase a_phase, const MFAMRCellData& a_mfdata) const;
428
435 EBAMRFluxData
436 alias(const phase::which_phase a_phase, const MFAMRFluxData& a_mfdata) const;
437
444 EBAMRIVData
445 alias(const phase::which_phase a_phase, const MFAMRIVData& a_mfdata) const;
446
456 void
457 computeGradient(EBAMRCellData& a_gradient,
458 const EBAMRCellData& a_phi,
459 const std::string& a_realm,
460 const phase::which_phase a_phase) const;
461
471 void
472 computeGradient(EBAMRFluxData& a_gradient,
473 const EBAMRCellData& a_phi,
474 const std::string& a_realm,
475 const phase::which_phase a_phase) const;
476
485 void
486 computeGradient(MFAMRCellData& a_gradient, const MFAMRCellData& a_phi, const std::string& a_realm) const;
487
496 void
497 computeGradient(MFAMRFluxData& a_gradient, const MFAMRCellData& a_phi, const std::string& a_realm) const;
498
508 void
509 allocate(EBAMRCellData& a_data,
510 const std::string& a_realm,
511 const phase::which_phase a_phase,
512 const int a_nComp,
513 const int a_ghost = -1) const;
514
525 void
526 allocate(LevelData<EBCellFAB>& a_data,
527 const std::string& a_realm,
528 const phase::which_phase a_phase,
529 const int a_level,
530 const int a_nComp,
531 const int a_ghost = -1) const;
532
542 void
543 allocate(EBAMRFluxData& a_data,
544 const std::string& a_realm,
545 const phase::which_phase a_phase,
546 const int a_nComp,
547 const int a_ghost = -1) const;
548
558 void
559 allocate(EBAMRIVData& a_data,
560 const std::string& a_realm,
561 const phase::which_phase a_phase,
562 const int a_nComp,
563 const int a_ghost = -1) const;
564
574 void
575 allocate(EBAMRIFData& a_data,
576 const std::string& a_realm,
577 const phase::which_phase a_phase,
578 const int a_nComp,
579 const int a_ghost = -1) const;
580
589 void
590 allocate(EBAMRBool& a_data, const std::string& a_realm, const int a_nComp, const int a_ghost = -1) const;
591
600 void
601 allocate(MFAMRCellData& a_data, const std::string& a_realm, const int a_nComp, const int a_ghost = -1) const;
602
611 void
612 allocate(MFAMRFluxData& a_data, const std::string& a_realm, const int a_nComp, const int a_ghost = -1) const;
613
622 void
623 allocate(MFAMRIVData& a_data, const std::string& a_realm, const int a_nComp, const int a_ghost = 0) const;
624
632 void
633 reallocate(EBAMRCellData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
634
642 void
643 reallocate(EBAMRFluxData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
644
652 void
653 reallocate(EBAMRIVData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
654
662 void
663 reallocate(EBAMRIFData& a_data, const phase::which_phase a_phase, const int a_lmin) const;
664
671 void
672 reallocate(EBAMRBool& a_data, const int a_lmin) const;
673
680 void
681 reallocate(MFAMRCellData& a_data, const int a_lmin) const;
682
689 void
690 reallocate(MFAMRFluxData& a_data, const int a_lmin) const;
691
698 void
699 reallocate(MFAMRIVData& a_data, const int a_lmin) const;
700
707 void
708 average(MFAMRCellData& a_data, const std::string& a_realm, const Average& a_average) const;
709
715 void
716 arithmeticAverage(MFAMRCellData& a_data, const std::string& a_realm) const;
717
723 void
724 harmonicAverage(MFAMRCellData& a_data, const std::string& a_realm) const;
725
731 void
732 conservativeAverage(MFAMRCellData& a_data, const std::string& a_realm) const;
733
740 void
741 average(MFAMRFluxData& a_data, const std::string& a_realm, const Average& a_average) const;
742
749 void
750 arithmeticAverage(MFAMRFluxData& a_data, const std::string& a_realm) const;
751
758 void
759 harmonicAverage(MFAMRFluxData& a_data, const std::string& a_realm) const;
760
766 void
767 conservativeAverage(MFAMRFluxData& a_data, const std::string& a_realm) const;
768
776 void
777 average(EBAMRCellData& a_data,
778 const std::string& a_realm,
779 const phase::which_phase a_phase,
780 const Average& a_average) const;
781
790 void
791 average(EBAMRCellData& a_data,
792 const std::string& a_realm,
793 const phase::which_phase a_phase,
794 const Average& a_average,
795 const Interval& a_variables) const;
796
803 void
804 arithmeticAverage(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
805
813 void
814 arithmeticAverage(EBAMRCellData& a_data,
815 const std::string& a_realm,
816 const phase::which_phase a_phase,
817 const Interval& a_variables) const;
818
825 void
826 harmonicAverage(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
827
835 void
836 harmonicAverage(EBAMRCellData& a_data,
837 const std::string& a_realm,
838 const phase::which_phase a_phase,
839 const Interval& a_variables) const;
840
847 void
848 conservativeAverage(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
849
857 void
858 conservativeAverage(EBAMRCellData& a_data,
859 const std::string& a_realm,
860 const phase::which_phase a_phase,
861 const Interval& a_variables) const;
862
870 void
871 average(EBAMRFluxData& a_data,
872 const std::string& a_realm,
873 const phase::which_phase a_phase,
874 const Average& a_average) const;
875
883 void
884 arithmeticAverage(EBAMRFluxData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
885
893 void
894 harmonicAverage(EBAMRFluxData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
895
902 void
903 conservativeAverage(EBAMRFluxData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
904
912 void
913 average(EBAMRIVData& a_data,
914 const std::string& a_realm,
915 const phase::which_phase a_phase,
916 const Average& a_average) const;
917
924 void
925 arithmeticAverage(EBAMRIVData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
926
933 void
934 harmonicAverage(EBAMRIVData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
935
942 void
943 conservativeAverage(EBAMRIVData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
944
957 template <class P, typename Ret, Ret (P::*MemberFunc)() const>
958 void
959 depositParticles(EBAMRCellData& a_meshData,
960 const std::string& a_realm,
961 const phase::which_phase& a_phase,
962 const ParticleContainer<P>& a_particles,
963 const DepositionType a_depositionType,
964 const CoarseFineDeposition a_coarseFineDeposition,
965 const bool a_forceIrregNGP = false);
966
976 template <class P, typename Ret, Ret (P::*MemberFunc)() const>
977 void
978 depositParticles(EBAMRIVData& a_meshData,
979 const std::string& a_realm,
980 const phase::which_phase& a_phase,
981 const ParticleContainer<P>& a_particles) const noexcept;
982
994 template <class P, class Ret, Ret (P::*MemberFunc)()>
995 void
996 interpolateParticles(ParticleContainer<P>& a_particles,
997 const std::string& a_realm,
998 const phase::which_phase& a_phase,
999 const EBAMRCellData& a_meshScalarField,
1000 const DepositionType a_interpType,
1001 const bool a_forceIrregNGP = false) const;
1002
1013 template <class P>
1014 void
1015 removeCoveredParticlesIF(ParticleContainer<P>& a_particles,
1016 const phase::which_phase& a_phase,
1017 const Real a_tolerance = 0.0) const;
1018
1029 template <class P>
1030 void
1031 removeCoveredParticlesDiscrete(ParticleContainer<P>& a_particles,
1032 const phase::which_phase& a_phase,
1033 const Real a_tolerance = 0.0) const;
1034
1043 template <class P>
1044 void
1045 removeCoveredParticlesVoxels(ParticleContainer<P>& a_particles, const phase::which_phase& a_phase) const;
1046
1058 template <class P>
1059 void
1060 transferCoveredParticlesIF(ParticleContainer<P>& a_particlesFrom,
1061 ParticleContainer<P>& a_particlesTo,
1062 const phase::which_phase& a_phase,
1063 const Real a_tolerance = 0.0) const;
1064
1076 template <class P>
1077 void
1078 transferCoveredParticlesDiscrete(ParticleContainer<P>& a_particlesFrom,
1079 ParticleContainer<P>& a_particlesTo,
1080 const phase::which_phase& a_phase,
1081 const Real a_tolerance = 0.0) const;
1082
1092 template <class P>
1093 void
1094 transferCoveredParticlesVoxels(ParticleContainer<P>& a_particlesFrom,
1095 ParticleContainer<P>& a_particlesTo,
1096 const phase::which_phase& a_phase) const;
1097
1105 template <class P>
1106 void
1107 transferIrregularParticles(
1108 ParticleContainer<P>& a_dstParticles,
1109 ParticleContainer<P>& a_srcParticles,
1110 const phase::which_phase a_phase,
1111 const std::function<void(P&)>& a_transferModifier = [](P&) -> void {
1112 return;
1113 }) const noexcept;
1114
1134 template <class P>
1135 void
1136 intersectParticlesRaycastIF(
1137 ParticleContainer<P>& a_activeParticles,
1138 ParticleContainer<P>& a_ebParticles,
1139 ParticleContainer<P>& a_domainParticles,
1140 const phase::which_phase a_phase,
1141 const Real a_tolerance,
1142 const bool a_deleteParticles,
1143 const std::function<void(P&)>& a_nonDeletionModifier = [](P&) -> void {
1144 return;
1145 }) const noexcept;
1146
1162 template <class P>
1163 void
1164 intersectParticlesBisectIF(
1165 ParticleContainer<P>& a_activeParticles,
1166 ParticleContainer<P>& a_ebParticles,
1167 ParticleContainer<P>& a_domainParticles,
1168 const phase::which_phase a_phase,
1169 const Real a_bisectionStep,
1170 const bool a_deleteParticles,
1171 const std::function<void(P&)>& a_nonDeletionModifier = [](P&) -> void {
1172 return;
1173 }) const noexcept;
1174
1181 void
1182 interpGhost(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
1183
1192 void
1193 interpGhost(LevelData<EBCellFAB>& a_fineData,
1194 const LevelData<EBCellFAB>& a_coarData,
1195 const int a_fineLevel,
1196 const std::string& a_realm,
1197 const phase::which_phase a_phase) const;
1198
1204 void
1205 interpGhost(MFAMRCellData& a_data, const std::string& a_realm) const;
1206
1212 void
1213 interpGhostPwl(MFAMRCellData& a_data, const std::string& a_realm) const;
1214
1221 void
1222 interpGhostPwl(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
1223
1231 void
1232 interpGhostMG(MFAMRCellData& a_data, const std::string& a_realm) const;
1233
1242 void
1243 interpGhostMG(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const;
1244
1256 void
1257 interpToNewGrids(MFAMRCellData& a_newData,
1258 const MFAMRCellData& a_oldData,
1259 const int a_lmin,
1260 const int a_oldFinestLevel,
1261 const int a_newFinestLevel,
1262 const EBCoarseToFineInterp::Type a_type);
1263
1276 void
1277 interpToNewGrids(EBAMRCellData& a_newData,
1278 const EBAMRCellData& a_oldData,
1279 const phase::which_phase a_phase,
1280 const int a_lmin,
1281 const int a_oldFinestLevel,
1282 const int a_newFinestLevel,
1283 const EBCoarseToFineInterp::Type a_type);
1284
1297 void
1298 interpToNewGrids(EBAMRIVData& a_newData,
1299 const EBAMRIVData& a_oldData,
1300 const phase::which_phase a_phase,
1301 const int a_lmin,
1302 const int a_oldFinestLevel,
1303 const int a_newFinestLevel,
1304 const EBCoarseToFineInterp::Type a_type);
1305
1313 template <class P>
1314 void
1315 remapToNewGrids(ParticleContainer<P>& a_particles, const int a_lmin, const int a_newFinestLevel) const noexcept;
1316
1324 void
1325 interpToCentroids(EBAMRCellData& a_data, const std::string& a_realm, const phase::which_phase a_phase) const noexcept;
1326
1334 void
1335 interpToCentroids(LevelData<EBCellFAB>& a_data,
1336 const std::string& a_realm,
1337 const phase::which_phase a_phase,
1338 const int a_level) const noexcept;
1339
1349 void
1350 interpToCentroids(EBCellFAB& a_centroidData,
1351 const EBCellFAB& a_cellData,
1352 const std::string& a_realm,
1353 const phase::which_phase a_phase,
1354 const int a_level,
1355 const DataIndex& a_din) const noexcept;
1356
1365 void
1366 interpToCentroids(EBAMRIVData& a_centroidData,
1367 const EBAMRCellData& a_cellData,
1368 const std::string& a_realm,
1369 const phase::which_phase a_phase) const noexcept;
1370
1379 void
1380 interpToCentroids(LevelData<BaseIVFAB<Real>>& a_centroidData,
1381 const LevelData<EBCellFAB>& a_cellData,
1382 const std::string& a_realm,
1383 const phase::which_phase a_phase,
1384 const int a_level) const noexcept;
1385
1393 void
1394 interpToEB(EBAMRIVData& a_centroidData,
1395 const EBAMRCellData& a_cellData,
1396 const std::string& a_realm,
1397 const phase::which_phase a_phase) const noexcept;
1398
1407 void
1408 interpToEB(LevelData<BaseIVFAB<Real>>& a_centroidData,
1409 const LevelData<EBCellFAB>& a_cellData,
1410 const std::string& a_realm,
1411 const phase::which_phase a_phase,
1412 const int a_level) const noexcept;
1413
1423 void
1424 interpToEB(BaseIVFAB<Real>& a_centroidData,
1425 const EBCellFAB& a_cellData,
1426 const std::string& a_realm,
1427 const phase::which_phase a_phase,
1428 const int a_level,
1429 const DataIndex& a_din) const noexcept;
1430
1440 void
1441 nonConservativeDivergence(EBAMRIVData& a_nonConsDivF,
1442 const EBAMRCellData& a_kappaDivF,
1443 const std::string& a_realm,
1444 const phase::which_phase& a_phase) const noexcept;
1445
1456 void
1457 nonConservativeDivergence(LevelData<BaseIVFAB<Real>>& a_nonConsDivF,
1458 const LevelData<EBCellFAB>& a_kappaDivF,
1459 const int& a_level,
1460 const std::string& a_realm,
1461 const phase::which_phase& a_phase) const noexcept;
1462
1467 void
1468 setMultifluidIndexSpace(const RefCountedPtr<MultiFluidIndexSpace>& a_multiFluidIndexSpace);
1469
1475 void
1476 setBaseImplicitFunction(const phase::which_phase a_phase, const RefCountedPtr<BaseIF>& a_baseIF);
1477
1481 void
1482 preRegrid();
1483
1490 void
1491 regridAmr(const Vector<IntVectSet>& a_tags, const int a_lmin, const int a_hardcap = -1);
1492
1500 void
1501 regridRealm(const std::string& a_realm,
1502 const Vector<Vector<int>>& a_procs,
1503 const Vector<Vector<Box>>& a_boxes,
1504 const int a_lmin);
1505
1509 void
1510 postRegrid();
1511
1518 void
1519 setGrids(const Vector<Vector<Box>>& a_boxes, const std::map<std::string, Vector<Vector<long int>>>& a_realmsAndLoads);
1520
1525 void
1526 regridOperators(const int a_lmin);
1527
1533 void
1534 regridOperators(const std::string& a_realm, const int a_lmin);
1535
1539 void
1540 sanityCheck() const;
1541
1546 int
1547 getFinestLevel() const;
1548
1553 int
1554 getMaxAmrDepth() const;
1555
1560 int
1561 getMaxSimulationDepth() const;
1562
1567 int
1568 getBlockingFactor() const;
1569
1574 int
1575 getMaxBoxSize() const;
1576
1581 int
1582 getBrBuffer() const;
1583
1588 int
1589 getMaxEbisBoxSize() const;
1590
1595 int
1596 getNumberOfGhostCells() const;
1597
1602 int
1603 getNumberOfEbGhostCells() const;
1604
1609 int
1610 getRedistributionRadius() const;
1611
1616 Real
1617 getFinestDx() const;
1618
1623 RealVect
1624 getProbLo() const;
1625
1630 RealVect
1631 getProbHi() const;
1632
1637 ProblemDomain
1638 getFinestDomain() const;
1639
1644 const Vector<Real>&
1645 getDx() const;
1646
1651 const Vector<int>&
1652 getRefinementRatios() const;
1653
1661 int
1662 getRefinementRatio(const int a_level1, const int a_level2) const;
1663
1669 const RefCountedPtr<BaseIF>&
1670 getBaseImplicitFunction(const phase::which_phase a_phase) const;
1671
1677 const Vector<DisjointBoxLayout>&
1678 getProxyGrids() const;
1679
1685 const Vector<DisjointBoxLayout>&
1686 getGrids(const std::string& a_realm) const;
1687
1692 const Vector<ProblemDomain>&
1693 getDomains() const;
1694
1702 const AMRMask&
1703 getMask(const std::string& a_mask, const int a_buffer, const std::string& a_realm) const;
1704
1710 const AMRMask&
1711 getValidCells(const std::string& a_realm) const;
1712
1718 const Vector<RefCountedPtr<LevelTiles>>&
1719 getLevelTiles(const std::string& a_realm) const;
1720
1727 const Vector<RefCountedPtr<EBLevelGrid>>&
1728 getEBLevelGrid(const std::string& a_realm, const phase::which_phase a_phase) const;
1729
1736 const Vector<RefCountedPtr<EBLevelGrid>>&
1737 getEBLevelGridCoFi(const std::string& a_realm, const phase::which_phase a_phase) const;
1738
1745 const Vector<EBISLayout>&
1746 getEBISLayout(const std::string& a_realm, const phase::which_phase a_phase) const;
1747
1753 const Vector<RefCountedPtr<MFLevelGrid>>&
1754 getMFLevelGrid(const std::string& a_realm) const;
1755
1764 Vector<RefCountedPtr<LayoutData<VoFIterator>>>&
1765 getVofIterator(const std::string& a_realm, const phase::which_phase a_phase) const;
1766
1776 Vector<RefCountedPtr<LayoutData<VoFIterator>>>&
1777 getMultiCutVofIterator(const std::string& a_realm, const phase::which_phase a_phase) const;
1778
1787 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>&
1788 getFaceIterator(const std::string& a_realm, const phase::which_phase a_phase) const;
1789
1798 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>&
1799 getFaceIteratorNoBoundary(const std::string& a_realm, const phase::which_phase a_phase) const;
1800
1812 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>&
1813 getFaceIteratorWithTangentialGhosts(const std::string& a_realm, const phase::which_phase a_phase) const;
1814
1825 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>&
1826 getMultiCutFaceIterator(const std::string& a_realm, const phase::which_phase a_phase) const;
1827
1834 const EBAMRFAB&
1835 getLevelset(const std::string& a_realm, const phase::which_phase a_phase) const;
1836
1843 const EBAMRCellData&
1844 getRegularCells(const std::string& a_realm, const phase::which_phase a_phase) const;
1845
1852 const EBAMRCellData&
1853 getCoveredCells(const std::string& a_realm, const phase::which_phase a_phase) const;
1854
1861 const EBAMRCellData&
1862 getNotCoveredCells(const std::string& a_realm, const phase::which_phase a_phase) const;
1863
1870 const EBAMRCellData&
1871 getIrregularCells(const std::string& a_realm, const phase::which_phase a_phase) const;
1872
1879 EBAMRParticleMesh&
1880 getParticleMesh(const std::string& a_realm, const phase::which_phase a_phase) const;
1881
1888 EBAMRSurfaceDeposition&
1889 getSurfaceDeposition(const std::string& a_realm, const phase::which_phase a_phase) const;
1890
1897 Vector<RefCountedPtr<EBCoarAve>>&
1898 getCoarseAverage(const std::string& a_realm, const phase::which_phase a_phase) const;
1899
1906 Vector<RefCountedPtr<EBMultigridInterpolator>>&
1907 getMultigridInterpolator(const std::string& a_realm, const phase::which_phase a_phase) const;
1908
1915 Vector<RefCountedPtr<EBCoarseToFineInterp>>&
1916 getFineInterp(const std::string& a_realm, const phase::which_phase a_phase) const;
1917
1924 Vector<RefCountedPtr<EBReflux>>&
1925 getFluxRegister(const std::string& a_realm, const phase::which_phase a_phase) const;
1926
1933 Vector<RefCountedPtr<EBFluxRedistribution>>&
1934 getRedistributionOp(const std::string& a_realm, const phase::which_phase a_phase) const;
1935
1936#ifdef CH_USE_PETSC
1941 const RefCountedPtr<PetscGrid>&
1942 getPetscGrid(const std::string& a_realm) const noexcept;
1943#endif
1944
1949 std::vector<std::string>
1950 getRealms() const;
1951
1956 BoxSorting
1957 getBoxSorting() const;
1958
1964 RefCountedPtr<EBIndexSpace>&
1965 getEBIndexSpace(const phase::which_phase a_phase)
1966 {
1967 return m_multifluidIndexSpace->getEBIndexSpace(a_phase);
1968 }
1969
1970protected:
1974 enum class GridGenerationMethod
1975 {
1976 BergerRigoutsous,
1977 Tiled,
1978 };
1979
1983 mutable std::map<std::string, RefCountedPtr<Realm>> m_realms;
1984
1988 std::map<phase::which_phase, RefCountedPtr<BaseIF>> m_baseif;
1989
1993 std::map<std::string, Vector<DisjointBoxLayout>> m_oldGrids;
1994
1999 std::map<std::string, Vector<Copier>> m_oldToNewCellCopiers;
2000
2005 std::map<std::string, Vector<Copier>> m_oldToNewEBCopiers;
2006
2011 std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validToValidRealmCopiers;
2012
2017 std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validToValidGhostRealmCopiers;
2018
2023 std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validGhostToValidRealmCopiers;
2024
2029 std::map<std::pair<std::string, std::string>, Vector<Copier>> m_validGhostToValidGhostRealmCopiers;
2030
2034 GridGenerationMethod m_gridGenerationMethod;
2035
2039 BoxSorting m_boxSort;
2040
2044 RefCountedPtr<MultiFluidIndexSpace> m_multifluidIndexSpace;
2045
2049 CellCentroidInterpolation::Type m_cellCentroidInterpolationType;
2050
2054 EBCentroidInterpolation::Type m_ebCentroidInterpolationType;
2055
2059 IntVect m_numCells;
2060
2064 Real m_fillRatioBR;
2065
2069 RealVect m_probLo;
2070
2074 RealVect m_probHi;
2075
2079 int m_refRatio;
2080
2084 int m_verbosity;
2085
2089 int m_finestLevel;
2090
2094 int m_oldFinestLevel;
2095
2099 int m_maxAmrDepth;
2100
2104 int m_maxSimulationDepth;
2105
2109 int m_maxBoxSize;
2110
2114 int m_maxEbisBoxSize;
2115
2119 int m_bufferSizeBR;
2120
2124 int m_blockingFactor;
2125
2129 int m_numEbGhostsCells;
2130
2134 int m_numGhostCells;
2135
2139 int m_numLsfGhostCells;
2140
2144 int m_multigridInterpOrder;
2145
2149 int m_multigridInterpRadius;
2150
2154 int m_multigridInterpWeight;
2155
2159 int m_redistributionRadius;
2160
2164 bool m_hasGrids;
2165
2169 bool m_hasRegridCopiers;
2170
2174 Vector<DisjointBoxLayout> m_grids;
2175
2179 Vector<ProblemDomain> m_domains;
2180
2184 Vector<int> m_refinementRatios;
2185
2189 Vector<Real> m_dx;
2190
2194 void
2195 defineRealms();
2196
2204 void
2205 buildGrids(const Vector<IntVectSet>& a_tags, const int a_lmin, const int a_hardcap = -1);
2206
2210 void
2211 buildCopiers();
2212
2216 void
2217 parseProbLoHiCorners();
2218
2222 void
2223 parseGridGeneration();
2224
2228 void
2229 parseVerbosity();
2230
2234 void
2235 parseCoarsestLevelNumCells();
2236
2240 void
2241 parseMaxAmrDepth();
2242
2246 void
2247 parseMaxSimulationDepth();
2248
2252 void
2253 parseRefinementRatios();
2254
2258 void
2259 parseBrBufferSize();
2260
2264 void
2265 parseBrFillRatio();
2266
2270 void
2271 parseMaxBoxSize();
2272
2277 void
2278 parseMaxEbisBoxSize();
2279
2283 void
2284 parseBlockingFactor();
2285
2290 void
2291 parseEbGhostCells();
2292
2296 void
2297 parseNumGhostCells();
2298
2302 void
2303 parseMultigridInterpolator();
2304
2308 void
2309 parseRedistributionRadius();
2310
2314 void
2315 parseCellCentroidInterpolation();
2316
2320 void
2321 parseEBCentroidInterpolation();
2322
2333 void
2334 computeGradient(LevelData<EBCellFAB>& a_gradient,
2335 const LevelData<EBCellFAB>& a_phi,
2336 const std::string& a_realm,
2337 const phase::which_phase a_phase,
2338 const int a_lvl) const;
2339};
2340
2341#include <CD_NamespaceFooter.H>
2342
2343#include <CD_AmrMeshImplem.H>
2344
2345#endif
CD_AmrMeshImplem.H
Implementation of CD_AmrMesh.H.
Average
Average
Various averaging methods.
Definition CD_Average.H:25
BoxSorting
BoxSorting
Enum for sorting boxes.
Definition CD_BoxSorting.H:22
CoarseFineDeposition
CoarseFineDeposition
Coarse-fine deposition types (see CD_EBAMRParticleMesh for how these are handled).
Definition CD_CoarseFineDeposition.H:27
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:24
DepositionType
DepositionType
Deposition types.
Definition CD_DepositionType.H:24
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_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_ParticleContainer.H
Declaration of a class for holding particles on an AMR hierarchy.
CD_Realm.H
Declaration of the Realm class.
AmrMesh
Class for handling spatial operations.
Definition CD_AmrMesh.H:43
AmrMesh::operator=
AmrMesh & operator=(const AmrMesh &a_other)=delete
Disallowed copy assignment.
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:606
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:739
AmrMesh::m_refRatio
int m_refRatio
Refinement ratio.
Definition CD_AmrMesh.H:2079
AmrMesh::getFaceIterator
Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > & getFaceIterator(const std::string &a_realm, const phase::which_phase a_phase) const
Get face iterators over cut-cell faces in each valid grid box (FaceStop::SurroundingWithBoundary).
Definition CD_AmrMesh.cpp:3252
AmrMesh::parseRuntimeOptions
void parseRuntimeOptions()
Parse runtime options.
Definition CD_AmrMesh.cpp:956
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:1999
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:3721
AmrMesh::allocatePointer
void allocatePointer(Vector< RefCountedPtr< T > > &a_data) const
Allocate pointer but not any memory blocks.
Definition CD_AmrMeshImplem.H:304
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:3521
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:3321
AmrMesh::getRealms
std::vector< std::string > getRealms() const
Get the name of all Realms.
Definition CD_AmrMesh.cpp:3841
AmrMesh::getProbHi
RealVect getProbHi() const
Get upper-right corner of computational domain.
Definition CD_AmrMesh.cpp:2967
AmrMesh::parseBrBufferSize
void parseBrBufferSize()
Parse buffer size for Berger-Rigoutsous grid algorithm.
Definition CD_AmrMesh.cpp:2653
AmrMesh::getFinestDx
Real getFinestDx() const
Get resolution on the finest grid level.
Definition CD_AmrMesh.cpp:3118
AmrMesh::registerRealm
void registerRealm(const std::string &a_realm)
Register a new realm.
Definition CD_AmrMesh.cpp:3708
AmrMesh::getBaseImplicitFunction
const RefCountedPtr< BaseIF > & getBaseImplicitFunction(const phase::which_phase a_phase) const
Get implicit function for a specific phase.
Definition CD_AmrMesh.cpp:3151
AmrMesh::AmrMesh
AmrMesh(const AmrMesh &&a_other)=delete
Disallowed move constructor.
AmrMesh::m_fillRatioBR
Real m_fillRatioBR
Fill ratio.
Definition CD_AmrMesh.H:2064
AmrMesh::getFinestDomain
ProblemDomain getFinestDomain() const
Get finest domain.
Definition CD_AmrMesh.cpp:3107
AmrMesh::parseCellCentroidInterpolation
void parseCellCentroidInterpolation()
Parse centroid interpolation stencils.
Definition CD_AmrMesh.cpp:2811
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:3538
AmrMesh::getLevelTiles
const Vector< RefCountedPtr< LevelTiles > > & getLevelTiles(const std::string &a_realm) const
Get the tiled space representation.
Definition CD_AmrMesh.cpp:3353
AmrMesh::m_oldFinestLevel
int m_oldFinestLevel
Finest level before a regrid.
Definition CD_AmrMesh.H:2094
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:23
AmrMesh::getProxyGrids
const Vector< DisjointBoxLayout > & getProxyGrids() const
Get the "proxy" grids in AmrMesh.
Definition CD_AmrMesh.cpp:3174
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:665
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:2029
AmrMesh::m_numCells
IntVect m_numCells
Coarsest box where we compute.
Definition CD_AmrMesh.H:2059
AmrMesh::postRegrid
void postRegrid()
Run post-regrid operations.
Definition CD_AmrMesh.cpp:1049
AmrMesh::parseMaxAmrDepth
void parseMaxAmrDepth()
Parse the maximum permitted AMR depth.
Definition CD_AmrMesh.cpp:2489
AmrMesh::m_probHi
RealVect m_probHi
Domain simulation corner.
Definition CD_AmrMesh.H:2074
AmrMesh::m_maxBoxSize
int m_maxBoxSize
Max box size.
Definition CD_AmrMesh.H:2109
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:3337
AmrMesh::m_grids
Vector< DisjointBoxLayout > m_grids
Grids.
Definition CD_AmrMesh.H:2174
AmrMesh::m_multifluidIndexSpace
RefCountedPtr< MultiFluidIndexSpace > m_multifluidIndexSpace
MultiFluidIndexSpace.
Definition CD_AmrMesh.H:2044
AmrMesh::m_probLo
RealVect m_probLo
Domain simulation corner.
Definition CD_AmrMesh.H:2069
AmrMesh::nonConservativeDivergence
void nonConservativeDivergence(EBAMRIVData &a_nonConsDivF, const EBAMRCellData &a_kappaDivF, const std::string &a_realm, const phase::which_phase &a_phase) const noexcept
Compute a non-conservative divergence.
Definition CD_AmrMesh.cpp:3641
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:363
AmrMesh::m_bufferSizeBR
int m_bufferSizeBR
Set buffer size.
Definition CD_AmrMesh.H:2119
AmrMesh::setCoarsestGrid
void setCoarsestGrid(const IntVect &a_nCells)
Set the coarsest grid cells.
Definition CD_AmrMesh.cpp:2448
AmrMesh::parseMaxSimulationDepth
void parseMaxSimulationDepth()
Set maximum simulation depth.
Definition CD_AmrMesh.cpp:2508
AmrMesh::getMaxSimulationDepth
int getMaxSimulationDepth() const
Get maximum permitted simulation depth.
Definition CD_AmrMesh.cpp:3000
AmrMesh::m_maxSimulationDepth
int m_maxSimulationDepth
Maximum allowed depth for simulation.
Definition CD_AmrMesh.H:2104
AmrMesh::setFinestLevel
void setFinestLevel(const int a_finestLevel)
Set the finest level.
Definition CD_AmrMesh.cpp:2565
AmrMesh::~AmrMesh
virtual ~AmrMesh()
Destructor.
AmrMesh::getGrids
const Vector< DisjointBoxLayout > & getGrids(const std::string &a_realm) const
Get the grids.
Definition CD_AmrMesh.cpp:3185
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:450
AmrMesh::m_multigridInterpRadius
int m_multigridInterpRadius
Multigrid interpolation radius.
Definition CD_AmrMesh.H:2149
AmrMesh::buildDomains
void buildDomains()
Build domains.
Definition CD_AmrMesh.cpp:992
AmrMesh::m_gridGenerationMethod
GridGenerationMethod m_gridGenerationMethod
Grid generation method.
Definition CD_AmrMesh.H:2034
AmrMesh::parseRefinementRatios
void parseRefinementRatios()
Parse refinement ratios.
Definition CD_AmrMesh.cpp:2527
AmrMesh::m_refinementRatios
Vector< int > m_refinementRatios
AMR resolutions.
Definition CD_AmrMesh.H:2184
AmrMesh::getFaceIteratorNoBoundary
Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > & getFaceIteratorNoBoundary(const std::string &a_realm, const phase::which_phase a_phase) const
Get face iterators over cut-cell faces in each valid grid box (FaceStop::SurroundingNoBoundary).
Definition CD_AmrMesh.cpp:3269
AmrMesh::getNotCoveredCells
const EBAMRCellData & getNotCoveredCells(const std::string &a_realm, const phase::which_phase a_phase) const
Get the non-covered-cell mask (1 in regular/irregular cells, 0 in covered cells) for a Realm and phas...
Definition CD_AmrMesh.cpp:3470
AmrMesh::m_verbosity
int m_verbosity
Verbosity.
Definition CD_AmrMesh.H:2084
AmrMesh::getMultiCutVofIterator
Vector< RefCountedPtr< LayoutData< VoFIterator > > > & getMultiCutVofIterator(const std::string &a_realm, const phase::which_phase a_phase) const
Get vof iterators restricted to multiply-cut cells for a Realm and phase.
Definition CD_AmrMesh.cpp:3235
AmrMesh::interpToCentroids
void interpToCentroids(EBAMRCellData &a_data, const std::string &a_realm, const phase::which_phase a_phase) const noexcept
Interpolate data to centroids on realm and phase.
Definition CD_AmrMesh.cpp:2270
AmrMesh::getIrregularCells
const EBAMRCellData & getIrregularCells(const std::string &a_realm, const phase::which_phase a_phase) const
Get the irregular-cell mask (value 1 in irregular cells, 0 elsewhere) for a Realm and phase.
Definition CD_AmrMesh.cpp:3487
AmrMesh::m_hasRegridCopiers
bool m_hasRegridCopiers
Has regrid copiers or not.
Definition CD_AmrMesh.H:2169
AmrMesh::parseMaxBoxSize
void parseMaxBoxSize()
Parse the maximum permitted box size.
Definition CD_AmrMesh.cpp:2615
AmrMesh::m_multigridInterpOrder
int m_multigridInterpOrder
Multigrid interpolation order.
Definition CD_AmrMesh.H:2144
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:2011
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:2023
AmrMesh::parseNumGhostCells
void parseNumGhostCells()
Parse the number of ghost cells.
Definition CD_AmrMesh.cpp:2747
AmrMesh::sanityCheck
void sanityCheck() const
Do a sanity check to make sure everything is set up correctly.
Definition CD_AmrMesh.cpp:2897
AmrMesh::m_numGhostCells
int m_numGhostCells
Number of ghost cells.
Definition CD_AmrMesh.H:2134
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:3606
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:1139
AmrMesh::parseVerbosity
void parseVerbosity()
Parse the verbosity for AmrMesh.
Definition CD_AmrMesh.cpp:2459
AmrMesh::GridGenerationMethod
GridGenerationMethod
Enum for having understandable notation for grid generation.
Definition CD_AmrMesh.H:1975
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:2029
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:516
AmrMesh::harmonicAverage
void harmonicAverage(MFAMRCellData &a_data, const std::string &a_realm) const
Harmonic coarsening of multifluid data.
Definition CD_AmrMesh.cpp:1534
AmrMesh::getEbCf
bool getEbCf() const
Check if mesh has an EBCF.
Definition CD_AmrMesh.cpp:3697
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:1154
AmrMesh::m_dx
Vector< Real > m_dx
Level resolutions.
Definition CD_AmrMesh.H:2189
AmrMesh::m_numEbGhostsCells
int m_numEbGhostsCells
Number of ghost cells to use for eb stuff.
Definition CD_AmrMesh.H:2129
AmrMesh::parseOptions
void parseOptions()
Parse options. Called during the constructor.
Definition CD_AmrMesh.cpp:922
AmrMesh::buildCopiers
void buildCopiers()
Build copiers for copying between realms.
Definition CD_AmrMesh.cpp:1280
AmrMesh::getProbLo
RealVect getProbLo() const
Get lower-left corner of computational domain.
Definition CD_AmrMesh.cpp:2956
AmrMesh::m_boxSort
BoxSorting m_boxSort
Box sorting.
Definition CD_AmrMesh.H:2039
AmrMesh::setMultifluidIndexSpace
void setMultifluidIndexSpace(const RefCountedPtr< MultiFluidIndexSpace > &a_multiFluidIndexSpace)
Sets multifluid index space.
Definition CD_AmrMesh.cpp:898
AmrMesh::m_maxAmrDepth
int m_maxAmrDepth
Maximum amr depth.
Definition CD_AmrMesh.H:2099
AmrMesh::defineRealms
void defineRealms()
Define Realms.
Definition CD_AmrMesh.cpp:3754
AmrMesh::m_finestLevel
int m_finestLevel
Finest level.
Definition CD_AmrMesh.H:2089
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:3386
AmrMesh::reallocate
void reallocate(EBAMRCellData &a_data, const phase::which_phase a_phase, const int a_lmin) const
Reallocate data.
Definition CD_AmrMesh.cpp:586
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:982
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:2005
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:384
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:3369
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:898
AmrMesh::m_baseif
std::map< phase::which_phase, RefCountedPtr< BaseIF > > m_baseif
Implicit functions.
Definition CD_AmrMesh.H:1988
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:3784
AmrMesh::getNumberOfGhostCells
int getNumberOfGhostCells() const
Get the default number of ghost cells.
Definition CD_AmrMesh.cpp:3011
AmrMesh::queryRealm
bool queryRealm(const std::string &a_realm) const
Query if a realm exists.
Definition CD_AmrMesh.cpp:3680
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:3420
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:1087
AmrMesh::operator=
AmrMesh & operator=(const AmrMesh &&a_other)=delete
Disallowed move assignment.
AmrMesh::parseCoarsestLevelNumCells
void parseCoarsestLevelNumCells()
Parse the coarsest domain grid.
Definition CD_AmrMesh.cpp:2472
AmrMesh::getEBIndexSpace
RefCountedPtr< EBIndexSpace > & getEBIndexSpace(const phase::which_phase a_phase)
Get EBIndexSpace corresponding to a particular phase.
Definition CD_AmrMesh.H:1965
AmrMesh::arithmeticAverage
void arithmeticAverage(MFAMRCellData &a_data, const std::string &a_realm) const
Arithmetic coarsening of multifluid data.
Definition CD_AmrMesh.cpp:1523
AmrMesh::parseRedistributionRadius
void parseRedistributionRadius()
Parse the default redistribution radius.
Definition CD_AmrMesh.cpp:2794
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:2578
AmrMesh::getFaceIteratorWithTangentialGhosts
Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > & getFaceIteratorWithTangentialGhosts(const std::string &a_realm, const phase::which_phase a_phase) const
Get per-direction face iterators covering cut-cell faces that reach into the tangential ghost layer.
Definition CD_AmrMesh.cpp:3286
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:2017
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:1345
AmrMesh::getDx
const Vector< Real > & getDx() const
Get spatial resolutions.
Definition CD_AmrMesh.cpp:3129
AmrMesh::getBrBuffer
int getBrBuffer() const
Return buffer for B-R mesh refinement algorithm.
Definition CD_AmrMesh.cpp:3066
AmrMesh::getMultiCutFaceIterator
Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > & getMultiCutFaceIterator(const std::string &a_realm, const phase::which_phase a_phase) const
Get face iterators restricted to multiply-cut faces in each valid grid box.
Definition CD_AmrMesh.cpp:3304
AmrMesh::getParticleMesh
EBAMRParticleMesh & getParticleMesh(const std::string &a_realm, const phase::which_phase a_phase) const
Get EBAMRParticleMesh operator.
Definition CD_AmrMesh.cpp:3504
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:3201
AmrMesh::getMaxBoxSize
int getMaxBoxSize() const
Get maximum permitted box size.
Definition CD_AmrMesh.cpp:3055
AmrMesh::m_realms
std::map< std::string, RefCountedPtr< Realm > > m_realms
These are all the Realms.
Definition CD_AmrMesh.H:1983
AmrMesh::getMaxEbisBoxSize
int getMaxEbisBoxSize() const
Get maximum box size for EBIS generation.
Definition CD_AmrMesh.cpp:3077
AmrMesh::AmrMesh
AmrMesh()
Default constructor.
Definition CD_AmrMesh.cpp:32
AmrMesh::parseEbGhostCells
void parseEbGhostCells()
Parse number of ghost cells for eb stuff.
Definition CD_AmrMesh.cpp:2730
AmrMesh::parseMultigridInterpolator
void parseMultigridInterpolator()
Parse settings for the multigrid interpolator.
Definition CD_AmrMesh.cpp:2769
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:46
AmrMesh::m_blockingFactor
int m_blockingFactor
Blocking factor.
Definition CD_AmrMesh.H:2124
AmrMesh::getRefinementRatio
int getRefinementRatio(const int a_level1, const int a_level2) const
Get refinement factor between two levels.
Definition CD_AmrMesh.cpp:3088
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:3738
AmrMesh::getRegularCells
const EBAMRCellData & getRegularCells(const std::string &a_realm, const phase::which_phase a_phase) const
Get the regular-cell mask (value 1 in regular cells, 0 elsewhere) for a Realm and phase.
Definition CD_AmrMesh.cpp:3436
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:3555
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 mesh data onto a particle position.
Definition CD_AmrMeshImplem.H:428
AmrMesh::parseBlockingFactor
void parseBlockingFactor()
Parse the B-R blocking factor. For tiled mesh refinement this parses the tile size.
Definition CD_AmrMesh.cpp:2714
AmrMesh::interpToEB
void interpToEB(EBAMRIVData &a_centroidData, const EBAMRCellData &a_cellData, const std::string &a_realm, const phase::which_phase a_phase) const noexcept
Interpolate data to EB centroids on realm and phase.
Definition CD_AmrMesh.cpp:2379
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:1964
AmrMesh::getBlockingFactor
int getBlockingFactor() const
Get blocking factor (i.e. the smallest possible box).
Definition CD_AmrMesh.cpp:3044
AmrMesh::m_cellCentroidInterpolationType
CellCentroidInterpolation::Type m_cellCentroidInterpolationType
Interpolation type for cell to cell centroid interpolation.
Definition CD_AmrMesh.H:2049
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:238
AmrMesh::parseEBCentroidInterpolation
void parseEBCentroidInterpolation()
Parse EB interpolation (or extrapolation) stencils.
Definition CD_AmrMesh.cpp:2854
AmrMesh::getMFLevelGrid
const Vector< RefCountedPtr< MFLevelGrid > > & getMFLevelGrid(const std::string &a_realm) const
Get EBISLayouts for a Realm.
Definition CD_AmrMesh.cpp:3403
AmrMesh::m_multigridInterpWeight
int m_multigridInterpWeight
Multigrid interpolation weights.
Definition CD_AmrMesh.H:2154
AmrMesh::m_oldGrids
std::map< std::string, Vector< DisjointBoxLayout > > m_oldGrids
Old grids.
Definition CD_AmrMesh.H:1993
AmrMesh::m_ebCentroidInterpolationType
EBCentroidInterpolation::Type m_ebCentroidInterpolationType
Interpolation type for cell to EB centroid interpolation
Definition CD_AmrMesh.H:2054
AmrMesh::getVofIterator
Vector< RefCountedPtr< LayoutData< VoFIterator > > > & getVofIterator(const std::string &a_realm, const phase::which_phase a_phase) const
Get VoF iterators covering all cut-cells (singly- and multiply-cut) for a Realm and phase.
Definition CD_AmrMesh.cpp:3218
AmrMesh::getCoveredCells
const EBAMRCellData & getCoveredCells(const std::string &a_realm, const phase::which_phase a_phase) const
Get the covered-cell mask (value 1 in covered cells, 0 elsewhere) for a Realm and phase.
Definition CD_AmrMesh.cpp:3453
AmrMesh::preRegrid
void preRegrid()
Run pre-regrid operations.
Definition CD_AmrMesh.cpp:1017
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:1483
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:3589
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:2549
AmrMesh::m_domains
Vector< ProblemDomain > m_domains
Problem domains.
Definition CD_AmrMesh.H:2179
AmrMesh::parseProbLoHiCorners
void parseProbLoHiCorners()
Parse the low/high corners of the computational domain.
Definition CD_AmrMesh.cpp:974
AmrMesh::parseGridGeneration
void parseGridGeneration()
Parse the grid generation algorithm.
Definition CD_AmrMesh.cpp:2672
AmrMesh::getFineInterp
Vector< RefCountedPtr< EBCoarseToFineInterp > > & getFineInterp(const std::string &a_realm, const phase::which_phase a_phase) const
Get interpolator.
Definition CD_AmrMesh.cpp:3572
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:2093
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:833
AmrMesh::parseMaxEbisBoxSize
void parseMaxEbisBoxSize()
Parse the maximum permitted box size.
Definition CD_AmrMesh.cpp:2634
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:210
AmrMesh::m_hasGrids
bool m_hasGrids
Has grids or not.
Definition CD_AmrMesh.H:2164
AmrMesh::deallocate
void deallocate(Vector< T * > &a_data) const
Deallocate data.
Definition CD_AmrMeshImplem.H:162
AmrMesh::getMaxAmrDepth
int getMaxAmrDepth() const
Get maximum permitted amr depth.
Definition CD_AmrMesh.cpp:2989
AmrMesh::conservativeAverage
void conservativeAverage(MFAMRCellData &a_data, const std::string &a_realm) const
Conservative coarsening of multifluid data.
Definition CD_AmrMesh.cpp:1545
AmrMesh::getRefinementRatios
const Vector< int > & getRefinementRatios() const
Get refinement ratios.
Definition CD_AmrMesh.cpp:3140
AmrMesh::regridOperators
void regridOperators(const int a_lmin)
Regrid AMR operators. This is done for all realms.
Definition CD_AmrMesh.cpp:1108
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:1874
AmrMesh::getRedistributionRadius
int getRedistributionRadius() const
Get default redistribution radius.
Definition CD_AmrMesh.cpp:3033
AmrMesh::getBoxSorting
BoxSorting getBoxSorting() const
Get box sorting method.
Definition CD_AmrMesh.cpp:3859
AmrMesh::m_maxEbisBoxSize
int m_maxEbisBoxSize
Maximum box size for EBIS generation.
Definition CD_AmrMesh.H:2114
AmrMesh::getNumberOfEbGhostCells
int getNumberOfEbGhostCells() const
Get number of ghost cells used for EB grid generation.
Definition CD_AmrMesh.cpp:3022
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:911
AmrMesh::getFinestLevel
int getFinestLevel() const
Get finest grid level.
Definition CD_AmrMesh.cpp:2978
AmrMesh::m_redistributionRadius
int m_redistributionRadius
Redistribution radius.
Definition CD_AmrMesh.H:2159
AmrMesh::getDomains
const Vector< ProblemDomain > & getDomains() const
Get domains.
Definition CD_AmrMesh.cpp:3163
AmrMesh::m_numLsfGhostCells
int m_numLsfGhostCells
Number of ghost cells to use when writing level-set to grid.
Definition CD_AmrMesh.H:2139
CellCentroidInterpolation::Type
Type
Supported interpolation types.
Definition CD_CellCentroidInterpolation.H:44
EBAMRParticleMesh
Class for handling particle-mesh operations with AMR.
Definition CD_EBAMRParticleMesh.H:53
EBAMRSurfaceDeposition
class for handling surface deposition of particles with EB and AMR.
Definition CD_EBAMRSurfaceDeposition.H:30
EBCentroidInterpolation::Type
Type
Supported interpolation types.
Definition CD_EBCentroidInterpolation.H:44
EBCoarseToFineInterp::Type
Type
Type of interpolation methods supported. PWC = Piecewise constant, ignoring the embedded boundary....
Definition CD_EBCoarseToFineInterp.H:43
TracerParticleSolver
Base class for a tracer particle solver. This solver can advance particles in a pre-defined velocity ...
Definition CD_TracerParticleSolver.H:38
TracerParticleSolver::TracerParticleSolver
TracerParticleSolver()
Default constructor.
Definition CD_TracerParticleSolverImplem.H:26
phase::which_phase
which_phase
Enumeration of supported phases.
Definition CD_MultiFluidIndexSpace.H:38
Generated by doxygen 1.9.8