Skip to content
Snippets Groups Projects
Commit 0b54ed95 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Bug fix: Truncate coarse correction!

parent 52fb127a
No related branches found
No related tags found
No related merge requests found
...@@ -76,7 +76,11 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem { ...@@ -76,7 +76,11 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
Linearization const &linearization) const { Linearization const &linearization) const {
// TODO: implement (not urgent) // TODO: implement (not urgent)
projected_v = v; projected_v = v;
}; for (size_t i = 0; i < v.size(); ++i)
for (int j = 0; j < block_size; ++j)
if (linearization.truncation[i][j])
projected_v[i][j] = 0;
}
double computeDampingParameter(VectorType const &u, double computeDampingParameter(VectorType const &u,
VectorType const &projected_v) const { VectorType const &projected_v) const {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment