![]() |
![]() |
![]() |
In this demo, we consider a problem derived from the Laplace equation, namely the Steklov-Poincaré problem. The problem is the following: Given a function g defined on the boundary of the domain. We look for the normal derivative of the solution of the Laplace equation whose trace on the boundary is g. The problem is solved by an integral representation by using P0 boundary elements. In the OFELI library contains the class SteklovPoincare2DBE for such a purpose.
Vect |
g.setNodeBC(1,"0."); g.setNodeBC(2,"1."); |
SteklovPoincare2DBE eq(ms,u); eq.setInput(BOUNDARY_CONDITION,g); |
eq.run(); cout << u; cout << "L2 solution norm: " << u.Norm(WNORM2) << endl; |
![]() |
![]() |
![]() |