#include <JacSubMatrix.h>

Collaboration diagram for ExpandableRowVector:

Classes
struct	ExpandableRowElement

Public Member Functions
	ExpandableRowVector (void)

	ExpandableRowVector (const integer n)

virtual	~ExpandableRowVector (void)

void	ReDim (const integer n)

integer	GetDim (void) const

void	Zero (void)

void	Reset (void)

void	Link (const integer i, const ExpandableRowVector *const xp, const integer rhs_block=1)

void	SetIdx (integer i, integer iidx)

void	Set (doublereal xx, integer i, integer iidx)

void	Set (Vec3 xx, integer i, integer iidx)

doublereal &	operator() (integer i)

const doublereal &	operator() (integer i) const

void	Set (doublereal xx, integer i)

void	Set (Vec3 xx, integer i)

void	Add (doublereal xx, integer i)

void	Sub (doublereal xx, integer i)

void	Add (SubVectorHandler &WorkVec, const doublereal c=1.) const

void	Sub (SubVectorHandler &WorkVec, const doublereal c=1.) const

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

void	Add (FullSubMatrixHandler &WM, const std::vector< integer > &eq, const std::vector< doublereal > &cc, const doublereal c=1.) const

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

void	Sub (FullSubMatrixHandler &WM, const std::vector< integer > &eq, const std::vector< doublereal > &cc, const doublereal c=1.) const

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

Private Attributes
std::vector< ExpandableRowElement >	v

Detailed Description

Definition at line 42 of file JacSubMatrix.h.

Constructor & Destructor Documentation

ExpandableRowVector::ExpandableRowVector ( void )

Definition at line 45 of file JacSubMatrix.cc.

45 {};

ExpandableRowVector::ExpandableRowVector ( const integer n )

Definition at line 46 of file JacSubMatrix.cc.

References ReDim().

                                                         {
         ReDim(n);
 }

Here is the call graph for this function:

ExpandableRowVector::~ExpandableRowVector ( void )

virtual

Definition at line 49 of file JacSubMatrix.cc.

49 {};

Member Function Documentation

void ExpandableRowVector::Add	(	doublereal	xx,
		integer	i
	)

Definition at line 119 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Brake::AssJac(), ModLugreFriction::AssJac(), DiscreteCoulombFriction::AssJac(), and AxialRotationJoint::AssJac().

                                                       {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::Add() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::Add() overflow");
         v[i - 1].x += xx;
 }

void ExpandableRowVector::Add	(	SubVectorHandler &	WorkVec,
		const doublereal	c = `1.`
	)		const

Definition at line 129 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, VectorHandler::IncCoef(), and v.

                                                                                  {
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 if (v[i].idx != 0) {
                         WorkVec.IncCoef(v[i].idx, c*v[i].x);
                 } else {
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Add() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Add(WorkVec, c*v[i].x);
                         }       
                 }
         }
 }

Here is the call graph for this function:

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

Definition at line 161 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, FullSubMatrixHandler::IncCoef(), and v.

                                   {
 //      std::cerr << "\t\t\tExpandableRowVector::Add" << std::endl;
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
 //              std::cerr << "\t\t\t\ti " << i << std::endl;
 
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 if (v[i].idx != 0) {
 //                      std::cerr << "Adding("<< eq << ", " << v[i].idx << ") = " <<
 //                      c*v[i].x << std::endl;
 //                      std::cerr << "\t\t\t\tdentro if" << std::endl;
                         WM.IncCoef(eq, v[i].idx, c*v[i].x);
                 } else {
 //                      std::cerr << "\t\t\t\tdentro else" << std::endl;
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
 //                              std::cerr << "\t\t\t\t\trhs_block " << rhs_block << std::endl;
 //                              std::cerr << "\t\t\t\t\t\trhs_block_ptr " << v[i].xm[rhs_block] << std::endl;
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Add() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Add(WM, eq, c*v[i].x);
                         }
                 }
         }
 //      std::cerr << "\t\t\t ---- ExpandableRowVector::Add" << std::endl;
 }

Here is the call graph for this function:

void ExpandableRowVector::Add	(	FullSubMatrixHandler &	WM,
		const std::vector< integer > &	eq,
		const std::vector< doublereal > &	cc,
		const doublereal	c = `1.`
	)		const

Definition at line 189 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, FullSubMatrixHandler::IncCoef(), and v.

 {
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 integer idx = v[i].idx;
                 if (idx != 0) {
                         doublereal d = c*v[i].x;
                         for (std::vector<integer>::size_type j = 0; j < eq.size(); j++) {
                                 WM.IncCoef(eq[j], idx, cc[j]*d);
                         }
                 } else {
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Add() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Add(WM, eq, cc, c*v[i].x);
                         }
                 }
         }
 }

Here is the call graph for this function:

integer ExpandableRowVector::GetDim ( void ) const

Definition at line 57 of file JacSubMatrix.cc.

References v.

                                               {
         return v.size();
 };

void ExpandableRowVector::Link	(	const integer	i,
		const ExpandableRowVector *const	xp,
		const integer	rhs_block = `1`
	)

Definition at line 68 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Brake::AssJac(), ModLugreFriction::AssJac(), DiscreteCoulombFriction::AssJac(), AxialRotationJoint::AssJac(), SimpleShapeCoefficient::dSh_c(), SimplePlaneHingeJointSh_c::dSh_c(), and ScrewJointSh_c::dSh_c().

                                                                                                            {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::Link() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::Link() overflow");
         ASSERTMSGBREAK(v[i - 1].idx == 0, "ExpandableRowVector::Link() fatal error");
         if (std::vector<ExpandableRowVector *>::size_type(rhs_block) > v[ i - 1].xm.size()) {
                 v[ i - 1].xm.resize(rhs_block);
         }
         v[i - 1].xm[rhs_block - 1] = xp;
 }

doublereal & ExpandableRowVector::operator() ( integer i )

Definition at line 106 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

 {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::() overflow");
         return v[i - 1].x;
 }

const doublereal & ExpandableRowVector::operator() ( integer i ) const

Definition at line 113 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

 {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::() overflow");
         return v[i - 1].x;
 }

void ExpandableRowVector::ReDim ( const integer n )

Definition at line 50 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Brake::AssJac(), PlaneHingeJoint::AssJac(), ModLugreFriction::AssJac(), DiscreteCoulombFriction::AssJac(), AxialRotationJoint::AssJac(), SimpleShapeCoefficient::dSh_c(), SimplePlaneHingeJointSh_c::dSh_c(), ScrewJointSh_c::dSh_c(), and ExpandableRowVector().

                                                {
         //we have to accept = 0, some elements do ReDim(0,0) (PointForceElement))
         ASSERTMSGBREAK(n>=0, "ExpandableRowVector:ReDim(), n shold be >= 0");
 //      std::cerr << "ExpandableRowVector::ReDim("  << n << ")" << std::endl;
 //      std::cerr << "----------------------" << std::endl;
         v.resize(n, er_Zero);
 }

void ExpandableRowVector::Reset ( void )

Definition at line 65 of file JacSubMatrix.cc.

References er_Zero, and v.

                                     {
         std::fill(v.begin(), v.end(), er_Zero);
 }

void ExpandableRowVector::Set	(	doublereal	xx,
		integer	i,
		integer	iidx
	)

Definition at line 95 of file JacSubMatrix.cc.

References SetIdx().

Referenced by Brake::AssJac(), PlaneHingeJoint::AssJac(), ModLugreFriction::AssJac(), DiscreteCoulombFriction::AssJac(), AxialRotationJoint::AssJac(), SimpleShapeCoefficient::dSh_c(), SimplePlaneHingeJointSh_c::dSh_c(), ScrewJointSh_c::dSh_c(), and Set().

                                                                     {
         Set(xx, i);
         SetIdx(i, iidx);
 }

Here is the call graph for this function:

void ExpandableRowVector::Set	(	Vec3	xx,
		integer	i,
		integer	iidx
	)

Definition at line 99 of file JacSubMatrix.cc.

References Set(), and SetIdx().

                                                               {
         Set(xx, i);
         SetIdx(i, iidx);
         SetIdx(i+1, iidx+1);
         SetIdx(i+2, iidx+2);
 }

Here is the call graph for this function:

void ExpandableRowVector::Set	(	doublereal	xx,
		integer	i
	)

Definition at line 77 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

                                                       {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::Set() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::Set() overflow");
         v[i - 1].x = xx;
 }

void ExpandableRowVector::Set	(	Vec3	xx,
		integer	i
	)

Definition at line 82 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

                                                 {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::Set() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i+3) <= v.size(), "ExpandableRowVector::Set() overflow");
         v[i - 1].x = xx(1);
         v[i - 0].x = xx(2);
         v[i + 1].x = xx(3);
 }

void ExpandableRowVector::SetIdx	(	integer	i,
		integer	iidx
	)

Definition at line 89 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Set().

                                                         {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::SetIdx() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::SetIdx() overflow");
 //      std::cerr << "ExpandableRowVector::SetIdx(" << i <<", " <<  iidx << ")" << std::endl;
         v[i - 1].idx = iidx;
 }

void ExpandableRowVector::Sub	(	doublereal	xx,
		integer	i
	)

Definition at line 124 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, and v.

Referenced by Brake::AssJac(), DiscreteCoulombFriction::AssJac(), and AxialRotationJoint::AssJac().

                                                       {
         ASSERTMSGBREAK(i > 0, "ExpandableRowVector::Sub() underflow");
         ASSERTMSGBREAK(std::vector<ExpandableRowVector>::size_type(i) <= v.size(), "ExpandableRowVector::Sub() overflow");
         v[i - 1].x -= xx;
 }

void ExpandableRowVector::Sub	(	SubVectorHandler &	WorkVec,
		const doublereal	c = `1.`
	)		const

Definition at line 145 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, VectorHandler::DecCoef(), and v.

                                                                                  {
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 if (v[i].idx != 0) {
                         WorkVec.DecCoef(v[i].idx, c*v[i].x);
                 } else {
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Sub() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Sub(WorkVec, c*v[i].x);
                         }
                 }
         }
 }

Here is the call graph for this function:

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

Definition at line 213 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, FullSubMatrixHandler::DecCoef(), and v.

                                   {
         for (std::vector<ExpandableRowElement>::size_type i = 0 ; i < v.size(); i++) {
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 if (v[i].idx != 0) {
                         WM.DecCoef(eq, v[i].idx, c*v[i].x);
                 } else {
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Sub() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Sub(WM, eq, c*v[i].x);
                         }
                 }
         }
 }

Here is the call graph for this function:

void ExpandableRowVector::Sub	(	FullSubMatrixHandler &	WM,
		const std::vector< integer > &	eq,
		const std::vector< doublereal > &	cc,
		const doublereal	c = `1.`
	)		const

Definition at line 231 of file JacSubMatrix.cc.

References ASSERTMSGBREAK, FullSubMatrixHandler::DecCoef(), and v.

 {
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 if (v[i].x == 0.) {
                         continue;
                 }
 
                 integer idx = v[i].idx;
                 if (idx != 0) {
                         doublereal d = c*v[i].x;
                         for (std::vector<integer>::size_type j = 0; j < eq.size(); j++) {
                                 WM.DecCoef(eq[j], idx, cc[j]*d);
                         }
                 } else {
                         for (std::vector<const ExpandableRowVector *>::size_type rhs_block = 0; rhs_block < v[i].xm.size(); rhs_block++) {
                                 ASSERTMSGBREAK(v[i].xm[rhs_block] != 0, "ExpandableRowVector::Add() null pointer to ExpandableRowVector");
                                 v[i].xm[rhs_block]->Sub(WM, eq, cc, c*v[i].x);
                         }
                 }
         }
 }

Here is the call graph for this function:

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

Definition at line 256 of file JacSubMatrix.cc.

References v.

Referenced by operator<<().

                                                                                 {
         out << "LocalDof: ";
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 if (v[i].idx != 0) {
                         out << sFill << std::setw(12) << v[i].idx;
                 } else {
                         out << sFill << std::setw(12) << "linked";
                 }               
         }
         out << std::endl;
         out << "   Value: ";
         for (std::vector<ExpandableRowElement>::size_type i = 0; i < v.size(); i++) {
                 out << sFill << std::setw(12) << v[i].x;
         }
         out << std::endl;
         return out;
 }

void ExpandableRowVector::Zero ( void )

Definition at line 60 of file JacSubMatrix.cc.

References v.

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

Member Data Documentation

std::vector<ExpandableRowElement> ExpandableRowVector::v

private

Definition at line 51 of file JacSubMatrix.h.

Referenced by Add(), GetDim(), Link(), operator()(), ReDim(), Reset(), Set(), SetIdx(), Sub(), Write(), and Zero().

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

Classes

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation