@@ -2462,7 +2462,6 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
24622462 /* --- Compute if marker is a straight line or plane. ---*/
24632463 ComputeSurf_Straightness (config, false );
24642464
2465-
24662465 /* --- Check which markers are symmetry or Euler wall and store in a list. ---*/
24672466 for (auto iMarker = 0u ; iMarker < nMarker; ++iMarker) {
24682467 if ((config->GetMarker_All_KindBC (iMarker) == SYMMETRY_PLANE) ||
@@ -2471,7 +2470,6 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
24712470 }
24722471 }
24732472
2474-
24752473 // first check Euler walls, we 'merge' all euler walls with other euler walls
24762474
24772475 /* --- Loop over all markers and find nodes on symmetry planes that are shared with other symmetries. ---*/
@@ -2499,9 +2497,8 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
24992497 const auto jVertex = nodes->GetVertex (iPoint, jMarker);
25002498 if (jVertex < 0 ) continue ;
25012499
2502-
2503- /* --- If both symmetries are curved, we sum the normals. ---*/
2504- if ((bound_is_straight[iMarker]) == true && (bound_is_straight[jMarker]) == true ) {
2500+ /* --- If both symmetries are curved, we sum the normals. ---*/
2501+ if ((bound_is_straight[iMarker]) == true && (bound_is_straight[jMarker]) == true ) {
25052502 isShared = true ;
25062503 }
25072504
@@ -2518,24 +2515,20 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
25182515 }
25192516
25202517 // averaging
2521- for (auto iDim = 0u ; iDim < nDim; iDim++)
2522- kNormal [iDim] = iNormal[iDim] + jNormal[iDim];
2518+ for (auto iDim = 0u ; iDim < nDim; iDim++) kNormal [iDim] = iNormal[iDim] + jNormal[iDim];
25232519 const su2double areak = GeometryToolbox::Norm (nDim, kNormal .data ());
25242520 for (auto iDim = 0u ; iDim < nDim; iDim++) kNormal [iDim] /= areak;
25252521
25262522 for (auto iDim = 0u ; iDim < nDim; iDim++) {
25272523 symmetryNormals[iMarker][iVertex][iDim] += kNormal [iDim];
25282524 symmetryNormals[jMarker][jVertex][iDim] += kNormal [iDim];
25292525 }
2530-
25312526 }
25322527
25332528 if (!isShared) continue ;
2534-
25352529 }
25362530 }
25372531
2538-
25392532 /* --- Now do Gramm-Schmidt Process for symmetries that are not both curved lines or planes ---*/
25402533
25412534 /* --- Loop over all markers and find nodes on symmetry planes that are shared with other symmetries. ---*/
@@ -2566,12 +2559,11 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
25662559 isShared = true ;
25672560
25682561 /* --- Both symmetries are curved, no further correction. ---*/
2569- if ((bound_is_straight[iMarker]) == false && (bound_is_straight[jMarker]) == false ) {
2562+ if ((bound_is_straight[iMarker]) == false && (bound_is_straight[jMarker]) == false ) {
25702563 isShared = false ;
25712564 continue ;
25722565 }
25732566
2574-
25752567 std::array<su2double, MAXNDIM> jNormal = {};
25762568 const auto it = symmetryNormals[jMarker].find (jVertex);
25772569
@@ -2599,7 +2591,6 @@ void CGeometry::ComputeModifiedSymmetryNormals(const CConfig* config) {
25992591 }
26002592 }
26012593 }
2602-
26032594}
26042595
26052596void CGeometry::ComputeSurf_Straightness (const CConfig* config, bool print_on_screen) {
0 commit comments