From d1cbcb59c4b583b8e3405aa02ee717b4fb1cfb03 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Sat, 9 Jul 2016 22:24:02 +0200
Subject: [PATCH] Use a true second-order function as the base of the
second-order vector field
It seemed strange to me to demonstrate a second-order tangent field on a
first-order base function. People may think that this is a necessary restriction.
Also, this patch includes a bugfix: The basis of the tangent space was hard-wired
to (0,1,0) and (0,0,1). But this only works if we are at the tangent space at (1,0,0).
---
test/interillustration.cc | 14 ++++++++------
1 file changed, 8 insertions(+), 6 deletions(-)
diff --git a/test/interillustration.cc b/test/interillustration.cc
index 17deabe4..89a01645 100644
--- a/test/interillustration.cc
+++ b/test/interillustration.cc
@@ -106,11 +106,13 @@ void interpolate(const LocalFEFunctionType& localGeodesicFEFunction,
{
FieldMatrix<double,3,3> derivative;
localGeodesicFEFunction.evaluateDerivativeOfValueWRTCoefficient(v.geometry().corner(0),
- partialDerivative, // derive wrt value at first triangle corner
+ partialDerivative, // select the Lagrange node
derivative);
- variation0[gridView.indexSet().index(v)] = derivative[1];
- variation1[gridView.indexSet().index(v)] = derivative[2];
+ Dune::FieldMatrix<double,2,3> basis = localGeodesicFEFunction.coefficient(partialDerivative).orthonormalFrame();
+
+ derivative.mtv(basis[0], variation0[gridView.indexSet().index(v)]);
+ derivative.mtv(basis[1], variation1[gridView.indexSet().index(v)]);
}
// sample a checkerboard pattern for nicer pictures
@@ -191,10 +193,10 @@ int main(int argc, char* argv[]) try
coefficients.resize(6);
coefficients[0] = TargetSpace(FieldVector<double,3>({1,0,0}));
- coefficients[1] = TargetSpace(FieldVector<double,3>({0.5,0.5,0}));
+ coefficients[1] = TargetSpace(FieldVector<double,3>({0.5,0.5,0.15}));
coefficients[2] = TargetSpace(FieldVector<double,3>({0,1,0}));
- coefficients[3] = TargetSpace(FieldVector<double,3>({0.5,0,0.5}));
- coefficients[4] = TargetSpace(FieldVector<double,3>({0,0.5,0.5}));
+ coefficients[3] = TargetSpace(FieldVector<double,3>({0.5,0.15,0.5}));
+ coefficients[4] = TargetSpace(FieldVector<double,3>({0.15,0.5,0.5}));
coefficients[5] = TargetSpace(FieldVector<double,3>({0,0,1}));
typedef LocalGeodesicFEFunction<dim, double, P2LocalFiniteElement<double,double,dim>, TargetSpace> P2LocalGFEFunctionType;
--
GitLab