Navigation

Main Page
Download
Support
Installation
Tutorial
Examples
Reference Manual
   Version 1.0.5new_solver
   Version 1.0.5dev
      Class Index
      File List
   Version 1.0.5b
   Version 1.0.4dev
   Version 1.0.4
Publications


Hosted by Get Ion Beam Simulator at SourceForge.net. Fast, secure and Free Open Source software downloads

gs_solver.hpp

Go to the documentation of this file.
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 GS_SOLVER_HPP
00044 #define GS_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 "crowmatrix.hpp"
00053 #include "mvector.hpp"
00054 
00055 
00068 class GSSolver : public Solver {
00069     double   _eps;
00070     uint32_t _imax;
00071     double   _w;
00072 
00073     double   _newton_Reps;  
00074     double   _newton_dXeps; 
00075     double   _newton_imax;  
00077     static void gauss_seidel_error( const std::string func, int status );
00078 
00079 public:
00080 
00083     GSSolver( double w = 1.66,
00084               double eps = 1.0e-6, 
00085               uint32_t imax = 10000, 
00086               double newton_Reps = 1.0e-5, 
00087               double newton_dXeps = 1.0e-6, 
00088               uint32_t newton_imax = 10 );
00089 
00092     ~GSSolver() {}
00093 
00097     virtual void solve( const Problem &p, Vector &X );
00098 
00105     virtual void reset( void );
00106 
00109     void set_newton_imax( uint32_t newton_imax ) {
00110         _newton_imax = newton_imax;
00111     }
00112 
00115     void set_newton_residual_eps( double newton_Reps ) {
00116         _newton_Reps = newton_Reps;
00117     }
00118 
00121     void set_newton_step_eps( double newton_dXeps ) {
00122         _newton_dXeps = newton_dXeps;
00123     }
00124 
00131     static bool gauss_seidel( const CRowMatrix &mat, const Vector &rhs, Vector &sol,
00132                               uint32_t &imax, double &eps, double w );
00133 };
00134 
00135 
00136 #endif
00137 
00138 
00139 
00140 
00141 
00142 


Reference manual for Ion Beam Simulator 1.0.5dev
Generated by Doxygen 1.7.1 on Mon Feb 6 2012 15:07:15.