modalforce.cc File Reference

#include "mbconfig.h"
#include <fstream>
#include <algorithm>
#include "dataman.h"
#include "modalforce.h"

Include dependency graph for modalforce.cc:

Go to the source code of this file.

Functions
Elem *	ReadModalForce (DataManager *pDM, MBDynParser &HP, unsigned int uLabel)

Function Documentation

Elem* ReadModalForce	(	DataManager *	pDM,
		MBDynParser &	HP,
		unsigned int	uLabel
	)

Definition at line 145 of file modalforce.cc.

References Elem::FORCE, DataManager::fReadOutput(), MBDynParser::GetDriveCaller(), HighParser::GetInt(), WithLabel::GetLabel(), IncludeParser::GetLineData(), Modal::GetModeList(), HighParser::GetReal(), HighParser::IsKeyWord(), Elem::JOINT, MBDYN_EXCEPT_ARGS, Modal::pGetModalNode(), DataManager::ReadElem(), SAFENEWWITHCONSTRUCTOR, and Modal::uGetNModes().

Referenced by ReadForce().

{
        const Modal *pModal = pDM->ReadElem<const Modal, const Joint, Elem::JOINT>(HP);

        std::vector<unsigned int> modeList;
        if (HP.IsKeyWord("list")) {
                const std::vector<unsigned int>& uModeList = pModal->GetModeList();
                std::vector<bool> gotIt(pModal->uGetNModes());
                for (integer i = 0; i < pModal->uGetNModes(); i++) {
                        gotIt[i] = false;
                }

                int iNumModes = HP.GetInt();
                if (iNumModes <= 0 || iNumModes > pModal->uGetNModes()) {
                        silent_cerr("ModalForce(" << uLabel << "): "
                                "illegal mode number " << iNumModes
                                << " at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                modeList.resize(iNumModes);

                for (int i = 0; i < iNumModes; i++) {
                        int iM = HP.GetInt();
                        if (iM <= 0) {
                                silent_cerr("ModalForce(" << uLabel << "): "
                                        "illegal mode " << iM
                                        << " at line " << HP.GetLineData()
                                        << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        std::vector<unsigned int>::const_iterator
                                iv = std::find(uModeList.begin(), uModeList.end(), (unsigned int)iM);
                        if (iv == uModeList.end()) {
                                silent_cerr("ModalForce(" << uLabel << "): "
                                        "mode " << iM << " not active "
                                        "in Modal(" << pModal->GetLabel() << ")"
                                        "at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                        int iModeIdx = iv - uModeList.begin();

                        if (gotIt[iModeIdx]) {
                                silent_cerr("ModalForce(" << uLabel << "): "
                                        "mode " << iModeIdx << " already set "
                                        "at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }

                        modeList[i] = iModeIdx + 1;
                        gotIt[iModeIdx] = true;
                }

        } else {
                modeList.resize(pModal->uGetNModes());
                for (integer i = 0; i < pModal->uGetNModes(); i++) {
                        modeList[i] = i + 1;
                }
        }

        std::vector<DriveCaller *> f(modeList.size());
        Mat3xN *Mt = 0;
        Mat3xN *Mr = 0;
        const StructNode *pNode = pModal->pGetModalNode();
        if (pNode) {
                Mt = new Mat3xN(modeList.size(), 0.);
                Mr = new Mat3xN(modeList.size(), 0.);
        }
        for (unsigned i = 0; i < f.size(); i++) {
                f[i] = HP.GetDriveCaller(false);
                if (f[i] == 0) {
                        silent_cerr("ModalForce(" << uLabel << "): "
                                "unable to read DriveCaller for mode #" << i + 1
                                << " (mode number " << modeList[i] << ") "
                                "at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                if (pNode && HP.IsKeyWord("resultant")) {
                        for (unsigned r = 1; r <= 3; r++) {
                                (*Mt)(r, i + 1) = HP.GetReal();
                        }

                        for (unsigned r = 1; r <= 3; r++) {
                                (*Mr)(r, i + 1) = HP.GetReal();
                        }
                }
        }

        flag fOut = pDM->fReadOutput(HP, Elem::FORCE);
        Elem *pEl = 0;
        SAFENEWWITHCONSTRUCTOR(pEl, ModalForce,
                ModalForce(uLabel, pModal, modeList, f, Mt, Mr, fOut));

        return pEl;
}

Here is the call graph for this function: