tboudreaux/SERiF
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

style(NonlinearPowerIntegrator): removed old comments

Browse Source
This commit is contained in:
Emily Boudreaux
2025-05-11 14:38:50 -04:00
parent 3d33839028
commit 8b9d46e996
1 changed files with 1 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/poly/utils/private/integrators.cpp
View File

@@ -70,10 +70,6 @@ namespace polyMFEMUtils {
} else { } else {
u_nl = std::pow(u_val, m_polytropicIndex); u_nl = std::pow(u_val, m_polytropicIndex);
} }
// const double u_safe = std::max(u_val, 0.0);
// const double u_nl = std::pow(u_safe, m_polytropicIndex);
//
// const double x2_u_nl = u_nl;
for (int i = 0; i < dof; i++){ for (int i = 0; i < dof; i++){
elvect(i) += shape(i) * u_nl * weight; elvect(i) += shape(i) * u_nl * weight;
@@ -121,10 +117,7 @@ namespace polyMFEMUtils {
} else { } else {
d_u_nl = m_polytropicIndex * std::pow(u_val, m_polytropicIndex - 1.0); d_u_nl = m_polytropicIndex * std::pow(u_val, m_polytropicIndex - 1.0);
} }
// const double u_safe = std::max(u_val, 0.0);
// const double d_u_nl = m_polytropicIndex * std::pow(u_safe, m_polytropicIndex - 1);
// const double x2_d_u_nl = d_u_nl;
for (int i = 0; i < dof; i++) { for (int i = 0; i < dof; i++) {
for (int j = 0; j < dof; j++) { for (int j = 0; j < dof; j++) {
elmat(i, j) += shape(i) * d_u_nl * shape(j) * weight; elmat(i, j) += shape(i) * d_u_nl * shape(j) * weight;