#include <matvecass.h>

Inheritance diagram for grad::GradientAssVec< Gradient< N_SIZE > >:

Collaboration diagram for grad::GradientAssVec< Gradient< N_SIZE > >:

Public Member Functions
	GradientAssVec (SparseSubMatrixHandler &mh, enum mode_t mode=RESET)

void	AddItem (integer iRow, const Gradient< N_SIZE > &g)

template<index_type N_rows>
void	AddItem (integer iFirstRow, const Vector< Gradient< N_SIZE >, N_rows > &v)

Static Public Member Functions
template<typename T >
static void	AssJac (T pElem, SparseSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, enum FunctionCall func, LocalDofMap pDofMap, enum mode_t mode=RESET)

template<typename T >
static void	InitialAssJac (T pElem, SparseSubMatrixHandler &WorkMat, const VectorHandler &XCurr, enum FunctionCall func, LocalDofMap pDofMap, enum mode_t mode=RESET)

Private Attributes
SparseSubMatrixHandler &	WorkMat

integer	iNextItem

Additional Inherited Members
Public Types inherited from grad::GradientAssVecBase
enum	mode_t { RESET, APPEND }

Detailed Description

template<index_type N_SIZE>
class grad::GradientAssVec< Gradient< N_SIZE > >

Definition at line 200 of file matvecass.h.

Constructor & Destructor Documentation

template<index_type N_SIZE>

grad::GradientAssVec< Gradient< N_SIZE > >::GradientAssVec	(	SparseSubMatrixHandler &	mh,
		enum mode_t	mode = `RESET`
	)

inline

Definition at line 202 of file matvecass.h.

References GRADIENT_ASSERT.

                 :WorkMat(mh) {
 
                 switch (mode) {
                 case RESET:
                         iNextItem = 1;
                         WorkMat.ResizeReset(iNextItem, 0);
                         WorkMat.PutItem(1, 1, 1, 0.); //FIXME: avoid SIGSEGV if the matrix is empty
                         break;
 
                 case APPEND:
                         iNextItem = WorkMat.iGetNumRows() + 1;
                         break;
 
                 default:
                         GRADIENT_ASSERT(0);
                 }
 
         }

Member Function Documentation

template<index_type N_SIZE>

void grad::GradientAssVec< Gradient< N_SIZE > >::AddItem	(	integer	iRow,
		const Gradient< N_SIZE > &	g
	)

inline

Definition at line 263 of file matvecass.h.

References grad::Gradient< N_SIZE >::dGetDerivativeLocal(), GRADIENT_ASSERT, grad::Gradient< N_SIZE >::iGetEndIndexLocal(), grad::Gradient< N_SIZE >::iGetGlobalDof(), and grad::Gradient< N_SIZE >::iGetStartIndexLocal().

                                                               {
                 for (index_type i = g.iGetStartIndexLocal(); i < g.iGetEndIndexLocal(); ++i) {
                         const double dCoef = g.dGetDerivativeLocal(i);
                 if (dCoef) {
                         const index_type iDofIndex = g.iGetGlobalDof(i);
                     WorkMat.Resize(iNextItem, 0);
                         GRADIENT_ASSERT(iNextItem <= WorkMat.iGetNumRows());
 #if GRADIENT_DEBUG > 0 && HAVE_FINITE == 1
                     if (!std::isfinite(dCoef)) {
                         silent_cerr("PutItem(" << iNextItem << ", " << iRow << ", " << iDofIndex << ", " << dCoef << ")" << std::endl);
                     }
 
                     GRADIENT_ASSERT(std::isfinite(dCoef));
 #endif
                         WorkMat.PutItem(iNextItem++, iRow, iDofIndex, dCoef);
                 }
         }
         }

Here is the call graph for this function:

template<index_type N_SIZE>

template<index_type N_rows>

void grad::GradientAssVec< Gradient< N_SIZE > >::AddItem	(	integer	iFirstRow,
		const Vector< Gradient< N_SIZE >, N_rows > &	v
	)

inline

Definition at line 283 of file matvecass.h.

                                                                                    {
                 // zero based index according to VectorHandler::Put(integer iRow, const Vec3& v)
 
                 for (integer i = 0; i < N_rows; ++i) {
                         AddItem(iFirstRow + i, v(i + 1));
                 }
         }

template<index_type N_SIZE>

template<typename T >

static void grad::GradientAssVec< Gradient< N_SIZE > >::AssJac	(	T *	pElem,
		SparseSubMatrixHandler &	WorkMat,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr,
		enum FunctionCall	func,
		LocalDofMap *	pDofMap,
		enum mode_t	mode = `RESET`
	)

inlinestatic

Definition at line 223 of file matvecass.h.

References grad::LocalDofMap::Reset().

                                                                      {
 
                 const GradientVectorHandler<Gradient<N_SIZE> > XCurr_grad(XCurr);
                 const GradientVectorHandler<Gradient<N_SIZE> > XPrimeCurr_grad(XPrimeCurr);
 
                 GradientAssVec WorkMat_grad(WorkMat, mode);
 
                 if (pDofMap) {
                         pDofMap->Reset(func);
                 }
 
                 pElem->AssRes(WorkMat_grad, dCoef, XCurr_grad, XPrimeCurr_grad, func);
         }

Here is the call graph for this function:

template<index_type N_SIZE>

template<typename T >

static void grad::GradientAssVec< Gradient< N_SIZE > >::InitialAssJac	(	T *	pElem,
		SparseSubMatrixHandler &	WorkMat,
		const VectorHandler &	XCurr,
		enum FunctionCall	func,
		LocalDofMap *	pDofMap,
		enum mode_t	mode = `RESET`
	)

inlinestatic

Definition at line 245 of file matvecass.h.

References grad::LocalDofMap::Reset().

                                                                                     {
 
                 const GradientVectorHandler<Gradient<N_SIZE> > XCurr_grad(XCurr);
 
                 GradientAssVec WorkMat_grad(WorkMat, mode);
 
                 if (pDofMap) {
                         pDofMap->Reset(func);
                 }
 
                 pElem->InitialAssRes(WorkMat_grad, XCurr_grad, func);
         }