chombo-discharge
Loading...
Searching...
No Matches
CD_DataOps.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
13#ifndef CD_DATAOPS_H
14#define CD_DATAOPS_H
15
16// Std includes
17#include <array>
18
19// Chombo includes
20#include <FaceIterator.H>
21#include <VoFIterator.H>
22
23// Our includes
24#include <CD_Average.H>
25#include <CD_EBAMRData.H>
26#include <CD_Decorations.H>
27#include <CD_MFInterfaceFAB.H>
28#include <CD_NamespaceHeader.H>
29
35class DataOps
36{
37public:
43 template <typename T>
44 static void
45 setValue(LevelData<MFInterfaceFAB<T>>& a_lhs, const T& a_value);
46
52 template <typename T>
53 static int
54 sgn(const T a_value);
55
64 static void
65 averageCellVelocityToFaceVelocity(EBAMRFluxData& a_faceData,
66 const EBAMRCellData& a_cellData,
67 const Vector<ProblemDomain>& a_domains,
68 const int a_tanGhosts,
69 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
70
79 static void
80 averageCellVelocityToFaceVelocity(LevelData<EBFluxFAB>& a_faceData,
81 const LevelData<EBCellFAB>& a_cellData,
82 const ProblemDomain& a_domain,
83 const int a_tanGhosts,
84 LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter);
85
96 static void
97 averageCellToFace(EBAMRFluxData& a_faceData,
98 const EBAMRCellData& a_cellData,
99 const Vector<ProblemDomain>& a_domains,
100 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
101
114 static void
115 averageCellToFace(EBAMRFluxData& a_faceData,
116 const EBAMRCellData& a_cellData,
117 const Vector<ProblemDomain>& a_domains,
118 const int a_tanGhosts,
119 const Interval& a_faceInterval,
120 const Interval& a_cellInterval,
121 const Average& a_average,
122 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
123
135 static void
136 averageCellToFace(LevelData<EBFluxFAB>& a_faceData,
137 const LevelData<EBCellFAB>& a_cellData,
138 const ProblemDomain& a_domain,
139 const int a_tanGhosts,
140 const Interval& a_faceInterval,
141 const Interval& a_cellInterval,
142 const Average& a_average,
143 LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter);
144
153 static void
154 averageFaceToCell(EBAMRCellData& a_cellData,
155 const EBAMRFluxData& a_faceData,
156 const Vector<ProblemDomain>& a_domains,
157 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
158
167 static void
168 averageFaceToCell(LevelData<EBCellFAB>& a_cellData,
169 const LevelData<EBFluxFAB>& a_fluxData,
170 const ProblemDomain& a_domain,
171 LayoutData<VoFIterator>& a_vofIter);
172
181 static void
182 axby(LevelData<EBCellFAB>& a_lhs,
183 const LevelData<EBCellFAB>& a_x,
184 const LevelData<EBCellFAB>& a_y,
185 const Real a_a,
186 const Real a_b) noexcept;
187
195 static void
196 compute(EBAMRCellData& a_data,
197 const std::function<Real(const Real a_cellValue)>& a_func,
198 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter) noexcept;
199
207 static void
208 compute(LevelData<EBCellFAB>& a_data,
209 const std::function<Real(const Real a_cellValue)>& a_func,
210 LayoutData<VoFIterator>& a_vofIter) noexcept;
211
219 static void
220 dotProduct(MFAMRCellData& a_result, const MFAMRCellData& a_data1, const MFAMRCellData& a_data2);
221
230 static void
231 dotProduct(LevelData<MFCellFAB>& a_result, const LevelData<MFCellFAB>& a_data1, const LevelData<MFCellFAB>& a_data2);
232
241 static void
242 dotProduct(EBAMRCellData& a_result,
243 const EBAMRCellData& a_data1,
244 const EBAMRCellData& a_data2,
245 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
246
255 static void
256 dotProduct(LevelData<EBCellFAB>& a_result,
257 const LevelData<EBCellFAB>& a_data1,
258 const LevelData<EBCellFAB>& a_data2,
259 LayoutData<VoFIterator>& a_vofIter);
260
270 static void
271 dotProduct(EBCellFAB& a_result,
272 const EBCellFAB& a_data1,
273 const EBCellFAB& a_data2,
274 const Box& a_box,
275 VoFIterator& a_vofIter);
276
287 static void
288 filterSmooth(EBAMRCellData& a_data, const Real a_alpha, const int a_stride, const bool a_zeroEB) noexcept;
289
300 static void
301 filterSmooth(LevelData<EBCellFAB>& a_data, const Real a_alpha, const int a_stride, const bool a_zeroEB) noexcept;
302
309 static void
310 incr(MFAMRCellData& a_lhs, const MFAMRCellData& a_rhs, const Real a_scale) noexcept;
311
318 static void
319 incr(LevelData<MFCellFAB>& a_lhs, const LevelData<MFCellFAB>& a_rhs, const Real a_scale) noexcept;
320
327 static void
328 incr(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs, const Real& a_scale) noexcept;
329
336 static void
337 incr(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const Real& a_scale) noexcept;
338
345 static void
346 incr(EBAMRFluxData& a_lhs, const EBAMRFluxData& a_rhs, const Real& a_scale);
347
354 static void
355 incr(LevelData<EBFluxFAB>& a_lhs, const LevelData<EBFluxFAB>& a_rhs, const Real& a_scale);
356
363 static void
364 incr(EBFluxFAB& a_lhs, const EBFluxFAB& a_rhs, const Real& a_scale);
365
373 static void
374 incr(EBAMRIVData& a_lhs, const EBAMRIVData& a_rhs, const Real& a_scale);
375
383 static void
384 incr(LevelData<BaseIVFAB<Real>>& a_lhs, const LevelData<BaseIVFAB<Real>>& a_rhs, const Real& a_scale);
385
393 static void
394 incr(EBAMRIVData& a_lhs,
395 const EBAMRIVData& a_rhs,
396 const Real& a_scale,
397 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
398
406 static void
407 incr(LevelData<BaseIVFAB<Real>>& a_lhs,
408 const LevelData<BaseIVFAB<Real>>& a_rhs,
409 const Real& a_scale,
410 LayoutData<VoFIterator>& a_vofIter);
411
418 static void
419 incr(EBAMRIFData& a_lhs, const EBAMRIFData& a_rhs, const Real& a_scale);
420
427 static void
428 incr(LevelData<DomainFluxIFFAB>& a_lhs, const LevelData<DomainFluxIFFAB>& a_rhs, const Real& a_scale);
429
438 static void
439 incr(EBAMRCellData& a_lhs,
440 const EBAMRIVData& a_rhs,
441 const Real a_scale,
442 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
443
452 static void
453 incr(LevelData<EBCellFAB>& a_lhs,
454 const LevelData<BaseIVFAB<Real>>& a_rhs,
455 const Real a_scale,
456 LayoutData<VoFIterator>& a_vofIter);
457
466 static void
467 incr(EBAMRIVData& a_lhs,
468 const EBAMRCellData& a_rhs,
469 const Real a_scale,
470 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
471
480 static void
481 incr(LevelData<BaseIVFAB<Real>>& a_lhs,
482 const LevelData<EBCellFAB>& a_rhs,
483 const Real a_scale,
484 LayoutData<VoFIterator>& a_vofIter);
485
491 static void
492 copy(MFAMRCellData& a_dst, const MFAMRCellData& a_src);
493
499 static void
500 copy(EBAMRCellData& a_dst, const EBAMRCellData& a_src);
501
507 static void
508 copy(EBAMRIVData& a_dst, const EBAMRIVData& a_src);
509
517 template <typename T>
518 static void
519 copy(EBAMRData<T>& a_dst, const EBAMRData<T>& a_src, const Interval a_dstInterv, const Interval a_srcInterv);
520
528 static void
529 divide(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs, const int a_lhsComp, const int a_rhsComp);
530
538 static void
539 divide(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs, const int a_lhsComp, const int a_rhsComp);
540
546 static void
547 divideByScalar(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs);
548
554 static void
555 divideByScalar(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs);
556
564 static void
565 divideFallback(EBAMRCellData& a_numerator,
566 const EBAMRCellData& a_denominator,
567 const Real a_fallback,
568 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
569
578 static void
579 divideFallback(LevelData<EBCellFAB>& a_numerator,
580 const LevelData<EBCellFAB>& a_denominator,
581 const Real a_fallback,
582 LayoutData<VoFIterator>& a_vofIter);
583
591 static void
592 divideFallback(EBAMRCellData& a_numerator,
593 const EBAMRCellData& a_denominator,
594 const EBAMRCellData& a_fallback,
595 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
596
605 static void
606 divideFallback(LevelData<EBCellFAB>& a_numerator,
607 const LevelData<EBCellFAB>& a_denominator,
608 const LevelData<EBCellFAB>& a_fallback,
609 LayoutData<VoFIterator>& a_vofIter);
610
618 static void
619 floor(EBAMRCellData& a_lhs, const Real a_value, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
620
628 static void
629 floor(LevelData<EBCellFAB>& a_lhs, const Real a_value, LayoutData<VoFIterator>& a_vofIter);
630
637 static void
638 floor(EBAMRIVData& a_lhs, const Real a_value, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
639
646 static void
647 floor(LevelData<BaseIVFAB<Real>>& a_lhs, const Real a_value, LayoutData<VoFIterator>& a_vofIter);
648
656 static void
657 roof(EBAMRCellData& a_lhs, const Real a_value, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
658
666 static void
667 roof(LevelData<EBCellFAB>& a_lhs, const Real a_value, LayoutData<VoFIterator>& a_vofIter);
668
675 static void
676 roof(EBAMRIVData& a_lhs, const Real a_value, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
677
684 static void
685 roof(LevelData<BaseIVFAB<Real>>& a_lhs, const Real a_value, LayoutData<VoFIterator>& a_vofIter);
686
694 static void
695 max(EBAMRCellData& a_data,
696 const EBAMRCellData& a_data1,
697 const EBAMRCellData& a_data2,
698 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
699
707 static void
708 max(LevelData<EBCellFAB>& a_data,
709 const LevelData<EBCellFAB>& a_data1,
710 const LevelData<EBCellFAB>& a_data2,
711 LayoutData<VoFIterator>& a_vofIter);
712
722 static void
723 getMaxMin(Real& max,
724 Real& min,
725 EBAMRCellData& a_data,
726 const int a_comp,
727 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
728
737 static void
738 getMaxMin(Real& a_max,
739 Real& a_min,
740 LevelData<EBCellFAB>& a_data,
741 const int a_comp,
742 LayoutData<VoFIterator>& a_vofIter);
743
753 static void
754 getMaxMin(Real& a_max,
755 Real& a_min,
756 EBAMRFluxData& a_data,
757 const int a_comp,
758 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter) noexcept;
759
768 static void
769 getMaxMin(Real& a_max,
770 Real& a_min,
771 LevelData<EBFluxFAB>& a_data,
772 const int a_comp,
773 LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter) noexcept;
774
782 static void
783 getMaxMin(Vector<Real>& a_max, Vector<Real>& a_min, Vector<EBAMRCellData>& a_data);
784
793 static void
794 getMaxMinNorm(Real& a_max,
795 Real& a_min,
796 EBAMRCellData& data,
797 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
798
806 static void
807 getMaxMinNorm(Real& a_max, Real& a_min, LevelData<EBCellFAB>& data, LayoutData<VoFIterator>& a_vofIter);
808
817 static void
818 getMaxMinNorm(Real& a_max,
819 Real& a_min,
820 EBAMRIVData& a_data,
821 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
822
830 static void
831 getMaxMinNorm(Real& a_max, Real& a_min, LevelData<BaseIVFAB<Real>>& a_data, LayoutData<VoFIterator>& a_vofIter);
832
838 static void
839 invert(EBAMRFluxData& a_data);
840
847 static void
848 invert(EBAMRFluxData& a_data, Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
849
856 static void
857 invert(LevelData<EBFluxFAB>& a_data, LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter);
858
864 static void
865 invert(EBFluxFAB& a_data);
866
875 static void
876 kappaSum(Real& a_mass,
877 const LevelData<EBCellFAB>& a_lhs,
878 const LevelData<EBCellFAB>& a_regularMask,
879 const int a_comp,
880 LayoutData<VoFIterator>& a_vofIter);
881
887 static void
888 kappaScale(EBAMRCellData& a_data, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter) noexcept;
889
895 static void
896 kappaScale(LevelData<EBCellFAB>& a_data, LayoutData<VoFIterator>& a_vofIter) noexcept;
897
904 static void
905 kappaScale(MFAMRCellData& a_data,
906 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase0,
907 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase1) noexcept;
908
915 static void
916 kappaScale(LevelData<MFCellFAB>& a_data,
917 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase0,
918 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase1) noexcept;
919
925 static void
926 volumeScale(EBAMRCellData& a_data, const Vector<Real>& a_dx);
927
934 static void
935 multiply(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs);
936
943 static void
944 multiply(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs);
945
952 static void
953 multiply(EBAMRFluxData& a_lhs, const EBAMRFluxData& a_rhs);
954
961 static void
962 multiply(LevelData<EBFluxFAB>& a_lhs, const LevelData<EBFluxFAB>& a_rhs);
963
971 static void
972 multiply(EBAMRIVData& a_lhs,
973 const EBAMRIVData& a_rhs,
974 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
975
983 static void
984 multiply(LevelData<BaseIVFAB<Real>>& a_lhs,
985 const LevelData<BaseIVFAB<Real>>& a_rhs,
986 LayoutData<VoFIterator>& a_vofIter);
987
994 static void
995 multiplyScalar(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs);
996
1003 static void
1004 multiplyScalar(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_rhs);
1005
1012 static void
1013 multiplyScalar(EBAMRIVData& a_lhs,
1014 const EBAMRIVData& a_rhs,
1015 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1016
1023 static void
1024 multiplyScalar(LevelData<BaseIVFAB<Real>>& a_lhs,
1025 const LevelData<BaseIVFAB<Real>>& a_rhs,
1026 LayoutData<VoFIterator>& a_vofIter);
1027
1038 static Real
1039 norm(const LevelData<EBCellFAB>& a_data,
1040 const LevelData<EBCellFAB>& a_regularMask,
1041 const int a_p,
1042 LayoutData<VoFIterator>& a_vofIter,
1043 const int a_comp = 0);
1044
1054 static void
1055 plus(EBAMRCellData& a_lhs, const EBAMRCellData& a_rhs, const int a_srcComp, const int a_dstComp, const int a_numComp);
1056
1066 static void
1067 plus(LevelData<EBCellFAB>& a_lhs,
1068 const LevelData<EBCellFAB>& a_rhs,
1069 const int a_srcComp,
1070 const int a_dstComp,
1071 const int a_numComp);
1072
1078 static void
1079 scale(MFAMRCellData& a_lhs, const Real& a_scale) noexcept;
1080
1086 static void
1087 scale(LevelData<MFCellFAB>& a_lhs, const Real& a_scale) noexcept;
1088
1094 static void
1095 scale(MFAMRFluxData& a_lhs, const Real& a_scale);
1096
1102 static void
1103 scale(LevelData<MFFluxFAB>& a_lhs, const Real& a_scale);
1104
1111 static void
1112 scale(EBAMRIVData& a_lhs, const Real& a_scale, const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1113
1119 static void
1120 scale(EBAMRCellData& a_lhs, const Real a_scale) noexcept;
1121
1127 static void
1128 scale(LevelData<EBCellFAB>& a_lhs, const Real a_scale) noexcept;
1129
1135 static void
1136 scale(EBAMRFluxData& a_lhs, const Real a_scale);
1137
1143 static void
1144 scale(LevelData<EBFluxFAB>& a_lhs, const Real a_scale);
1145
1152 static void
1153 scale(LevelData<BaseIVFAB<Real>>& a_lhs, const Real& a_scale, LayoutData<VoFIterator>& a_vofIter);
1154
1162 static void
1163 setCoveredValue(EBAMRCellData& a_lhs, const EBAMRCellData& a_coveredMask, const int a_comp, const Real a_value);
1164
1172 static void
1173 setCoveredValue(LevelData<EBCellFAB>& a_lhs,
1174 const LevelData<EBCellFAB>& a_coveredMask,
1175 const int a_comp,
1176 const Real a_value);
1177
1184 static void
1185 setCoveredValue(EBAMRCellData& a_lhs, const EBAMRCellData& a_coveredMask, const Real a_value);
1186
1193 static void
1194 setCoveredValue(LevelData<EBCellFAB>& a_lhs, const LevelData<EBCellFAB>& a_coveredMask, const Real a_value);
1195
1202 static void
1203 setCoveredValue(EBAMRFluxData& a_lhs, const int a_comp, const Real a_value) noexcept;
1204
1211 static void
1212 setCoveredValue(LevelData<EBFluxFAB>& a_lhs, const int a_comp, const Real a_value) noexcept;
1213
1219 static void
1220 setCoveredValue(EBAMRFluxData& a_lhs, const Real a_value) noexcept;
1221
1227 static void
1228 setCoveredValue(LevelData<EBFluxFAB>& a_lhs, const Real a_value) noexcept;
1229
1236 static void
1237 setInvalidValue(EBAMRCellData& a_lhs, const Vector<int>& a_refRat, const Real a_value);
1238
1249 static void
1250 setValue(MFAMRCellData& a_lhs,
1251 const std::function<Real(const RealVect)>& a_function,
1252 const RealVect& a_probLo,
1253 const Vector<Real>& a_dx,
1254 const int a_comp,
1255 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase0,
1256 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase1);
1257
1268 static void
1269 setValue(LevelData<MFCellFAB>& a_lhs,
1270 const std::function<Real(const RealVect)>& a_function,
1271 const RealVect a_probLo,
1272 const Real a_dx,
1273 const int a_comp,
1274 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase0,
1275 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase1);
1276
1286 static void
1287 setValue(EBAMRCellData& a_lhs,
1288 const std::function<Real(const RealVect)>& a_function,
1289 const RealVect& a_probLo,
1290 const Vector<Real>& a_dx,
1291 const int a_comp,
1292 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1293
1303 static void
1304 setValue(LevelData<EBCellFAB>& a_lhs,
1305 const std::function<Real(const RealVect)>& a_function,
1306 const RealVect a_probLo,
1307 const Real a_dx,
1308 const int a_comp,
1309 LayoutData<VoFIterator>& a_vofIter);
1310
1320 static void
1321 setValue(EBAMRFluxData& a_lhs,
1322 const std::function<Real(const RealVect)>& a_function,
1323 const RealVect& a_probLo,
1324 const Vector<Real>& a_dx,
1325 const int a_comp,
1326 Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
1327
1337 static void
1338 setValue(LevelData<EBFluxFAB>& a_lhs,
1339 const std::function<Real(const RealVect)>& a_function,
1340 const RealVect a_probLo,
1341 const Real a_dx,
1342 const int a_comp,
1343 LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter);
1344
1354 static void
1355 setValue(EBAMRIVData& a_lhs,
1356 const std::function<Real(const RealVect)>& a_function,
1357 const RealVect& a_probLo,
1358 const Vector<Real>& a_dx,
1359 const int a_comp,
1360 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1361
1371 static void
1372 setValue(LevelData<BaseIVFAB<Real>>& a_lhs,
1373 const std::function<Real(const RealVect)>& a_function,
1374 const RealVect a_probLo,
1375 const Real a_dx,
1376 const int a_comp,
1377 LayoutData<VoFIterator>& a_vofIter);
1378
1387 static void
1388 setValue(EBAMRCellData& a_lhs,
1389 const std::function<RealVect(const RealVect)>& a_function,
1390 const RealVect& a_probLo,
1391 const Vector<Real>& a_dx,
1392 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1393
1402 static void
1403 setValue(LevelData<EBCellFAB>& a_lhs,
1404 const std::function<RealVect(const RealVect)>& a_function,
1405 const RealVect a_probLo,
1406 const Real a_dx,
1407 LayoutData<VoFIterator>& a_vofIter);
1408
1414 static void
1415 setValue(EBAMRCellData& a_data, const Real& a_value);
1416
1423 static void
1424 setValue(EBAMRCellData& a_lhs, const Real a_value, const int a_comp);
1425
1432 static void
1433 setValue(LevelData<EBCellFAB>& a_lhs, const Real a_value, const int a_comp);
1434
1440 static void
1441 setValue(LevelData<EBCellFAB>& a_lhs, const Real a_value);
1442
1448 static void
1449 setValue(LevelData<EBFluxFAB>& a_lhs, const Real a_value);
1450
1456 static void
1457 setValue(LevelData<BaseIVFAB<Real>>& a_lhs, const Real a_value);
1458
1464 static void
1465 setValue(EBAMRFluxData& a_data, const Real& a_value);
1466
1472 static void
1473 setValue(EBAMRIVData& a_data, const Real& a_value);
1474
1480 static void
1481 setValue(MFAMRCellData& a_lhs, const Real& a_value) noexcept;
1482
1488 static void
1489 setValue(LevelData<MFCellFAB>& a_lhs, const Real& a_value) noexcept;
1490
1496 static void
1497 setValue(MFAMRFluxData& a_lhs, const Real& a_value);
1498
1504 static void
1505 setValue(LevelData<MFFluxFAB>& a_lhs, const Real& a_value);
1506
1512 static void
1513 setValue(MFAMRIVData& a_lhs, const Real& a_value);
1514
1520 static void
1521 setValue(LevelData<MFBaseIVFAB>& a_lhs, const Real& a_value);
1522
1528 static void
1529 setValue(EBAMRIFData& a_lhs, const Real a_value);
1530
1536 static void
1537 setValue(LevelData<DomainFluxIFFAB>& a_lhs, const Real a_value);
1538
1543 static void
1544 sum(Real& a_value);
1545
1552 static void
1553 squareRoot(EBAMRFluxData& a_lhs, Vector<RefCountedPtr<LayoutData<std::array<FaceIterator, SpaceDim>>>>& a_faceIter);
1554
1561 static void
1562 squareRoot(LevelData<EBFluxFAB>& a_lhs, LayoutData<std::array<FaceIterator, SpaceDim>>& a_faceIter);
1563
1571 static void
1572 squareRoot(MFAMRCellData& a_lhs,
1573 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase0,
1574 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIterPhase1);
1575
1583 static void
1584 squareRoot(LevelData<MFCellFAB>& a_lhs,
1585 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase0,
1586 const RefCountedPtr<LayoutData<VoFIterator>>& a_vofIterPhase1);
1587
1596 static void
1597 vectorLength(EBAMRCellData& a_lhs,
1598 const EBAMRCellData& a_rhs,
1599 const EBAMRCellData& a_notCovered,
1600 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1601
1610 static void
1611 vectorLength(LevelData<EBCellFAB>& a_lhs,
1612 const LevelData<EBCellFAB>& a_rhs,
1613 const LevelData<EBCellFAB>& a_notCovered,
1614 LayoutData<VoFIterator>& a_vofIter);
1615
1625 static void
1626 vectorLength(EBCellFAB& a_lhs,
1627 const EBCellFAB& a_rhs,
1628 const EBCellFAB& a_notCovered,
1629 const Box& box,
1630 VoFIterator& a_vofIter);
1631
1640 static void
1641 vectorLength2(EBAMRCellData& a_lhs,
1642 const EBAMRCellData& a_rhs,
1643 const EBAMRCellData& a_notCovered,
1644 const Vector<RefCountedPtr<LayoutData<VoFIterator>>>& a_vofIter);
1645
1654 static void
1655 vectorLength2(LevelData<EBCellFAB>& a_lhs,
1656 const LevelData<EBCellFAB>& a_rhs,
1657 const LevelData<EBCellFAB>& a_notCovered,
1658 LayoutData<VoFIterator>& a_vofIter);
1659
1669 static void
1670 vectorLength2(EBCellFAB& a_lhs,
1671 const EBCellFAB& a_rhs,
1672 const EBCellFAB& a_notCovered,
1673 const Box& box,
1674 VoFIterator& a_vofIter);
1675
1683 static void
1684 computeMinValidBox(RealVect& a_lo, RealVect& a_hi, const RealVect& a_normal, const RealVect& a_centroid);
1685
1693 static bool
1694 allCornersInsideEb(const Vector<RealVect>& a_corners, const RealVect& a_normal, const RealVect& a_centroid);
1695
1702 static void
1703 shiftCorners(Vector<RealVect>& a_corners, const RealVect& a_distance);
1704};
1705
1706#include <CD_NamespaceFooter.H>
1707
1708#include <CD_DataOpsImplem.H>
1709
1710#endif
CD_Average.H
Declaration of averaging methods.
Average
Average
Various averaging methods.
Definition CD_Average.H:25
CD_DataOpsImplem.H
Implementation of CD_DataOps.H.
CD_Decorations.H
Declaration of various decoration macros.
CD_EBAMRData.H
Class for holding data across EBAMR hierarchies.
CD_MFInterfaceFAB.H
Declaration of a class that allocates irregular data holders over the interface between two phases.
DataOps
Agglomeration of useful data operations.
Definition CD_DataOps.H:36
DataOps::scale
static void scale(MFAMRCellData &a_lhs, const Real &a_scale) noexcept
Scale data by factor.
Definition CD_DataOps.cpp:2503
DataOps::floor
static void floor(EBAMRCellData &a_lhs, const Real a_value, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Floor values in data holder. This sets all values below a_value to a_value.
Definition CD_DataOps.cpp:1465
DataOps::getMaxMin
static void getMaxMin(Real &max, Real &min, EBAMRCellData &a_data, const int a_comp, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Get maximum and minimum value of specified component.
Definition CD_DataOps.cpp:1711
DataOps::averageFaceToCell
static void averageFaceToCell(EBAMRCellData &a_cellData, const EBAMRFluxData &a_faceData, const Vector< ProblemDomain > &a_domains, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Average all components of face centered data to cell centers.
Definition CD_DataOps.cpp:362
DataOps::volumeScale
static void volumeScale(EBAMRCellData &a_data, const Vector< Real > &a_dx)
Scale data by dx^SpaceDim.
Definition CD_DataOps.cpp:2236
DataOps::vectorLength2
static void vectorLength2(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs, const EBAMRCellData &a_notCovered, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Compute the squared vector length of a data holder. Sets a_lhs = |a_rhs|^2 where a_rhs contains Space...
Definition CD_DataOps.cpp:3594
DataOps::sgn
static int sgn(const T a_value)
Sign function. Returns +/- if the value is > 0 or < 0.
Definition CD_DataOpsImplem.H:43
DataOps::getMaxMinNorm
static void getMaxMinNorm(Real &a_max, Real &a_min, EBAMRCellData &data, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Get maximum and minimum value of normed data.
Definition CD_DataOps.cpp:1879
DataOps::filterSmooth
static void filterSmooth(EBAMRCellData &a_data, const Real a_alpha, const int a_stride, const bool a_zeroEB) noexcept
Apply a convolved filter phi = alpha * phi_i + 0.5*(1-alpha) * [phi_(i+s) + phi_(i-s)] in each direct...
Definition CD_DataOps.cpp:679
DataOps::divide
static void divide(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs, const int a_lhsComp, const int a_rhsComp)
Divide one data holder by another. Does specified components.
Definition CD_DataOps.cpp:1243
DataOps::roof
static void roof(EBAMRCellData &a_lhs, const Real a_value, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Roof values in data holder. This sets all values above a_value to a_value.
Definition CD_DataOps.cpp:1556
DataOps::setInvalidValue
static void setInvalidValue(EBAMRCellData &a_lhs, const Vector< int > &a_refRat, const Real a_value)
Set value in cells that are covered by a finer grid.
Definition CD_DataOps.cpp:2771
DataOps::shiftCorners
static void shiftCorners(Vector< RealVect > &a_corners, const RealVect &a_distance)
Shift points by a specified distance.
Definition CD_DataOps.cpp:3771
DataOps::divideByScalar
static void divideByScalar(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs)
Divide all components in a data holder by a scalar value.
Definition CD_DataOps.cpp:1279
DataOps::multiply
static void multiply(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs)
Multiply data holder by another data holder.
Definition CD_DataOps.cpp:2246
DataOps::invert
static void invert(EBAMRFluxData &a_data)
Invert data.
DataOps::computeMinValidBox
static void computeMinValidBox(RealVect &a_lo, RealVect &a_hi, const RealVect &a_normal, const RealVect &a_centroid)
Compute the tightest possible valid box around a cut-cell volume.
Definition CD_DataOps.cpp:3687
DataOps::invert
static void invert(EBFluxFAB &a_data)
Invert data.
DataOps::sum
static void sum(Real &a_value)
Compute the sum of the input value.
Definition CD_DataOps.cpp:3367
DataOps::vectorLength
static void vectorLength(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs, const EBAMRCellData &a_notCovered, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Compute the vector length of a data holder. Sets a_lhs = |a_rhs| where a_rhs contains SpaceDim compon...
Definition CD_DataOps.cpp:3500
DataOps::incr
static void incr(MFAMRCellData &a_lhs, const MFAMRCellData &a_rhs, const Real a_scale) noexcept
Function which increments data in the form a_lhs = a_lhs + a_rhs*a_scale for all components.
Definition CD_DataOps.cpp:820
DataOps::norm
static Real norm(const LevelData< EBCellFAB > &a_data, const LevelData< EBCellFAB > &a_regularMask, const int a_p, LayoutData< VoFIterator > &a_vofIter, const int a_comp=0)
Compute specified norm of data.
Definition CD_DataOps.cpp:2418
DataOps::kappaScale
static void kappaScale(EBAMRCellData &a_data, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter) noexcept
Scale data by volume fraction.
Definition CD_DataOps.cpp:2142
DataOps::setValue
static void setValue(LevelData< MFInterfaceFAB< T > > &a_lhs, const T &a_value)
Set value in an MFInterfaceFAB data holder.
Definition CD_DataOpsImplem.H:24
DataOps::allCornersInsideEb
static bool allCornersInsideEb(const Vector< RealVect > &a_corners, const RealVect &a_normal, const RealVect &a_centroid)
Check if all corners/points are inside an EB.
Definition CD_DataOps.cpp:3754
DataOps::divideFallback
static void divideFallback(EBAMRCellData &a_numerator, const EBAMRCellData &a_denominator, const Real a_fallback, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Divide data. If the denominator is zero, set the value to a fallback option.
Definition CD_DataOps.cpp:1387
DataOps::dotProduct
static void dotProduct(MFAMRCellData &a_result, const MFAMRCellData &a_data1, const MFAMRCellData &a_data2)
Compote the cell-wise dot product between two data holders.
Definition CD_DataOps.cpp:537
DataOps::squareRoot
static void squareRoot(EBAMRFluxData &a_lhs, Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > &a_faceIter)
Compute the square root of the input data.
Definition CD_DataOps.cpp:3384
DataOps::setCoveredValue
static void setCoveredValue(EBAMRCellData &a_lhs, const EBAMRCellData &a_coveredMask, const int a_comp, const Real a_value)
Set value in covered cells. Does specified component.
Definition CD_DataOps.cpp:2655
DataOps::max
static void max(EBAMRCellData &a_data, const EBAMRCellData &a_data1, const EBAMRCellData &a_data2, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter)
Make a_data = max(a_data1, a_data2).
Definition CD_DataOps.cpp:1645
DataOps::kappaSum
static void kappaSum(Real &a_mass, const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_regularMask, const int a_comp, LayoutData< VoFIterator > &a_vofIter)
Compute the kappa-weighted sum of the input data. Assumes one component.
Definition CD_DataOps.cpp:2087
DataOps::compute
static void compute(EBAMRCellData &a_data, const std::function< Real(const Real a_cellValue)> &a_func, const Vector< RefCountedPtr< LayoutData< VoFIterator > > > &a_vofIter) noexcept
Compute a new value given the old cell value.
Definition CD_DataOps.cpp:482
DataOps::plus
static void plus(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs, const int a_srcComp, const int a_dstComp, const int a_numComp)
General addition operator for adding together data. The user can choose which components to add.
Definition CD_DataOps.cpp:885
DataOps::copy
static void copy(MFAMRCellData &a_dst, const MFAMRCellData &a_src)
Copy data from one data holder to another.
Definition CD_DataOps.cpp:1201
DataOps::averageCellVelocityToFaceVelocity
static void averageCellVelocityToFaceVelocity(EBAMRFluxData &a_faceData, const EBAMRCellData &a_cellData, const Vector< ProblemDomain > &a_domains, const int a_tanGhosts, Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > &a_faceIter)
Routine which computes the average of a cell-centered quantity on faces for the normal component only...
Definition CD_DataOps.cpp:28
DataOps::averageCellToFace
static void averageCellToFace(EBAMRFluxData &a_faceData, const EBAMRCellData &a_cellData, const Vector< ProblemDomain > &a_domains, Vector< RefCountedPtr< LayoutData< std::array< FaceIterator, SpaceDim > > > > &a_faceIter)
Average all components of the cell-centered data to faces (arithmetic, no tangential ghost faces).
Definition CD_DataOps.cpp:148
DataOps::multiplyScalar
static void multiplyScalar(EBAMRCellData &a_lhs, const EBAMRCellData &a_rhs)
Multiply data holder by another data holder.
Definition CD_DataOps.cpp:2341
DataOps::axby
static void axby(LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_x, const LevelData< EBCellFAB > &a_y, const Real a_a, const Real a_b) noexcept
Compute lhs = a * x + b * y.
Definition CD_DataOps.cpp:462
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
Generated by doxygen 1.9.8