fsi/sfc-amroc/TubeCJBurnFlaps/src/ShellEnforceBCFunctor.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// 
//                           Fehmi Cirak, Ralf Deiterding
//                        California Institute of Technology
//                           (C) 2005 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef SHELLENFORCEBCFUNCTOR_H
#define SHELLENFORCEBCFUNCTOR_H

#include "shells/driverCC/SVertexFunctors.h"
#include "shells/driverCC/PrescribeVarFunctor.h"

#include <functional>


namespace shells {
    class ShellEnforceBCFunctor;
}


class shells::ShellEnforceBCFunctor :
    public std::unary_function<shells::SVertexS*, void> {
    
public:
    ShellEnforceBCFunctor(double minz, double tubeLength, double torsionAngle):
        _minz(minz+1.e-8), _maxz(tubeLength-1.e-8), _torsionAngle(torsionAngle){}
    
    ~ShellEnforceBCFunctor() {}

    void operator()(shells::SVertex *vtx) const {
        SVertexCoordinate extractC;
        SVertexCoordinate::DataType *coor = extractC(vtx);      
        const double xcoor = coor[0];
        const double ycoor = coor[1];
        const double zcoor = coor[2];
            
        // set the velocities and accelerations to zero at the boundaries
        if ((zcoor<_minz)||(zcoor>_maxz)) {         
            SVertexVelocity extractV;
            SVertexAcceleration extractA;

            SVertexVelocity::DataType *vel = extractV(vtx);
            SVertexAcceleration::DataType *acc = extractA(vtx);
            for (unsigned i=0; i<SVertexVelocity::numVar; ++i) {
                vel[i] = 0.0;
                acc[i] = 0.0;
            }
        }       
        return;
    }
    
private:    
    const double     _minz;
    const double     _maxz;    
    const double     _torsionAngle;
};
#endif

---