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;
};