Skip to content
Snippets Groups Projects
Commit d91d2cf2 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

Remove old hack around the ambiguity in the interface between 3-vectors and skew-symmetric matrices

[[Imported from SVN: r8861]]
parent aaf4108c
No related branches found
No related tags found
No related merge requests found
......@@ -107,15 +107,7 @@ void TargetSpaceRiemannianTRSolver<TargetSpace>::solve()
// ////////////////////////////////////////////////////
TargetSpace newIterate = x_;
#if 0 // out-commented until the Rotation class can distinguish skew-symmetric matrices from three-vectors
newIterate = TargetSpace::exp(newIterate, corr[0]);
#else
Dune::FieldMatrix<double,TargetSpace::TangentVector::dimension,TargetSpace::EmbeddedTangentVector::dimension> B = x_.orthonormalFrame();
Dune::FieldVector<double,TargetSpace::EmbeddedTangentVector::dimension> embeddedCorr(0);
B.mtv(corr[0], embeddedCorr);
newIterate = TargetSpace::exp(newIterate, embeddedCorr);
#endif
/** \todo Don't always recompute oldEnergy */
double oldEnergy = assembler_->value(x_);
......
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