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

fix(PolytropeOperator): uncommented code that should never have been commented in GetGradient

Browse Source
This commit is contained in:
Emily Boudreaux
2025-05-11 15:09:22 -04:00
parent bee8f693c1
commit 14eb21bd31
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/poly/utils/private/operator.cpp
View File

@@ -267,9 +267,10 @@ mfem::Operator& PolytropeOperator::GetGradient(const mfem::Vector &x) const {
// auto *gradPtr = &grad;
// updatePreconditioner(grad);
// // TODO: Save the state vector somewhere so that I can inspect how the actual solution
// // to theta and phi are evolving
// m_jacobian->SetBlock(0, 0, &grad);
// TODO: Save the state vector somewhere so that I can inspect how the actual solution
// to theta and phi are evolving
m_jacobian->SetBlock(0, 0, &grad);
// // TODO: Save the jacobian here
// std::cout << "Writing out jacobian to file..." << std::endl;
// std::vector<mfem::Operator*> ops;
@@ -280,6 +281,7 @@ mfem::Operator& PolytropeOperator::GetGradient(const mfem::Vector &x) const {
// saveBlockFormToBinary(ops, {{0, 0}, {0, 1}, {1, 0}, {1, 1}}, {false, false, false, false}, "jacobian_" + std::to_string(s_newtonStepGrad) + ".bin");
// s_newtonStepGrad++;
// std::cout << "Done writing out jacobian to file." << std::endl;
return *m_jacobian;
}