shells/fragmentation/ShellManagerFragmented.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// 
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2005 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef SHELLMANAGERFRAGMENTED_H
#define SHELLMANAGERFRAGMENTED_H
#include <string>
#include <vector>


namespace shells {
    class MShell;
    class ELCCoder;
    class CheckpointData; 
    class ShellManagerFragmented;
}

namespace fragment {
    class CCheckpointData;
}


class shells::ShellManagerFragmented {
    //
    // Facet class for serial fragmenting shells 
    // Provides a basic interface (with limited functionality) for shells library    
    // 
    
public:
    ShellManagerFragmented(const std::string& controlFileName);
    virtual ~ShellManagerFragmented();
    
    enum PressureLocation {
        element = 0,           // apply pressure to the elements
        vertex = 1             // apply pressure to the vertices 
    };

    // initialization 
    virtual void computeMassPrepareAdvance();

    // time integration primitives
    virtual void predictAndEnforceBC();
    virtual void internalExternalForces();
    virtual void correct();

    // time integration primitives combined
    virtual void advance();

    virtual void printData(bool flush=false, unsigned myRank=0);    

    virtual void checkPointing(unsigned fileIDFirst=0, unsigned fileIDSecond=0);
    virtual void restart(unsigned fileIDFirst=0, unsigned fileIDSecond=0);
    
    // decode shell data structure
    virtual void decode(double **coordinates,
                        double **velocities,
                        int **globalNodeIDs,
                        int *numNodes,
                        int **connectivity,
                        int *numElements);

    // encode pressures
    virtual void encodePressure(double* pressures, unsigned size, PressureLocation loc); 

    // print interface mesh and pressure for debugging purposes
    void printIFaceMeshPressure(unsigned myRank=0);

    // accessors
    virtual double stableTimeStep() {
        return _timeStepSize;
    }
        
// copy and assignment constructors - not implemented
private:
    ShellManagerFragmented(const ShellManagerFragmented &);
    const ShellManagerFragmented & operator=(const ShellManagerFragmented &);
    
// accessors for derived classes
protected:
    shells::MShell*& mShell(){return _mshell;} 
    int getCurrentStepNum(){return _currentStepNum;}
    void incrementCurrentTimeAndStep(){
        _currentStepNum++;
        _currentTime += _timeStepSize;
    }

    void setTimeStep(double timeStepSize){_timeStepSize=timeStepSize;}
    double getTimeStep(){return _timeStepSize;}
    
    double getCurrentTime(){return _currentTime;}
    double Thickness() {return _thickness;}

private:
    shells::MShell                   *_mshell;

    shells::CheckpointData           *_checkPointerS;
    fragment::CCheckpointData        *_checkPointerC;

    // coder/decoder used for eulerian lagrangian coupling
    shells::ELCCoder                 *_elc;

    // vector of vertex/element pressures - temporary vector
    std::vector<double>     _pressures;
    PressureLocation        _pType;                  

    // control variables
    unsigned int    _numSurfIntgPoints;
    unsigned int    _numThickIntgPoints;
    unsigned int    _printAtEach;

    unsigned int    _currentStepNum;
    double          _currentTime;

    double          _newmarkGamma; 
    double          _newmarkDamping;
    double          _timeStepSize;
    double          _thickness;
    
    std::string     _outputFileName;
    std::string     _checkPointFileName;
};

#endif

---