← Back to branch summary

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

« back to all changes in this revision

Viewing changes to nfem/ptst/field_multi_reinterpolate_tst.cc

  • 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
http://bugs.debian.org/667356
http://bugs.debian.org/661689
* 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:
nfem/ptst/field_multi_reinterpolate_tst.cc
 
1
///
 
2
/// This file is part of Rheolef.
 
3
///
 
4
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
 
5
///
 
6
/// Rheolef is free software; you can redistribute it and/or modify
 
7
/// it under the terms of the GNU General Public License as published by
 
8
/// the Free Software Foundation; either version 2 of the License, or
 
9
/// (at your option) any later version.
 
10
///
 
11
/// Rheolef is distributed in the hope that it will be useful,
 
12
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
/// GNU General Public License for more details.
 
15
///
 
16
/// You should have received a copy of the GNU General Public License
 
17
/// along with Rheolef; if not, write to the Free Software
 
18
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
19
///
 
20
/// =========================================================================
 
21
//  given u1h on a mesh omega1, re-interpolate it on a mesh omega2 as u2h
 
22
//
 
23
// usage: prog omega1 omega2
 
24
//
 
25
#include "rheolef.h"
 
26
using namespace rheolef;
 
27
using namespace std;
 
28
point u(const point& x) { return point(-x[0]+x[1]+x[2],
 
29
                                        x[0]-x[1]+x[2],
 
30
                                        x[0]+x[1]-x[2]); }
 
31
int main(int argc, char**argv) {
 
32
  rheolef::environment rheolef (argc,argv);
 
33
  Float tol = 1e-7;
 
34
  geo omega1 (argv[1]);
 
35
  geo omega2 (argv[2]);
 
36
  space V1h (omega1, "P1", "vector");
 
37
  space V2h (omega2, "P1", "vector");
 
38
  field u1h     = interpolate (V1h, u);
 
39
  field u2h     = interpolate (V2h, u1h);
 
40
  field u2h_bis = interpolate (V2h, u);
 
41
  field eh = u2h - u2h_bis;
 
42
  Float err = sqrt(dot(eh,eh));
 
43
  derr << "err = " << err << endl;
 
44
  if (argc > 3) dout << eh;
 
45
  return (err < tol) ? 0 : 1;
 
46
}