shells/driverCC/ShellManagerBasic.h

Go to the documentation of this file.

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


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


class shells::ShellManagerBasic {
    //
    // Facet class for serial non-fragmenting shells 
    // Provides a basic interface (with limited functionality) for shells library    
    // 
    
public:
    ShellManagerBasic(const std::string& controlFileName);
    virtual ~ShellManagerBasic();
    
    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 predict();
    virtual void predictAndEnforceBC();
    virtual void internalExternalForces();
    virtual void correct();

    // time integration primitives combined
    virtual void advance();
    
    // print shell mesh and data
    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:
    ShellManagerBasic(const ShellManagerBasic &);
    const ShellManagerBasic & operator=(const ShellManagerBasic &);
    
// 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           *_checkPointer;

    // 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    _bcType[3];

    unsigned int    _currentStepNum;
    double          _currentTime;

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

#endif

---