NearestReporterCoordinatesFunction.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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

#pragma once

#include "OptimizationFunction.h"
#include "ReporterInterface.h"
class NearestReporterCoordinatesFunction : public OptimizationFunction, public ReporterInterface
{
public:
  static InputParameters validParams();

  NearestReporterCoordinatesFunction(const InputParameters & parameters);

  using Function::value;
  virtual Real value(Real t, const Point & p) const override;
  virtual RealGradient gradient(Real t, const Point & p) const override;
  virtual Real timeDerivative(Real t, const Point & p) const override;
  virtual std::vector<Real> parameterGradient(Real t, const Point & p) const override;

protected:
  void buildCoordinateMapping() const;

  std::array<std::pair<Real, std::size_t>, 2> findNearestPoint(Real t, const Point & p) const;

  const std::vector<Real> & _coordx;
  const std::vector<Real> & _coordy;
  const std::vector<Real> & _coordz;
  const std::vector<Real> & _coordt;
  const std::vector<Real> & _values;

  mutable std::size_t _nval;
  mutable std::map<Point, std::vector<std::pair<Real, std::size_t>>> _coord_mapping;

private:
  const std::vector<Real> _empty_vec = {};
};