From: Wolfgang Bangerth Date: Fri, 8 Jul 2016 16:42:12 +0000 (-0500) Subject: Ensure that charts do not map from higher into lower dimensional spaces. X-Git-Tag: v8.5.0-rc1~884^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3743d91a7512984598dcacd903db003df6eddecf;p=dealii.git Ensure that charts do not map from higher into lower dimensional spaces. --- 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