Skip to content
Snippets Groups Projects
Commit 2d50f53a authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

abort when requesting Lagrangenodes for d>1, p!=2. Not implemented yet

[[Imported from SVN: r8868]]
parent 12cd2be9
No related branches found
No related tags found
No related merge requests found
...@@ -47,6 +47,7 @@ std::vector<FieldVector<double,domainDim> > lagrangeNodes(const GeometryType& ty ...@@ -47,6 +47,7 @@ std::vector<FieldVector<double,domainDim> > lagrangeNodes(const GeometryType& ty
return result; return result;
} }
assert(elementOrder==2);
PQkLocalFiniteElementCache<double,double,domainDim,elementOrder> feCache; PQkLocalFiniteElementCache<double,double,domainDim,elementOrder> feCache;
result.resize(feCache.get(type).localBasis().size()); result.resize(feCache.get(type).localBasis().size());
...@@ -153,7 +154,7 @@ void test(const GeometryType& element) ...@@ -153,7 +154,7 @@ void test(const GeometryType& element)
LocalGeodesicFEFunction<domainDim,double,LocalFiniteElement,TargetSpace> f(feCache.get(element),corners); LocalGeodesicFEFunction<domainDim,double,LocalFiniteElement,TargetSpace> f(feCache.get(element),corners);
static const int highElementOrder = 2; static const int highElementOrder = 3;
testNestedness<domainDim,TargetSpace,highElementOrder>(f); testNestedness<domainDim,TargetSpace,highElementOrder>(f);
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment