PODReducedBasisSurrogate.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 "SurrogateModel.h"

class PODReducedBasisSurrogate : public SurrogateModel
{
public:
  static InputParameters validParams();

  PODReducedBasisSurrogate(const InputParameters & parameters);

  void evaluateSolution(const std::vector<Real> & params);

  void evaluateSolution(const std::vector<Real> & params,
                        DenseVector<Real> & inp_vector,
                        std::string var_name);

  const std::vector<DenseVector<Real>> & getApproximateSolution() const { return _approx_solution; }

  Real getNodalQoI(std::string var_name, unsigned int qoi_type) const;

protected:
  void initializeReducedSystem();

  void initializeApproximateSolution();

  void solveReducedSystem(const std::vector<Real> & params);

  void reconstructApproximateSolution();

  void reconstructApproximateSolution(DenseVector<Real> & inp_vector, std::string var_name);

  std::vector<std::string> _change_rank;

  std::vector<unsigned int> _new_ranks;

  std::vector<unsigned int> _final_ranks;

  std::vector<unsigned int> _comulative_ranks;

  const std::vector<std::string> & _var_names;

  const std::vector<std::string> & _tag_types;

  const std::vector<std::vector<DenseVector<Real>>> & _base;

  const std::vector<DenseMatrix<Real>> & _red_operators;

  DenseVector<Real> _coeffs;

  DenseMatrix<Real> _sys_mx;

  DenseVector<Real> _rhs;

  std::vector<DenseVector<Real>> _approx_solution;

  const Real _penalty;

  bool _initialized;
};