NavierStokesProblem.C File Reference

Go to the source code of this file.

Functions
	registerMooseObject ("NavierStokesApp", NavierStokesProblem)
	_commute_lsc (getParam< bool >("commute_lsc"))
	_mass_matrix (getParam< TagName >("mass_matrix"))
	_L_matrix (getParam< TagName >("L_matrix"))
	_have_mass_matrix (!_mass_matrix.empty())
	_have_L_matrix (!_L_matrix.empty())
	_pressure_mass_matrix_as_pre (getParam< bool >("use_pressure_mass_matrix"))
	_schur_fs_index (getParam< std::vector< unsigned int >>("schur_fs_index"))
PetscErrorCode	navierStokesKSPPreSolve (KSP root_ksp, Vec, Vec, void *context)

Function Documentation

_commute_lsc()

_commute_lsc

(

getParam< bool >

"commute_lsc"

)

Referenced by _schur_fs_index().

_have_L_matrix()

_have_L_matrix

(

!_L_matrix.

empty()

)

Referenced by _schur_fs_index().

_have_mass_matrix()

_have_mass_matrix

(

!_mass_matrix.

empty()

)

Referenced by _schur_fs_index().

_L_matrix()

_L_matrix

(

getParam< TagName >

"L_matrix"

)

_mass_matrix()

_mass_matrix

(

getParam< TagName >

"mass_matrix"

)

_pressure_mass_matrix_as_pre()

_pressure_mass_matrix_as_pre

(

getParam< bool >

"use_pressure_mass_matrix"

)

Referenced by _schur_fs_index().

_schur_fs_index()

_schur_fs_index

(

getParam< std::vector< unsigned int >>

"schur_fs_index"

)

Definition at line 62 of file NavierStokesProblem.C.

{
  if (_commute_lsc)
  {
    if (!_have_mass_matrix)
      paramError("mass_matrix",
                 "A pressure mass matrix must be provided if we are commuting the LSC commutator.");
    if (!_have_L_matrix)
      paramError("L_matrix",
                 "A matrix corresponding to the viscous component of the momentum equation must be "
                 "provided if we are commuting the LSC commutator.");
  }
  else if (_have_L_matrix)
    paramError("L_matrix",
               "If not commuting the LSC commutator, then the 'L_matrix' should not be provided "
               "because it will not be used. For Elman LSC preconditioning, L will be computed "
               "automatically using system matrix data (e.g. the off-diagonal blocks in the "
               "velocity-pressure system).");

  if (_pressure_mass_matrix_as_pre && !_have_mass_matrix)
    paramError("mass_matrix",
               "If 'use_pressure_mass_matrix', then a pressure 'mass_matrix' must be provided");
}

navierStokesKSPPreSolve()

PetscErrorCode navierStokesKSPPreSolve	(	KSP	root_ksp,
		Vec	,
		Vec	,
		void *	context
	)

Definition at line 364 of file NavierStokesProblem.C.

Referenced by NavierStokesProblem::initPetscOutputAndSomeSolverSettings().

{
  PetscFunctionBegin;

  auto * ns_problem = static_cast<NavierStokesProblem *>(context);
  ns_problem->clearIndexSets();
  auto schur_ksp = ns_problem->findSchurKSP(root_ksp, 0);
  ns_problem->setupLSCMatrices(schur_ksp);

  PetscFunctionReturn(PETSC_SUCCESS);
}

registerMooseObject()

registerMooseObject	(	"NavierStokesApp"	,
		NavierStokesProblem
	)