Commit 50b2c118 authored by Praetorius, Simon's avatar Praetorius, Simon
Browse files

some cleanup of convergence test

parent 7f8a2779
...@@ -20,7 +20,7 @@ ...@@ -20,7 +20,7 @@
#define STR_HELPER(s) #s #define STR_HELPER(s) #s
const int order = 4; const int order = 4;
const int quad_order = order+5; const int quad_order = order+10;
const int num_levels = 5; const int num_levels = 5;
// Test Hausdorf-distance // Test Hausdorf-distance
...@@ -94,7 +94,7 @@ typename Grid::ctype curvature_error(Grid const& grid, Projection const& project ...@@ -94,7 +94,7 @@ typename Grid::ctype curvature_error(Grid const& grid, Projection const& project
{ {
using GlobalCoordinate = typename Grid::template Codim<0>::Entity::Geometry::GlobalCoordinate; using GlobalCoordinate = typename Grid::template Codim<0>::Entity::Geometry::GlobalCoordinate;
using QuadProvider = Dune::QuadratureRules<typename Grid::ctype, 2>; using QuadProvider = Dune::QuadratureRules<typename Grid::ctype, 2>;
using FiniteElementType = Dune::PkLocalFiniteElement<typename Grid::ctype, typename Grid::ctype, 2, order+2>; using FiniteElementType = Dune::PkLocalFiniteElement<typename Grid::ctype, typename Grid::ctype, 2, order+4>;
FiniteElementType fe; FiniteElementType fe;
using LBTraits = typename FiniteElementType::Traits::LocalBasisType::Traits; using LBTraits = typename FiniteElementType::Traits::LocalBasisType::Traits;
...@@ -152,33 +152,6 @@ typename Grid::ctype edge_length(Grid const& grid) ...@@ -152,33 +152,6 @@ typename Grid::ctype edge_length(Grid const& grid)
return h; return h;
} }
// struct Sphere
// {
// struct DSphere
// {
// template <class K, int n>
// Dune::FieldMatrix<K,n,n> operator() (Dune::FieldVector<K,n> const& x) const
// {
// K nrmInv = K(1)/x.two_norm();
// Dune::FieldMatrix<K,n,n> out;
// for (int i = 0; i < n; ++i)
// for (int j = 0; j < n; ++j)
// out[i][j] = nrmInv * (i == j ? K(1) : K(0)) + x[i] * x[j];
// return out;
// }
// };
// template <class K, int n>
// Dune::FieldVector<K,n> operator() (Dune::FieldVector<K,n> const& x) const
// {
// return x / x.two_norm();
// }
// friend DSphere derivative(Sphere)
// {
// return DSphere{};
// }
// };
int main(int argc, char** argv) int main(int argc, char** argv)
{ {
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment