From 9a02ce7ccba2c45eef1a657925af7e32f1247f9b Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 8 Jul 2013 13:45:44 +0000
Subject: [PATCH] Remove extra transformations from the geodesic update
These were only there to hack our way around a bug in the
interface of the Rotation class. This bug has been fixed
quite a while ago.
[[Imported from SVN: r9293]]
---
dune/gfe/riemanniantrsolver.cc | 12 ------------
1 file changed, 12 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 6c1bae0e..12203ab2 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -396,20 +396,8 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
// ////////////////////////////////////////////////////
SolutionType newIterate = x_;
-#if 0 // out-commented until the Rotation class can distinguish skew-symmetric matrices from three-vectors
for (int j=0; j<newIterate.size(); j++)
newIterate[j] = TargetSpace::exp(newIterate[j], corr[j]);
-#else
- //std::cout << "embedded correction:\n";
- for (size_t j=0; j<newIterate.size(); j++) {
- Dune::FieldMatrix<double,TargetSpace::TangentVector::dimension,TargetSpace::EmbeddedTangentVector::dimension> B = x_[j].orthonormalFrame();
- Dune::FieldVector<double,TargetSpace::EmbeddedTangentVector::dimension> embeddedCorr(0);
- //std::cout << "B[" << j << "]:\n" << B << std::endl;
- B.mtv(corr[j], embeddedCorr);
- newIterate[j] = TargetSpace::exp(newIterate[j], embeddedCorr);
- //std::cout << embeddedCorr << " " << newIterate[j] << std::endl;
- }
-#endif
double energy = assembler_->computeEnergy(newIterate);
--
GitLab