mosolov_augmented_lagrangian2.icc

Go to the documentation of this file.

mosolov_augmented_lagrangian::mosolov_augmented_lagrangian()
: Bi(0), n(1), r(1), tol(1e-10), max_iter(1000000), 
  Xh(), Th(), lh(), a(), b(), inv_mt(), pa()
{}
void mosolov_augmented_lagrangian::reset (geo omega, string approx) {
  Xh = space (omega, approx);
  Xh.block ("boundary");
  string grad_approx = "P" + to_string(Xh.degree()-1) + "d";
  Th = space (omega, grad_approx, "vector");
  trial u (Xh), sigma(Th);
  test  v (Xh), tau  (Th);
  lh = integrate(2*v);
  a = integrate (dot(grad(u),grad(v)));
  b = integrate (dot(grad(u),tau));
  integrate_option iopt;
  iopt.invert = true;
  inv_mt = integrate(dot(sigma,tau), iopt);
  pa = problem (a);
}
void
mosolov_augmented_lagrangian::initial (field& sigma_h, field& uh) const {
  uh = field(Xh);
  uh ["boundary"] = 0;
  pa.solve (lh, uh);
  test tau (Th);
  field mt_grad_uh = integrate(dot(grad(uh),tau)); 
  sigma_h = inv_mt*mt_grad_uh;
}
void mosolov_augmented_lagrangian::put (odiststream& out, 
    field& sigma_h, field& uh) const
{
  out << catchmark("Bi")    << Bi << endl
      << catchmark("n")     << n  << endl
      << catchmark("r")     << r  << endl
      << catchmark("sigma") << sigma_h
      << catchmark("u")     << uh;
}