POD.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 "ParallelSolutionStorage.h"
#include "MooseTypes.h"
#include "libmesh/petsc_vector.h"
#include "libmesh/libmesh_config.h"

#include <slepcsvd.h>

namespace StochasticTools
{
class POD
{
public:
  POD(const ParallelSolutionStorage * const parallel_storage,
      const std::string & extra_slepc_options,
      const Parallel::Communicator & comm);
  void computePOD(const VariableName & vname,
                  std::vector<DenseVector<Real>> & left_basis_functions,
                  std::vector<DenseVector<Real>> & right_basis_functions,
                  std::vector<Real> & singular_values,
                  const dof_id_type num_modes,
                  const Real energy) const;

private:
  dof_id_type determineNumberOfModes(const std::vector<Real> & singular_values,
                                     const dof_id_type num_modes_compute,
                                     const Real energy) const;

#if !PETSC_VERSION_LESS_THAN(3, 14, 0)
  const ParallelSolutionStorage * const _parallel_storage;
  const std::string & _extra_slepc_options;
  const Parallel::Communicator & _communicator;
#endif
};
}