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 _BiCGStabSmoother_H_ 00030 #define _BiCGStabSmoother_H_ 00031 00032 #include <iostream> 00033 #include <math.h> 00034 #include <assert.h> 00035 #include <stdlib.h> 00036 #include <REAL.H> 00037 #include <IntVect.H> 00038 #include <Box.H> 00039 #include <DisjointBoxLayout.H> 00040 #include <LevelData.H> 00041 #include <FArrayBox.H> 00042 #include <ProblemDomain.H> 00043 #include "BaseBottomSmoother.H" 00044 00046 00049 class BiCGStabSmoother : public BaseBottomSmoother 00050 { 00051 00052 public: 00053 00055 00056 00058 BiCGStabSmoother(); 00059 00061 virtual ~BiCGStabSmoother(); 00062 00064 00066 virtual BaseBottomSmoother* new_bottomSmoother() const; 00067 00069 00071 void setMaxIter(int a_max_iter); 00072 00074 00076 void setSolverTol(Real a_solver_tol); 00077 00079 00082 void setVerbose(bool a_verbose); 00083 00084 00086 00090 virtual void doBottomSmooth(LevelData<FArrayBox>& a_phi, 00091 const LevelData<FArrayBox>& a_rhs, 00092 LevelOp* a_levelop_ptr); 00093 00094 protected: 00095 00096 int m_maxIter; 00097 00098 Real m_solver_tol; 00099 00100 // small parameter to check for division by zero 00101 Real m_small; 00102 00103 // default is false 00104 bool m_verbose; 00105 00106 }; 00107 00108 #endif 00109 00110 00111 00112