00001 00005 /* Copyright (c) 2005-2011 Taneli Kalvas. All rights reserved. 00006 * 00007 * You can redistribute this software and/or modify it under the terms 00008 * of the GNU General Public License as published by the Free Software 00009 * Foundation; either version 2 of the License, or (at your option) 00010 * any later version. 00011 * 00012 * This library is distributed in the hope that it will be useful, but 00013 * WITHOUT ANY WARRANTY; without even the implied warranty of 00014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00015 * General Public License for more details. 00016 * 00017 * You should have received a copy of the GNU General Public License 00018 * along with this library (file "COPYING" included in the package); 00019 * if not, write to the Free Software Foundation, Inc., 51 Franklin 00020 * Street, Fifth Floor, Boston, MA 02110-1301 USA 00021 * 00022 * If you have questions about your rights to use or distribute this 00023 * software, please contact Berkeley Lab's Technology Transfer 00024 * Department at TTD@lbl.gov. Other questions, comments and bug 00025 * reports should be sent directly to the author via email at 00026 * taneli.kalvas@jyu.fi. 00027 * 00028 * NOTICE. This software was developed under partial funding from the 00029 * U.S. Department of Energy. As such, the U.S. Government has been 00030 * granted for itself and others acting on its behalf a paid-up, 00031 * nonexclusive, irrevocable, worldwide license in the Software to 00032 * reproduce, prepare derivative works, and perform publicly and 00033 * display publicly. Beginning five (5) years after the date 00034 * permission to assert copyright is obtained from the U.S. Department 00035 * of Energy, and subject to any subsequent five (5) year renewals, 00036 * the U.S. Government is granted for itself and others acting on its 00037 * behalf a paid-up, nonexclusive, irrevocable, worldwide license in 00038 * the Software to reproduce, prepare derivative works, distribute 00039 * copies to the public, perform publicly and display publicly, and to 00040 * permit others to do so. 00041 */ 00042 00043 #ifndef BICGSTAB_SOLVER_HPP 00044 #define BICGSTAB_SOLVER_HPP 1 00045 00046 00047 #include <iostream> 00048 #include "solver.hpp" 00049 #include "geometry.hpp" 00050 #include "problem.hpp" 00051 #include "matrix.hpp" 00052 #include "mvector.hpp" 00053 00054 00062 class BiCGSTABSolver : public Solver { 00063 double _eps; 00064 uint32_t _imax; 00066 bool _gnewton; 00067 double _newton_Reps; 00068 double _newton_dXeps; 00069 uint32_t _newton_imax; 00071 public: 00072 00075 BiCGSTABSolver( double eps = 1.0e-6, uint32_t imax = 10000, 00076 double newton_Reps = 1.0e-5, double newton_dXeps = 1.0e-6, 00077 uint32_t newton_imax = 10, bool gnewton = false ); 00078 00081 ~BiCGSTABSolver() {} 00082 00086 virtual void solve( const Problem &p, Vector &X ); 00087 00094 virtual void reset( void ); 00095 00098 void set_gnewton( bool enable ) { 00099 _gnewton = enable; 00100 } 00101 00104 void set_eps( double eps ) { 00105 _eps = eps; 00106 } 00107 00110 void set_imax( uint32_t imax ) { 00111 _imax = imax; 00112 } 00113 00116 void set_newton_imax( uint32_t newton_imax ) { 00117 _newton_imax = newton_imax; 00118 } 00119 00122 void set_newton_residual_eps( double newton_Reps ) { 00123 _newton_Reps = newton_Reps; 00124 } 00125 00128 void set_newton_step_eps( double newton_dXeps ) { 00129 _newton_dXeps = newton_dXeps; 00130 } 00131 00132 }; 00133 00134 00135 #endif 00136 00137 00138 00139 00140 00141 00142 00143 00144 00145 00146 00147 00148 00149 00150 00151 00152 00153 00154 00155 00156