Skip to content

Commit

Permalink
Changing IntegrationValues::evaluateEverything to not recompute Jacob…
Browse files Browse the repository at this point in the history
…ian values at every request
  • Loading branch information
seanofthemillers committed Apr 4, 2024
1 parent fbdc0df commit 19ab932
Showing 1 changed file with 18 additions and 17 deletions.
35 changes: 18 additions & 17 deletions packages/panzer/disc-fe/src/Panzer_IntegrationValues2.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -1793,43 +1793,44 @@ evaluateEverything()
const bool is_cv = (int_rule->getType() == ID::CV_VOLUME) or (int_rule->getType() == ID::CV_SIDE) or (int_rule->getType() == ID::CV_BOUNDARY);
const bool is_side = int_rule->isSide();

// This will force all values to be re-evaluated
resetArrays();

// Base cubature stuff
if(is_cv){
getCubaturePoints(true,true);
getCubaturePointsRef(true,true);
getCubaturePoints(true);
getCubaturePointsRef(true);
} else {
if(not is_surface){
getUniformCubaturePointsRef(true,true);
getUniformCubatureWeightsRef(true,true);
getUniformCubaturePointsRef(true);
getUniformCubatureWeightsRef(true);
if(is_side)
getUniformSideCubaturePointsRef(true,true);
getUniformSideCubaturePointsRef(true);
}
getCubaturePointsRef(true,true);
getCubaturePoints(true,true);
getCubaturePointsRef(true);
getCubaturePoints(true);
}

// Measure stuff
getJacobian(true,true);
getJacobianDeterminant(true,true);
getJacobianInverse(true,true);
getJacobian(true);
getJacobianDeterminant(true);
getJacobianInverse(true);
if(int_rule->cv_type == "side")
getWeightedNormals(true,true);
getWeightedNormals(true);
else
getWeightedMeasure(true,true);
getWeightedMeasure(true);

// Surface stuff
if(is_surface){
getSurfaceNormals(true,true);
getSurfaceRotationMatrices(true,true);
getSurfaceNormals(true);
getSurfaceRotationMatrices(true);
}

// Stabilization stuff
if(not (is_surface or is_cv)){
getContravarientMatrix(true,true);
getCovarientMatrix(true,true);
getNormContravarientMatrix(true,true);
getContravarientMatrix(true);
getCovarientMatrix(true);
getNormContravarientMatrix(true);
}
}

Expand Down

0 comments on commit 19ab932

Please sign in to comment.