Commit 4f0a5026 authored by Oliver Sander's avatar Oliver Sander
Browse files

[bugfix] The directionalSubDiff method was simply wrong

Another hurray for unit tests!
parent 766a697e
......@@ -118,10 +118,8 @@ namespace Dune {
subdifferential[0] += coefficients_[row] * (u[row]*v[row]) / norm;
subdifferential[1] += coefficients_[row] * (u[row]*v[row]) / norm;
} else {
for (int j=0; j<blockSize; j++) {
subdifferential[0] -= coefficients_[row] * v[row][j];
subdifferential[1] += coefficients_[row] * v[row][j];
}
subdifferential[0] -= coefficients_[row] * v[row].two_norm();
subdifferential[1] += coefficients_[row] * v[row].two_norm();
}
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment