From d4fcd896a9cb52ca5e47948f4818c4803b563056 Mon Sep 17 00:00:00 2001
From: Jonathan Youett <youett@math.fu-berlin.de>
Date: Fri, 2 Mar 2018 14:57:23 +0100
Subject: [PATCH] Use correct field_type
---
dune/gfe/geodesicfeassembler.hh | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/geodesicfeassembler.hh b/dune/gfe/geodesicfeassembler.hh
index 36640199..7be9133a 100644
--- a/dune/gfe/geodesicfeassembler.hh
+++ b/dune/gfe/geodesicfeassembler.hh
@@ -15,6 +15,7 @@
template <class Basis, class TargetSpace>
class GeodesicFEAssembler {
+ using field_type = typename TargetSpace::field_type;
typedef typename Basis::GridView GridView;
using LocalStiffness = LocalGeodesicFEStiffness<Basis, TargetSpace>;
@@ -78,7 +79,7 @@ public:
* anyway to compute the Riemannian Hessian.
*/
virtual void assembleGradientAndHessian(const std::vector<TargetSpace>& sol,
- Dune::BlockVector<Dune::FieldVector<double, blocksize> >& gradient,
+ Dune::BlockVector<Dune::FieldVector<field_type, blocksize> >& gradient,
Dune::BCRSMatrix<MatrixBlock>& hessian,
bool computeOccupationPattern=true) const;
@@ -134,7 +135,7 @@ getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
template <class Basis, class TargetSpace>
void GeodesicFEAssembler<Basis,TargetSpace>::
assembleGradientAndHessian(const std::vector<TargetSpace>& sol,
- Dune::BlockVector<Dune::FieldVector<double, blocksize> >& gradient,
+ Dune::BlockVector<Dune::FieldVector<field_type, blocksize> > &gradient,
Dune::BCRSMatrix<MatrixBlock>& hessian,
bool computeOccupationPattern) const
{
--
GitLab