dirichlet-nh.cc

The Poisson problem with non-homogeneous Dirichlet boundary conditions.

The Poisson problem with non-homogeneous Dirichlet boundary conditions

#include "rheolef.h"
using namespace rheolef;
using namespace std;
#include "cosinusprod_laplace.h"
int main(int argc, char**argv) {
  environment rheolef(argc, argv);
  geo omega (argv[1]);
  size_t d = omega.dimension();
  space Xh (omega, argv[2]);
  Xh.block ("boundary");
  trial u (Xh); test v (Xh);
  form  a  = integrate (dot(grad(u),grad(v)));
  field lh = integrate (f(d)*v);
  field uh (Xh);
  space Wh (omega["boundary"], argv[2]);
  uh ["boundary"] = lazy_interpolate(Wh, g(d));
  problem p (a);
  p.solve (lh, uh);
  dout << uh;
}