Function.cc

Go to the documentation of this file.

#include <utl/Function.h>

#include <utl/ExpressionParser.h>
#include <utl/AugerUnits.h>
#include <utl/SymbolTable.h>

#include <iostream>
using namespace std;


namespace utl {
  
  Function::Function(const std::string& function,
                     const std::vector<std::string>& vars)
    : fFunction(function) {
    SymbolTable fVariables;
    for (const auto &v : vars) {
      fVariableNames.push_back(v);
      fVariables[v] = 0;
    }    
    fFunc = FunctionParser(fVariables);
  }

  Function::Function(const std::string& function,
                     const std::string& var)
    : fFunction(function) { 
    fVariableNames.push_back(var);
    SymbolTable fVariables;
    fVariables[var] = 0;
    fFunc = FunctionParser(fVariables);
  }

  double
  Function::operator()(const double& var)
    const
  {
    if (fVariableNames.size() != 1) {
      ostringstream err;
      err << "Number of input variables " << 1 << " is different from number of function parameters " << fVariableNames.size() << " for function \'" << fFunction << "\'";
      ERROR(err);
      throw XMLParseException(err.str());
    }
    fFunc.SetVariable(fVariableNames[0], var);
    return fFunc.Evaluate(fFunction);
  }
  
  double
  Function::operator()(const std::vector<double>& vars)
    const
  {
    if (vars.size() != fVariableNames.size()) {
      ostringstream err;
      err << "Number of input variables " << vars.size() << " is different from number of function parameters " << fVariableNames.size() << " for function \'" << fFunction << "\'";
      ERROR(err);
      throw XMLParseException(err.str());
    }
    unsigned int i = 0;
    for (const auto &v : vars) {
      fFunc.SetVariable(fVariableNames[i], v);
      ++i;
    }
    return fFunc.Evaluate(fFunction);
  }
  
}