SegregatedSolverUtils.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

// Libmesh includes
#include "MooseTypes.h"
#include "MooseMesh.h"
#include "libmesh/petsc_vector.h"
#include "libmesh/petsc_matrix.h"

namespace NS
{
namespace FV
{
void relaxMatrix(SparseMatrix<Number> & matrix_in,
                 const Real relaxation_parameter,
                 NumericVector<Number> & diff_diagonal);

void relaxRightHandSide(NumericVector<Number> & rhs_in,
                        const NumericVector<Number> & solution_in,
                        const NumericVector<Number> & diff_diagonal);

void relaxSolutionUpdate(NumericVector<Number> & vec_new,
                         const NumericVector<Number> & vec_old,
                         const Real relaxation_factor);

void limitSolutionUpdate(NumericVector<Number> & solution,
                         const Real min_limit = std::numeric_limits<Real>::epsilon(),
                         const Real max_limit = 1e10);

Real computeNormalizationFactor(const NumericVector<Number> & solution,
                                const SparseMatrix<Number> & mat,
                                const NumericVector<Number> & rhs);

void constrainSystem(SparseMatrix<Number> & mx,
                     NumericVector<Number> & rhs,
                     const Real desired_value,
                     const dof_id_type dof_id);

dof_id_type findPointDoFID(const MooseVariableFieldBase & variable,
                           const MooseMesh & mesh,
                           const Point & point);

bool converged(const std::vector<std::pair<unsigned int, Real>> & residuals,
               const std::vector<Real> & abs_tolerances);
} // End FV namespace
} // End Moose namespace