polymerConductivity.py

import os
import modena
from modena import ForwardMappingModel, BackwardMappingModel, SurrogateModel, CFunction
import modena.Strategy as Strategy
from fireworks.user_objects.firetasks.script_task import FireTaskBase, ScriptTask
from fireworks import Firework, Workflow, FWAction
from fireworks.utilities.fw_utilities import explicit_serialize
from blessings import Terminal
from jinja2 import Template


term = Terminal()


f_polymer_thermal_conductivity = CFunction(
    Ccode='''
#include "modena.h"
#include "math.h"

void thermal_conductivity
(
    const modena_model_t* model,
    const double* inputs,
    double *outputs
)
{
    {% block variables %}{% endblock %}

    const double a = parameters[0];
    const double b = parameters[1];

    outputs[0] = (a*T)+b;
}
''',
    # These are global bounds for the function
    inputs={
        'T': {'min': 273, 'max': 550},
    },
    outputs={
        'polymer_thermal_conductivity': {
            'min': 0, 'max': +9e99, 'argPos': 0
        },
    },
    parameters={
        'param0': {'min': -9e99, 'max': +9e99, 'argPos': 0},
        'param1': {'min': -9e99, 'max': +9e99, 'argPos': 1},
    },
)


m_polymer_thermal_conductivity = ForwardMappingModel(
    _id='polymer_thermal_conductivity',
    surrogateFunction=f_polymer_thermal_conductivity,
    substituteModels=[],
    parameters=[0.198e-3, 151.08e-3],
    inputs={
        'T': {'min': 273, 'max': 450},
    },
    outputs={
        'polymer_thermal_conductivity': {'min': 0, 'max': +9e99},
    },
)