Skip to content
Snippets Groups Projects
Commit bffe0873 authored by Boris Martin's avatar Boris Martin
Browse files

further doc

parent 8b26c10d
No related branches found
No related tags found
1 merge request!17New tutorial case: overlapping DDM (Dirichlet) for Laplace equation
......@@ -6,8 +6,6 @@
#include <gmshfem/Message.h>
using namespace gmshddm;
using namespace gmshddm::common;
using namespace gmshddm::domain;
......@@ -78,11 +76,13 @@ int main(int argc, char **argv)
formulation.addInterfaceField(g);
// First domain
// Laplacian(u) = 0 on the domain, and u matches g(1,0) on the interface
formulation(0).integral(grad(dof(u(0))), grad(tf(u(0))), omega(0), gauss);
formulation(0).integral((-g(1, 0)), tf(u(0)), interface2, gauss);
formulation(0).integral(dof(u(0)), tf(u(0)), interface2, gauss);
// First interface (What 0 sends to 1)
// Compute g(0,1) such that g(0,1) + g(1,0) = 2 u(0)
formulation(0, 1).integral(-dof(g(0, 1)), tf(g(0, 1)), sigma(0, 1), gauss);
formulation(0, 1).integral(-(g(1, 0)), tf(g(0, 1)), sigma(0, 1), gauss);
formulation(0, 1).integral((2*u(0)), tf(g(0, 1)), sigma(0, 1), gauss);
......@@ -108,7 +108,6 @@ int main(int argc, char **argv)
if(isItMySubdomain(i)) {
save(u(i), omega(i), "u_" + std::to_string(i));
save(u(i) - trueSolution, omega(i), "err_" + std::to_string(i));
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment