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

#ifndef _LINE3_BASIS_BASSO_H_
#define _LINE3_BASIS_BASSO_H_

#include "basso.h"
#include "Basis.h"

namespace basso {

class Line3Basis : public Basis {

public:
        
        // constructors
        Line3Basis() { }
        
        // destructors
        virtual ~Line3Basis() { }
        
        // member functions
        virtual BasisType       Type() const { return kLINE3; }
        virtual BasisShape      Shape() const { return kLINE; }
        virtual int                     NumNodes() const { return 3; }     
        virtual int                     NumEdges() const { return 0; }      
        virtual int                     NumFaces() const { return 0; }      
        virtual int                     Dimension() const { return 1; }
        virtual int                     Order() const { return 2; }
        
        virtual void            ParentCoord( Array< Point > &xi )       const;
        virtual Point           Centroid() const { return Point( 0.0 ); }
        virtual void            NumNodesOnFaces( iArray &nn_face )      const;
        virtual void            NumNodesOnEdges( iArray &nn_edge )      const;
        virtual void            NodesOnFace ( int f, iArray &face_nodeids )             const;
        virtual void            NodesOnEdge ( int e, iArray &edge_nodeids )                     const;
        virtual void            FaceBasisType( Array<BasisType> &face_basis )   const;
        virtual void            EdgeBasisType( Array< BasisType > &edge_basis ) const;
        
        virtual void            Quadrature( int order, QuadratureRule &qrule ) const;
        
        virtual void            Na( const Point &p, nArray &shapefunct ) const { shape_line3(shapefunct,p.x()); }
        virtual void            DNa( const Point &p, nMatrix &grad_shape )      const { dshape_line3(grad_shape,p.x()); }
        virtual void            Na( nArray &shapefunct )    const { Na( Centroid(), shapefunct ); }
        virtual void            DNa( nMatrix &grad_shape ) const { DNa( Centroid(), grad_shape ); }

private:
        
};

void Line3Basis::ParentCoord( Array< Point > &xi ) const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( xi.size()!=NumNodes() )
                xi.resize(NumNodes());
#endif
        xi[0]=Point(-1.0);
        xi[1]=Point(1.0);
}

void    Line3Basis::NumNodesOnFaces( iArray &nn_face )  const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( nn_face.size()!=NumEdges() )
                nn_face.resize( NumEdges() );
#endif
}

void    Line3Basis::NumNodesOnEdges( iArray &nn_edge )  const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( nn_edge.size()!=NumEdges() )
                nn_edge.resize( NumEdges() );
#endif 
}

void Line3Basis::NodesOnFace ( int e, iArray &face_nodeids ) const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( face_nodeids.size()!=0 )
                face_nodeids.resize(0);
#endif
}

void Line3Basis::NodesOnEdge ( int e, iArray &edge_nodeids ) const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( edge_nodeids.size()!=0 )
                edge_nodeids.resize(0);
#endif
}

void Line3Basis::FaceBasisType( Array< BasisType > &face_basis ) const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( face_basis.size()!=NumFaces() )
                face_basis.resize(NumFaces());
#endif
}

void Line3Basis::EdgeBasisType( Array< BasisType > &edge_basis ) const
{
#ifdef ALLOW_DYNAMIC_RESIZE 
        if ( edge_basis.size()!=NumEdges() )
                edge_basis.resize(NumEdges());
#endif
}

void Line3Basis::Quadrature( int order, QuadratureRule &qrule ) const
{
        
        int npts;
        npts=(order+1)/2;

        if ( vect_size(qrule)!=npts )
                resize(qrule,npts);

        quadrature_gauss1d( qrule );
        
}

} // end namespace

#endif

---