← Back to branch summary

~ubuntu-branches/ubuntu/trusty/rheolef/trusty-proposed

« back to all changes in this revision

Viewing changes to skit/tst/solve_gmres_tst.cc

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2010-06-12 09:08:59 UTC
  • Revision ID: james.westby@ubuntu.com-20100612090859-8gpm2gc7j3ab43et
Tags: upstream-5.89
Import upstream version 5.89

Show diffs side-by-side

added added

removed removed

Lines of Context:
skit/tst/solve_gmres_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
//
 
22
// test gmres solver
 
23
//
 
24
# include "rheolef/skit.h"
 
25
# include "rheolef/gmres.h"
 
26
using namespace std;
 
27
 
 
28
extern
 
29
void
 
30
check (
 
31
    const char* method, 
 
32
    int max_iter, 
 
33
    const csr<Float>& a, 
 
34
    const vec<Float>& x, 
 
35
    const vec<Float>& b);
 
36
 
 
37
int solve_gmres (const csr<Float>& a, const vec<Float>& b)
 
38
{
 
39
    vec<Float> x(a.ncol());
 
40
    x = 0; 
 
41
    int max_iter = 100; 
 
42
    Float tol = 1e-12;
 
43
    int m = 6;
 
44
    dns<Float> H(m+1,m+1);
 
45
    int status = gmres(a, x, b, EYE, H, m, max_iter, tol);
 
46
    check("gmres",max_iter,a,x,b);
 
47
    return status;
 
48
}
 
49