diff --git a/dune/gfe/targetspacertrsolver.cc b/dune/gfe/targetspacertrsolver.cc index d59f290067b0be7bba9456d1a340e3da1d4983c4..e05bcbf546fa7abb5d83c57dff61f901a381423c 100644 --- a/dune/gfe/targetspacertrsolver.cc +++ b/dune/gfe/targetspacertrsolver.cc @@ -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_);