Function Interpolate

kratos/python_scripts/point_output_process.py:231–251

(variable, entity, sf_values, historical_value)

    return header

def Interpolate(variable, entity, sf_values, historical_value):
    if type(entity) == KratosMultiphysics.Node:
        if historical_value:
            return entity.GetSolutionStepValue(variable)
        else:
            return entity.GetValue(variable)
    else: # entity is element or condition
        nodes = entity.GetNodes()
        # Initializing 'value' like this, i don't need to know its type
        # => this way it works both for scalar and array3 variables
        if historical_value:
            value = nodes[0].GetSolutionStepValue(variable) * sf_values[0]
            for n,c in zip(nodes[1:], sf_values[1:]):
                value = value + c * n.GetSolutionStepValue(variable)
        else:
            value = nodes[0].GetValue(variable) * sf_values[0]
            for n,c in zip(nodes[1:], sf_values[1:]):
                value = value + c * n.GetValue(variable)


        return value

def IsArrayVariable(var):
    return type(var) == KratosMultiphysics.Array1DVariable3