module-udunits.cc File Reference

#include "mbconfig.h"
#include "mathp.h"
#include "dataman.h"
#include <udunits.h>

Include dependency graph for module-udunits.cc:

Go to the source code of this file.

Classes
class	UDUnitsNameSpace

Functions
static int	unit_convert (const MathParser::MathArgs &args)
int	module_init (const char *module_name, void *pdm, void *php)
	This function registers our user defined element for the math parser. More...

Function Documentation

int module_init	(	const char *	module_name,
		void *	pdm,
		void *	php
	)

This function registers our user defined element for the math parser.

It is called when the "module load" statement appears in the input file.

Definition at line 226 of file module-udunits.cc.

References IncludeParser::GetFileName(), HighParser::GetMathParser(), HighParser::IsArg(), path, pHP, and MathParser::RegisterNameSpace().

{
        MBDynParser     *pHP = (MBDynParser *)php;

        const char *path = 0;
        if (pHP->IsArg()) {
                path = pHP->GetFileName();
        }

        /* registers unit conversion namespace */
        MathParser::NameSpace *pNS = new UDUnitsNameSpace(path);
        int rc = pHP->GetMathParser().RegisterNameSpace(pNS);
        if (rc != 0) {
                delete pNS;
        }
        return rc;
}

Here is the call graph for this function:

static int unit_convert ( const MathParser::MathArgs &  args )

static

Definition at line 62 of file module-udunits.cc.

References MathParser::AF_OPTIONAL_NON_PRESENT, ASSERT, MathParser::AT_REAL, MathParser::AT_STRING, MathParser::MathArg_t::IsFlag(), and MBDYN_EXCEPT_ARGS.

Referenced by UDUnitsNameSpace::UDUnitsNameSpace().

{
        ASSERT(args.size() == 1 + 2 + 1);
        ASSERT(args[0]->Type() == MathParser::AT_REAL);
        ASSERT(args[1]->Type() == MathParser::AT_STRING);
        ASSERT(args[2]->Type() == MathParser::AT_STRING);
        ASSERT(args[3]->Type() == MathParser::AT_REAL);

        MathParser::MathArgReal_t *out = dynamic_cast<MathParser::MathArgReal_t *>(args[0]);
        ASSERT(out != 0);

        MathParser::MathArgString_t *arg1 = dynamic_cast<MathParser::MathArgString_t *>(args[1]);
        ASSERT(arg1 != 0);

        MathParser::MathArgString_t *arg2 = dynamic_cast<MathParser::MathArgString_t *>(args[2]);
        ASSERT(arg2 != 0);

        MathParser::MathArgReal_t *arg3 = dynamic_cast<MathParser::MathArgReal_t *>(args[3]);
        ASSERT(arg3 != 0);

        int rc;

        utUnit  u_from;
        rc = utScan((*arg1)().c_str(), &u_from);
        switch (rc) {
        case 0:
                break;

        default:
                /* TODO: handle specific errors */
                silent_cerr("module-udunits: utScan could not interpret "
                        "unit \"" << (*arg1)() << "\" in first arg"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        utUnit  u_to;
        rc = utScan((*arg2)().c_str(), &u_to);
        switch (rc) {
        case 0:
                break;

        default:
                silent_cerr("module-udunits: utScan could not interpret "
                        "unit \"" << (*arg2)() << "\" in second arg"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        double slope, intercept;
        rc = utConvert(&u_from, &u_to, &slope, &intercept);
        switch (rc) {
        case 0:
                break;

        default:
                silent_cerr("module-udunits: utConvert failed (rc=" << rc << ")" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        if (arg3->IsFlag(MathParser::AF_OPTIONAL_NON_PRESENT)) {
                if (intercept != 0) {
                        silent_cerr("module-udunits: conversion between \"" << (*arg1)() << "\" "
                                "and \"" << (*arg2)() << "\" "
                                "has non-zero intercept" << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                *out = slope;

        } else {
                *out = (*arg3)()*slope + intercept;
        }

        return 0;
}

Here is the call graph for this function: