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

Dead code

parent 7e286c1e
No related branches found
No related tags found
No related merge requests found
...@@ -49,26 +49,6 @@ template <int dim> class EllipticEnergy { ...@@ -49,26 +49,6 @@ template <int dim> class EllipticEnergy {
} }
} }
// returns false if the direction is tangential
void descentDirectionNew(SmallVector const &x, SmallVector &ret) const {
SmallVector pg;
gradient(x, pg);
SmallVector mg;
gradient(x, mg);
double const pgx = pg * x;
double const mgx = mg * x;
if (pgx >= 0 && mgx >= 0) {
ret = pg;
ret *= -1;
} else if (pgx <= 0 && mgx <= 0) {
ret = mg;
ret *= -1;
} else {
// We know that no descent can be expected in the tangential direction
ret = SmallVector(0);
}
}
bool descentDirection(SmallVector const &x, SmallVector &ret) const { bool descentDirection(SmallVector const &x, SmallVector &ret) const {
// Check the squared norm rather than each component because // Check the squared norm rather than each component because
// complementaryProjection() divides by it // complementaryProjection() divides by it
......
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