FunctionEvaluation.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <map>
#include <memory>
#include <span>
#include <string>
#include <vector>
 
#include "ExprtkUtils.h"
#include "SpatialPosition.h"
 
namespace ParameterLib
{

class FunctionEvaluation
{
public:
    using Expression = exprtk::expression<double>;
    using SymbolTable = exprtk::symbol_table<double>;

    FunctionEvaluation(
        int num_threads,
        std::vector<std::string> const& variables,
        std::vector<std::string> const& expression_strings,
        std::map<std::string,
                 std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
            curves);

    std::vector<double> evaluate(SpatialPosition const& pos, double t) const;

    void evaluate(SpatialPosition const& pos, double t,
                  std::span<double> result) const;

    int getNumberOfComponents() const;
 
    bool isTimeDependent() const { return true; }
 
private:
    struct PerThreadData
    {
        PerThreadData(
            std::vector<std::string> const& variables,
            bool spatial_position_is_required,
            std::vector<std::string> const& used_curve_names,
            std::map<
                std::string,
                std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
                curves,
            std::vector<std::string> const& expression_strings);
 
        PerThreadData(PerThreadData&& other);
        PerThreadData& operator=(PerThreadData&& other) noexcept = delete;
        PerThreadData(PerThreadData const&) = delete;
        PerThreadData& operator=(PerThreadData const&) = delete;

        std::map<std::string, CurveWrapper> curve_wrappers;
        SymbolTable symbol_table;
        std::vector<Expression> value_expressions;
        SymbolTableCache symbol_table_cache;
    };

    std::vector<PerThreadData> per_thread_data_;
};
 
}  // namespace ParameterLib