bicgstab_solver.hpp
Go to the documentation of this file.
1 
5 /* Copyright (c) 2005-2011 Taneli Kalvas. All rights reserved.
6  *
7  * You can redistribute this software and/or modify it under the terms
8  * of the GNU General Public License as published by the Free Software
9  * Foundation; either version 2 of the License, or (at your option)
10  * any later version.
11  *
12  * This library is distributed in the hope that it will be useful, but
13  * WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  * General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License
18  * along with this library (file "COPYING" included in the package);
19  * if not, write to the Free Software Foundation, Inc., 51 Franklin
20  * Street, Fifth Floor, Boston, MA 02110-1301 USA
21  *
22  * If you have questions about your rights to use or distribute this
23  * software, please contact Berkeley Lab's Technology Transfer
24  * Department at TTD@lbl.gov. Other questions, comments and bug
25  * reports should be sent directly to the author via email at
26  * taneli.kalvas@jyu.fi.
27  *
28  * NOTICE. This software was developed under partial funding from the
29  * U.S. Department of Energy. As such, the U.S. Government has been
30  * granted for itself and others acting on its behalf a paid-up,
31  * nonexclusive, irrevocable, worldwide license in the Software to
32  * reproduce, prepare derivative works, and perform publicly and
33  * display publicly. Beginning five (5) years after the date
34  * permission to assert copyright is obtained from the U.S. Department
35  * of Energy, and subject to any subsequent five (5) year renewals,
36  * the U.S. Government is granted for itself and others acting on its
37  * behalf a paid-up, nonexclusive, irrevocable, worldwide license in
38  * the Software to reproduce, prepare derivative works, distribute
39  * copies to the public, perform publicly and display publicly, and to
40  * permit others to do so.
41  */
42 
43 #ifndef BICGSTAB_SOLVER_HPP
44 #define BICGSTAB_SOLVER_HPP 1
45 
46 
47 #include <iostream>
48 #include "solver.hpp"
49 #include "geometry.hpp"
50 #include "problem.hpp"
51 #include "matrix.hpp"
52 #include "mvector.hpp"
53 
54 
62 class BiCGSTABSolver : public Solver {
63  double _eps;
64  uint32_t _imax;
66  bool _gnewton;
67  double _newton_Reps;
68  double _newton_dXeps;
69  uint32_t _newton_imax;
71 public:
72 
75  BiCGSTABSolver( double eps = 1.0e-6, uint32_t imax = 10000,
76  double newton_Reps = 1.0e-5, double newton_dXeps = 1.0e-6,
77  uint32_t newton_imax = 10, bool gnewton = false );
78 
82 
86  virtual void solve( const Problem &p, Vector &X );
87 
94  virtual void reset( void );
95 
98  void set_gnewton( bool enable ) {
99  _gnewton = enable;
100  }
101 
104  void set_eps( double eps ) {
105  _eps = eps;
106  }
107 
110  void set_imax( uint32_t imax ) {
111  _imax = imax;
112  }
113 
116  void set_newton_imax( uint32_t newton_imax ) {
117  _newton_imax = newton_imax;
118  }
119 
122  void set_newton_residual_eps( double newton_Reps ) {
123  _newton_Reps = newton_Reps;
124  }
125 
128  void set_newton_step_eps( double newton_dXeps ) {
129  _newton_dXeps = newton_dXeps;
130  }
131 
132 };
133 
134 
135 #endif
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156