00001
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
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