Skip to content

Commit 15e128a

Browse files
WallyMaierWallyMaier
committed
replace area+= with geometry toolboxes. Will need to check include statements
1 parent 987b825 commit 15e128a

38 files changed

Lines changed: 131 additions & 378 deletions

‎Common/src/geometry/CGeometry.cpp‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -2687,10 +2687,7 @@ void CGeometry::ComputeSurf_Straightness(CConfig *config,
26872687
UnitNormal = Normal;
26882688

26892689
/*--- Compute unit normal. ---*/
2690-
Area = 0.0;
2691-
for (iDim = 0; iDim < nDim; iDim++)
2692-
Area += Normal[iDim]*Normal[iDim];
2693-
Area = sqrt(Area);
2690+
Area = GeometryToolbox::Norm(nDim, Normal);
26942691

26952692
/*--- Negate for outward convention. ---*/
26962693
for (iDim = 0; iDim < nDim; iDim++)

‎Common/src/geometry/CMultiGridGeometry.cpp‎

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1135,8 +1135,7 @@ void CMultiGridGeometry::SetControlVolume(CConfig *config, CGeometry *fine_grid,
11351135

11361136
for (iEdge = 0; iEdge < nEdge; iEdge++) {
11371137
const auto NormalFace = edges->GetNormal(iEdge);
1138-
Area = 0.0; for (iDim = 0; iDim < nDim; iDim++) Area += NormalFace[iDim]*NormalFace[iDim];
1139-
Area = sqrt(Area);
1138+
Area = GeometryToolbox::Norm(nDim, NormalFace);
11401139
if (Area == 0.0) {
11411140
su2double DefaultNormal[3] = {EPS*EPS};
11421141
edges->SetNormal(iEdge, DefaultNormal);
@@ -1181,8 +1180,7 @@ void CMultiGridGeometry::SetBoundControlVolume(CConfig *config, CGeometry *fine_
11811180
for (iMarker = 0; iMarker < nMarker; iMarker ++)
11821181
for (iVertex = 0; iVertex < nVertex[iMarker]; iVertex++) {
11831182
NormalFace = vertex[iMarker][iVertex]->GetNormal();
1184-
Area = 0.0; for (iDim = 0; iDim < nDim; iDim++) Area += NormalFace[iDim]*NormalFace[iDim];
1185-
Area = sqrt(Area);
1183+
Area = GeometryToolbox::Norm(nDim, NormalFace);
11861184
if (Area == 0.0) for (iDim = 0; iDim < nDim; iDim++) NormalFace[iDim] = EPS*EPS;
11871185
}
11881186

‎Common/src/geometry/CPhysicalGeometry.cpp‎

Lines changed: 5 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -4590,8 +4590,8 @@ void CPhysicalGeometry::SetPositive_ZArea(CConfig *config) {
45904590
if (axisymmetric) AxiFactor = 2.0*PI_NUMBER*nodes->GetCoord(iPoint, 1);
45914591
else AxiFactor = 1.0;
45924592

4593-
if (nDim == 2) WettedArea += AxiFactor * sqrt (Normal[0]*Normal[0] + Normal[1]*Normal[1]);
4594-
if (nDim == 3) WettedArea += sqrt (Normal[0]*Normal[0] + Normal[1]*Normal[1] + Normal[2]*Normal[2]);
4593+
if (nDim == 2) WettedArea = AxiFactor * GeometryToolbox::Norm(nDim, Normal);
4594+
if (nDim == 3) WettedArea = GeometryToolbox::Norm(nDim, Normal);
45954595

45964596
if (Normal[0] < 0) PositiveXArea -= Normal[0];
45974597
if (Normal[1] < 0) PositiveYArea -= Normal[1];
@@ -5633,10 +5633,8 @@ void CPhysicalGeometry::SetTurboVertex(CConfig *config, unsigned short val_iZone
56335633
/*--- compute the face area associated with the vertex ---*/
56345634
vertex[iMarker][iVertex]->GetNormal(NormalArea);
56355635
for (iDim = 0; iDim < nDim; iDim++) NormalArea[iDim] = -NormalArea[iDim];
5636-
Area = 0.0;
5637-
for (iDim = 0; iDim < nDim; iDim++)
5638-
Area += NormalArea[iDim]*NormalArea[iDim];
5639-
Area = sqrt(Area);
5636+
Area = GeometryToolbox::Norm(nDim, Normal);
5637+
56405638
for (iDim = 0; iDim < nDim; iDim++) NormalArea[iDim] /= Area;
56415639
/*--- store all the all the info into the auxiliary containers ---*/
56425640
disordered[jSpan][nVertexSpanHalo[jSpan]] = iPoint;
@@ -6814,8 +6812,7 @@ void CPhysicalGeometry::SetBoundControlVolume(CConfig *config, unsigned short ac
68146812
for (iMarker = 0; iMarker < nMarker; iMarker ++)
68156813
for (iVertex = 0; iVertex < nVertex[iMarker]; iVertex++) {
68166814
NormalFace = vertex[iMarker][iVertex]->GetNormal();
6817-
Area = 0.0; for (iDim = 0; iDim < nDim; iDim++) Area += NormalFace[iDim]*NormalFace[iDim];
6818-
Area = sqrt(Area);
6815+
Area = GeometryToolbox::Norm(nDim, NormalFace);
68196816
if (Area == 0.0) for (iDim = 0; iDim < nDim; iDim++) NormalFace[iDim] = EPS*EPS;
68206817
}
68216818

‎Common/src/interface_interpolation/CSlidingMesh.cpp‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -451,10 +451,7 @@ void CSlidingMesh::SetTransferCoeff(const CConfig* const* config) {
451451
target_geometry->vertex[markTarget][iVertex]->GetNormal(Normal);
452452

453453
/*--- The value of Area computed here includes also portion of boundary belonging to different marker ---*/
454-
Area = 0.0;
455-
for (iDim = 0; iDim < nDim; iDim++)
456-
Area += Normal[iDim]*Normal[iDim];
457-
Area = sqrt(Area);
454+
Area = GeometryToolbox::Norm(nDim, Normal);
458455

459456
for (iDim = 0; iDim < nDim; iDim++)
460457
Normal[iDim] /= Area;

‎SU2_CFD/include/solvers/CFVMFlowSolverBase.inl‎

Lines changed: 4 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -706,9 +706,7 @@ su2double CFVMFlowSolverBase<V, R>::GetInletAtVertex(su2double* val_inlet, unsig
706706
/*-- Compute boundary face area for this vertex. ---*/
707707

708708
geometry->vertex[iMarker][iVertex]->GetNormal(Normal);
709-
Area = 0.0;
710-
for (iDim = 0; iDim < nDim; iDim++) Area += Normal[iDim] * Normal[iDim];
711-
Area = sqrt(Area);
709+
Area = GeometryToolbox::Norm(nDim, Normal);
712710

713711
/*--- Access and store the inlet variables for this vertex. ---*/
714712

@@ -845,9 +843,7 @@ void CFVMFlowSolverBase<V, R>::BC_Sym_Plane(CGeometry* geometry, CSolver** solve
845843

846844
/*--- Compute unit normal, to be used for unit tangential, projected velocity and velocity
847845
component gradients. ---*/
848-
Area = 0.0;
849-
for (iDim = 0; iDim < nDim; iDim++) Area += Normal[iDim] * Normal[iDim];
850-
Area = sqrt(Area);
846+
Area = GeometryToolbox::Norm(nDim, Normal);
851847

852848
for (iDim = 0; iDim < nDim; iDim++) UnitNormal[iDim] = -Normal[iDim] / Area;
853849

@@ -1794,11 +1790,7 @@ void CFVMFlowSolverBase<V, FlowRegime>::Momentum_Forces(const CGeometry* geometr
17941790
Normal = geometry->vertex[iMarker][iVertex]->GetNormal();
17951791
Coord = geometry->nodes->GetCoord(iPoint);
17961792
Density = nodes->GetDensity(iPoint);
1797-
1798-
Area = 0.0;
1799-
for (iDim = 0; iDim < nDim; iDim++) Area += Normal[iDim] * Normal[iDim];
1800-
Area = sqrt(Area);
1801-
1793+
Area = GeometryToolbox::Norm(nDim, Normal);
18021794
MassFlow = 0.0;
18031795
su2double Velocity[MAXNDIM] = {0.0}, MomentDist[MAXNDIM] = {0.0};
18041796
for (iDim = 0; iDim < nDim; iDim++) {
@@ -2180,10 +2172,7 @@ void CFVMFlowSolverBase<V, FlowRegime>::Friction_Forces(const CGeometry* geometr
21802172
Grad_PrimVar = nodes->GetGradient_Primitive(iPoint);
21812173
}
21822174

2183-
Area = 0.0;
2184-
for (iDim = 0; iDim < nDim; iDim++) Area += Normal[iDim] * Normal[iDim];
2185-
Area = sqrt(Area);
2186-
2175+
Area = GeometryToolbox::Norm(nDim, Normal);
21872176
for (iDim = 0; iDim < nDim; iDim++) {
21882177
UnitNormal[iDim] = Normal[iDim] / Area;
21892178
}

‎SU2_CFD/src/interfaces/cht/CConjugateHeatInterface.cpp‎

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -75,11 +75,9 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
7575

7676
for (iDim = 0; iDim < nDim; iDim++) {
7777
Edge_Vector[iDim] = Coord_Normal[iDim] - Coord[iDim];
78-
dist2 += Edge_Vector[iDim]*Edge_Vector[iDim];
79-
Area += Normal[iDim]*Normal[iDim];
8078
}
81-
dist = sqrt(dist2);
82-
Area = sqrt(Area);
79+
dist = GeometryToolbox::Norm(nDim, Edge_Vector);
80+
Area = GeometryToolbox::Norm(nDim, Normal);
8381

8482
/*--- Retrieve temperature solution and its gradient ---*/
8583

‎SU2_CFD/src/numerics/NEMO/NEMO_diffusion.cpp‎

Lines changed: 2 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -90,10 +90,7 @@ CNumerics::ResidualType<> CAvgGrad_NEMO::ComputeResidual(const CConfig *config)
9090
unsigned short iSpecies, iVar, iDim;
9191

9292
/*--- Normalized normal vector ---*/
93-
Area = 0;
94-
for (iDim = 0; iDim < nDim; iDim++)
95-
Area += Normal[iDim]*Normal[iDim];
96-
Area = sqrt(Area);
93+
Area = GeometryToolbox::Norm(nDim, Normal);
9794

9895
for (iDim = 0; iDim < nDim; iDim++)
9996
UnitNormal[iDim] = Normal[iDim]/Area;
@@ -238,10 +235,7 @@ CNumerics::ResidualType<> CAvgGradCorrected_NEMO::ComputeResidual(const CConfig
238235
su2double dist_ij_2;
239236

240237
/*--- Normalized normal vector ---*/
241-
Area = 0;
242-
for (iDim = 0; iDim < nDim; iDim++)
243-
Area += Normal[iDim]*Normal[iDim];
244-
Area = sqrt(Area);
238+
Area = GeometryToolbox::Norm(nDim, Normal);
245239

246240
for (iDim = 0; iDim < nDim; iDim++)
247241
UnitNormal[iDim] = Normal[iDim]/Area;

‎SU2_CFD/src/numerics/NEMO/convection/ausm.cpp‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -70,10 +70,7 @@ CNumerics::ResidualType<> CUpwAUSM_NEMO::ComputeResidual(const CConfig *config)
7070
e_ve_i, e_ve_j, mL, mR, mLP, mRM, mF, pLP, pRM, pF, Phi;
7171

7272
/*--- Compute geometric quantities ---*/
73-
Area = 0;
74-
for (iDim = 0; iDim < nDim; iDim++)
75-
Area += Normal[iDim]*Normal[iDim];
76-
Area = sqrt(Area);
73+
Area = GeometryToolbox::Norm(nDim, Normal);
7774

7875
for (iDim = 0; iDim < nDim; iDim++)
7976
UnitNormal[iDim] = Normal[iDim]/Area;

‎SU2_CFD/src/numerics/NEMO/convection/ausmplusup2.cpp‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -109,10 +109,7 @@ CNumerics::ResidualType<> CUpwAUSMPLUSUP2_NEMO::ComputeResidual(const CConfig *c
109109
e_ve_i, e_ve_j, mL, mR, mLP, mRM, mF, pLP, pRM, pF, Phi, sq_veli, sq_velj;
110110

111111
/*--- Face area ---*/
112-
Area = 0.0;
113-
for (iDim = 0; iDim < nDim; iDim++)
114-
Area += Normal[iDim]*Normal[iDim];
115-
Area = sqrt(Area);
112+
Area = GeometryToolbox::Norm(nDim, Normal);
116113

117114
/*-- Unit Normal ---*/
118115
for (iDim = 0; iDim < nDim; iDim++)

‎SU2_CFD/src/numerics/NEMO/convection/ausmpwplus.cpp‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -75,10 +75,7 @@ CNumerics::ResidualType<> CUpwAUSMPWplus_NEMO::ComputeResidual(const CConfig *co
7575
Flux[iVar] = 0.0;
7676

7777
/*--- Calculate geometric quantities ---*/
78-
Area = 0.0;
79-
for (iDim = 0; iDim < nDim; iDim++)
80-
Area += Normal[iDim]*Normal[iDim];
81-
Area = sqrt(Area);
78+
Area = GeometryToolbox::Norm(nDim, Normal);
8279

8380
for (iDim = 0; iDim < nDim; iDim++)
8481
UnitNormal[iDim] = Normal[iDim]/Area;

0 commit comments

Comments
 (0)