ScalarFunctionsImpl.cc File Reference

#include "mbconfig.h"
#include <cmath>
#include <typeinfo>
#include "myassert.h"
#include "ScalarFunctionsImpl.h"
#include "interp.h"
#include "mbpar.h"
#include "dataman.h"

Include dependency graph for ScalarFunctionsImpl.cc:

Go to the source code of this file.

Classes
struct	ptr_cast_helper< T1, T2 >
struct	ptr_cast_helper< T1 *, T2 >
struct	ptr_cast_helper< T1 *const, T2 >
struct	ptr_cast_helper< const T1 *, T2 >
struct	ptr_cast_helper< const T1 *const, T2 >
struct	ConstSFR
struct	LinearSFR
struct	PowSFR
struct	LogSFR
struct	ExpSFR
struct	CubicSplineSFR
struct	MultiLinearSFR
struct	ChebychevSFR
struct	SumSFR
struct	SubSFR
struct	MulSFR
struct	DivSFR
struct	SFWordSetType
class	ScalarFunctionDriveCaller
struct	ScalarFunctionDCR
class	ScalarFunctionIsotropicCL< T, Tder >
class	ScalarFunctionIsotropicCL< doublereal, doublereal >
struct	ScalarFunctionIsotropicCLR< T, Tder >
class	ScalarFunctionOrthotropicCL< T, Tder >
struct	ScalarFunctionOrthotropicCLR< T, Tder >

Typedefs
typedef std::map< std::string, const ScalarFunctionRead *, ltstrcase >	SFReadType

Functions
template<class T1 , class T2 >
T1	ptr_cast (T2 &arg)
const BasicScalarFunction *const	ParseScalarFunction (MBDynParser &HP, DataManager *const pDM)
bool	SetSF (const std::string &s, const ScalarFunctionRead *rf)
void	InitSF (void)
void	DestroySF (void)

Variables
static SFReadType	SFRead
static SFWordSetType	SFWordSet
static unsigned	done = 0

Typedef Documentation

typedef std::map<std::string, const ScalarFunctionRead *, ltstrcase> SFReadType

Definition at line 1000 of file ScalarFunctionsImpl.cc.

Function Documentation

void DestroySF

(

void

)

Definition at line 1527 of file ScalarFunctionsImpl.cc.

References done, MBDYN_EXCEPT_ARGS, and SFRead.

Referenced by MBDynParser::~MBDynParser().

{
        if (::done == 0) {
                silent_cerr("DestroySF() called once too many" << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        if (--done > 0) {
                return;
        }

        for (SFReadType::iterator i = SFRead.begin(); i != SFRead.end(); ++i) {
                delete i->second;
        }
        SFRead.clear();
}

void InitSF

(

void

)

Definition at line 1409 of file ScalarFunctionsImpl.cc.

References ASSERT, done, MBDYN_EXCEPT_ARGS, SetCL1D(), SetCL3D(), SetCL6D(), SetDriveCallerData(), and SetSF().

Referenced by MBDynParser::MBDynParser().

{
        if (::done++ > 0) {
                return;
        }

        bool b;

        b = SetSF("const", new ConstSFR);
        ASSERT(b);
        b = SetSF("linear", new LinearSFR);
        ASSERT(b);
        b = SetSF("pow", new PowSFR);
        ASSERT(b);
        b = SetSF("log", new LogSFR);
        ASSERT(b);
        b = SetSF("exp", new ExpSFR);
        ASSERT(b);
        b = SetSF("sum", new SumSFR);
        ASSERT(b);
        b = SetSF("sub", new SubSFR);
        ASSERT(b);
        b = SetSF("mul", new MulSFR);
        ASSERT(b);
        b = SetSF("div", new DivSFR);
        ASSERT(b);
        b = SetSF("cubic" "spline", new CubicSplineSFR);
        ASSERT(b);
        b = SetSF("multilinear", new MultiLinearSFR);
        ASSERT(b);
        b = SetSF("chebychev", new ChebychevSFR);
        ASSERT(b);

        /* this is about initializing the scalar function drive */
        ScalarFunctionDCR *rf = new ScalarFunctionDCR;
        if (!SetDriveCallerData("scalar" "function", rf)) {
                delete rf;

                silent_cerr("unable to register scalar function drive caller"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        /* this is about initializing the scalar function constitutive law(s) */
        ConstitutiveLawRead<doublereal, doublereal> *rf1D
                = new ScalarFunctionIsotropicCLR<doublereal, doublereal>;
        if (!SetCL1D("scalar" "function" "elastic" "isotropic", rf1D)) {
                delete rf1D;

                silent_cerr("unable to register scalar function isotropic 1D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        rf1D = new ScalarFunctionIsotropicCLR<doublereal, doublereal>;
        if (!SetCL1D("scalar" "function" "elastic" "orthotropic", rf1D)) {
                delete rf1D;

                silent_cerr("unable to register scalar function orthotropic 1D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        rf1D = new ScalarFunctionIsotropicCLR<doublereal, doublereal>;
        if (!SetCL1D("scalar" "function" "elastic", rf1D)) {
                delete rf1D;

                silent_cerr("unable to register scalar function 1D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        ConstitutiveLawRead<Vec3, Mat3x3> *rf3D = new ScalarFunctionIsotropicCLR<Vec3, Mat3x3>;
        if (!SetCL3D("scalar" "function" "elastic" "isotropic", rf3D)) {
                delete rf3D;

                silent_cerr("unable to register scalar function isotropic 3D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        rf3D = new ScalarFunctionOrthotropicCLR<Vec3, Mat3x3>;
        if (!SetCL3D("scalar" "function" "elastic" "orthotropic", rf3D)) {
                delete rf3D;

                silent_cerr("unable to register scalar function orthotropic 3D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        ConstitutiveLawRead<Vec6, Mat6x6> *rf6D = new ScalarFunctionIsotropicCLR<Vec6, Mat6x6>;
        if (!SetCL6D("scalar" "function" "elastic" "isotropic", rf6D)) {
                delete rf6D;

                silent_cerr("unable to register scalar function isotropic 6D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        rf6D = new ScalarFunctionOrthotropicCLR<Vec6, Mat6x6>;
        if (!SetCL6D("scalar" "function" "elastic" "orthotropic", rf6D)) {
                delete rf6D;

                silent_cerr("unable to register scalar function orthotropic 6D constitutive law"
                        << std::endl);

                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
}

Here is the call graph for this function:

const BasicScalarFunction* const ParseScalarFunction	(	MBDynParser &	HP,
		DataManager *const	pDM
	)

Definition at line 1011 of file ScalarFunctionsImpl.cc.

References func(), IncludeParser::GetLineData(), MBDynParser::GetScalarFunction(), HighParser::GetStringWithDelims(), HighParser::IsWord(), MBDYN_EXCEPT_ARGS, MBDynParser::SetScalarFunction(), and SFRead.

Referenced by MBDynParser::GetScalarFunction(), ParseFriction(), ScalarFuncGR::Read(), NLSFViscoElasticCLR< T, Tder, Typ >::Read(), NLPViscoElasticCLR< T, Tder, Typ >::Read(), SumSFR::Read(), SubSFR::Read(), MulSFR::Read(), DivSFR::Read(), ScalarFunctionDCR::Read(), ScalarFunctionIsotropicCLR< T, Tder >::Read(), ScalarFunctionOrthotropicCLR< T, Tder >::Read(), and MBDynParser::ScalarFunction_int().

{
        std::string func_name(HP.GetStringWithDelims());
        
        const BasicScalarFunction *sf = HP.GetScalarFunction(func_name);
        if (sf == 0) {
                const char *s = HP.IsWord(SFWordSet);
                if (s == 0) {
                        s = "const";
                }

                SFReadType::iterator func = SFRead.find(std::string(s));
                if (func == SFRead.end()) {
                        silent_cerr("unknown scalar function type \"" << s << "\" "
                                "for function \"" << func_name << "\" "
                                "at line " << HP.GetLineData() << std::endl);
                        throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                try {
                        sf = func->second->Read(pDM, HP);
                } catch (...) {
                        silent_cerr("Unable to parse "
                                "ScalarFunction(\"" << func_name << "\") "
                                "at line " << HP.GetLineData() << std::endl);
                        throw;
                }
                if (!HP.SetScalarFunction(func_name, sf)) {
                        silent_cerr("scalar function \"" << func_name << "\" "
                                "already defined at line " << HP.GetLineData() << std::endl);
                        throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

        } else if (HP.IsWord(SFWordSet)) {
                silent_cerr("Error: redefinition of "
                        "\"" << func_name << "\" scalar function "
                        "at line " << HP.GetLineData() << std::endl);
                throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
        
        return sf;
}

Here is the call graph for this function:

template<class T1 , class T2 >

T1 ptr_cast

(

T2 &

arg

)

Definition at line 126 of file ScalarFunctionsImpl.cc.

References ptr_cast_helper< T1, T2 >::cast().

Referenced by ScalarFunctionDriveCaller::dGetP().

                     {
        return ptr_cast_helper<T1,T2>::cast(arg);
}

Here is the call graph for this function:

bool SetSF	(	const std::string &	s,
		const ScalarFunctionRead *	rf
	)

Definition at line 1399 of file ScalarFunctionsImpl.cc.

References SFRead.

Referenced by InitSF(), mbdyn_octave_set(), and module_init().

{
        pedantic_cout("registering scalar function \"" << s << "\""
                << std::endl );
        return SFRead.insert(SFReadType::value_type(s, rf)).second;
}

Variable Documentation

unsigned done = 0

static

Definition at line 1406 of file ScalarFunctionsImpl.cc.

Referenced by DestroySF(), and InitSF().

SFReadType SFRead

static

Definition at line 1001 of file ScalarFunctionsImpl.cc.

Referenced by DestroySF(), SFWordSetType::IsWord(), ParseScalarFunction(), and SetSF().

SFWordSetType SFWordSet

static

Definition at line 1008 of file ScalarFunctionsImpl.cc.