#include <JacSubMatrix.h>

Collaboration diagram for ExpandableMatrix:

Classes
struct	ExpandableColBlock

Public Member Functions
	ExpandableMatrix (void)

	ExpandableMatrix (const integer n, const integer m)

virtual	~ExpandableMatrix (void)

void	ReDim (const integer n, const integer m)

void	SetBlockDim (const integer block, const integer ncols)

void	Zero (void)

integer	GetNRows () const

integer	GetNBlocks () const

integer	GetBlockNCols (const integer block) const

void	Reset (void)

void	Link (const integer block, const ExpandableMatrix *const xp)

void	Link (const integer block, const ExpandableRowVector *const xp)

void	SetBlockIdx (integer block, integer iidx)

void	Set (doublereal xx, integer eq, integer block, integer block_col=1)

void	Set (Vec3 xx, integer eq, integer block, integer block_col=1)

void	Set (Mat3x3 xx, integer eq, integer block, integer block_col=1)

void	Add (doublereal xx, integer eq, integer block, integer block_col=1)

void	Sub (doublereal xx, integer eq, integer block, integer block_col=1)

void	Add (Vec3 xx, integer eq, integer block, integer block_col=1)

void	Sub (Vec3 xx, integer eq, integer block, integer block_col=1)

void	Add (Mat3x3 xx, integer eq, integer block, integer block_col=1)

void	Sub (Mat3x3 xx, integer eq, integer block, integer block_col=1)

void	Add (FullSubMatrixHandler &WM, const integer eq, const doublereal c=1.) const

void	Sub (FullSubMatrixHandler &WM, const integer eq, const doublereal c=1.) const

std::ostream &	Write (std::ostream &out, const char *sFill="") const

Private Member Functions
ExpandableMatrix &	operator= (const ExpandableMatrix &)

	ExpandableMatrix (const ExpandableMatrix &)

Private Attributes
std::vector< ExpandableColBlock >	v

Detailed Description

Definition at line 87 of file JacSubMatrix.h.

Constructor & Destructor Documentation

ExpandableMatrix::ExpandableMatrix ( const ExpandableMatrix & )

private

ExpandableMatrix::ExpandableMatrix ( void )

Definition at line 282 of file JacSubMatrix.cc.

282 {};

ExpandableMatrix::ExpandableMatrix	(	const integer	n,
		const integer	m
	)

Definition at line 283 of file JacSubMatrix.cc.

References ReDim().

                                                                    {
         ReDim(n, m);
 }

Here is the call graph for this function:

ExpandableMatrix::~ExpandableMatrix ( void )

virtual

Definition at line 286 of file JacSubMatrix.cc.

286 {};

Member Function Documentation

void ExpandableMatrix::Add	(	doublereal	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 405 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetBlockNCols(), GetNBlocks(), GetNRows(), and v.

Referenced by Add(), and PlaneHingeJoint::AssJac().

                                                                                       {
         ASSERTMSGBREAK(eq > 0, "ExpandableMatrix::Add() underflow");
         ASSERTMSGBREAK(block > 0, "ExpandableMatrix::Add() underflow");
         ASSERTMSGBREAK(block_col > 0, "ExpandableMatrix::Add() underflow");
         ASSERTMSGBREAK(eq <= GetNRows(), "ExpandableRowVector::Add() overflow");
         ASSERTMSGBREAK(block <= GetNBlocks(), "ExpandableRowVector::Add() overflow");
         ASSERTMSGBREAK(block_col <= GetBlockNCols(block), "ExpandableRowVector::Add() overflow");
         v[block - 1].rows[eq - 1].Add(xx, block_col);
 }

Here is the call graph for this function:

void ExpandableMatrix::Add	(	Vec3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 423 of file JacSubMatrix.cc.

References Add().

                                                                                 {
         for (integer i = 0; i <3; i++) {
                 Add(xx(i + 1), eq + i, block, block_col);
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::Add	(	Mat3x3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 428 of file JacSubMatrix.cc.

References Add().

                                                                                   {
         for (integer i = 0; i <3; i++) {
                 for (integer ii = 0; ii <3; ii++) {
                         Add(xx(i + 1, ii + 1), eq + i, block, block_col + ii);
                 }
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::Add	(	FullSubMatrixHandler &	WM,
		const integer	eq,
		const doublereal	c = `1.`
	)		const

Definition at line 475 of file JacSubMatrix.cc.

References GetNBlocks(), and v.

                                   {
 //      std::cerr << "ExpandableMatrix::Add" << std::endl;
 //              std::cerr << "\teq " << eq << std::endl;
 //              std::cerr << "\tc " << c << std::endl;
         for (integer block = 0; block < GetNBlocks(); block++) {
 //              std::cerr << "\t\tblock " << block << std::endl;
                 v[block].Add(WM, eq, c);
         }
 }

Here is the call graph for this function:

integer ExpandableMatrix::GetBlockNCols ( const integer block ) const

Definition at line 325 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Add(), Link(), Set(), and Sub().

                                                                  {
         ASSERTMSGBREAK(block >= 0, "ExpandableMatrix:GetBlockNCols(), "
                 "block shold be >= 0");
         ASSERTMSGBREAK((unsigned long)block < v.size(), "ExpandableMatrix:GetBlockNCols(), "
                 "block shold be < nblocks");
         return v[block].GetBlockNCols();
 }

integer ExpandableMatrix::GetNBlocks ( ) const

Definition at line 322 of file JacSubMatrix.cc.

References v.

Referenced by Add(), ExpandableMatrix::ExpandableColBlock::Link(), Set(), SetBlockDim(), SetBlockIdx(), Sub(), and Write().

                                            {
         return v.size();
 };

integer ExpandableMatrix::GetNRows ( ) const

Definition at line 315 of file JacSubMatrix.cc.

References v.

Referenced by Add(), Link(), Set(), and Sub().

                                          {
         if (v.begin() != v.end() ) {
                 return v[0].GetBlockNRows(); 
         } else {
                 return 0;
         }
 };

void ExpandableMatrix::Link	(	const integer	block,
		const ExpandableMatrix *const	xp
	)

Definition at line 340 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetBlockNCols(), GetNRows(), and v.

Referenced by PlaneHingeJoint::AssJac().

                                                                             {
         ASSERTMSGBREAK(i > 0, "ExpandableMatrix::Link() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableColBlock>::size_type(i) <= v.size(), "ExpandableMatrix::Link() overflow");
         //FIXME
         //ASSERTMSGBREAK(v[i - 1].idx == 0, "ExpandableMatrix::Link() fatal error");
         ASSERTMSGBREAK(v[i - 1].GetBlockNCols() == xp->GetNRows(), 
                 "ExpandableMatrix::Link() dimension mismatch");
         v[i - 1].Link(xp);
 }

Here is the call graph for this function:

void ExpandableMatrix::Link	(	const integer	block,
		const ExpandableRowVector *const	xp
	)

Definition at line 350 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetBlockNCols(), and v.

                                                                                {
         ASSERTMSGBREAK(i > 0, "ExpandableMatrix::Link() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableColBlock>::size_type(i) <= v.size(), "ExpandableMatrix::Link() overflow");
         //FIXME
         //ASSERTMSGBREAK(v[i - 1].idx == 0, "ExpandableMatrix::Link() fatal error");
         ASSERTMSGBREAK(v[i - 1].GetBlockNCols() == 1, 
                 "ExpandableMatrix::Link() dimension mismatch");
         v[i - 1].Link(xp);
 }

Here is the call graph for this function:

ExpandableMatrix& ExpandableMatrix::operator= ( const ExpandableMatrix & )

private

void ExpandableMatrix::ReDim	(	const integer	n,
		const integer	m
	)

Definition at line 287 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by PlaneHingeJoint::AssJac(), and ExpandableMatrix().

                                                              {
         //we have to accept = 0, some elements do ReDim(0,0) (PointForceElement))
         ASSERTMSGBREAK(n>=0, "ExpandableMatrix:ReDim(), n shold be >= 0");
         ASSERTMSGBREAK(m>=0, "ExpandableMatrix:ReDim(), m shold be >= 0");
         v.resize(m, ecb_Zero);
         for (std::vector<ExpandableColBlock>::iterator i = v.begin(); i != v.end(); ++i) {
                 i->ReDim(n);
         }
 }

void ExpandableMatrix::Reset ( void )

Definition at line 332 of file JacSubMatrix.cc.

References v.

                                  {
         for (std::vector<ExpandableColBlock>::iterator i = v.begin();
                 i != v.end(); ++i)
         {
                 i->Reset();
         }
 }

void ExpandableMatrix::Set	(	doublereal	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 360 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetBlockNCols(), GetNBlocks(), GetNRows(), and v.

Referenced by PlaneHingeJoint::AssJac(), and Set().

                                                                                       {
         ASSERTMSGBREAK(eq > 0, "ExpandableMatrix::Set() underflow");
         ASSERTMSGBREAK(block > 0, "ExpandableMatrix::Set() underflow");
         ASSERTMSGBREAK(block_col > 0, "ExpandableMatrix::Set() underflow");
         ASSERTMSGBREAK(eq <= GetNRows(), "ExpandableRowVector::Set() overflow");
         ASSERTMSGBREAK(block <= GetNBlocks(), "ExpandableRowVector::Set() overflow");
         ASSERTMSGBREAK(block_col <= GetBlockNCols(block), "ExpandableRowVector::Set() overflow");
         v[block - 1].rows[eq - 1].Set(xx, block_col);
 }

Here is the call graph for this function:

void ExpandableMatrix::Set	(	Vec3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 369 of file JacSubMatrix.cc.

References Set().

                                                                                 {
         for (integer i = 0; i <3; i++) {
                 Set(xx(i + 1), eq + i, block, block_col);
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::Set	(	Mat3x3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 374 of file JacSubMatrix.cc.

References Set().

                                                                                   {
         for (integer i = 0; i <3; i++) {
                 for (integer ii = 0; ii <3; ii++) {
                         Set(xx(i + 1, ii + 1), eq + i, block, block_col + ii);
                 }
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::SetBlockDim	(	const integer	block,
		const integer	ncols
	)

Definition at line 296 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetNBlocks(), and v.

Referenced by PlaneHingeJoint::AssJac().

                                                                            {
         ASSERTMSGBREAK(block>0, "ExpandableMatrix:SetBlockDim(), block shold be > 0");
         ASSERTMSGBREAK(block <= GetNBlocks(), "ExpandableMatrix:SetBlockDim(), block shold be "
                 "<= GetNBlocks()");
 //      std::cerr << "ExpandableMatrix::SetBlockDim(" << block <<", " << ncols << ")" << std::endl;
         int pippo = 0;
         for (std::vector<ExpandableRowVector>::iterator i = v[block-1].rows.begin();
                 i != v[block-1].rows.end(); ++i)
         {
                         pippo++;
 //                      std::cerr << "\t" << pippo << std::endl;
                         i->ReDim(ncols);
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::SetBlockIdx	(	integer	block,
		integer	iidx
	)

Definition at line 381 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetNBlocks(), and v.

                                                               {
         ASSERTMSGBREAK(block > 0, "ExpandableMatrix::SetColIdx() underflow");
         ASSERTMSGBREAK(block <= GetNBlocks(), "ExpandableMatrix::SetColIdx() overflow");
 //      std::cerr << "ExpandableMatrix::SetBlockIdx" << std::endl;
         v[block - 1].SetColIdx(iidx);
 }

Here is the call graph for this function:

void ExpandableMatrix::Sub	(	doublereal	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 414 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, GetBlockNCols(), GetNBlocks(), GetNRows(), and v.

Referenced by PlaneHingeJoint::AssJac(), and Sub().

                                                                                       {
         ASSERTMSGBREAK(eq > 0, "ExpandableMatrix::Sub() underflow");
         ASSERTMSGBREAK(block > 0, "ExpandableMatrix::Sub() underflow");
         ASSERTMSGBREAK(block_col > 0, "ExpandableMatrix::Sub() underflow");
         ASSERTMSGBREAK(eq <= GetNRows(), "ExpandableRowVector::Sub() overflow");
         ASSERTMSGBREAK(block <= GetNBlocks(), "ExpandableRowVector::Sub() overflow");
         ASSERTMSGBREAK(block_col <= GetBlockNCols(block), "ExpandableRowVector::Sub() overflow");
         v[block - 1].rows[eq - 1].Sub(xx, block_col);
 }

Here is the call graph for this function:

void ExpandableMatrix::Sub	(	Vec3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 435 of file JacSubMatrix.cc.

References Sub().

                                                                                 {
         for (integer i = 0; i <3; i++) {
                 Sub(xx(i + 1), eq + i, block, block_col);
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::Sub	(	Mat3x3	xx,
		integer	eq,
		integer	block,
		integer	block_col = `1`
	)

Definition at line 440 of file JacSubMatrix.cc.

References Sub().

                                                                                   {
         for (integer i = 0; i <3; i++) {
                 for (integer ii = 0; ii <3; ii++) {
                         Sub(xx(i + 1, ii + 1), eq + i, block, block_col + ii);
                 }
         }
 }

Here is the call graph for this function:

void ExpandableMatrix::Sub	(	FullSubMatrixHandler &	WM,
		const integer	eq,
		const doublereal	c = `1.`
	)		const

Definition at line 508 of file JacSubMatrix.cc.

References GetNBlocks(), and v.

                                   {
         for (integer block = 0; block < GetNBlocks(); block++) {
                 v[block].Sub(WM, eq, c);
         }
 }

Here is the call graph for this function:

std::ostream & ExpandableMatrix::Write	(	std::ostream &	out,
		const char *	sFill = `""`
	)		const

Definition at line 538 of file JacSubMatrix.cc.

References GetNBlocks(), and v.

Referenced by operator<<().

                                                                              {
         for (integer block = 0; block < GetNBlocks(); block++) {
                 out << "Block: " <<  block << std::endl;
                 v[block].Write(out, "\t");
         }
         return out;
 }

Here is the call graph for this function:

void ExpandableMatrix::Zero ( void )

Definition at line 310 of file JacSubMatrix.cc.

References v.

                                 {
         for (std::vector<ExpandableColBlock>::iterator i = v.begin(); i != v.end(); ++i) {
                 i->Zero();
         }
 }

Member Data Documentation

std::vector<ExpandableColBlock> ExpandableMatrix::v

private

Definition at line 89 of file JacSubMatrix.h.

Referenced by Add(), GetBlockNCols(), GetNBlocks(), GetNRows(), ExpandableMatrix::ExpandableColBlock::Link(), Link(), ReDim(), Reset(), Set(), SetBlockDim(), SetBlockIdx(), Sub(), Write(), and Zero().

The documentation for this class was generated from the following files:

Classes

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation