#include <modalext.h>

Inheritance diagram for ModalExt:

Collaboration diagram for ModalExt:

Public Member Functions
	ModalExt (unsigned int uL, DataManager pDM, const Modal pmodal, const StructNode pnode, bool bOutputAccelerations, ExtFileHandlerBase pEFH, ExtModalForceBase *pEMF, bool bSendAfterPredict, int iCoupling, ExtModalForceBase::BitMask bm, flag fOut)

virtual	~ModalExt (void)

virtual Force::Type	GetForceType (void) const

void	WorkSpaceDim (integer piNumRows, integer piNumCols) const

SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	Output (OutputHandler &OH) const

virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const

Public Member Functions inherited from Elem
	Elem (unsigned int uL, flag fOut)

virtual	~Elem (void)

virtual unsigned int	iGetNumDof (void) const

virtual std::ostream &	DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const

virtual void	DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const

virtual std::ostream &	DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const

virtual void	DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const

virtual DofOrder::Order	GetDofType (unsigned int) const

virtual void	AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual bool	bInverseDynamics (void) const

void	SetInverseDynamicsFlags (unsigned uIDF)

unsigned	GetInverseDynamicsFlags (void) const

bool	bIsErgonomy (void) const

bool	bIsRightHandSide (void) const

virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)

virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)

virtual int	GetNumConnectedNodes (void) const

Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")

virtual	~WithLabel (void)

void	PutLabel (unsigned int uL)

void	PutName (const std::string &sN)

unsigned int	GetLabel (void) const

const std::string &	GetName (void) const

Public Member Functions inherited from SimulationEntity
	SimulationEntity (void)

virtual	~SimulationEntity (void)

virtual bool	bIsValidIndex (unsigned int i) const

virtual DofOrder::Order	GetEqType (unsigned int i) const

virtual Hint *	ParseHint (DataManager pDM, const char s) const

virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const

virtual void	DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)

virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)

virtual unsigned int	iGetNumPrivData (void) const

virtual unsigned int	iGetPrivDataIdx (const char *s) const

virtual doublereal	dGetPrivData (unsigned int i) const

virtual std::ostream &	OutputAppend (std::ostream &out) const

virtual void	ReadInitialState (MBDynParser &HP)

Public Member Functions inherited from ToBeOutput
	ToBeOutput (flag fOut=fDefaultOut)

virtual	~ToBeOutput (void)

virtual void	OutputPrepare (OutputHandler &OH)

virtual void	Output (OutputHandler &OH, const VectorHandler &X, const VectorHandler &XP) const

virtual flag	fToBeOutput (void) const

virtual bool	bToBeOutput (void) const

virtual void	SetOutputFlag (flag f=flag(1))

Public Member Functions inherited from ExtForce
	ExtForce (unsigned int uL, DataManager pDM, ExtFileHandlerBase pEFH, bool bSendAfterPredict, int iCoupling, flag fOut)

virtual	~ExtForce (void)

virtual void	SetValue (DataManager pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints h=0)

virtual void	Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)

virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)

virtual void	InitialWorkSpaceDim (integer piNumRows, integer piNumCols) const

virtual VariableSubMatrixHandler &	InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)

virtual SubVectorHandler &	InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)

Public Member Functions inherited from Force
	Force (unsigned int uL, flag fOut)

virtual	~Force (void)

virtual Elem::Type	GetElemType (void) const

virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal, const VectorHandler &, const VectorHandler &)

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

virtual unsigned int	iGetInitialNumDof (void) const

Public Member Functions inherited from InitialAssemblyElem
	InitialAssemblyElem (unsigned int uL, flag fOut)

virtual	~InitialAssemblyElem (void)

Public Member Functions inherited from SubjectToInitialAssembly
	SubjectToInitialAssembly (void)

virtual	~SubjectToInitialAssembly (void)

Protected Member Functions
bool	Prepare (ExtFileHandlerBase *pEFH)

void	Send (ExtFileHandlerBase *pEFH, ExtFileHandlerBase::SendWhen when)

void	Recv (ExtFileHandlerBase *pEFH)

Protected Member Functions inherited from ExtForce
void	Send (ExtFileHandlerBase::SendWhen when)

void	Recv (void)

Protected Attributes
const Modal *	pModal

const StructNode *	pNode

bool	bOutputAccelerations

ExtModalForceBase *	pEMF

unsigned	uFlags

Vec3	F

Vec3	M

std::vector< doublereal >	f

std::vector< doublereal >	q

std::vector< doublereal >	qP

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

Protected Attributes inherited from ToBeOutput
flag	fOutput

Protected Attributes inherited from ExtForce
Converged	c

ExtFileHandlerBase *	pEFH

bool	bSendAfterPredict

int	iCoupling

int	iCouplingCounter

bool	bFirstSend

bool	bFirstRecv

Additional Inherited Members
Public Types inherited from Elem
enum	Type { UNKNOWN = -1, AIRPROPERTIES = 0, INDUCEDVELOCITY, AUTOMATICSTRUCTURAL, GRAVITY, BODY, JOINT, JOINT_REGULARIZATION, BEAM, PLATE, FORCE, INERTIA, ELECTRICBULK, ELECTRIC, THERMAL, HYDRAULIC, BULK, LOADABLE, DRIVEN, EXTERNAL, AEROMODAL, AERODYNAMIC, GENEL, SOCKETSTREAM_OUTPUT, RTAI_OUTPUT = SOCKETSTREAM_OUTPUT, LASTELEMTYPE }

Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints

Public Types inherited from ToBeOutput
enum	{ OUTPUT = 0x1U, OUTPUT_MASK = 0xFU, OUTPUT_PRIVATE = 0x10U, OUTPUT_PRIVATE_MASK = ~OUTPUT_MASK }

Public Types inherited from ExtForce
enum	{ COUPLING_NONE = -2, COUPLING_STAGGERED = -1, COUPLING_LOOSE = 0, COUPLING_TIGHT = 1 }

Public Types inherited from Force
enum	Type { UNKNOWN = -1, ABSTRACTFORCE = 0, ABSTRACTINTERNALFORCE, ABSOLUTEDISPFORCE, ABSOLUTEINTERNALDISPFORCE, ABSOLUTEFORCE, FOLLOWERFORCE, ABSOLUTECOUPLE, FOLLOWERCOUPLE, ABSOLUTEINTERNALFORCE, FOLLOWERINTERNALFORCE, ABSOLUTEINTERNALCOUPLE, FOLLOWERINTERNALCOUPLE, TOTALFORCE, TOTALINTERNALFORCE, EXTERNALSTRUCTURAL, MODALFORCE, EXTERNALMODAL, LASTFORCETYPE }

Detailed Description

Definition at line 144 of file modalext.h.

Constructor & Destructor Documentation

ModalExt::ModalExt	(	unsigned int	uL,
		DataManager *	pDM,
		const Modal *	pmodal,
		const StructNode *	pnode,
		bool	bOutputAccelerations,
		ExtFileHandlerBase *	pEFH,
		ExtModalForceBase *	pEMF,
		bool	bSendAfterPredict,
		int	iCoupling,
		ExtModalForceBase::BitMask	bm,
		flag	fOut
	)

Definition at line 333 of file modalext.cc.

References ASSERT, ExtModalForceBase::EMF_ALL, ExtModalForceBase::EMF_MODAL, ExtModalForceBase::EMF_RIGID, f, pModal, pNode, q, qP, uFlags, and Modal::uGetNModes().

 : Elem(uL, fOut),
 ExtForce(uL, pDM, pEFH, bSendAfterPredict, iCoupling, fOut),
 pModal(pmodal),
 pNode(pnode),
 bOutputAccelerations(bOutputAccelerations),
 pEMF(pEMF),
 uFlags(ExtModalForceBase::EMF_NONE),
 F(Zero3),
 M(Zero3)
 {
         ASSERT((pModal != 0) || (pNode != 0));
         ASSERT(pEMF != 0);
 
         ASSERT((bm & ExtModalForceBase::EMF_ALL) != 0);
         ASSERT((bm & ~ExtModalForceBase::EMF_ALL) == 0);
 
         if (bm & ExtModalForceBase::EMF_MODAL) {
                 f.resize(pModal->uGetNModes());
 
                 // Temporary?
                 q.resize(pModal->uGetNModes());
                 qP.resize(pModal->uGetNModes());
 
                 ASSERT(pModal->uGetNModes() > 0);
                 uFlags |= ExtModalForceBase::EMF_MODAL;
 
         } else {
                 f.resize(0);
                 q.resize(0);
                 qP.resize(0);
         }
 
         if (bm & ExtModalForceBase::EMF_RIGID) {
                 ASSERT(pNode != 0);
                 uFlags |= ExtModalForceBase::EMF_RIGID;
         }
 }

Here is the call graph for this function:

ModalExt::~ModalExt ( void )

virtual

Definition at line 382 of file modalext.cc.

References pEMF, and SAFEDELETE.

 {
         if (pEMF) {
                 SAFEDELETE(pEMF);
         }
 }

Member Function Documentation

SubVectorHandler & ModalExt::AssRes	(	SubVectorHandler &	WorkVec,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr
	)

virtual

Implements Elem.

Definition at line 507 of file modalext.cc.

References ExtForce::COUPLING_NONE, ExtModalForceBase::EMF_MODAL, ExtModalForceBase::EMF_RIGID, F, f, ExtForce::iCoupling, StructDispNode::iGetFirstMomentumIndex(), Modal::iGetModalIndex(), M, pModal, pNode, VectorHandler::Put(), SubVectorHandler::PutItem(), SubVectorHandler::PutRowIndex(), ExtForce::Recv(), VectorHandler::Resize(), VectorHandler::ResizeReset(), uFlags, Modal::uGetNModes(), and WorkSpaceDim().

 {
         ExtForce::Recv();
 
         if (iCoupling == COUPLING_NONE) {
                 WorkVec.Resize(0);
                 return WorkVec;
         }
 
         integer iR, iC;
         WorkSpaceDim(&iR, &iC);
         WorkVec.ResizeReset(iR);
 
         integer iIdx = 1;
         if (uFlags & ExtModalForceBase::EMF_RIGID) {
                 integer iFirstIndex = pNode->iGetFirstMomentumIndex();
 
                 for (integer iCnt = 1; iCnt <= 6; iCnt++) {
                         WorkVec.PutRowIndex(iCnt, iFirstIndex + iCnt);
                 }
 
                 WorkVec.Put(1, F);
                 WorkVec.Put(4, M);
         }
 
         if (pNode) {
                 iIdx += 6;
         }
 
         if (uFlags & ExtModalForceBase::EMF_MODAL) {
                 integer iModalIndex = pModal->iGetModalIndex() + pModal->uGetNModes() + 1;
                 for (integer iMode = 0; iMode < pModal->uGetNModes(); iMode++) {
                         WorkVec.PutItem(iIdx + iMode, iModalIndex + iMode, f[iMode]);
                 }
         }
 
         return WorkVec;
 }

Here is the call graph for this function:

void ModalExt::GetConnectedNodes ( std::vector< const Node * > & connectedNodes ) const

virtual

Reimplemented from Elem.

Definition at line 578 of file modalext.cc.

References pNode.

 {
         if (pNode) {
                 connectedNodes.resize(1);
                 connectedNodes[0] = pNode;
         } else {
                 connectedNodes.resize(0);
         }
 }

Force::Type ModalExt::GetForceType ( void ) const

virtual

Implements Force.

Definition at line 488 of file modalext.cc.

References Force::EXTERNALMODAL.

 {
         return Force::EXTERNALMODAL;
 }

void ModalExt::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 550 of file modalext.cc.

References ToBeOutput::bToBeOutput(), ExtModalForceBase::EMF_MODAL, ExtModalForceBase::EMF_RIGID, F, f, OutputHandler::FORCES, OutputHandler::Forces(), WithLabel::GetLabel(), Modal::GetModeList(), M, pModal, pNode, uFlags, and OutputHandler::UseText().

 {
         if (bToBeOutput()) {
                 if (OH.UseText(OutputHandler::FORCES)) {
                         std::ostream& out = OH.Forces();
 
                         if (uFlags & ExtModalForceBase::EMF_RIGID) {
                                 out << GetLabel() << "#" << pNode->GetLabel()
                                         << " " << F << " " << M
                                         << std::endl;
                         }
 
                         if (uFlags & ExtModalForceBase::EMF_MODAL) {
                                 unsigned cnt = 0;
                                 const std::vector<unsigned int>& ML(pModal->GetModeList());
                                 for (std::vector<doublereal>::const_iterator i = f.begin(); i != f.end(); ++i) {
                                         out << GetLabel() << '.' << ML[cnt]
                                                 << " " << *i << std::endl;
                                         cnt++;
                                 }
                         }
                 }
 
                 /* TODO: NetCDF */
         }
 }

Here is the call graph for this function:

bool ModalExt::Prepare ( ExtFileHandlerBase * pEFH )

protectedvirtual

Implements ExtForce.

Definition at line 390 of file modalext.cc.

References ExtModalForceBase::EMF_RIGID, WithLabel::GetLabel(), pEMF, pModal, ExtModalForceBase::Prepare(), uFlags, and Modal::uGetNModes().

 {
         return pEMF->Prepare(pEFH, GetLabel(),
                 uFlags & ExtModalForceBase::EMF_RIGID,
                 pModal ? pModal->uGetNModes() : 0);
 }

Here is the call graph for this function:

void ModalExt::Recv ( ExtFileHandlerBase * pEFH )

protectedvirtual

Implements ExtForce.

Definition at line 467 of file modalext.cc.

References ExtModalForceBase::EMF_ERR, F, f, WithLabel::GetLabel(), M, MBDYN_EXCEPT_ARGS, pEMF, ExtModalForceBase::Recv(), uFlags, and WithLabel::uLabel.

 {
         unsigned uLabel = 0;
         unsigned uOutFlags = pEMF->Recv(pEFH, uFlags, uLabel, F, M, f);
 
         if (uOutFlags & ExtModalForceBase::EMF_ERR) {
                 silent_cerr("ModalExt(" << GetLabel() << "): "
                         "error while reading modal force data" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         if (uFlags != uOutFlags) {
                 silent_cerr("ModalExt(" << GetLabel() << "): "
                         "error while reading modal force data "
                         "(" << uOutFlags << "!=" << uFlags << ")"
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 }

Here is the call graph for this function:

void ModalExt::Send	(	ExtFileHandlerBase *	pEFH,
		ExtFileHandlerBase::SendWhen	when
	)

protectedvirtual

Implements ExtForce.

Definition at line 401 of file modalext.cc.

References a, bOutputAccelerations, ExtModalForceBase::EMF_MODAL, ExtModalForceBase::EMF_RIGID, Eye3, Modal::GetA(), Modal::GetB(), Modal::GetBP(), WithLabel::GetLabel(), StructNode::GetRCurr(), StructDispNode::GetVCurr(), StructNode::GetWCurr(), StructNode::GetWPCurr(), StructDispNode::GetXCurr(), StructDispNode::GetXPPCurr(), pEMF, Modal::pGetModalNode(), pModal, pNode, q, qP, R, ExtModalForceBase::Send(), Mat3x3::Transpose(), uFlags, Modal::uGetNModes(), and Zero3.

 {
         Vec3 x;
         Mat3x3 R;
         Vec3 v;
         Vec3 w;
 
         if (uFlags & ExtModalForceBase::EMF_RIGID) {
                 x = pNode->GetXCurr();
                 R = pNode->GetRCurr();
                 v = pNode->GetVCurr();
                 w = pNode->GetWCurr();
 
         } else {
                 x = Zero3;
                 R = Eye3;
                 v = Zero3;
                 w = Zero3;
         }
 
         // Temporary?
         if (uFlags & ExtModalForceBase::EMF_MODAL) {
                 const VecN& a = pModal->GetA();
                 const VecN& b = pModal->GetB();
                 for (unsigned i = 0; i < q.size(); i++) {
                         q[i] = a(i + 1);
                         qP[i] = b(i + 1);
                 }
         }
 
         pEMF->Send(pEFH, uFlags, GetLabel(), x, R, v, w, q, qP);
 
 #if 0
         if (uFlags & ExtModalForceBase::EMF_RIGID) {
                 const StructNode *pNode = pModal->pGetModalNode();
 
                 Mat3x3 RT = pNode->GetRCurr().Transpose();
                 outf << pNode->GetLabel()
                         << " " << pNode->GetXCurr()
                         << " " << pNode->GetRCurr()
                         << " " << RT*pNode->GetVCurr()
                         << " " << RT*pNode->GetWCurr();
                 if (bOutputAccelerations) {
                         outf
                                 << " " << pNode->GetXPPCurr()
                                 << " " << pNode->GetWPCurr();
                 }
                 outf << std::endl;
         }
 
         if (uFlags & ExtModalForceBase::EMF_MODAL) {
                 const VecN& a = pModal->GetA();
                 const VecN& b = pModal->GetB();
                 const VecN& bP = pModal->GetBP();
                 for (integer i = 1; i < pModal->uGetNModes(); i++) {
                         outf << a(i) << " " << b(i);
                         if (bOutputAccelerations) {
                                 outf << " " << bP(i);
                         }
                         outf << std::endl;
                 }
         }
 #endif
 }

Here is the call graph for this function:

void ModalExt::WorkSpaceDim	(	integer *	piNumRows,
		integer *	piNumCols
	)		const

virtual

Implements Elem.

Definition at line 494 of file modalext.cc.

References ExtForce::COUPLING_NONE, ExtForce::iCoupling, pModal, pNode, and Modal::uGetNModes().

Referenced by AssRes().

 {
         if (iCoupling == COUPLING_NONE) {
                 *piNumRows = 0;
                 *piNumCols = 0;
 
         } else {
                 *piNumRows = (pNode ? 6 : 0) + (pModal ? pModal->uGetNModes() : 0);
                 *piNumCols = 1;
         }
 }