simulation.h

Go to the documentation of this file.

#ifndef SIMULATION_H
#define SIMULATION_H

class BezierMesh;
class BoundaryMesh;
class DataStore;
class PersistantStore;
#include "util.h" // LinearData, ControlPoint

class Simulation {
private:
    /* Define copy and assignment constructors as private to prevent their use*/
    Simulation( const Simulation &o);
    const Simulation& operator=( const Simulation &o);


    void initialize_structures();

    const LinearData& get_data(const ControlPoint&) const;

    PersistantStore *persistant_store;
    DataStore *data_store; 
    BoundaryMesh *bdry_mesh; 
    BezierMesh *bezier_mesh; 
public:

    Simulation( char *filename );
    Simulation();

    BezierMesh* get_bezier_mesh();
    BoundaryMesh* get_boundary_mesh();
    DataStore* get_data_store();

    // Methods "move_" put the new position into the data string

    // 5 arguments, t, vx, vy, x1, y1
    void move_(double time_step, unsigned vx_pos, unsigned vy_pos, unsigned newx_pos, unsigned newy_pos);
    // 7 arguments, t, v0x, v0y, v1x, v1y, x1, y1
    void move_(double time_step, unsigned vx1_pos, unsigned vy1_pos, unsigned vx0_pos, unsigned vy0_pos,  unsigned newx_pos, unsigned newy_pos);
    
    // Shorthands, assume y index is always x+1
    // 3 arguments, t, vxy, xy1
    void move_(double time_step, unsigned vx_pos, unsigned newx_pos);

    // 4 arguments, t, v1xy, v0xy, xy1
    void move_(double time_step, unsigned vx1_pos, unsigned vx0_pos, unsigned newx_pos);


    // Methods "move" stores position in XposIndex and swaps the new position into the geometry
    unsigned XposIndex;
    unsigned YposIndex;
    bool positions_set;

    void set_positions(unsigned xPosI,unsigned yPosI);
    
    // 3 arguments, t, vx, vy
    void move(double time_step, unsigned vx_pos, unsigned vy_pos);
    // 5 arguments, t, v0x, v0y, v1x, v1y
    void move(double time_step, unsigned vx1_pos, unsigned vy1_pos, unsigned vx0_pos, unsigned vy0_pos);
    // Shorthands, assume y index is always x+1
    // 2 arguments, t, vxy
    void move(double time_step, unsigned vx_pos);

    // 3 arguments, t, vx, vy
    double move_safe(double time_step, unsigned vx_pos, unsigned vy_pos);
    // 5 arguments, t, v0x, v0y, v1x, v1y
    double move_safe(double time_step, unsigned vx1_pos, unsigned vy1_pos, unsigned vx0_pos, unsigned vy0_pos);
    // Shorthands, assume y index is always x+1
    // 2 arguments, t, vxy
    double move_safe(double time_step, unsigned vx_pos);

    void store_state(unsigned vIndexFrom, unsigned vIndexTo, unsigned xIndexTo);
    void restore_state(unsigned vIndexFrom, unsigned vIndexTo, unsigned xIndexFrom);

    /* Read in a text file, and use it to initialize bdry_mesh and bezier_mesh
     * never read in more than one file to a simulation.  There are 4 seperate
     * files for text file I/O.  see meshio.h
     */
    int read_from_text_file( char *filename );

    /* Same as the text file version, except this function loads from a binary
     * file.  This preserves state much better between runs as floating point
     * roundoff is minimized.
     *
     * this function requires filename end in ".geo"
     */
    int read_from_binary_file( char *filename );

    /* write the current state of the meshes to a text file, using the 4
     * filenames specified in meshio.C
     */
    int to_text_file( char *filename );

    /* write the current state of the meshes to a binary file.  This is much more
     * accurate as it doesn't cause as much roundoff of floating point values.
     * filename must end in ".geo"
     */
    int to_binary_file( char *filename);

    /* write the current state to an EPS file.  This is done according to the default
     * EPSWriter values.
     * epsfile should end in ".eps"
     */
    int to_eps(char *epsfile);

    /* Print out data for both bezier and boundary meshes */
    void print();

    /* Print out size statistics for both meshes */
    void print_statistics();



};

#endif /* SIMULATION_H */

---