meshvectorfield.hpp
Go to the documentation of this file.
00001 00005 /* Copyright (c) 2005-2011 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 MESHVECTORFIELD_HPP 00044 #define MESHVECTORFIELD_HPP 1 00045 00046 00047 #include "transformation.hpp" 00048 #include "vectorfield.hpp" 00049 #include "mesh.hpp" 00050 #include "types.hpp" 00051 00052 00077 class MeshVectorField : public VectorField, public Mesh { 00078 00079 field_extrpl_e _extrpl[6]; 00080 double *_F[3]; 00086 Transformation _T; 00087 Transformation _Tinv; 00089 void check_definition(); 00090 00091 void convert_3d_to_3d( const MeshVectorField &fin ); 00092 void convert_cyl_to_3d( const MeshVectorField &fin ); 00093 00094 public: 00095 00096 00104 MeshVectorField(); 00105 00112 MeshVectorField( const Mesh &m, const bool fout[3] ); 00113 00124 MeshVectorField( geom_mode_e geom_mode, const bool fout[3], Int3D size, 00125 Vec3D origo, double h ); 00126 00152 MeshVectorField( geom_mode_e geom_mode, const bool fout[3], double xscale, 00153 double fscale, const std::string &filename ); 00154 00178 MeshVectorField( geom_mode_e geom_mode, const bool fout[3], Int3D size, 00179 Vec3D origo, double h, const MeshVectorField &fin ); 00180 00183 MeshVectorField( const MeshVectorField &f ); 00184 00187 MeshVectorField( std::istream &is ); 00188 00191 virtual ~MeshVectorField(); 00192 00212 void set_extrapolation( const field_extrpl_e extrpl[6] ) { 00213 memcpy( _extrpl, extrpl, 6*sizeof(field_extrpl_e) ); 00214 } 00215 00218 void reset_transformation( void ); 00219 00222 void set_transformation( const Transformation &T ); 00223 00226 void translate( const Vec3D &dx ); 00227 00230 void scale( const Vec3D &sx ); 00231 00236 void rotate_x( double a ); 00237 00244 void rotate_y( double a ); 00245 00250 void rotate_z( double a ); 00251 00254 void clear(); 00255 00262 void reset( geom_mode_e geom_mode, const bool fout[3], Int3D size, 00263 Vec3D origo, double h ); 00264 00268 void get_minmax( double &min, double &max ) const; 00269 00275 void get_minmax( Vec3D &min, Vec3D &max ) const; 00276 00279 void get_defined_components( bool fout[3] ) const; 00280 00283 MeshVectorField &operator=( const MeshVectorField &f ); 00284 00287 MeshVectorField &operator+=( const MeshVectorField &f ); 00288 00291 MeshVectorField &operator*=( double x ); 00292 00295 MeshVectorField &operator/=( double x ); 00296 00299 const Vec3D operator()( int32_t i ) const; 00300 00303 const Vec3D operator()( int32_t i, int32_t j ) const; 00304 00307 const Vec3D operator()( int32_t i, int32_t j, int32_t k ) const; 00308 00311 void set( int32_t i, const Vec3D &v ); 00312 00315 void set( int32_t i, int32_t j, const Vec3D &v ); 00316 00319 void set( int32_t i, int32_t j, int32_t k, const Vec3D &v ); 00320 00324 virtual const Vec3D operator()( Vec3D x ) const; 00325 00328 void save( const std::string &filename ) const; 00329 00332 void save( std::ostream &os ) const; 00333 00336 void debug_print( std::ostream &os ) const; 00337 }; 00338 00339 00340 #endif 00341