chombo-discharge
Loading...
Searching...
No Matches
CD_EBHelmholtzOp.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
12#ifndef CD_EBHelmholtzOp_H
13#define CD_EBHelmholtzOp_H
14
15// Std includes
16#include <map>
17
18// Chombo includes
19#include <BaseEBBC.H>
20#include <LevelTGA.H>
21#include <BaseIVFAB.H>
22#include <VCAggStencil.H>
23
24// Our includes
25#include <CD_Timer.H>
26#include <CD_Location.H>
27#include <CD_EBMultigridInterpolator.H>
28#include <CD_EBCoarAve.H>
29#include <CD_EBReflux.H>
30#include <CD_EBMGRestrict.H>
31#include <CD_EBMGProlong.H>
32#include <CD_EBHelmholtzEBBC.H>
33#include <CD_EBHelmholtzDomainBC.H>
34#include <CD_NamespaceHeader.H>
35
40class EBHelmholtzOp : public LevelTGAHelmOp<LevelData<EBCellFAB>, EBFluxFAB>
41{
42public:
46 enum class Smoother
47 {
48 NoRelax,
49 PointJacobi,
50 GauSaiRedBlack,
51 GauSaiMultiColor,
52 };
53
57 EBHelmholtzOp() = delete;
58
63 EBHelmholtzOp(const EBHelmholtzOp& a_other) = delete;
64
69 EBHelmholtzOp(const EBHelmholtzOp&& a_other) = delete;
70
102 EBHelmholtzOp(const Location::Cell a_dataLocation,
103 const EBLevelGrid& a_eblgFine,
104 const EBLevelGrid& a_eblg,
105 const EBLevelGrid& a_eblgCoFi,
106 const EBLevelGrid& a_eblgCoar,
107 const EBLevelGrid& a_eblgCoarMG,
108 const RefCountedPtr<LevelData<BaseFab<bool>>>& a_validCells,
109 const RefCountedPtr<EBMultigridInterpolator>& a_interpolator,
110 const RefCountedPtr<EBReflux>& a_fluxReg,
111 const RefCountedPtr<EBCoarAve>& a_coarAve,
112 const RefCountedPtr<EBHelmholtzDomainBC>& a_domainBC,
113 const RefCountedPtr<EBHelmholtzEBBC>& a_ebBC,
114 const RealVect& a_probLo,
115 const Real& a_dx,
116 const int& a_refToFine,
117 const int& a_refToCoar,
118 const bool& a_hasFine,
119 const bool& a_hasCoar,
120 const bool& a_hasMGObjects,
121 const Real& a_alpha,
122 const Real& a_beta,
123 const RefCountedPtr<LevelData<EBCellFAB>>& a_Acoef,
124 const RefCountedPtr<LevelData<EBFluxFAB>>& a_Bcoef,
125 const RefCountedPtr<LevelData<BaseIVFAB<Real>>>& a_BcoIrreg,
126 const IntVect& a_ghostCellsPhi,
127 const IntVect& a_ghostCellsRHS,
128 const Smoother& a_smoother,
129 const Real& a_relaxFactor);
130
134 virtual ~EBHelmholtzOp();
135
140 EBHelmholtzOp&
141 operator=(const EBHelmholtzOp& a_oper) = delete;
142
147 EBHelmholtzOp&
148 operator=(const EBHelmholtzOp&& a_oper) = delete;
149
153 void
154 turnOffCFInterp();
155
159 void
160 turnOnCFInterp();
161
166 void
167 turnOffExchange();
168
173 void
174 turnOnExchange();
175
180 void
181 turnOffCoarsening();
182
187 void
188 turnOnCoarsening();
189
197 void
198 setAcoAndBco(const RefCountedPtr<LevelData<EBCellFAB>>& a_Acoef,
199 const RefCountedPtr<LevelData<EBFluxFAB>>& a_Bcoef,
200 const RefCountedPtr<LevelData<BaseIVFAB<Real>>>& a_BcoefIrreg);
201
206 const RefCountedPtr<LevelData<EBCellFAB>>&
207 getAcoef();
208
213 const RefCountedPtr<LevelData<EBFluxFAB>>&
214 getBcoef();
215
220 const RefCountedPtr<LevelData<BaseIVFAB<Real>>>&
221 getBcoefIrreg();
222
229 void
230 coarsenCell(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_phiFine);
231
239 void
240 coarsenFlux(LevelData<EBFluxFAB>& a_flux, const LevelData<EBFluxFAB>& a_fineFlux);
241
253 void
254 pointJacobiKernel(EBCellFAB& a_Lcorr,
255 EBCellFAB& a_corr,
256 const EBCellFAB& a_resid,
257 const EBCellFAB& a_Acoef,
258 const EBFluxFAB& a_Bcoef,
259 const BaseIVFAB<Real>& a_BcoefIrreg,
260 const Box& a_cellBox,
261 const DataIndex& a_dit) const noexcept;
262
275 void
276 gauSaiRedBlackKernel(EBCellFAB& a_Lcorr,
277 EBCellFAB& a_corr,
278 const EBCellFAB& a_resid,
279 const EBCellFAB& a_Acoef,
280 const EBFluxFAB& a_Bcoef,
281 const BaseIVFAB<Real>& a_BcoefIrreg,
282 const Box& a_cellBox,
283 const DataIndex& a_dit,
284 const int& a_redBlack) const noexcept;
285
298 void
299 gauSaiMultiColorKernel(EBCellFAB& a_Lcorr,
300 EBCellFAB& a_corr,
301 const EBCellFAB& a_resid,
302 const EBCellFAB& a_Acoef,
303 const EBFluxFAB& a_Bcoef,
304 const BaseIVFAB<Real>& a_BcoefIrreg,
305 const Box& a_cellBox,
306 const DataIndex& a_dit,
307 const IntVect& a_color) const noexcept;
308
316 void
317 residual(LevelData<EBCellFAB>& a_residual,
318 const LevelData<EBCellFAB>& a_phi,
319 const LevelData<EBCellFAB>& a_rhs,
320 const bool a_homogeneousPhysBc) override;
321
328 void
329 preCond(LevelData<EBCellFAB>& a_corr, const LevelData<EBCellFAB>& a_residual) override final;
330
337 void
338 interpolateCF(LevelData<EBCellFAB>& a_phiFine, const LevelData<EBCellFAB>* a_phiCoar, const bool a_homogeneousCFBC);
339
344 void
345 homogeneousCFInterp(LevelData<EBCellFAB>& a_phi);
346
352 void
353 inhomogeneousCFInterp(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_phiCoar);
354
362 void
363 gauSaiMultiColor(LevelData<EBCellFAB>& a_phi,
364 const LevelData<EBCellFAB>& a_Lphi,
365 const LevelData<EBCellFAB>& a_rhs,
366 const IntVect& a_color) const;
367
377 void
378 applyOp(LevelData<EBCellFAB>& a_Lphi,
379 const LevelData<EBCellFAB>& a_phi,
380 const LevelData<EBCellFAB>* const a_phiCoar,
381 const bool a_homogeneousPhysBC,
382 const bool a_homogeneousCFBC);
383
391 void
392 applyOp(LevelData<EBCellFAB>& a_Lphi, const LevelData<EBCellFAB>& a_phi, bool a_homogeneousPhysBc) override final;
393
402 void
403 applyOp(EBCellFAB& a_Lphi,
404 EBCellFAB& a_phi,
405 const EBCellFAB& a_Acoef,
406 const EBFluxFAB& a_Bcoef,
407 const BaseIVFAB<Real>& a_BcoefIrreg,
408 const Box& a_cellBox,
409 const DataIndex& a_dit,
410 const bool a_homogeneousPhysBC) const noexcept;
411
420 void
421 applyOpRegular(EBCellFAB& a_Lphi,
422 EBCellFAB& a_phi,
423 const EBCellFAB& a_Acoef,
424 const EBFluxFAB& a_Bcoef,
425 const BaseIVFAB<Real>& a_BcoefIrreg,
426 const Box& a_cellBox,
427 const DataIndex& a_dit,
428 const bool a_homogeneousPhysBC) const noexcept;
429
440 void
441 applyDomainFlux(EBCellFAB& a_phi,
442 const EBFluxFAB& a_Bcoef,
443 const Box& a_cellBox,
444 const DataIndex& a_dit,
445 const bool a_homogeneousPhysBc) const noexcept;
446
455 void
456 fillDomainFlux(EBFluxFAB& a_flux, const EBCellFAB& a_phi, const Box& a_cellBox, const DataIndex& a_dit);
457
466 void
467 applyOpIrregular(EBCellFAB& a_Lphi,
468 const EBCellFAB& a_phi,
469 const EBCellFAB& a_Acoef,
470 const EBFluxFAB& a_Bcoef,
471 const BaseIVFAB<Real>& a_BcoefIrreg,
472 const BaseIVFAB<Real>& a_alphaDiagWeight,
473 const Box& a_cellBox,
474 const DataIndex& a_dit,
475 const bool a_homogeneousPhysBC) const noexcept;
476
482 void
483 create(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
484
489 void
490 assign(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
491
500 void
501 assignCopier(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const Copier& a_copier) override final;
502
508 void
509 assignLocal(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
510
517 void
518 buildCopier(Copier& a_copier, const LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override;
519
523 Real
524 dotProduct(const LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
525
532 void
533 incr(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const Real a_scale) override final;
534
543 void
544 axby(LevelData<EBCellFAB>& a_lhs,
545 const LevelData<EBCellFAB>& a_x,
546 const LevelData<EBCellFAB>& a_y,
547 const Real a_a,
548 const Real a_b) override final;
549
555 void
556 scale(LevelData<EBCellFAB>& a_lhs, const Real& a_scale) override final;
557
563 void
564 scaleLocal(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) const noexcept;
565
571 Real
572 norm(const LevelData<EBCellFAB>& a_rhs, const int a_order) override final;
573
578 void
579 setToZero(LevelData<EBCellFAB>& a_lhs) override final;
580
587 void
588 createCoarser(LevelData<EBCellFAB>& a_coarse, const LevelData<EBCellFAB>& a_fine, bool a_ghosted) override final;
589
596 void
597 relax(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, int a_iterations) override final;
598
605 void
606 restrictResidual(LevelData<EBCellFAB>& a_resCoar,
607 LevelData<EBCellFAB>& a_phi,
608 const LevelData<EBCellFAB>& a_rhs) override final;
609
615 void
616 prolongIncrement(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_correctCoarse) override final;
617
621 int
622 refToCoarser() override final;
623
634 void
635 AMROperator(LevelData<EBCellFAB>& a_Lphi,
636 const LevelData<EBCellFAB>& a_phiFine,
637 const LevelData<EBCellFAB>& a_phi,
638 const LevelData<EBCellFAB>& a_phiCoar,
639 const bool a_homogeneousPhysBC,
640 AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
641
653 void
654 refluxFreeAMROperator(LevelData<EBCellFAB>& a_Lphi,
655 const LevelData<EBCellFAB>& a_phiFine,
656 const LevelData<EBCellFAB>& a_phi,
657 const LevelData<EBCellFAB>& a_phiCoar,
658 const bool a_homogeneousPhysBC,
659 AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp);
660
671 void
672 AMRResidual(LevelData<EBCellFAB>& a_residual,
673 const LevelData<EBCellFAB>& a_phiFine,
674 const LevelData<EBCellFAB>& a_phi,
675 const LevelData<EBCellFAB>& a_phiCoar,
676 const LevelData<EBCellFAB>& a_rhs,
677 bool a_homogeneousPhysBC,
678 AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
679
688 void
689 AMRResidualNF(LevelData<EBCellFAB>& a_residual,
690 const LevelData<EBCellFAB>& a_phi,
691 const LevelData<EBCellFAB>& a_phiCoar,
692 const LevelData<EBCellFAB>& a_rhs,
693 bool a_homogeneousPhysBC) override final;
694
704 void
705 AMRResidualNC(LevelData<EBCellFAB>& a_residual,
706 const LevelData<EBCellFAB>& a_phiFine,
707 const LevelData<EBCellFAB>& a_phi,
708 const LevelData<EBCellFAB>& a_rhs,
709 bool a_homogeneousPhysBC,
710 AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
711
720 void
721 AMROperatorNF(LevelData<EBCellFAB>& a_Lphi,
722 const LevelData<EBCellFAB>& a_phi,
723 const LevelData<EBCellFAB>& a_phiCoar,
724 bool a_homogeneousPhysBC) override final;
725
734 void
735 AMROperatorNC(LevelData<EBCellFAB>& a_Lphi,
736 const LevelData<EBCellFAB>& a_phi,
737 const LevelData<EBCellFAB>& a_phiCoar,
738 bool a_homogeneousPhysBC,
739 AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
740
749 void
750 AMRRestrict(LevelData<EBCellFAB>& a_residualCoarse,
751 const LevelData<EBCellFAB>& a_residual,
752 const LevelData<EBCellFAB>& a_correction,
753 const LevelData<EBCellFAB>& a_coarseCorrection,
754 bool a_skip_res) override final;
755
761 void
762 AMRProlong(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_coarseCorrection) override final;
763
770 void
771 AMRUpdateResidual(LevelData<EBCellFAB>& a_residual,
772 const LevelData<EBCellFAB>& a_correction,
773 const LevelData<EBCellFAB>& a_coarseCorrection) override final;
774
780 void
781 setAlphaAndBeta(const Real& a_alpha, const Real& a_beta) override final;
782
789 void
790 createCoarsened(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const int& a_refRat) override final;
791
797 void
798 diagonalScale(LevelData<EBCellFAB>& a_rhs, bool a_kappaWeighted) override final;
799
804 void
805 divideByIdentityCoef(LevelData<EBCellFAB>& a_rhs) override final;
806
810 void
811 applyOpNoBoundary(LevelData<EBCellFAB>& a_ans, const LevelData<EBCellFAB>& a_phi) override final;
812
817 void
818 fillGrad(const LevelData<EBCellFAB>& a_phi) override final;
819
828 void
829 getFlux(EBFluxFAB& a_flux,
830 const LevelData<EBCellFAB>& a_data,
831 const Box& a_grid,
832 const DataIndex& a_dit,
833 Real a_scale) override final;
834
838 void
839 allocateFlux() const noexcept;
840
844 void
845 deallocateFlux() const noexcept;
846
850 LevelData<EBFluxFAB>&
851 getFlux() const;
852
856 const LevelData<EBCellFAB>&
857 getRelaxationCoeff() const noexcept;
858
859protected:
863 static constexpr int m_nComp = 1;
864
868 static constexpr int m_comp = 0;
869
873 Timer m_timer;
874
878 Interval m_interval;
879
883 Smoother m_smoother;
884
888 Location::Cell m_dataLocation;
889
893 bool m_refluxFree;
894
898 bool m_profile;
899
903 bool m_hasMGObjects;
904
908 bool m_hasFine;
909
913 bool m_hasCoar;
914
918 int m_refToCoar;
919
923 int m_refToFine;
924
928 int m_numSmoothPreCond;
929
934 bool m_doInterpCF;
935
940 bool m_doExchange;
941
946 bool m_doCoarsen;
947
951 IntVect m_ghostPhi;
952
956 IntVect m_ghostRhs;
957
961 Real m_alpha;
962
966 Real m_beta;
967
971 Real m_dx;
972
976 Real m_relaxFactor;
977
981 RealVect m_vecDx;
982
986 RealVect m_probLo;
987
991 Copier m_exchangeCopier;
992
996 Copier m_exchangeCopierFine;
997
1001 EBLevelGrid m_eblgFine;
1002
1006 EBLevelGrid m_eblg;
1007
1011 EBLevelGrid m_eblgCoFi;
1012
1016 EBLevelGrid m_eblgCoar;
1017
1021 EBLevelGrid m_eblgCoarMG;
1022
1026 EBMGRestrict m_restrictOp;
1027
1031 EBMGRestrict m_restrictOpMG;
1032
1036 EBMGProlong m_prolongOp;
1037
1041 EBMGProlong m_prolongOpMG;
1042
1046 RefCountedPtr<LevelData<BaseFab<bool>>> m_validCells;
1047
1051 RefCountedPtr<EBHelmholtzDomainBC> m_domainBc;
1052
1056 RefCountedPtr<EBHelmholtzEBBC> m_ebBc;
1057
1061 RefCountedPtr<EBMultigridInterpolator> m_interpolator;
1062
1066 RefCountedPtr<EBReflux> m_fluxReg;
1067
1071 RefCountedPtr<EBCoarAve> m_coarAve;
1072
1076 RefCountedPtr<LevelData<EBCellFAB>> m_Acoef;
1077
1081 RefCountedPtr<LevelData<EBFluxFAB>> m_Bcoef;
1082
1086 RefCountedPtr<LevelData<BaseIVFAB<Real>>> m_BcoefIrreg;
1087
1091 LevelData<EBCellFAB> m_relCoef;
1092
1096 mutable LevelData<EBFluxFAB>* m_flux;
1097
1101 LayoutData<BaseIFFAB<Real>> m_interpolant[SpaceDim];
1102
1106 LayoutData<BaseIFFAB<FaceStencil>> m_interpStencil[SpaceDim];
1107
1111 LayoutData<BaseIFFAB<VoFStencil>> m_centroidFluxStencil[SpaceDim];
1112
1118 LayoutData<BaseIVFAB<VoFStencil>> m_relaxStencils;
1119
1124 LayoutData<RefCountedPtr<VCAggStencil>> m_aggRelaxStencil;
1125
1129 LayoutData<BaseIVFAB<Real>> m_alphaDiagWeight;
1130
1134 LayoutData<BaseIVFAB<Real>> m_betaDiagWeight;
1135
1139 mutable LayoutData<VoFIterator> m_vofIterIrreg;
1140
1144 mutable LayoutData<VoFIterator> m_vofIterMulti;
1145
1149 mutable LayoutData<VoFIterator> m_vofIterStenc;
1150
1154 mutable LayoutData<VoFIterator> m_vofIterDomLo[SpaceDim];
1155
1159 mutable LayoutData<VoFIterator> m_vofIterDomHi[SpaceDim];
1160
1164 std::map<std::pair<int, Side::LoHiSide>, Box> m_sideBox;
1165
1169 Vector<IntVect> m_colors;
1170
1177 void
1178 relaxPointJacobi(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1179
1186 void
1187 relaxGSRedBlack(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1188
1195 void
1196 relaxGSMultiColor(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1197
1201 void
1202 computeDiagWeight();
1203
1207 void
1208 computeRelaxationCoefficient();
1209
1213 void
1214 makeAggStencil();
1215
1219 void
1220 defineStencils();
1221
1229 VoFStencil
1230 getFaceCenterFluxStencil(const FaceIndex& a_face, const DataIndex& a_dit) const;
1231
1239 VoFStencil
1240 getFaceCentroidFluxStencil(const FaceIndex& a_face, const DataIndex& a_dit) const;
1241
1250 void
1251 computeFlux(EBFaceFAB& a_fluxCentroid,
1252 const EBCellFAB& a_phi,
1253 const Box& a_cellBox,
1254 const DataIndex& a_dit,
1255 const int a_dir);
1256
1265 void
1266 computeFaceCenteredFlux(EBFaceFAB& a_fluxCenter,
1267 const EBCellFAB& a_phi,
1268 const Box& a_cellBox,
1269 const DataIndex& a_dit,
1270 const int a_dir);
1271
1280 void
1281 computeFaceCentroidFlux(EBFaceFAB& a_flux,
1282 const EBCellFAB& a_phi,
1283 const Box& a_cellBox,
1284 const DataIndex& a_dit,
1285 const int a_dir);
1286
1292 void
1293 computeFlux(const LevelData<EBCellFAB>& a_phi);
1294
1302 void
1303 reflux(LevelData<EBCellFAB>& a_Lphi,
1304 const LevelData<EBCellFAB>& a_phiFine,
1305 const LevelData<EBCellFAB>& a_phi,
1306 AMRLevelOp<LevelData<EBCellFAB>>& a_finerOp);
1307};
1308
1309#include <CD_NamespaceFooter.H>
1310
1311#endif
CD_EBCoarAve.H
Declaration of conservative coarsening utility.
CD_EBMGProlong.H
Declaration of a prolongation operator for EB geometric multigrid.
CD_EBMGRestrict.H
Declaration of a restriction operator for EB geometric multigrid.
CD_EBMultigridInterpolator.H
Declaration of a class that can interpolate more ghost cells near the coarse-fine boundary near the E...
CD_EBReflux.H
Declaration of a class which can reflux over the coarse-fine interface.
CD_Location.H
Declaration of cell positions.
CD_Timer.H
Implementation of CD_Timer.H.
EBCoarAve
Class which replaces data at coarse level of refinement with average at fine level of refinement.
Definition CD_EBCoarAve.H:31
EBHelmholtzDomainBC
Base class for passing domain boundary conditions into EBHelmholtzOp.
Definition CD_EBHelmholtzDomainBC.H:29
EBHelmholtzEBBC
Base class for passing EB boundary conditions into EBHelmholtzOp.
Definition CD_EBHelmholtzEBBC.H:29
EBHelmholtzOp
Helmholtz operator for equations like alpha*a(x)*phi(x) + beta*div(b(x)*grad(phi(x))) = rho.
Definition CD_EBHelmholtzOp.H:41
EBHelmholtzOp::getFaceCentroidFluxStencil
VoFStencil getFaceCentroidFluxStencil(const FaceIndex &a_face, const DataIndex &a_dit) const
Get the face-centroid flux stencil.
Definition CD_EBHelmholtzOp.cpp:2127
EBHelmholtzOp::relaxGSMultiColor
void relaxGSMultiColor(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, const int a_iterations)
Multi-colored gauss-seidel relaxation.
Definition CD_EBHelmholtzOp.cpp:1829
EBHelmholtzOp::setAcoAndBco
void setAcoAndBco(const RefCountedPtr< LevelData< EBCellFAB > > &a_Acoef, const RefCountedPtr< LevelData< EBFluxFAB > > &a_Bcoef, const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &a_BcoefIrreg)
Update with new A and B coefficients.
Definition CD_EBHelmholtzOp.cpp:165
EBHelmholtzOp::reflux
void reflux(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > &a_finerOp)
Reflux algorithm.
Definition CD_EBHelmholtzOp.cpp:2304
EBHelmholtzOp::getFlux
LevelData< EBFluxFAB > & getFlux() const
Returns m_flux. This is used in refluxing routines.
Definition CD_EBHelmholtzOp.cpp:268
EBHelmholtzOp::m_nComp
static constexpr int m_nComp
Number of components that we solve for (always one..)
Definition CD_EBHelmholtzOp.H:863
EBHelmholtzOp::AMRResidualNC
void AMRResidualNC(LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousPhysBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) override final
Compute AMR residual on coarsest.
Definition CD_EBHelmholtzOp.cpp:1069
EBHelmholtzOp::m_coarAve
RefCountedPtr< EBCoarAve > m_coarAve
Conservative coarsener.
Definition CD_EBHelmholtzOp.H:1071
EBHelmholtzOp::axby
void axby(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_x, const LevelData< EBCellFAB > &a_y, const Real a_a, const Real a_b) override final
Set a_lhs = a*x + b*y.
Definition CD_EBHelmholtzOp.cpp:643
EBHelmholtzOp::create
void create(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Create data which clones the layout of the other.
Definition CD_EBHelmholtzOp.cpp:536
EBHelmholtzOp::gauSaiMultiColor
void gauSaiMultiColor(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_rhs, const IntVect &a_color) const
Multi-colored Gauss-Seidel kernel. Public because MFHelmholtzOp may want to use use.
EBHelmholtzOp::m_comp
static constexpr int m_comp
Component that we solve for.
Definition CD_EBHelmholtzOp.H:868
EBHelmholtzOp::turnOnCFInterp
void turnOnCFInterp()
Turn on BCs.
Definition CD_EBHelmholtzOp.cpp:212
EBHelmholtzOp::m_eblgCoFi
EBLevelGrid m_eblgCoFi
Coarsened of m_eblg.
Definition CD_EBHelmholtzOp.H:1011
EBHelmholtzOp::getBcoefIrreg
const RefCountedPtr< LevelData< BaseIVFAB< Real > > > & getBcoefIrreg()
Get the Helmholtz B-coefficient on the EB.
Definition CD_EBHelmholtzOp.cpp:193
EBHelmholtzOp::AMRProlong
void AMRProlong(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection) override final
Prolongation onto AMR level.
Definition CD_EBHelmholtzOp.cpp:1114
EBHelmholtzOp::m_smoother
Smoother m_smoother
Relaxation method.
Definition CD_EBHelmholtzOp.H:883
EBHelmholtzOp::applyOpRegular
void applyOpRegular(EBCellFAB &a_Lphi, EBCellFAB &a_phi, const EBCellFAB &a_Acoef, const EBFluxFAB &a_Bcoef, const BaseIVFAB< Real > &a_BcoefIrreg, const Box &a_cellBox, const DataIndex &a_dit, const bool a_homogeneousPhysBC) const noexcept
Apply operator in regular cells.
Definition CD_EBHelmholtzOp.cpp:1224
EBHelmholtzOp::m_relaxStencils
LayoutData< BaseIVFAB< VoFStencil > > m_relaxStencils
Operator stencils in irregular cells (and ones that border irregular cells if using a centroid discre...
Definition CD_EBHelmholtzOp.H:1118
EBHelmholtzOp::m_hasCoar
bool m_hasCoar
True if there's a coarser level.
Definition CD_EBHelmholtzOp.H:913
EBHelmholtzOp::AMRResidualNF
void AMRResidualNF(LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousPhysBC) override final
Compute AMR residual on finest AMR level.
Definition CD_EBHelmholtzOp.cpp:1054
EBHelmholtzOp::computeFaceCentroidFlux
void computeFaceCentroidFlux(EBFaceFAB &a_flux, const EBCellFAB &a_phi, const Box &a_cellBox, const DataIndex &a_dit, const int a_dir)
Compute face-centered fluxes.
Definition CD_EBHelmholtzOp.cpp:2250
EBHelmholtzOp::applyDomainFlux
void applyDomainFlux(EBCellFAB &a_phi, const EBFluxFAB &a_Bcoef, const Box &a_cellBox, const DataIndex &a_dit, const bool a_homogeneousPhysBc) const noexcept
Apply domain flux.
Definition CD_EBHelmholtzOp.cpp:1278
EBHelmholtzOp::applyOp
void applyOp(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoar, const bool a_homogeneousPhysBC, const bool a_homogeneousCFBC)
Apply operator on this level. This is a more general version which can turn on/off homogeneous and CF...
Definition CD_EBHelmholtzOp.cpp:1150
EBHelmholtzOp::allocateFlux
void allocateFlux() const noexcept
Allocate m_flux.
Definition CD_EBHelmholtzOp.cpp:252
EBHelmholtzOp::computeDiagWeight
void computeDiagWeight()
Calculate the weight of the diagonal term.
Definition CD_EBHelmholtzOp.cpp:1942
EBHelmholtzOp::scale
void scale(LevelData< EBCellFAB > &a_lhs, const Real &a_scale) override final
Scale data. Returns a_lhs = a_lhs*a_scale.
Definition CD_EBHelmholtzOp.cpp:655
EBHelmholtzOp::EBHelmholtzOp
EBHelmholtzOp()=delete
Disallowed default constructor.
EBHelmholtzOp::~EBHelmholtzOp
virtual ~EBHelmholtzOp()
Dtor.
Definition CD_EBHelmholtzOp.cpp:198
EBHelmholtzOp::m_exchangeCopierFine
Copier m_exchangeCopierFine
Pre-built exchange copier.
Definition CD_EBHelmholtzOp.H:996
EBHelmholtzOp::relax
void relax(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, int a_iterations) override final
Relaxation method. This does smoothing for the system L(correction) = residual.
Definition CD_EBHelmholtzOp.cpp:1620
EBHelmholtzOp::turnOffCoarsening
void turnOffCoarsening()
Turn off coarsening operation.
Definition CD_EBHelmholtzOp.cpp:236
EBHelmholtzOp::relaxPointJacobi
void relaxPointJacobi(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, const int a_iterations)
Jacobi relaxation.
Definition CD_EBHelmholtzOp.cpp:1654
EBHelmholtzOp::m_ebBc
RefCountedPtr< EBHelmholtzEBBC > m_ebBc
Domain bc object.
Definition CD_EBHelmholtzOp.H:1056
EBHelmholtzOp::makeAggStencil
void makeAggStencil()
Compute aggregated stencils.
Definition CD_EBHelmholtzOp.cpp:2057
EBHelmholtzOp::homogeneousCFInterp
void homogeneousCFInterp(LevelData< EBCellFAB > &a_phi)
Do homogeneous coarse-fine interpolation.
Definition CD_EBHelmholtzOp.cpp:1569
EBHelmholtzOp::m_eblg
EBLevelGrid m_eblg
Grid.
Definition CD_EBHelmholtzOp.H:1006
EBHelmholtzOp::m_aggRelaxStencil
LayoutData< RefCountedPtr< VCAggStencil > > m_aggRelaxStencil
For making irregular stencil applications go faster.
Definition CD_EBHelmholtzOp.H:1124
EBHelmholtzOp::m_restrictOp
EBMGRestrict m_restrictOp
Restriction operator for AMR levels.
Definition CD_EBHelmholtzOp.H:1026
EBHelmholtzOp::pointJacobiKernel
void pointJacobiKernel(EBCellFAB &a_Lcorr, EBCellFAB &a_corr, const EBCellFAB &a_resid, const EBCellFAB &a_Acoef, const EBFluxFAB &a_Bcoef, const BaseIVFAB< Real > &a_BcoefIrreg, const Box &a_cellBox, const DataIndex &a_dit) const noexcept
Point Jacobi kernel.
Definition CD_EBHelmholtzOp.cpp:1694
EBHelmholtzOp::m_refluxFree
bool m_refluxFree
Use reflux-free formulation or not.
Definition CD_EBHelmholtzOp.H:893
EBHelmholtzOp::m_domainBc
RefCountedPtr< EBHelmholtzDomainBC > m_domainBc
Domain bc object.
Definition CD_EBHelmholtzOp.H:1051
EBHelmholtzOp::dotProduct
Real dotProduct(const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Compute the dot product??
Definition CD_EBHelmholtzOp.cpp:568
EBHelmholtzOp::m_fluxReg
RefCountedPtr< EBReflux > m_fluxReg
Flux register.
Definition CD_EBHelmholtzOp.H:1066
EBHelmholtzOp::preCond
void preCond(LevelData< EBCellFAB > &a_corr, const LevelData< EBCellFAB > &a_residual) override final
Precondition system before bottom solve.
Definition CD_EBHelmholtzOp.cpp:525
EBHelmholtzOp::m_numSmoothPreCond
int m_numSmoothPreCond
Number of smoothings in the preconditioner.
Definition CD_EBHelmholtzOp.H:928
EBHelmholtzOp::m_doInterpCF
bool m_doInterpCF
Do coarse-fine interpolation or not.
Definition CD_EBHelmholtzOp.H:934
EBHelmholtzOp::m_refToFine
int m_refToFine
Refinement factor to fine level.
Definition CD_EBHelmholtzOp.H:923
EBHelmholtzOp::Smoother
Smoother
Relaxation method for the operators.
Definition CD_EBHelmholtzOp.H:47
EBHelmholtzOp::m_flux
LevelData< EBFluxFAB > * m_flux
For holding fluxes.
Definition CD_EBHelmholtzOp.H:1096
EBHelmholtzOp::m_hasMGObjects
bool m_hasMGObjects
True if there is a multigrid level below this operator.
Definition CD_EBHelmholtzOp.H:903
EBHelmholtzOp::m_dx
Real m_dx
Grid resolution;.
Definition CD_EBHelmholtzOp.H:971
EBHelmholtzOp::m_profile
bool m_profile
Profile the operator.
Definition CD_EBHelmholtzOp.H:898
EBHelmholtzOp::turnOffExchange
void turnOffExchange()
Turn off exchange operation.
Definition CD_EBHelmholtzOp.cpp:220
EBHelmholtzOp::relaxGSRedBlack
void relaxGSRedBlack(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, const int a_iterations)
Jacobi relaxation.
Definition CD_EBHelmholtzOp.cpp:1720
EBHelmholtzOp::m_alphaDiagWeight
LayoutData< BaseIVFAB< Real > > m_alphaDiagWeight
Weights of diagonal alpha terms.
Definition CD_EBHelmholtzOp.H:1129
EBHelmholtzOp::fillGrad
void fillGrad(const LevelData< EBCellFAB > &a_phi) override final
Not called, I think.
Definition CD_EBHelmholtzOp.cpp:1549
EBHelmholtzOp::getFaceCenterFluxStencil
VoFStencil getFaceCenterFluxStencil(const FaceIndex &a_face, const DataIndex &a_dit) const
Get the face-centered flux stencil.
Definition CD_EBHelmholtzOp.cpp:2108
EBHelmholtzOp::restrictResidual
void restrictResidual(LevelData< EBCellFAB > &a_resCoar, LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) override final
Restrict residual onto coarse level.
Definition CD_EBHelmholtzOp.cpp:768
EBHelmholtzOp::divideByIdentityCoef
void divideByIdentityCoef(LevelData< EBCellFAB > &a_rhs) override final
Divide by the a-coefficient.
Definition CD_EBHelmholtzOp.cpp:1523
EBHelmholtzOp::assign
void assign(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Assign data.
Definition CD_EBHelmholtzOp.cpp:544
EBHelmholtzOp::m_probLo
RealVect m_probLo
Lower-left corner of domain.
Definition CD_EBHelmholtzOp.H:986
EBHelmholtzOp::m_interpolator
RefCountedPtr< EBMultigridInterpolator > m_interpolator
Interpolator object.
Definition CD_EBHelmholtzOp.H:1061
EBHelmholtzOp::computeRelaxationCoefficient
void computeRelaxationCoefficient()
Calculate relaxation coefficient.
Definition CD_EBHelmholtzOp.cpp:1997
EBHelmholtzOp::computeFaceCenteredFlux
void computeFaceCenteredFlux(EBFaceFAB &a_fluxCenter, const EBCellFAB &a_phi, const Box &a_cellBox, const DataIndex &a_dit, const int a_dir)
Compute face-centered fluxes.
Definition CD_EBHelmholtzOp.cpp:2222
EBHelmholtzOp::AMROperatorNF
void AMROperatorNF(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, bool a_homogeneousPhysBC) override final
Apply the AMR operator, i.e. compute L(phi) in an AMR context, assuming no finer levels.
Definition CD_EBHelmholtzOp.cpp:1012
EBHelmholtzOp::m_prolongOp
EBMGProlong m_prolongOp
Prolongation operator for AMR levels.
Definition CD_EBHelmholtzOp.H:1036
EBHelmholtzOp::refluxFreeAMROperator
void refluxFreeAMROperator(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, const bool a_homogeneousPhysBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
Apply the AMR operator, i.e. compute L(phi) in an AMR context.
Definition CD_EBHelmholtzOp.cpp:840
EBHelmholtzOp::AMROperatorNC
void AMROperatorNC(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, bool a_homogeneousPhysBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) override final
Apply the AMR operator, i.e. compute L(phi) in an AMR context, assuming no coarser AMR levels.
Definition CD_EBHelmholtzOp.cpp:1023
EBHelmholtzOp::m_ghostRhs
IntVect m_ghostRhs
Ghost cells for rhs (note, the operator rhs)
Definition CD_EBHelmholtzOp.H:956
EBHelmholtzOp::refToCoarser
int refToCoarser() override final
Return coarsening factor to coarser level (1 if there is no coarser level);.
Definition CD_EBHelmholtzOp.cpp:793
EBHelmholtzOp::m_centroidFluxStencil
LayoutData< BaseIFFAB< VoFStencil > > m_centroidFluxStencil[SpaceDim]
Face centroid flux stencil. Defined on all faces connecting one or more irregular vofs.
Definition CD_EBHelmholtzOp.H:1111
EBHelmholtzOp::operator=
EBHelmholtzOp & operator=(const EBHelmholtzOp &a_oper)=delete
No copy assigment allowed.
EBHelmholtzOp::m_interval
Interval m_interval
Interval.
Definition CD_EBHelmholtzOp.H:878
EBHelmholtzOp::applyOpNoBoundary
void applyOpNoBoundary(LevelData< EBCellFAB > &a_ans, const LevelData< EBCellFAB > &a_phi) override final
Apply operator but turn off all BCs.
Definition CD_EBHelmholtzOp.cpp:1539
EBHelmholtzOp::assignLocal
void assignLocal(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Local assignment function.
Definition CD_EBHelmholtzOp.cpp:560
EBHelmholtzOp::m_doExchange
bool m_doExchange
Turn on/off exchange operation.
Definition CD_EBHelmholtzOp.H:940
EBHelmholtzOp::m_prolongOpMG
EBMGProlong m_prolongOpMG
Prolongation operator for MG levels.
Definition CD_EBHelmholtzOp.H:1041
EBHelmholtzOp::norm
Real norm(const LevelData< EBCellFAB > &a_rhs, const int a_order) override final
Compute norm of data.
Definition CD_EBHelmholtzOp.cpp:684
EBHelmholtzOp::AMRResidual
void AMRResidual(LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousPhysBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) override final
Compute residual on this level. AMR version.
Definition CD_EBHelmholtzOp.cpp:1036
EBHelmholtzOp::m_interpStencil
LayoutData< BaseIFFAB< FaceStencil > > m_interpStencil[SpaceDim]
Face centroid interpolation stencil.
Definition CD_EBHelmholtzOp.H:1106
EBHelmholtzOp::residual
void residual(LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, const bool a_homogeneousPhysBc) override
Compute residual on this level.
Definition CD_EBHelmholtzOp.cpp:512
EBHelmholtzOp::fillDomainFlux
void fillDomainFlux(EBFluxFAB &a_flux, const EBCellFAB &a_phi, const Box &a_cellBox, const DataIndex &a_dit)
Fill domain flux. This fills the flux on the domain face using centered differencing ala applyDomainF...
Definition CD_EBHelmholtzOp.cpp:1365
EBHelmholtzOp::EBHelmholtzOp
EBHelmholtzOp(const EBHelmholtzOp &&a_other)=delete
Disallowed move constructor.
EBHelmholtzOp::m_Bcoef
RefCountedPtr< LevelData< EBFluxFAB > > m_Bcoef
B-coefficient in Helmholtz equation.
Definition CD_EBHelmholtzOp.H:1081
EBHelmholtzOp::deallocateFlux
void deallocateFlux() const noexcept
Deallocate m_flux.
Definition CD_EBHelmholtzOp.cpp:260
EBHelmholtzOp::m_timer
Timer m_timer
Timer so user can profile.
Definition CD_EBHelmholtzOp.H:873
EBHelmholtzOp::m_dataLocation
Location::Cell m_dataLocation
Data centering.
Definition CD_EBHelmholtzOp.H:888
EBHelmholtzOp::createCoarser
void createCoarser(LevelData< EBCellFAB > &a_coarse, const LevelData< EBCellFAB > &a_fine, bool a_ghosted) override final
Create coarsened data.
Definition CD_EBHelmholtzOp.cpp:746
EBHelmholtzOp::m_eblgCoar
EBLevelGrid m_eblgCoar
Coarse level grid (if the operator has a coarse level)
Definition CD_EBHelmholtzOp.H:1016
EBHelmholtzOp::m_ghostPhi
IntVect m_ghostPhi
Ghost cells for phi.
Definition CD_EBHelmholtzOp.H:951
EBHelmholtzOp::computeFlux
void computeFlux(EBFaceFAB &a_fluxCentroid, const EBCellFAB &a_phi, const Box &a_cellBox, const DataIndex &a_dit, const int a_dir)
Compute centroid fluxes in a direction.
Definition CD_EBHelmholtzOp.cpp:2203
EBHelmholtzOp::m_eblgFine
EBLevelGrid m_eblgFine
Fine level grid (if the operator has a fine level)
Definition CD_EBHelmholtzOp.H:1001
EBHelmholtzOp::turnOffCFInterp
void turnOffCFInterp()
Turn off BCs.
Definition CD_EBHelmholtzOp.cpp:204
EBHelmholtzOp::m_relaxFactor
Real m_relaxFactor
Successive over-relaxation factor (Must be >= 1).
Definition CD_EBHelmholtzOp.H:976
EBHelmholtzOp::m_refToCoar
int m_refToCoar
Refinement factor to coarse level.
Definition CD_EBHelmholtzOp.H:918
EBHelmholtzOp::m_alpha
Real m_alpha
Alpha-coefficient.
Definition CD_EBHelmholtzOp.H:961
EBHelmholtzOp::m_relCoef
LevelData< EBCellFAB > m_relCoef
Relaxation coefficient.
Definition CD_EBHelmholtzOp.H:1091
EBHelmholtzOp::diagonalScale
void diagonalScale(LevelData< EBCellFAB > &a_rhs, bool a_kappaWeighted) override final
Do diagonal scaling.
Definition CD_EBHelmholtzOp.cpp:1501
EBHelmholtzOp::setToZero
void setToZero(LevelData< EBCellFAB > &a_lhs) override final
Set data to zero.
Definition CD_EBHelmholtzOp.cpp:738
EBHelmholtzOp::setAlphaAndBeta
void setAlphaAndBeta(const Real &a_alpha, const Real &a_beta) override final
Set alpha coefficient and beta coefficient (can change as diffusion solvers progress)
Definition CD_EBHelmholtzOp.cpp:498
EBHelmholtzOp::m_validCells
RefCountedPtr< LevelData< BaseFab< bool > > > m_validCells
Valid grid cells (might be nullptr on MG levels)
Definition CD_EBHelmholtzOp.H:1046
EBHelmholtzOp::m_exchangeCopier
Copier m_exchangeCopier
Pre-built exchange copier.
Definition CD_EBHelmholtzOp.H:991
EBHelmholtzOp::interpolateCF
void interpolateCF(LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > *a_phiCoar, const bool a_homogeneousCFBC)
Apply coarse-fine boundary conditions.
Definition CD_EBHelmholtzOp.cpp:1593
EBHelmholtzOp::coarsenCell
void coarsenCell(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiFine)
Coarsen data from fine to coar level.
Definition CD_EBHelmholtzOp.cpp:2342
EBHelmholtzOp::gauSaiRedBlackKernel
void gauSaiRedBlackKernel(EBCellFAB &a_Lcorr, EBCellFAB &a_corr, const EBCellFAB &a_resid, const EBCellFAB &a_Acoef, const EBFluxFAB &a_Bcoef, const BaseIVFAB< Real > &a_BcoefIrreg, const Box &a_cellBox, const DataIndex &a_dit, const int &a_redBlack) const noexcept
Red-black Gauss-Seidel kernel.
Definition CD_EBHelmholtzOp.cpp:1767
EBHelmholtzOp::m_interpolant
LayoutData< BaseIFFAB< Real > > m_interpolant[SpaceDim]
Interpolant for when we want centroid fluxes.
Definition CD_EBHelmholtzOp.H:1101
EBHelmholtzOp::turnOnExchange
void turnOnExchange()
Turn on exchange operation.
Definition CD_EBHelmholtzOp.cpp:228
EBHelmholtzOp::inhomogeneousCFInterp
void inhomogeneousCFInterp(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar)
Inhomogeneous coarse-fine interpolation.
Definition CD_EBHelmholtzOp.cpp:1581
EBHelmholtzOp::coarsenFlux
void coarsenFlux(LevelData< EBFluxFAB > &a_flux, const LevelData< EBFluxFAB > &a_fineFlux)
Coarsen fluxes on the fine level onto this level.
Definition CD_EBHelmholtzOp.cpp:2350
EBHelmholtzOp::m_hasFine
bool m_hasFine
True if there's a finer level.
Definition CD_EBHelmholtzOp.H:908
EBHelmholtzOp::m_doCoarsen
bool m_doCoarsen
Turn on/off exchange operation.
Definition CD_EBHelmholtzOp.H:946
EBHelmholtzOp::assignCopier
void assignCopier(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs, const Copier &a_copier) override final
Assign lhs.
Definition CD_EBHelmholtzOp.cpp:552
EBHelmholtzOp::gauSaiMultiColorKernel
void gauSaiMultiColorKernel(EBCellFAB &a_Lcorr, EBCellFAB &a_corr, const EBCellFAB &a_resid, const EBCellFAB &a_Acoef, const EBFluxFAB &a_Bcoef, const BaseIVFAB< Real > &a_BcoefIrreg, const Box &a_cellBox, const DataIndex &a_dit, const IntVect &a_color) const noexcept
Multi-color Gauss-Seidel kernel.
Definition CD_EBHelmholtzOp.cpp:1875
EBHelmholtzOp::m_sideBox
std::map< std::pair< int, Side::LoHiSide >, Box > m_sideBox
Domain boxes on each side.
Definition CD_EBHelmholtzOp.H:1164
EBHelmholtzOp::m_vofIterStenc
LayoutData< VoFIterator > m_vofIterStenc
VoFIterator which iterates over all cells that are 1) a cut-cell or 2) borders a cut-cell.
Definition CD_EBHelmholtzOp.H:1149
EBHelmholtzOp::scaleLocal
void scaleLocal(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) const noexcept
Local scaling function. Multiplies the left-hand side by the right-hand side.
Definition CD_EBHelmholtzOp.cpp:663
EBHelmholtzOp::getAcoef
const RefCountedPtr< LevelData< EBCellFAB > > & getAcoef()
Get the Helmholtz A-coefficient on cell centers.
Definition CD_EBHelmholtzOp.cpp:181
EBHelmholtzOp::m_vecDx
RealVect m_vecDx
Vector resolution.
Definition CD_EBHelmholtzOp.H:981
EBHelmholtzOp::incr
void incr(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs, const Real a_scale) override final
Increment operator.
Definition CD_EBHelmholtzOp.cpp:635
EBHelmholtzOp::m_vofIterMulti
LayoutData< VoFIterator > m_vofIterMulti
VoFIterator for "multi-cells".
Definition CD_EBHelmholtzOp.H:1144
EBHelmholtzOp::m_betaDiagWeight
LayoutData< BaseIVFAB< Real > > m_betaDiagWeight
Weights of diagonal beta terms.
Definition CD_EBHelmholtzOp.H:1134
EBHelmholtzOp::EBHelmholtzOp
EBHelmholtzOp(const EBHelmholtzOp &a_other)=delete
Disallowed copy constructor.
EBHelmholtzOp::m_colors
Vector< IntVect > m_colors
"Colors" for the relaxation methods
Definition CD_EBHelmholtzOp.H:1169
EBHelmholtzOp::getBcoef
const RefCountedPtr< LevelData< EBFluxFAB > > & getBcoef()
Get the Helmholtz B-coefficient on faces.
Definition CD_EBHelmholtzOp.cpp:187
EBHelmholtzOp::m_Acoef
RefCountedPtr< LevelData< EBCellFAB > > m_Acoef
A-coefficient in Helmholtz equation.
Definition CD_EBHelmholtzOp.H:1076
EBHelmholtzOp::turnOnCoarsening
void turnOnCoarsening()
Turn on coarsening operation.
Definition CD_EBHelmholtzOp.cpp:244
EBHelmholtzOp::m_vofIterDomLo
LayoutData< VoFIterator > m_vofIterDomLo[SpaceDim]
VoF iterators for lo domain side.
Definition CD_EBHelmholtzOp.H:1154
EBHelmholtzOp::m_vofIterIrreg
LayoutData< VoFIterator > m_vofIterIrreg
VoFIterator for irregular cells.
Definition CD_EBHelmholtzOp.H:1139
EBHelmholtzOp::AMROperator
void AMROperator(LevelData< EBCellFAB > &a_Lphi, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar, const bool a_homogeneousPhysBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) override final
Apply the AMR operator, i.e. compute L(phi) in an AMR context.
Definition CD_EBHelmholtzOp.cpp:799
EBHelmholtzOp::AMRRestrict
void AMRRestrict(LevelData< EBCellFAB > &a_residualCoarse, const LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, bool a_skip_res) override final
Restrict residual.
Definition CD_EBHelmholtzOp.cpp:1086
EBHelmholtzOp::m_restrictOpMG
EBMGRestrict m_restrictOpMG
Restriction operator if this is an MG level.
Definition CD_EBHelmholtzOp.H:1031
EBHelmholtzOp::createCoarsened
void createCoarsened(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs, const int &a_refRat) override final
Create coarsening of data holder.
Definition CD_EBHelmholtzOp.cpp:756
EBHelmholtzOp::m_beta
Real m_beta
Beta-coefficient.
Definition CD_EBHelmholtzOp.H:966
EBHelmholtzOp::operator=
EBHelmholtzOp & operator=(const EBHelmholtzOp &&a_oper)=delete
No move assigment allowed.
EBHelmholtzOp::applyOpIrregular
void applyOpIrregular(EBCellFAB &a_Lphi, const EBCellFAB &a_phi, const EBCellFAB &a_Acoef, const EBFluxFAB &a_Bcoef, const BaseIVFAB< Real > &a_BcoefIrreg, const BaseIVFAB< Real > &a_alphaDiagWeight, const Box &a_cellBox, const DataIndex &a_dit, const bool a_homogeneousPhysBC) const noexcept
Apply operator in irregular cells.
Definition CD_EBHelmholtzOp.cpp:1419
EBHelmholtzOp::AMRUpdateResidual
void AMRUpdateResidual(LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection) override final
Update AMR residual.
Definition CD_EBHelmholtzOp.cpp:1122
EBHelmholtzOp::m_vofIterDomHi
LayoutData< VoFIterator > m_vofIterDomHi[SpaceDim]
VoF iterators for hi domain side.
Definition CD_EBHelmholtzOp.H:1159
EBHelmholtzOp::buildCopier
void buildCopier(Copier &a_copier, const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override
Build copier.
Definition CD_EBHelmholtzOp.cpp:2358
EBHelmholtzOp::defineStencils
void defineStencils()
Define stencils.
Definition CD_EBHelmholtzOp.cpp:284
EBHelmholtzOp::prolongIncrement
void prolongIncrement(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_correctCoarse) override final
Prolongation method.
Definition CD_EBHelmholtzOp.cpp:785
EBHelmholtzOp::getRelaxationCoeff
const LevelData< EBCellFAB > & getRelaxationCoeff() const noexcept
Returns m_relCoef. Used by MFHelmholtzOp.
Definition CD_EBHelmholtzOp.cpp:276
EBHelmholtzOp::m_eblgCoarMG
EBLevelGrid m_eblgCoarMG
Coarser grids (multigrid level)
Definition CD_EBHelmholtzOp.H:1021
EBHelmholtzOp::m_BcoefIrreg
RefCountedPtr< LevelData< BaseIVFAB< Real > > > m_BcoefIrreg
B-coefficient in Helmholtz equation, but on EB faces.
Definition CD_EBHelmholtzOp.H:1086
EBMGProlong
Class for prolongation of multigrid residual onto a finer grid.
Definition CD_EBMGProlong.H:30
EBMGRestrict
Class for restricting multigrid residual onto a coarser grid.
Definition CD_EBMGRestrict.H:30
EBMultigridInterpolator
Multigrid interpolator class.
Definition CD_EBMGLeastSquaresInterpolator.H:48
EBReflux
Class which can do refluxing across a coarse-fine interface.
Definition CD_EBReflux.H:37
Timer
Class which is used for run-time monitoring of events.
Definition CD_Timer.H:31
TracerParticleSolver
Base class for a tracer particle solver. This solver can advance particles in a pre-defined velocity ...
Definition CD_TracerParticleSolver.H:37
Location
Namespace for encapsulating various data centerings.
Definition CD_Location.H:24
Location::Cell
Cell
Enum for distinguishing between cell locations.
Definition CD_Location.H:30
Generated by doxygen 1.9.8