Skip to content
Snippets Groups Projects
Commit eaa5fa0f authored by Dr. Felix Tobias Schindler's avatar Dr. Felix Tobias Schindler
Browse files

use new integrands

parent e722bf20
No related branches found
No related tags found
No related merge requests found
......@@ -27,9 +27,9 @@
#include <dune/gdt/functionals/vector-based.hh>
#include <dune/gdt/local/functionals/integrals.hh>
#include <dune/gdt/local/bilinear-forms/integrals.hh>
#include <dune/gdt/local/integrands/elliptic-ipdg.hh>
#include <dune/gdt/local/integrands/laplace.hh>
#include <dune/gdt/local/integrands/laplace-ipdg.hh>
#include <dune/gdt/local/integrands/ipdg.hh>
#include <dune/gdt/local/integrands/elliptic.hh>
#include <dune/gdt/local/integrands/product.hh>
#include <dune/gdt/local/integrands/conversion.hh>
#include <dune/gdt/norms.hh>
......@@ -265,10 +265,10 @@ protected:
(space_type_.size() >= 2 && space_type_.substr(0, 2) == "cg") ? Stencil::element
: Stencil::element_and_intersection));
// - volume term
lhs_op->append(LocalElementIntegralBilinearForm<E>(LocalEllipticIntegrand<E>(one, this->diffusion())));
lhs_op->append(LocalElementIntegralBilinearForm<E>(LocalLaplaceIntegrand<E>(this->diffusion())));
// - inner faces
lhs_op->append(LocalIntersectionIntegralBilinearForm<I>(
LocalEllipticIPDGIntegrands::InnerCoupling<I>(1., this->diffusion(), this->diffusion())),
LocalLaplaceIPDGIntegrands::InnerCoupling<I>(1., this->diffusion(), this->diffusion())),
{},
XT::Grid::ApplyOn::InnerIntersectionsOnce<GV>());
lhs_op->append(
......@@ -279,7 +279,7 @@ protected:
// - Dirichlet faces
lhs_op->append(
LocalIntersectionIntegralBilinearForm<I>(
LocalEllipticIPDGIntegrands::DirichletCoupling<I>(1., this->diffusion())),
LocalLaplaceIPDGIntegrands::DirichletCoupling<I>(1., this->diffusion())),
{},
XT::Grid::ApplyOn::CustomBoundaryIntersections<GV>(this->boundary_info(), new XT::Grid::DirichletBoundary()));
lhs_op->append(
......
......@@ -24,7 +24,7 @@
#include <dune/xt/grid/type_traits.hh>
#include <dune/gdt/local/bilinear-forms/integrals.hh>
#include <dune/gdt/local/integrands/elliptic.hh>
#include <dune/gdt/local/integrands/laplace.hh>
#include <dune/gdt/local/integrands/product.hh>
#include <dune/gdt/spaces/interface.hh>
......@@ -43,9 +43,9 @@ double estimate_inverse_inequality_constant(const SpaceInterface<GV, r>& space)
basis->bind(element);
const double h = XT::Grid::diameter(element);
auto H1_product_matrix = XT::LA::convert_to<XT::LA::CommonDenseMatrix<double>>(
LocalElementIntegralBilinearForm<E, r>(LocalEllipticIntegrand<E, r>(1.)).apply2(*basis, *basis));
LocalElementIntegralBilinearForm<E, r>(LocalLaplaceIntegrand<E, r>()).apply2(*basis, *basis));
auto L2_product_matrix = XT::LA::convert_to<XT::LA::CommonDenseMatrix<double>>(
LocalElementIntegralBilinearForm<E, r>(LocalElementProductIntegrand<E, r>(1.)).apply2(*basis, *basis));
LocalElementIntegralBilinearForm<E, r>(LocalElementProductIntegrand<E, r>()).apply2(*basis, *basis));
auto evs =
XT::LA::make_generalized_eigen_solver(H1_product_matrix,
L2_product_matrix,
......
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