00001 /* _______ __ 00002 / ___/ / ___ __ _ / / ___ 00003 / /__/ _ \/ _ \/ ' \/ _ \/ _ \ 00004 \___/_//_/\___/_/_/_/_.__/\___/ 00005 */ 00006 // 00007 // This software is copyright (C) by the Lawrence Berkeley 00008 // National Laboratory. Permission is granted to reproduce 00009 // this software for non-commercial purposes provided that 00010 // this notice is left intact. 00011 // 00012 // It is acknowledged that the U.S. Government has rights to 00013 // this software under Contract DE-AC03-765F00098 between 00014 // the U.S. Department of Energy and the University of 00015 // California. 00016 // 00017 // This software is provided as a professional and academic 00018 // contribution for joint exchange. Thus it is experimental, 00019 // is provided ``as is'', with no warranties of any kind 00020 // whatsoever, no support, no promise of updates, or printed 00021 // documentation. By using this software, you acknowledge 00022 // that the Lawrence Berkeley National Laboratory and 00023 // Regents of the University of California shall have no 00024 // liability with respect to the infringement of other 00025 // copyrights by any part of this software. 00026 // 00027 // DFMartin, Sun, May 5, 2002 00028 00029 #ifndef _GENBICGSTABSMOOTHER_H_ 00030 #define _GENBICGSTABSMOOTHER_H_ 00031 00032 #include <iostream> 00033 #include <math.h> 00034 #include <assert.h> 00035 #include <stdlib.h> 00036 00037 #include "REAL.H" 00038 00039 #include "GenSolverOp.H" 00040 #include "GenBaseBottomSmoother.H" 00041 00043 00046 template <class T> class GenBiCGStabSmoother : public GenBaseBottomSmoother<T> 00047 { 00048 public: 00050 GenBiCGStabSmoother(); 00051 00053 virtual ~GenBiCGStabSmoother(); 00054 00056 00058 virtual GenBiCGStabSmoother<T>* newBottomSmoother() const; 00059 00061 00063 void setMaxIter(int a_maxIter); 00064 00066 00068 void setSolverTol(Real a_solverTol); 00069 00071 00075 void setEnableRestart(bool a_enableRestart); 00076 00078 00081 void setVerbose(bool a_verbose); 00082 00084 00088 virtual void doBottomSmooth(T& a_phi, 00089 const T& a_rhs, 00090 GenSolverOp<T>* a_operatorPtr); 00091 00092 protected: 00093 int m_maxIter; 00094 00095 Real m_solverTol; 00096 00097 // small parameter to check for division by zero 00098 Real m_small; 00099 00100 // small parameter to check for size of correction (to prevent hanging) 00101 Real m_convergeSmall; 00102 00103 // if true, enable restarting capability 00104 bool m_enableRestart; 00105 00106 // default is false 00107 bool m_verbose; 00108 }; 00109 00110 #include "GenBiCGStabSmootherImplem.H" 00111 00112 #endif