StrutContent.py

import os
import modena
from modena import CFunction, ForwardMappingModel


f = CFunction(
    Ccode='''
#include "modena.h"
#include "math.h"
#define MAX(a,b) ((a) > (b) ? a : b)
void strutContent
(
    const modena_model_t* model,
    const double* inputs,
    double *outputs
)
{
    {% block variables %}{% endblock %}

    const double a = parameters[0];
    const double b = parameters[1];
    const double c = parameters[2];
    const double pi = 3.14159;

    outputs[0] = MAX(atan(a*rho + b)/pi*2,0.0);
}
''',
    # These are global bounds for the function
    inputs={
        'rho': { 'min': 0, 'max': 1e5},
    },
    outputs={
        'fs': { 'min': 0, 'max': 1, 'argPos': 0 },
    },
    parameters={
        'param0': { 'min': -9e99, 'max': +9e99, 'argPos': 0 },
        'param1': { 'min': -9e99, 'max': +9e99, 'argPos': 1 },
        'param2': { 'min': -9e99, 'max': +9e99, 'argPos': 2 },
    },
)


m = ForwardMappingModel(
    _id='strutContent',
    surrogateFunction=f,
    substituteModels= [ ],
    parameters=[0.06115509, -0.72513392,  1.],
)