← Back to branch summary

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

« 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: 2012-05-14 14:02:09 UTC
  • mfrom: (1.1.6)
  • Revision ID: package-import@ubuntu.com-20120514140209-dzbdlidkotyflf9e
Tags: 6.1-1
http://bugs.debian.org/671996
[ Pierre Saramito ]
* New upstream release 6.1 (minor changes):
  - support arbitrarily polynomial order Pk
  - source code supports g++-4.7 (closes: #671996)

[ Sylvestre Ledru ]
* update of the watch file

Show diffs side-by-side

added added

removed removed

Lines of Context:
doc/pexamples/neumann-laplace.cc
12
12
  form m (Xh, Xh, "mass");
13
13
  form a (Xh, Xh, "grad_grad");
14
14
  field b = m*field(Xh,1);
15
 
  field lh = riesz(Xh, f) + riesz(Xh, "boundary", g);
 
15
  field lh = riesz(Xh, f) + riesz(Xh, g, "boundary");
16
16
  csr<Float> A = {{  a.uu(),     b.u()}, 
17
 
                  {trans(b.u()),  0 }};
18
 
  vec<Float> B =  {  lh.u(),      0 };
19
 
  solver sa (A);
 
17
                  {trans(b.u()),  0   }};
 
18
  vec<Float> B =  {  lh.u(),      0   };
 
19
  A.set_symmetry(true);
 
20
  solver sa = ldlt(A);
20
21
  vec<Float> U = sa.solve (B);
21
22
  field uh(Xh);
22
23
  uh.set_u() = U [range(0,uh.u().size())];