#include <shape_impl.h>

Inheritance diagram for PiecewiseLinearShape1D:

Collaboration diagram for PiecewiseLinearShape1D:

Public Member Functions
	PiecewiseLinearShape1D (int n, doublereal x, doublereal v)

	~PiecewiseLinearShape1D (void)

doublereal	dGet (doublereal d, doublereal=0.) const

std::ostream &	Restart (std::ostream &out) const

Public Member Functions inherited from Shape1D
virtual	~Shape1D (void)

Public Member Functions inherited from Shape
virtual	~Shape (void)

Protected Attributes
int	nPoints

doublereal *	pdX

doublereal *	pdV

Detailed Description

Definition at line 76 of file shape_impl.h.

Constructor & Destructor Documentation

PiecewiseLinearShape1D::PiecewiseLinearShape1D	(	int	n,
		doublereal *	x,
		doublereal *	v
	)

Definition at line 171 of file shape.cc.

References ASSERT, nPoints, pdV, and pdX.

 : nPoints(n), pdX(x), pdV(v)
 {
         ASSERT(nPoints > 0);
         ASSERT(pdX != NULL);
         ASSERT(pdV != NULL);
 }

PiecewiseLinearShape1D::~PiecewiseLinearShape1D ( void )

Definition at line 180 of file shape.cc.

References pdV, pdX, and SAFEDELETEARR.

 {
         SAFEDELETEARR(pdX);
         SAFEDELETEARR(pdV);
 }

Member Function Documentation

doublereal PiecewiseLinearShape1D::dGet	(	doublereal	d,
		doublereal	= `0.`
	)		const

virtual

Implements Shape.

Definition at line 187 of file shape.cc.

References nPoints, pdV, and pdX.

 {
         if (d <= pdX[0]) {
                 return pdV[0];
         }
 
         for (int i = 1; i < nPoints; i++) {
                 if (d <= pdX[i]) {
                         doublereal dl = pdX[i] - pdX[i-1];
                         doublereal dw1 = pdX[i] - d;
                         doublereal dw2 = d - pdX[i - 1];
                         return (pdV[i]*dw2 + pdV[i - 1]*dw1)/dl;
                 }
         }
 
         return pdV[nPoints - 1];
 }

std::ostream & PiecewiseLinearShape1D::Restart ( std::ostream & out ) const

virtual

Implements Shape.

Definition at line 206 of file shape.cc.

References nPoints, pdV, and pdX.

 {
         out << "piecewise linear, " << nPoints;
 
         for (int i = 0; i < nPoints; i++) {
                 out << ", " << pdX[i] << ", " << pdV[i];
         }
 
         return out;
 }