yield_slip.h

The yield slip problem – class header.

The yield slip problem – class header

class yield_slip {
public:
  typedef field value_type;
  typedef Float float_type;
  yield_slip (Float S, Float n, Float Cf, Float r,
    const geo& omega, const geo& boundary, string approx = "P1");
  field residue          (const field& beta_h) const;
  void update_derivative (const field& beta_h) const;
  field derivative_solve      (const field& mrh) const;
  field derivative_trans_mult (const field& mrh) const;
  Float space_norm      (const field&) const;
  Float dual_space_norm (const field&) const;
  field initial () const;
  void post (const field& beta_h, field& uh, field& lambda_h) const;
protected:
  Float S, n, Cf, r;
  geo boundary;
  space Xh, Wh, Yh;
  field lh, mkh;
  form m, mb, a, b;
  mutable form c1;
  problem pmb, pa;
  mutable problem pA;
};
#include "yield_slip1.icc"
#include "yield_slip2.icc"