From 3743d91a7512984598dcacd903db003df6eddecf Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Fri, 8 Jul 2016 11:42:12 -0500 Subject: [PATCH] Ensure that charts do not map from higher into lower dimensional spaces. --- tests/manifold/chart_manifold_05_embedded.cc | 14 ++++++------- tests/manifold/chart_manifold_06_embedded.cc | 22 +++++++++----------- 2 files changed, 17 insertions(+), 19 deletions(-) diff --git a/tests/manifold/chart_manifold_05_embedded.cc b/tests/manifold/chart_manifold_05_embedded.cc index 552854fd5e..a208dd2440 100644 --- a/tests/manifold/chart_manifold_05_embedded.cc +++ b/tests/manifold/chart_manifold_05_embedded.cc @@ -22,14 +22,14 @@ template -class MyFlatManifold : public ChartManifold +class MyFlatManifold : public ChartManifold { public: virtual - Point + Point pull_back(const Point &space_point) const { - Point p; + Point p; for (unsigned int d=0; d - push_forward(const Point &chart_point) const + push_forward(const Point &chart_point) const { Point p; for (unsigned int d=0; d - push_forward_gradient(const Point &chart_point) const + DerivativeForm<1,spacedim,spacedim> + push_forward_gradient(const Point &chart_point) const { - DerivativeForm<1,spacedim+1,spacedim> x; + DerivativeForm<1,spacedim,spacedim> x; for (unsigned int d=0; d @@ -23,19 +21,19 @@ template -class MyFlatManifold : public ChartManifold +class MyFlatManifold : public ChartManifold { public: - MyFlatManifold (const Tensor<1,spacedim+1> &periodicity) + MyFlatManifold (const Tensor<1,spacedim> &periodicity) : - ChartManifold(periodicity) + ChartManifold(periodicity) {} virtual - Point + Point pull_back(const Point &space_point) const { - Point p; + Point p; for (unsigned int d=0; d - push_forward(const Point &chart_point) const + push_forward(const Point &chart_point) const { Point p; for (unsigned int d=0; d - push_forward_gradient(const Point &chart_point) const + DerivativeForm<1,spacedim,spacedim> + push_forward_gradient(const Point &chart_point) const { - DerivativeForm<1,spacedim+1,spacedim> x; + DerivativeForm<1,spacedim,spacedim> x; for (unsigned int d=0; d