laplace_s.cc

Go to the documentation of this file.

#include "rheolef.h"
using namespace rheolef;
using namespace std;
#include "torus.icc"
int main (int argc, char**argv) {
  environment rheolef (argc, argv);
  geo gamma (argv[1]);
  size_t d = gamma.dimension();
  space Wh (gamma, argv[2]);
  trial u (Wh); test v (Wh);
  form  a = integrate (dot(grad_s(u),grad_s(v)));
  field b = integrate(v);
  field lh = integrate (f(d)*v);
  form  A = {{    a,     b },
             { trans(b), 0 }};
  field Bh = {   lh,     0 };
  field Uh (Bh.get_space(), 0);
  A.set_symmetry(true);
  problem pa (A);
  pa.solve (Bh, Uh);
  dout << Uh[0];
}