IBSimu
1.0.4
|
00001 00005 /* Copyright (c) 2005-2010 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 SCALARFIELD_HPP 00044 #define SCALARFIELD_HPP 1 00045 00046 00047 #include <iostream> 00048 #include "geometry.hpp" 00049 #include "vec3d.hpp" 00050 00051 00066 class ScalarField { 00067 geom_mode_e _geom_mode; 00068 Int3D _size; 00069 Vec3D _origo; 00070 Vec3D _max; 00071 double _h; 00072 double _div_h; 00073 double *_F; 00075 public: 00076 00077 00080 ScalarField() : _geom_mode(MODE_3D), _h(1.0), _div_h(1.0), _F(0) {} 00081 00088 ScalarField( const Geometry &g ); 00089 00095 ScalarField( geom_mode_e geom_mode, Int3D size, Vec3D origo, double h ); 00096 00099 ScalarField( std::istream &s ); 00100 00103 ScalarField( const ScalarField &f ); 00104 00107 ~ScalarField(); 00108 00111 inline geom_mode_e geom_mode( void ) const { return( _geom_mode ); } 00112 00115 inline Int3D size( void ) const { return( _size ); } 00116 00119 inline int32_t size( int i ) const { return( _size[i] ); } 00120 00123 inline int32_t nodecount( void ) const { return( _size[0]*_size[1]*_size[2] ); } 00124 00127 inline Vec3D origo( void ) const { return( _origo ); } 00128 00131 inline double origo( int i ) const { return( _origo[i] ); } 00132 00136 Vec3D max( void ) const { return( _max ); } 00137 00141 double max( int i ) const { return( _max[i] ); } 00142 00145 inline double h( void ) const { return( _h ); } 00146 00149 inline double div_h( void ) const { return( _div_h ); } 00150 00153 void clear(); 00154 00160 void reset( geom_mode_e geom_mode, Int3D size, Vec3D origo, double h ); 00161 00170 void get_minmax( double &min, double &max ) const; 00171 00179 void epot_get_minmax( const Geometry &g, double &min, double &max ) const; 00180 00183 ScalarField &operator=( const ScalarField &f ); 00184 00187 ScalarField &operator+=( const ScalarField &f ); 00188 00191 ScalarField &operator*=( double x ); 00192 00195 ScalarField &operator/=( double x ); 00196 00199 inline double &operator()( int32_t i ) { 00200 return( _F[i] ); } 00201 00204 inline double &operator()( int32_t i, int32_t j ) { 00205 return( _F[i + j*_size[0]] ); } 00206 00209 inline double &operator()( int32_t i, int32_t j, int32_t k ) { 00210 return( _F[i + (j + k*_size[1])*_size(0)] ); } 00211 00214 inline const double &operator()( int32_t i ) const { 00215 return( _F[i] ); } 00216 00219 inline const double &operator()( int32_t i, int32_t j ) const { 00220 return( _F[i + j*_size[0]] ); } 00221 00224 inline const double &operator()( int32_t i, int32_t j, int32_t k ) const { 00225 return( _F[i + (j + k*_size[1])*_size(0)] ); } 00226 00234 double operator()( Vec3D x ) const; 00235 00238 void save( std::ostream &s ) const; 00239 00242 void debug_print( void ) const; 00243 }; 00244 00245 00246 #endif 00247 00248 00249 00250 00251 00252 00253 00254 00255 00256 00257 00258 00259 00260 00261 00262 00263 00264