~ubuntu-branches/ubuntu/saucy/rheolef/saucy-proposed

« back to all changes in this revision

Viewing changes to doc/pexamples/elasticity_solve.icc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito
  • Date: 2012-04-06 09:12:21 UTC
  • mfrom: (1.1.5)
  • Revision ID: package-import@ubuntu.com-20120406091221-m58me99p1nxqui49
Tags: 6.0-1
* New upstream release 6.0 (major changes):
  - massively distributed and parallel support
  - full FEM characteristic method (Lagrange-Gakerkin method) support
  - enhanced users documentation 
  - source code supports g++-4.7 (closes: #667356)
* debian/control: dependencies for MPI distributed solvers added
* debian/rules: build commands simplified
* debian/librheolef-dev.install: man1/* to man9/* added
* debian/changelog: package description rewritted (closes: #661689)

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
field elasticity_solve(const space& Xh, Float lambda) {
 
2
  size_t d = Xh.get_geo().dimension();
 
3
  point f (0,0,0);
 
4
  f[d-1] = -1;
 
5
  field lh = riesz (Xh, f);
 
6
  form m  (Xh, Xh, "mass");
 
7
  form a1 (Xh, Xh, "div_div");
 
8
  form a2 (Xh, Xh, "2D_D");
 
9
  form a = lambda*a1 + a2;
 
10
  solver sa (a.uu());
 
11
  field uh (Xh, 0);
 
12
  uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());
 
13
  return uh;
 
14
}