PODMapping.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 "VariableMappingBase.h"
#include "ParallelSolutionStorage.h"
#include "UserObjectInterface.h"

#include <slepcsvd.h>
#include "POD.h"

class PODMapping : public VariableMappingBase, public UserObjectInterface
{
public:
  static InputParameters validParams();
  PODMapping(const InputParameters & parameters);

  virtual void buildMapping(const VariableName & vname) override;

  void map(const VariableName & vname,
           const DenseVector<Real> & full_order_vector,
           std::vector<Real> & reduced_order_vector) const override;

  void map(const VariableName & vname,
           const unsigned int global_sample_i,
           std::vector<Real> & reduced_order_vector) const override;

  void inverse_map(const VariableName & vname,
                   const std::vector<Real> & reduced_order_vector,
                   DenseVector<Real> & full_order_vector) const override;

  const std::vector<DenseVector<Real>> & leftBasis(const VariableName & vname);

  const std::vector<DenseVector<Real>> & rightBasis(const VariableName & vname);

  const std::vector<Real> & singularValues(const VariableName & vname);

  const DenseVector<Real> & leftBasisFunction(const VariableName & vname,
                                              const unsigned int base_i);

  const DenseVector<Real> & rightBasisFunction(const VariableName & vname,
                                               const unsigned int base_i);

protected:
  const std::vector<dof_id_type> _num_modes;

  const std::vector<Real> & _energy_threshold;

  std::map<VariableName, std::vector<DenseVector<Real>>> & _left_basis_functions;

  std::map<VariableName, std::vector<DenseVector<Real>>> & _right_basis_functions;

  std::map<VariableName, std::vector<Real>> & _singular_values;

  const std::string & _extra_slepc_options;

private:
  const ParallelSolutionStorage * const _parallel_storage;

  const StochasticTools::POD _pod;
};