gravity.h File Reference

#include "elem.h"
#include "tpldrive.h"

Include dependency graph for gravity.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	Gravity
class	UniformGravity
class	CentralGravity
class	GravityOwner
class	ElemGravityOwner

Functions
Elem *	ReadGravity (DataManager *pDM, MBDynParser &HP)

Function Documentation

Elem* ReadGravity	(	DataManager *	pDM,
		MBDynParser &	HP
	)

Definition at line 237 of file gravity.cc.

References AbsRefFrame, DataManager::fReadOutput(), IncludeParser::GetLineData(), MBDynParser::GetPosAbs(), HighParser::GetReal(), Elem::GRAVITY, HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, ReadDC3D(), SAFENEWWITHCONSTRUCTOR, and Zero3.

Referenced by DataManager::ReadElems().

{
        Elem *pE = 0;

        if (HP.IsKeyWord("central")) {
                Vec3 X0 = ::Zero3;
                if (HP.IsKeyWord("origin")) {
                        X0 = HP.GetPosAbs(::AbsRefFrame);
                }

                if (!HP.IsKeyWord("mass")) {
                        silent_cerr("Gravity: \"mass\" keyword expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dM = HP.GetReal();
                if (dM <= std::numeric_limits<doublereal>::epsilon()) {
                        silent_cerr("Gravity: mass " << dM << " is negative or too small at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                if (!HP.IsKeyWord("G")) {
                        silent_cerr("Gravity: \"G\" keyword expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                // 6.673 84 x 10-11 m3 kg-1 s-2, http://physics.nist.gov/cuu/Constants/index.html
                doublereal dG = 6.67384e-11;
                if (!HP.IsKeyWord("si")) {
                        dG = HP.GetReal();
                        if (dG <= std::numeric_limits<doublereal>::epsilon()) {
                                silent_cerr("Gravity: gravity constant " << dG << " is negative or too small at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                }

                SAFENEWWITHCONSTRUCTOR(pE, CentralGravity, CentralGravity(X0, dM, dG, 0));

        } else {
                if (!HP.IsKeyWord("uniform")) {
                        silent_cerr("Gravity: "
                                "<type> expected "
                                "at line " << HP.GetLineData() << "; "
                                "assuming \"uniform\""
                                << std::endl);
                }

                TplDriveCaller<Vec3>* pDC = ReadDC3D(pDM, HP);

                flag fOut = pDM->fReadOutput(HP, Elem::GRAVITY);

                SAFENEWWITHCONSTRUCTOR(pE, UniformGravity, UniformGravity(pDC, fOut));
        }

        return pE;
}

Here is the call graph for this function: