← Back to branch summary

~ubuntu-branches/ubuntu/wily/rheolef/wily

« back to all changes in this revision

Viewing changes to doc/pexamples/neumann-laplace.cc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito, Pierre Saramito, Sylvestre Ledru
  • Date: 2013-04-30 11:29:51 UTC
  • mfrom: (1.2.1) (11.1.1 experimental)
  • Revision ID: package-import@ubuntu.com-20130430112951-g69nuqin77t75dek
Tags: 6.4-1
[ Pierre Saramito ]
* New upstream release 6.4 (major changes):
  - new powerful c++ expressions for FEM specifications
  - configure script improved (portability increased)
  - minor bug fixes
* control:
  - all architectures are considered (armel & armhf are back)
  - no more circular-dependency librheolef-dev <--> rheolef
* rules: add security compiler flags (hardening and -W,-l,relro)
* watch: include an updated version, thanks to B. Martens <bartm@debian.org>
* copyright: various file format fixes

[ Sylvestre Ledru ]
*  Standards-Version updated to 3.9.4

Show diffs side-by-side

added added

removed removed

Lines of Context:
doc/pexamples/neumann-laplace.cc
9
9
  geo omega (argv[1]);
10
10
  d = omega.dimension();
11
11
  space Xh (omega, argv[2]);
12
 
  form m (Xh, Xh, "mass");
13
 
  form a (Xh, Xh, "grad_grad");
14
 
  field b = m*field(Xh,1);
15
 
  field lh = riesz(Xh, f) + riesz(Xh, g, "boundary");
 
12
  trial u (Xh); test v (Xh);
 
13
  form  m  = integrate (u*v);
 
14
  form  a  = integrate (dot(grad(u),grad(v)));
 
15
  field b  = m*field(Xh,1);
 
16
  field lh = integrate (f*v) + integrate ("boundary", g*v);
16
17
  csr<Float> A = {{  a.uu(),     b.u()}, 
17
18
                  {trans(b.u()),  0   }};
18
19
  vec<Float> B =  {  lh.u(),      0   };