LinearSum.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "LinearSum.h"

registerMooseObject("StochasticToolsApp", LinearSum);

InputParameters
LinearSum::validParams()
{
  InputParameters params = LikelihoodFunctionBase::validParams();
  params.addClassDescription(
      "LinearSum function evaluating the model goodness against experiments.");
  return params;
}

LinearSum::LinearSum(const InputParameters & parameters) : Gaussian(parameters) {}

Real
LinearSum::function(const std::vector<Real> & exp,
                    const std::vector<Real> & model,
                    const Real & noise)
{
  mooseAssert(exp.size() == model.size(), "LinearSum: Vectors must be the same size.");
  mooseAssert(noise != 0.0, "LinearSum: Noise must not be zero.");
  Real result = 0.0;
  for (unsigned int i = 0; i < exp.size(); ++i)
  {
    const Real diff = exp[i] - model[i];
    result += diff * diff;
  }
  return result / (noise * noise);
}

Real
LinearSum::function(const std::vector<Real> & x) const
{
  return function(_exp_values, x, _noise);
}