*** Wartungsfenster jeden ersten Mittwoch vormittag im Monat ***

Skip to content
Snippets Groups Projects
Commit c29bc9e2 authored by Innerberger, Michael's avatar Innerberger, Michael
Browse files

Fix bug in Helmholtz linear form (missing coefficient)

parent 7c032c6c
No related branches found
No related tags found
No related merge requests found
...@@ -116,7 +116,7 @@ function [lf, perturbation] = setUpHelmholtzProblem(fes, uold, coeff) ...@@ -116,7 +116,7 @@ function [lf, perturbation] = setUpHelmholtzProblem(fes, uold, coeff)
lf = LinearForm(fes); lf = LinearForm(fes);
perturbation = BilinearForm(fes); perturbation = BilinearForm(fes);
lf.f = CompositeFunction(@(u) 1 + u, uold); lf.f = CompositeFunction(@(u) 1 + coeff^2*u, uold);
lf.qrf = QuadratureRule.ofOrder(2); lf.qrf = QuadratureRule.ofOrder(2);
perturbation.c = Constant(fes.mesh, -coeff^2); perturbation.c = Constant(fes.mesh, -coeff^2);
perturbation.qrc = QuadratureRule.ofOrder(2); perturbation.qrc = QuadratureRule.ofOrder(2);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment