@@ -566,8 +566,10 @@ CMeshFEM::CMeshFEM(CGeometry* geometry, CConfig* config) {
566566 the elements with constant and non-constant Jacobians are
567567 considered the same. */
568568 if (JacConstant) {
569- const auto orderExactStraight = static_cast <unsigned short >(ceil (nPolySol * config->GetQuadrature_Factor_Straight ()));
570- const auto orderExactCurved = static_cast <unsigned short >(ceil (nPolySol * config->GetQuadrature_Factor_Curved ()));
569+ const auto orderExactStraight =
570+ static_cast <unsigned short >(ceil (nPolySol * config->GetQuadrature_Factor_Straight ()));
571+ const auto orderExactCurved =
572+ static_cast <unsigned short >(ceil (nPolySol * config->GetQuadrature_Factor_Curved ()));
571573 if (orderExactStraight == orderExactCurved) JacConstant = false ;
572574 }
573575
@@ -2455,12 +2457,16 @@ void CMeshFEM_DG::CreateFaces(CConfig* config) {
24552457 is set to false. Hence it is only needed to carry out this check for faces
24562458 with a constant Jacobian. This is done to reduce the number of standard elements. */
24572459 if (thisFace.JacFaceIsConsideredConstant ) {
2458- auto orderExactStraight = static_cast <unsigned short >(ceil (thisFace.nPolyGrid0 * config->GetQuadrature_Factor_Straight ()));
2459- auto orderExactCurved = static_cast <unsigned short >(ceil (thisFace.nPolyGrid0 * config->GetQuadrature_Factor_Curved ()));
2460+ auto orderExactStraight =
2461+ static_cast <unsigned short >(ceil (thisFace.nPolyGrid0 * config->GetQuadrature_Factor_Straight ()));
2462+ auto orderExactCurved =
2463+ static_cast <unsigned short >(ceil (thisFace.nPolyGrid0 * config->GetQuadrature_Factor_Curved ()));
24602464
24612465 if (orderExactStraight == orderExactCurved) {
2462- orderExactStraight = static_cast <unsigned short >(ceil (thisFace.nPolySol0 * config->GetQuadrature_Factor_Straight ()));
2463- orderExactCurved = static_cast <unsigned short >(ceil (thisFace.nPolySol0 * config->GetQuadrature_Factor_Curved ()));
2466+ orderExactStraight =
2467+ static_cast <unsigned short >(ceil (thisFace.nPolySol0 * config->GetQuadrature_Factor_Straight ()));
2468+ orderExactCurved =
2469+ static_cast <unsigned short >(ceil (thisFace.nPolySol0 * config->GetQuadrature_Factor_Curved ()));
24642470 if (orderExactStraight == orderExactCurved) thisFace.JacFaceIsConsideredConstant = false ;
24652471 }
24662472 }
0 commit comments