/Users/jack/Code/basso_dev/inc/PointLoadBC.h

Go to the documentation of this file.

#ifndef _POINT_LOAD_BC_BASSO_H_
#define _POINT_LOAD_BC_BASSO_H_

// std includes
#include <list>
#include <set>

// Basso includes
#include "basso.h"
#include "NeumannBC.h"

namespace basso {


class PointLoadBC : public NeumannBC {

public:
        
        // PUBLIC DEFS
        
        // CONSTRUCTORS 
        PointLoadBC( const set<int> &nids, const Array<Dof> &dofs, const nArray &n, Numeric sf=1.0  ) 
                : NeumannBC( dofs, sf )  
        { 
                if ( n.size()!=dofs.size() )
                        WarningMessage("PointLoadBC::PointLoadBC","incorrect dimsnsion for load direction vector");
                fNodeIds = &nids; 
                mDirection=n; 
        }
        
        // CONSTRUCTORS 
        PointLoadBC( const set<int> &nids, Dof &ldof, Numeric sf  ) 
                : NeumannBC( ldof, sf ), mDirection(1)  
        { 
                fNodeIds = &nids; 
                mDirection[0]=1.0; 
        }
        
        virtual ~PointLoadBC() {  }
        
        // ACCESSORS
        
        // MEMBER FUNCTIONS
        virtual void ContributeForce( const DofMap &dofmap, nArray &f );
        
protected:

protected:
        const set<int> *fNodeIds;
        
        nArray mDirection;
        
};

void PointLoadBC::ContributeForce( const DofMap &dofmap, nArray &f ) 
{
        set<int>::const_iterator nItr;
        Array<Dof>::const_iterator dofItr;
        for ( nItr=fNodeIds->begin(); nItr!=fNodeIds->end(); ++nItr ) 
        {
                int i=0;
                for ( dofItr=this->mDofs.begin(); dofItr!=this->mDofs.end(); ++dofItr, ++i )
                {
                        int gdof = dofmap.GlobalDof( *nItr, *dofItr );
                        f[gdof] += (this->mScaleFactor)*mDirection[i];
                }
        }
}


} // end namespace

#endif

---