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 00028 #ifndef _Deriv_Stencil_H_ 00029 #define _Deriv_Stencil_H_ 00030 00031 #include <iostream> 00032 #include <math.h> 00033 #include <assert.h> 00034 #include <stdlib.h> 00035 #include "REAL.H" 00036 #include "IntVect.H" 00037 #include "Vector.H" 00038 00040 00051 class DerivStencil 00052 { 00053 public: 00055 void define(); 00056 00058 00060 bool isDefined() const 00061 {return isdefined;} 00062 00064 00067 DerivStencil(); 00068 00070 00073 DerivStencil(const DerivStencil& a_dsin); 00074 00076 ~DerivStencil(); 00077 00079 00082 void clear(); 00083 00085 00088 int size() const; 00089 00091 00094 const IntVect& getIndex(int a_ivec) const; 00095 00097 00100 const Real& getWeight(int a_ivec) const; 00101 00103 00108 void accumulate(const IntVect& a_iv, Real a_weight); 00109 00111 00114 const DerivStencil& operator=(const DerivStencil& a_dsin); 00115 00117 00121 const DerivStencil& operator*=(Real a_facin); 00122 00124 00131 const DerivStencil& operator/=(Real a_denom); 00132 00134 00138 const DerivStencil& operator+=(Real a_facin); 00139 00141 00145 const DerivStencil& operator-=(Real a_facin); 00146 00147 protected: 00149 Vector<IntVect> m_vectIV; 00151 Vector<Real> m_vectWgt; 00153 bool isdefined; 00154 }; 00155 00156 00157 00158 #endif