Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
yield_slip2.icc
Go to the documentation of this file.
1
25#include "poisson_robin.icc"
27 field mryh(Yh, 0.);
28 mryh[1] = -mrh;
29 field delta_yh(Yh);
30 pA.solve(mryh, delta_yh);
31 return delta_yh[1];
32}
34 field rh (Wh);
35 pmb.solve (mrh, rh);
36 field rhs = b.trans_mult(rh);
37 field delta_vh (Xh, 0.);
38 pa.solve (rhs, delta_vh);
39 field mgh = b*delta_vh + c1*rh;
40 field gh (Wh);
41 pmb.solve (mgh, gh);
42 return gh;
43}
45 return sqrt (mb(rh,rh));
46}
48 field rh (Wh,0);
49 pmb.solve (mrh, rh);
50 return sqrt (dual (mrh,rh));
51}
54 return (Cf+r)*uh["boundary"];
55}
56void yield_slip::post (const field& beta_h, field& uh, field& lambda_h) const {
57 field rhs = lh - b.trans_mult(beta_h);
58 uh = field (Xh, 0.);
59 pa.solve (rhs, uh);
60 lambda_h = beta_h - r*uh["boundary"];
61}
field gh(Float epsilon, Float t, const field &uh, const test &v)
see the Float page for the full documentation
see the field page for the full documentation
field derivative_trans_mult(const field &mrh) const
problem pa
Definition yield_slip.h:46
Float dual_space_norm(const field &) const
problem pA
Definition yield_slip.h:47
Float space_norm(const field &) const
field derivative_solve(const field &mrh) const
problem pmb
Definition yield_slip.h:46
field initial() const
void post(const field &beta_h, field &uh, field &lambda_h) const
The Poisson problem with Robin boundary condition – solver function.
field poisson_robin(Float Cf, const geo &boundary, const field &lh)