chombo-discharge
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 
40 class EBHelmholtzOp : public LevelTGAHelmOp<LevelData<EBCellFAB>, EBFluxFAB>
41 {
42 public:
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 
100  EBHelmholtzOp(const Location::Cell a_dataLocation,
101  const EBLevelGrid& a_eblgFine,
102  const EBLevelGrid& a_eblg,
103  const EBLevelGrid& a_eblgCoFi,
104  const EBLevelGrid& a_eblgCoar,
105  const EBLevelGrid& a_eblgCoarMG,
106  const RefCountedPtr<EBMultigridInterpolator>& a_interpolator,
107  const RefCountedPtr<EBReflux>& a_fluxReg,
108  const RefCountedPtr<EBCoarAve>& a_coarAve,
109  const RefCountedPtr<EBHelmholtzDomainBC>& a_domainBC,
110  const RefCountedPtr<EBHelmholtzEBBC>& a_ebBC,
111  const RealVect& a_probLo,
112  const Real& a_dx,
113  const int& a_refToFine,
114  const int& a_refToCoar,
115  const bool& a_hasFine,
116  const bool& a_hasCoar,
117  const bool& a_hasMGObjects,
118  const Real& a_alpha,
119  const Real& a_beta,
120  const RefCountedPtr<LevelData<EBCellFAB>>& a_Acoef,
121  const RefCountedPtr<LevelData<EBFluxFAB>>& a_Bcoef,
122  const RefCountedPtr<LevelData<BaseIVFAB<Real>>>& a_BcoIrreg,
123  const IntVect& a_ghostCellsPhi,
124  const IntVect& a_ghostCellsRHS,
125  const Smoother& a_smoother);
126 
130  virtual ~EBHelmholtzOp();
131 
136  EBHelmholtzOp&
137  operator=(const EBHelmholtzOp& a_oper) = delete;
138 
143  EBHelmholtzOp&
144  operator=(const EBHelmholtzOp&& a_oper) = delete;
145 
149  void
150  turnOffCFInterp();
151 
155  void
156  turnOnCFInterp();
157 
162  void
163  turnOffExchange();
164 
169  void
170  turnOnExchange();
171 
176  void
177  turnOffCoarsening();
178 
183  void
184  turnOnCoarsening();
185 
193  void
194  setAcoAndBco(const RefCountedPtr<LevelData<EBCellFAB>>& a_Acoef,
195  const RefCountedPtr<LevelData<EBFluxFAB>>& a_Bcoef,
196  const RefCountedPtr<LevelData<BaseIVFAB<Real>>>& a_BcoefIrreg);
197 
202  const RefCountedPtr<LevelData<EBCellFAB>>&
203  getAcoef();
204 
209  const RefCountedPtr<LevelData<EBFluxFAB>>&
210  getBcoef();
211 
216  const RefCountedPtr<LevelData<BaseIVFAB<Real>>>&
217  getBcoefIrreg();
218 
225  void
226  coarsenCell(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_phiFine);
227 
235  void
236  coarsenFlux(LevelData<EBFluxFAB>& a_flux, const LevelData<EBFluxFAB>& a_fineFlux);
237 
249  void
250  pointJacobiKernel(EBCellFAB& a_Lcorr,
251  EBCellFAB& a_corr,
252  const EBCellFAB& a_resid,
253  const EBCellFAB& a_Acoef,
254  const EBFluxFAB& a_Bcoef,
255  const BaseIVFAB<Real>& a_BcoefIrreg,
256  const Box& a_cellBox,
257  const DataIndex& a_dit) const noexcept;
258 
271  void
272  gauSaiRedBlackKernel(EBCellFAB& a_Lcorr,
273  EBCellFAB& a_corr,
274  const EBCellFAB& a_resid,
275  const EBCellFAB& a_Acoef,
276  const EBFluxFAB& a_Bcoef,
277  const BaseIVFAB<Real>& a_BcoefIrreg,
278  const Box& a_cellBox,
279  const DataIndex& a_dit,
280  const int& a_redBlack) const noexcept;
281 
294  void
295  gauSaiMultiColorKernel(EBCellFAB& a_Lcorr,
296  EBCellFAB& a_corr,
297  const EBCellFAB& a_resid,
298  const EBCellFAB& a_Acoef,
299  const EBFluxFAB& a_Bcoef,
300  const BaseIVFAB<Real>& a_BcoefIrreg,
301  const Box& a_cellBox,
302  const DataIndex& a_dit,
303  const IntVect& a_color) const noexcept;
304 
312  void
313  residual(LevelData<EBCellFAB>& a_residual,
314  const LevelData<EBCellFAB>& a_phi,
315  const LevelData<EBCellFAB>& a_rhs,
316  const bool a_homogeneousPhysBc) override;
317 
324  void
325  preCond(LevelData<EBCellFAB>& a_corr, const LevelData<EBCellFAB>& a_residual) override final;
326 
333  void
334  interpolateCF(LevelData<EBCellFAB>& a_phiFine, const LevelData<EBCellFAB>* a_phiCoar, const bool a_homogeneousCFBC);
335 
340  void
341  homogeneousCFInterp(LevelData<EBCellFAB>& a_phi);
342 
348  void
349  inhomogeneousCFInterp(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_phiCoar);
350 
358  void
359  gauSaiMultiColor(LevelData<EBCellFAB>& a_phi,
360  const LevelData<EBCellFAB>& a_Lphi,
361  const LevelData<EBCellFAB>& a_rhs,
362  const IntVect& a_color) const;
363 
373  void
374  applyOp(LevelData<EBCellFAB>& a_Lphi,
375  const LevelData<EBCellFAB>& a_phi,
376  const LevelData<EBCellFAB>* const a_phiCoar,
377  const bool a_homogeneousPhysBC,
378  const bool a_homogeneousCFBC);
379 
387  void
388  applyOp(LevelData<EBCellFAB>& a_Lphi, const LevelData<EBCellFAB>& a_phi, bool a_homogeneousPhysBc) override final;
389 
398  void
399  applyOp(EBCellFAB& a_Lphi,
400  EBCellFAB& a_phi,
401  const EBCellFAB& a_Acoef,
402  const EBFluxFAB& a_Bcoef,
403  const BaseIVFAB<Real>& a_BcoefIrreg,
404  const Box& a_cellBox,
405  const DataIndex& a_dit,
406  const bool a_homogeneousPhysBC) const noexcept;
407 
416  void
417  applyOpRegular(EBCellFAB& a_Lphi,
418  EBCellFAB& a_phi,
419  const EBCellFAB& a_Acoef,
420  const EBFluxFAB& a_Bcoef,
421  const BaseIVFAB<Real>& a_BcoefIrreg,
422  const Box& a_cellBox,
423  const DataIndex& a_dit,
424  const bool a_homogeneousPhysBC) const noexcept;
425 
436  void
437  applyDomainFlux(EBCellFAB& a_phi,
438  const EBFluxFAB& a_Bcoef,
439  const Box& a_cellBox,
440  const DataIndex& a_dit,
441  const bool a_homogeneousPhysBc) const noexcept;
442 
451  void
452  fillDomainFlux(EBFluxFAB& a_flux, const EBCellFAB& a_phi, const Box& a_cellBox, const DataIndex& a_dit);
453 
462  void
463  applyOpIrregular(EBCellFAB& a_Lphi,
464  const EBCellFAB& a_phi,
465  const EBCellFAB& a_Acoef,
466  const EBFluxFAB& a_Bcoef,
467  const BaseIVFAB<Real>& a_BcoefIrreg,
468  const BaseIVFAB<Real>& a_alphaDiagWeight,
469  const Box& a_cellBox,
470  const DataIndex& a_dit,
471  const bool a_homogeneousPhysBC) const noexcept;
472 
478  void
479  create(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
480 
485  void
486  assign(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
487 
496  void
497  assignCopier(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const Copier& a_copier) override final;
498 
504  void
505  assignLocal(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
506 
513  void
514  buildCopier(Copier& a_copier, const LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override;
515 
519  Real
520  dotProduct(const LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs) override final;
521 
528  void
529  incr(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const Real a_scale) override final;
530 
539  void
540  axby(LevelData<EBCellFAB>& a_lhs,
541  const LevelData<EBCellFAB>& a_x,
542  const LevelData<EBCellFAB>& a_y,
543  const Real a_a,
544  const Real a_b) override final;
545 
551  void
552  scale(LevelData<EBCellFAB>& a_lhs, const Real& a_scale) override final;
553 
559  Real
560  norm(const LevelData<EBCellFAB>& a_rhs, const int a_order) override final;
561 
566  void
567  setToZero(LevelData<EBCellFAB>& a_lhs) override final;
568 
575  void
576  createCoarser(LevelData<EBCellFAB>& a_coarse, const LevelData<EBCellFAB>& a_fine, bool a_ghosted) override final;
577 
584  void
585  relax(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, int a_iterations) override final;
586 
593  void
594  restrictResidual(LevelData<EBCellFAB>& a_resCoar,
595  LevelData<EBCellFAB>& a_phi,
596  const LevelData<EBCellFAB>& a_rhs) override final;
597 
603  void
604  prolongIncrement(LevelData<EBCellFAB>& a_phi, const LevelData<EBCellFAB>& a_correctCoarse) override final;
605 
609  int
610  refToCoarser() override final;
611 
622  void
623  AMROperator(LevelData<EBCellFAB>& a_Lphi,
624  const LevelData<EBCellFAB>& a_phiFine,
625  const LevelData<EBCellFAB>& a_phi,
626  const LevelData<EBCellFAB>& a_phiCoar,
627  const bool a_homogeneousPhysBC,
628  AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
629 
641  void
642  refluxFreeAMROperator(LevelData<EBCellFAB>& a_Lphi,
643  const LevelData<EBCellFAB>& a_phiFine,
644  const LevelData<EBCellFAB>& a_phi,
645  const LevelData<EBCellFAB>& a_phiCoar,
646  const bool a_homogeneousPhysBC,
647  AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp);
648 
659  void
660  AMRResidual(LevelData<EBCellFAB>& a_residual,
661  const LevelData<EBCellFAB>& a_phiFine,
662  const LevelData<EBCellFAB>& a_phi,
663  const LevelData<EBCellFAB>& a_phiCoar,
664  const LevelData<EBCellFAB>& a_rhs,
665  bool a_homogeneousPhysBC,
666  AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
667 
676  void
677  AMRResidualNF(LevelData<EBCellFAB>& a_residual,
678  const LevelData<EBCellFAB>& a_phi,
679  const LevelData<EBCellFAB>& a_phiCoar,
680  const LevelData<EBCellFAB>& a_rhs,
681  bool a_homogeneousPhysBC) override final;
682 
692  void
693  AMRResidualNC(LevelData<EBCellFAB>& a_residual,
694  const LevelData<EBCellFAB>& a_phiFine,
695  const LevelData<EBCellFAB>& a_phi,
696  const LevelData<EBCellFAB>& a_rhs,
697  bool a_homogeneousPhysBC,
698  AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
699 
708  void
709  AMROperatorNF(LevelData<EBCellFAB>& a_Lphi,
710  const LevelData<EBCellFAB>& a_phi,
711  const LevelData<EBCellFAB>& a_phiCoar,
712  bool a_homogeneousPhysBC) override final;
713 
722  void
723  AMROperatorNC(LevelData<EBCellFAB>& a_Lphi,
724  const LevelData<EBCellFAB>& a_phi,
725  const LevelData<EBCellFAB>& a_phiCoar,
726  bool a_homogeneousPhysBC,
727  AMRLevelOp<LevelData<EBCellFAB>>* a_finerOp) override final;
728 
737  void
738  AMRRestrict(LevelData<EBCellFAB>& a_residualCoarse,
739  const LevelData<EBCellFAB>& a_residual,
740  const LevelData<EBCellFAB>& a_correction,
741  const LevelData<EBCellFAB>& a_coarseCorrection,
742  bool a_skip_res) override final;
743 
749  void
750  AMRProlong(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_coarseCorrection) override final;
751 
758  void
759  AMRUpdateResidual(LevelData<EBCellFAB>& a_residual,
760  const LevelData<EBCellFAB>& a_correction,
761  const LevelData<EBCellFAB>& a_coarseCorrection) override final;
762 
768  void
769  setAlphaAndBeta(const Real& a_alpha, const Real& a_beta) override final;
770 
777  void
778  createCoarsened(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const int& a_refRat) override final;
779 
785  void
786  diagonalScale(LevelData<EBCellFAB>& a_rhs, bool a_kappaWeighted) override final;
787 
792  void
793  divideByIdentityCoef(LevelData<EBCellFAB>& a_rhs) override final;
794 
798  void
799  applyOpNoBoundary(LevelData<EBCellFAB>& a_ans, const LevelData<EBCellFAB>& a_phi) override final;
800 
805  void
806  fillGrad(const LevelData<EBCellFAB>& a_phi) override final;
807 
816  void
817  getFlux(EBFluxFAB& a_flux,
818  const LevelData<EBCellFAB>& a_data,
819  const Box& a_grid,
820  const DataIndex& a_dit,
821  Real a_scale) override final;
822 
826  void
827  allocateFlux() const noexcept;
828 
832  void
833  deallocateFlux() const noexcept;
834 
838  LevelData<EBFluxFAB>&
839  getFlux() const;
840 
841 protected:
845  static constexpr int m_nComp = 1;
846 
850  static constexpr int m_comp = 0;
851 
855  Timer m_timer;
856 
860  Interval m_interval;
861 
865  Smoother m_smoother;
866 
870  Location::Cell m_dataLocation;
871 
875  bool m_refluxFree;
876 
880  bool m_profile;
881 
885  bool m_hasMGObjects;
886 
890  bool m_hasFine;
891 
895  bool m_hasCoar;
896 
900  int m_refToCoar;
901 
905  int m_refToFine;
906 
911  bool m_doInterpCF;
912 
917  bool m_doExchange;
918 
923  bool m_doCoarsen;
924 
928  IntVect m_ghostPhi;
929 
933  IntVect m_ghostRhs;
934 
938  Real m_alpha;
939 
943  Real m_beta;
944 
948  Real m_dx;
949 
953  RealVect m_vecDx;
954 
958  RealVect m_probLo;
959 
963  Copier m_exchangeCopier;
964 
968  Copier m_exchangeCopierFine;
969 
973  EBLevelGrid m_eblgFine;
974 
978  EBLevelGrid m_eblg;
979 
983  EBLevelGrid m_eblgCoFi;
984 
988  EBLevelGrid m_eblgCoar;
989 
993  EBLevelGrid m_eblgCoarMG;
994 
998  EBMGRestrict m_restrictOp;
999 
1003  EBMGRestrict m_restrictOpMG;
1004 
1008  EBMGProlong m_prolongOp;
1009 
1013  EBMGProlong m_prolongOpMG;
1014 
1018  RefCountedPtr<EBHelmholtzDomainBC> m_domainBc;
1019 
1023  RefCountedPtr<EBHelmholtzEBBC> m_ebBc;
1024 
1028  RefCountedPtr<EBMultigridInterpolator> m_interpolator;
1029 
1033  RefCountedPtr<EBReflux> m_fluxReg;
1034 
1038  RefCountedPtr<EBCoarAve> m_coarAve;
1039 
1043  RefCountedPtr<LevelData<EBCellFAB>> m_Acoef;
1044 
1048  RefCountedPtr<LevelData<EBFluxFAB>> m_Bcoef;
1049 
1053  RefCountedPtr<LevelData<BaseIVFAB<Real>>> m_BcoefIrreg;
1054 
1058  LevelData<EBCellFAB> m_relCoef;
1059 
1063  mutable LevelData<EBFluxFAB>* m_flux;
1064 
1068  LayoutData<BaseIFFAB<Real>> m_interpolant[SpaceDim];
1069 
1073  LayoutData<BaseIFFAB<FaceStencil>> m_interpStencil[SpaceDim];
1074 
1078  LayoutData<BaseIFFAB<VoFStencil>> m_centroidFluxStencil[SpaceDim];
1079 
1085  LayoutData<BaseIVFAB<VoFStencil>> m_relaxStencils;
1086 
1091  LayoutData<RefCountedPtr<VCAggStencil>> m_aggRelaxStencil;
1092 
1096  LayoutData<BaseIVFAB<Real>> m_alphaDiagWeight;
1097 
1101  LayoutData<BaseIVFAB<Real>> m_betaDiagWeight;
1102 
1106  mutable LayoutData<VoFIterator> m_vofIterIrreg;
1107 
1111  mutable LayoutData<VoFIterator> m_vofIterMulti;
1112 
1116  mutable LayoutData<VoFIterator> m_vofIterStenc;
1117 
1121  mutable LayoutData<VoFIterator> m_vofIterDomLo[SpaceDim];
1122 
1126  mutable LayoutData<VoFIterator> m_vofIterDomHi[SpaceDim];
1127 
1131  std::map<std::pair<int, Side::LoHiSide>, Box> m_sideBox;
1132 
1136  Vector<IntVect> m_colors;
1137 
1144  void
1145  relaxPointJacobi(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1146 
1153  void
1154  relaxGSRedBlack(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1155 
1162  void
1163  relaxGSMultiColor(LevelData<EBCellFAB>& a_correction, const LevelData<EBCellFAB>& a_residual, const int a_iterations);
1164 
1168  void
1169  computeDiagWeight();
1170 
1174  void
1175  computeRelaxationCoefficient();
1176 
1180  void
1181  makeAggStencil();
1182 
1186  void
1187  defineStencils();
1188 
1196  VoFStencil
1197  getFaceCenterFluxStencil(const FaceIndex& a_face, const DataIndex& a_dit) const;
1198 
1206  VoFStencil
1207  getFaceCentroidFluxStencil(const FaceIndex& a_face, const DataIndex& a_dit) const;
1208 
1217  void
1218  computeFlux(EBFaceFAB& a_fluxCentroid,
1219  const EBCellFAB& a_phi,
1220  const Box& a_cellBox,
1221  const DataIndex& a_dit,
1222  const int a_dir);
1223 
1232  void
1233  computeFaceCenteredFlux(EBFaceFAB& a_fluxCenter,
1234  const EBCellFAB& a_phi,
1235  const Box& a_cellBox,
1236  const DataIndex& a_dit,
1237  const int a_dir);
1238 
1247  void
1248  computeFaceCentroidFlux(EBFaceFAB& a_flux,
1249  const EBCellFAB& a_phi,
1250  const Box& a_cellBox,
1251  const DataIndex& a_dit,
1252  const int a_dir);
1253 
1259  void
1260  computeFlux(const LevelData<EBCellFAB>& a_phi);
1261 
1269  void
1270  reflux(LevelData<EBCellFAB>& a_Lphi,
1271  const LevelData<EBCellFAB>& a_phiFine,
1272  const LevelData<EBCellFAB>& a_phi,
1273  AMRLevelOp<LevelData<EBCellFAB>>& a_finerOp);
1274 };
1275 
1276 #include <CD_NamespaceFooter.H>
1277 
1278 #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:2091
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:1789
EBHelmholtzOp::getFlux
LevelData< EBFluxFAB > & getFlux() const
Returns m_flux. This is used in refluxing routines.
Definition: CD_EBHelmholtzOp.cpp:261
EBHelmholtzOp::m_nComp
static constexpr int m_nComp
Number of components that we solve for (always one..)
Definition: CD_EBHelmholtzOp.H:845
EBHelmholtzOp::m_coarAve
RefCountedPtr< EBCoarAve > m_coarAve
Conservative coarsener.
Definition: CD_EBHelmholtzOp.H:1038
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:158
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:625
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:518
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:850
EBHelmholtzOp::turnOnCFInterp
void turnOnCFInterp()
Turn on BCs.
Definition: CD_EBHelmholtzOp.cpp:205
EBHelmholtzOp::m_eblgCoFi
EBLevelGrid m_eblgCoFi
Coarsened of m_eblg.
Definition: CD_EBHelmholtzOp.H:983
EBHelmholtzOp::getBcoefIrreg
const RefCountedPtr< LevelData< BaseIVFAB< Real > > > & getBcoefIrreg()
Get the Helmholtz B-coefficient on the EB.
Definition: CD_EBHelmholtzOp.cpp:186
EBHelmholtzOp::AMRProlong
void AMRProlong(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection) override final
Prolongation onto AMR level.
Definition: CD_EBHelmholtzOp.cpp:1075
EBHelmholtzOp::m_smoother
Smoother m_smoother
Relaxation method.
Definition: CD_EBHelmholtzOp.H:865
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:1185
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:1085
EBHelmholtzOp::operator=
EBHelmholtzOp & operator=(const EBHelmholtzOp &&a_oper)=delete
No move assigment allowed.
EBHelmholtzOp::m_hasCoar
bool m_hasCoar
True if there's a coarser level.
Definition: CD_EBHelmholtzOp.H:895
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:1015
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:2214
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:1239
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:1111
EBHelmholtzOp::allocateFlux
void allocateFlux() const noexcept
Allocate m_flux.
Definition: CD_EBHelmholtzOp.cpp:245
EBHelmholtzOp::computeDiagWeight
void computeDiagWeight()
Calculate the weight of the diagonal term.
Definition: CD_EBHelmholtzOp.cpp:1902
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:637
EBHelmholtzOp::EBHelmholtzOp
EBHelmholtzOp()=delete
Disallowed default constructor.
EBHelmholtzOp::~EBHelmholtzOp
virtual ~EBHelmholtzOp()
Dtor.
Definition: CD_EBHelmholtzOp.cpp:191
EBHelmholtzOp::m_exchangeCopierFine
Copier m_exchangeCopierFine
Pre-built exchange copier.
Definition: CD_EBHelmholtzOp.H:968
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:1581
EBHelmholtzOp::turnOffCoarsening
void turnOffCoarsening()
Turn off coarsening operation.
Definition: CD_EBHelmholtzOp.cpp:229
EBHelmholtzOp::relaxPointJacobi
void relaxPointJacobi(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, const int a_iterations)
Jacobi relaxation.
Definition: CD_EBHelmholtzOp.cpp:1615
EBHelmholtzOp::m_ebBc
RefCountedPtr< EBHelmholtzEBBC > m_ebBc
Domain bc object.
Definition: CD_EBHelmholtzOp.H:1023
EBHelmholtzOp::makeAggStencil
void makeAggStencil()
Compute aggregated stencils.
Definition: CD_EBHelmholtzOp.cpp:2021
EBHelmholtzOp::homogeneousCFInterp
void homogeneousCFInterp(LevelData< EBCellFAB > &a_phi)
Do homogeneous coarse-fine interpolation.
Definition: CD_EBHelmholtzOp.cpp:1530
EBHelmholtzOp::m_eblg
EBLevelGrid m_eblg
Grid.
Definition: CD_EBHelmholtzOp.H:978
EBHelmholtzOp::m_aggRelaxStencil
LayoutData< RefCountedPtr< VCAggStencil > > m_aggRelaxStencil
For making irregular stencil applications go faster.
Definition: CD_EBHelmholtzOp.H:1091
EBHelmholtzOp::m_restrictOp
EBMGRestrict m_restrictOp
Restriction operator for AMR levels.
Definition: CD_EBHelmholtzOp.H:998
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:801
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:1655
EBHelmholtzOp::m_refluxFree
bool m_refluxFree
Use reflux-free formulation or not.
Definition: CD_EBHelmholtzOp.H:875
EBHelmholtzOp::m_domainBc
RefCountedPtr< EBHelmholtzDomainBC > m_domainBc
Domain bc object.
Definition: CD_EBHelmholtzOp.H:1018
EBHelmholtzOp::dotProduct
Real dotProduct(const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Compute the dot product??
Definition: CD_EBHelmholtzOp.cpp:550
EBHelmholtzOp::m_fluxReg
RefCountedPtr< EBReflux > m_fluxReg
Flux register.
Definition: CD_EBHelmholtzOp.H:1033
EBHelmholtzOp::preCond
void preCond(LevelData< EBCellFAB > &a_corr, const LevelData< EBCellFAB > &a_residual) override final
Precondition system before bottom solve.
Definition: CD_EBHelmholtzOp.cpp:510
EBHelmholtzOp::m_doInterpCF
bool m_doInterpCF
Do coarse-fine interpolation or not.
Definition: CD_EBHelmholtzOp.H:911
EBHelmholtzOp::m_refToFine
int m_refToFine
Refinement factor to fine level.
Definition: CD_EBHelmholtzOp.H:905
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:1063
EBHelmholtzOp::m_hasMGObjects
bool m_hasMGObjects
True if there is a multigrid level below this operator.
Definition: CD_EBHelmholtzOp.H:885
EBHelmholtzOp::m_dx
Real m_dx
Grid resolution;.
Definition: CD_EBHelmholtzOp.H:948
EBHelmholtzOp::m_profile
bool m_profile
Profile the operator.
Definition: CD_EBHelmholtzOp.H:880
EBHelmholtzOp::turnOffExchange
void turnOffExchange()
Turn off exchange operation.
Definition: CD_EBHelmholtzOp.cpp:213
EBHelmholtzOp::relaxGSRedBlack
void relaxGSRedBlack(LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_residual, const int a_iterations)
Jacobi relaxation.
Definition: CD_EBHelmholtzOp.cpp:1681
EBHelmholtzOp::m_alphaDiagWeight
LayoutData< BaseIVFAB< Real > > m_alphaDiagWeight
Weights of diagonal alpha terms.
Definition: CD_EBHelmholtzOp.H:1096
EBHelmholtzOp::fillGrad
void fillGrad(const LevelData< EBCellFAB > &a_phi) override final
Not called, I think.
Definition: CD_EBHelmholtzOp.cpp:1510
EBHelmholtzOp::getFaceCenterFluxStencil
VoFStencil getFaceCenterFluxStencil(const FaceIndex &a_face, const DataIndex &a_dit) const
Get the face-centered flux stencil.
Definition: CD_EBHelmholtzOp.cpp:2072
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:729
EBHelmholtzOp::divideByIdentityCoef
void divideByIdentityCoef(LevelData< EBCellFAB > &a_rhs) override final
Divide by the a-coefficient.
Definition: CD_EBHelmholtzOp.cpp:1484
EBHelmholtzOp::assign
void assign(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Assign data.
Definition: CD_EBHelmholtzOp.cpp:526
EBHelmholtzOp::m_probLo
RealVect m_probLo
Lower-left corner of domain.
Definition: CD_EBHelmholtzOp.H:958
EBHelmholtzOp::m_interpolator
RefCountedPtr< EBMultigridInterpolator > m_interpolator
Interpolator object.
Definition: CD_EBHelmholtzOp.H:1028
EBHelmholtzOp::computeRelaxationCoefficient
void computeRelaxationCoefficient()
Calculate relaxation coefficient.
Definition: CD_EBHelmholtzOp.cpp:1957
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:2186
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:973
EBHelmholtzOp::m_prolongOp
EBMGProlong m_prolongOp
Prolongation operator for AMR levels.
Definition: CD_EBHelmholtzOp.H:1008
EBHelmholtzOp::m_ghostRhs
IntVect m_ghostRhs
Ghost cells for rhs (note, the operator rhs)
Definition: CD_EBHelmholtzOp.H:933
EBHelmholtzOp::refToCoarser
int refToCoarser() override final
Return coarsening factor to coarser level (1 if there is no coarser level);.
Definition: CD_EBHelmholtzOp.cpp:754
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:1078
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:760
EBHelmholtzOp::m_interval
Interval m_interval
Interval.
Definition: CD_EBHelmholtzOp.H:860
EBHelmholtzOp::operator=
EBHelmholtzOp & operator=(const EBHelmholtzOp &a_oper)=delete
No copy assigment allowed.
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:1500
EBHelmholtzOp::assignLocal
void assignLocal(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override final
Local assignment function.
Definition: CD_EBHelmholtzOp.cpp:542
EBHelmholtzOp::m_doExchange
bool m_doExchange
Turn on/off exchange operation.
Definition: CD_EBHelmholtzOp.H:917
EBHelmholtzOp::m_prolongOpMG
EBMGProlong m_prolongOpMG
Prolongation operator for MG levels.
Definition: CD_EBHelmholtzOp.H:1013
EBHelmholtzOp::norm
Real norm(const LevelData< EBCellFAB > &a_rhs, const int a_order) override final
Compute norm of data.
Definition: CD_EBHelmholtzOp.cpp:645
EBHelmholtzOp::m_interpStencil
LayoutData< BaseIFFAB< FaceStencil > > m_interpStencil[SpaceDim]
Face centroid interpolation stencil.
Definition: CD_EBHelmholtzOp.H:1073
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:497
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:1326
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:1048
EBHelmholtzOp::deallocateFlux
void deallocateFlux() const noexcept
Deallocate m_flux.
Definition: CD_EBHelmholtzOp.cpp:253
EBHelmholtzOp::m_timer
Timer m_timer
Timer so user can profile.
Definition: CD_EBHelmholtzOp.H:855
EBHelmholtzOp::m_dataLocation
Location::Cell m_dataLocation
Data centering.
Definition: CD_EBHelmholtzOp.H:870
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:707
EBHelmholtzOp::m_eblgCoar
EBLevelGrid m_eblgCoar
Coarse level grid (if the operator has a coarse level)
Definition: CD_EBHelmholtzOp.H:988
EBHelmholtzOp::m_ghostPhi
IntVect m_ghostPhi
Ghost cells for phi.
Definition: CD_EBHelmholtzOp.H:928
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:2167
EBHelmholtzOp::m_eblgFine
EBLevelGrid m_eblgFine
Fine level grid (if the operator has a fine level)
Definition: CD_EBHelmholtzOp.H:973
EBHelmholtzOp::turnOffCFInterp
void turnOffCFInterp()
Turn off BCs.
Definition: CD_EBHelmholtzOp.cpp:197
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:2268
EBHelmholtzOp::m_refToCoar
int m_refToCoar
Refinement factor to coarse level.
Definition: CD_EBHelmholtzOp.H:900
EBHelmholtzOp::m_alpha
Real m_alpha
Alpha-coefficient.
Definition: CD_EBHelmholtzOp.H:938
EBHelmholtzOp::m_relCoef
LevelData< EBCellFAB > m_relCoef
Relaxation coefficient.
Definition: CD_EBHelmholtzOp.H:1058
EBHelmholtzOp::diagonalScale
void diagonalScale(LevelData< EBCellFAB > &a_rhs, bool a_kappaWeighted) override final
Do diagonal scaling.
Definition: CD_EBHelmholtzOp.cpp:1462
EBHelmholtzOp::setToZero
void setToZero(LevelData< EBCellFAB > &a_lhs) override final
Set data to zero.
Definition: CD_EBHelmholtzOp.cpp:699
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:483
EBHelmholtzOp::m_exchangeCopier
Copier m_exchangeCopier
Pre-built exchange copier.
Definition: CD_EBHelmholtzOp.H:963
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:1554
EBHelmholtzOp::coarsenCell
void coarsenCell(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiFine)
Coarsen data from fine to coar level.
Definition: CD_EBHelmholtzOp.cpp:2306
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:1728
EBHelmholtzOp::m_interpolant
LayoutData< BaseIFFAB< Real > > m_interpolant[SpaceDim]
Interpolant for when we want centroid fluxes.
Definition: CD_EBHelmholtzOp.H:1068
EBHelmholtzOp::turnOnExchange
void turnOnExchange()
Turn on exchange operation.
Definition: CD_EBHelmholtzOp.cpp:221
EBHelmholtzOp::inhomogeneousCFInterp
void inhomogeneousCFInterp(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoar)
Inhomogeneous coarse-fine interpolation.
Definition: CD_EBHelmholtzOp.cpp:1542
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:2314
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:1030
EBHelmholtzOp::m_hasFine
bool m_hasFine
True if there's a finer level.
Definition: CD_EBHelmholtzOp.H:890
EBHelmholtzOp::m_doCoarsen
bool m_doCoarsen
Turn on/off exchange operation.
Definition: CD_EBHelmholtzOp.H:923
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:984
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:534
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:1835
EBHelmholtzOp::m_sideBox
std::map< std::pair< int, Side::LoHiSide >, Box > m_sideBox
Domain boxes on each side.
Definition: CD_EBHelmholtzOp.H:1131
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:1116
EBHelmholtzOp::getAcoef
const RefCountedPtr< LevelData< EBCellFAB > > & getAcoef()
Get the Helmholtz A-coefficient on cell centers.
Definition: CD_EBHelmholtzOp.cpp:174
EBHelmholtzOp::m_vecDx
RealVect m_vecDx
Vector resolution.
Definition: CD_EBHelmholtzOp.H:953
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:617
EBHelmholtzOp::m_vofIterMulti
LayoutData< VoFIterator > m_vofIterMulti
VoFIterator for "multi-cells".
Definition: CD_EBHelmholtzOp.H:1111
EBHelmholtzOp::m_betaDiagWeight
LayoutData< BaseIVFAB< Real > > m_betaDiagWeight
Weights of diagonal beta terms.
Definition: CD_EBHelmholtzOp.H:1101
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:1136
EBHelmholtzOp::getBcoef
const RefCountedPtr< LevelData< EBFluxFAB > > & getBcoef()
Get the Helmholtz B-coefficient on faces.
Definition: CD_EBHelmholtzOp.cpp:180
EBHelmholtzOp::m_Acoef
RefCountedPtr< LevelData< EBCellFAB > > m_Acoef
A-coefficient in Helmholtz equation.
Definition: CD_EBHelmholtzOp.H:1043
EBHelmholtzOp::turnOnCoarsening
void turnOnCoarsening()
Turn on coarsening operation.
Definition: CD_EBHelmholtzOp.cpp:237
EBHelmholtzOp::m_vofIterDomLo
LayoutData< VoFIterator > m_vofIterDomLo[SpaceDim]
VoF iterators for lo domain side.
Definition: CD_EBHelmholtzOp.H:1121
EBHelmholtzOp::m_vofIterIrreg
LayoutData< VoFIterator > m_vofIterIrreg
VoFIterator for irregular cells.
Definition: CD_EBHelmholtzOp.H:1106
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:1047
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:997
EBHelmholtzOp::m_restrictOpMG
EBMGRestrict m_restrictOpMG
Restriction operator if this is an MG level.
Definition: CD_EBHelmholtzOp.H:1003
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:717
EBHelmholtzOp::m_beta
Real m_beta
Beta-coefficient.
Definition: CD_EBHelmholtzOp.H:943
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:1380
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:1083
EBHelmholtzOp::m_vofIterDomHi
LayoutData< VoFIterator > m_vofIterDomHi[SpaceDim]
VoF iterators for hi domain side.
Definition: CD_EBHelmholtzOp.H:1126
EBHelmholtzOp::buildCopier
void buildCopier(Copier &a_copier, const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) override
Build copier.
Definition: CD_EBHelmholtzOp.cpp:2322
EBHelmholtzOp::defineStencils
void defineStencils()
Define stencils.
Definition: CD_EBHelmholtzOp.cpp:269
EBHelmholtzOp::prolongIncrement
void prolongIncrement(LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_correctCoarse) override final
Prolongation method.
Definition: CD_EBHelmholtzOp.cpp:746
EBHelmholtzOp::m_eblgCoarMG
EBLevelGrid m_eblgCoarMG
Coarser grids (multigrid level)
Definition: CD_EBHelmholtzOp.H:993
EBHelmholtzOp::m_BcoefIrreg
RefCountedPtr< LevelData< BaseIVFAB< Real > > > m_BcoefIrreg
B-coefficient in Helmholtz equation, but on EB faces.
Definition: CD_EBHelmholtzOp.H:1053
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
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.1