From b1a168094de913b9d9b532a96eb8e3c968295db8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Mon, 16 Apr 2018 17:51:03 +0200
Subject: [PATCH] HarmonicEnergyLocalStiffness needs a LocalInterpolationRule
nowaday
---
src/gradient-flow.cc | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/gradient-flow.cc b/src/gradient-flow.cc
index c630fb81..187b9c34 100644
--- a/src/gradient-flow.cc
+++ b/src/gradient-flow.cc
@@ -210,7 +210,8 @@ int main (int argc, char *argv[]) try
auto l2DistanceSquaredEnergy = std::make_shared<L2DistanceSquaredEnergy<FEBasis, ATargetSpace> >();
std::vector<std::shared_ptr<LocalGeodesicFEStiffness<FEBasis,ATargetSpace> > > addends(2);
- addends[0] = std::make_shared<HarmonicEnergyLocalStiffness<FEBasis, ATargetSpace> >();
+ using GeodesicInterpolationRule = LocalGeodesicFEFunction<dim, double, FEBasis::LocalView::Tree::FiniteElement, ATargetSpace>;
+ addends[0] = std::make_shared<HarmonicEnergyLocalStiffness<FEBasis, GeodesicInterpolationRule, ATargetSpace> >();
addends[1] = l2DistanceSquaredEnergy;
std::vector<double> weights = {1.0, 1.0/(2*timeStepSize)};
--
GitLab