← Back to branch summary

~ubuntu-branches/ubuntu/raring/rheolef/raring-proposed

« back to all changes in this revision

Viewing changes to nfem/tst/tensor_eig_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/tst/tensor_eig_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
 
// usage: tensor_eig_tst | octave -q
22
 
#include "rheolef.h"
23
 
using namespace rheolef;
24
 
using namespace std;
25
 
 
26
 
int eig_check (const tensor& a, size_t dim) {
27
 
    tensor q;
28
 
    point d = a.eig (q,dim);
29
 
    tensor err_a = a - q*diag(d)*trans(q);
30
 
    Float err = norm(err_a);
31
 
    cout << setprecision(15) << matlab
32
 
         << "a = " ; a.put(cout,dim); cout << endl
33
 
         << "d = " ; diag(d).put(cout,dim); cout << endl
34
 
         << "q = " ; q.put(cout,dim); cout << endl
35
 
         << "err = norm(a-q*d*q')" << endl
36
 
         << "%err =" << err << endl
37
 
        ;
38
 
    return (err < 1e-10) ? 0 : 1;
39
 
}
40
 
int eig2x2_tst () {
41
 
    tensor a;
42
 
    a(0,0) = 1;
43
 
    a(1,0) = 2;
44
 
    a(0,1) = 2;
45
 
    a(1,1) = 3;
46
 
    return eig_check (a, 2);
47
 
}
48
 
int eig3x3_tst () {
49
 
    tensor a;
50
 
    a(0,0) = 1;
51
 
    a(1,0) = 2;
52
 
    a(2,0) = 3;
53
 
    a(0,1) = 2;
54
 
    a(1,1) = 3;
55
 
    a(2,1) = 4;
56
 
    a(0,2) = 3;
57
 
    a(1,2) = 4;
58
 
    a(2,2) = 5;
59
 
    return eig_check (a, 3);
60
 
}
61
 
int main(int argc, char**argv) {
62
 
  int status = 0;
63
 
  status |= eig2x2_tst() ;
64
 
  status |= eig3x3_tst() ;
65
 
  return status;
66
 
}