StaticCondensationFieldSplitPreconditioner.C

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
 
#include "StaticCondensationFieldSplitPreconditioner.h"
 
// MOOSE includes
#include "PetscSupport.h"
#include "NonlinearSystemBase.h"
#include "PetscDMMoose.h"
#include "Split.h"
 
#include "libmesh/static_condensation.h"
#include "libmesh/static_condensation_dof_map.h"
#include "libmesh/petsc_linear_solver.h"
 
registerMooseObjectAliased("MooseApp", StaticCondensationFieldSplitPreconditioner, "SCFSP");
 
InputParameters
StaticCondensationFieldSplitPreconditioner::validParams()
{
  return FieldSplitPreconditionerTempl<MooseStaticCondensationPreconditioner>::validParams();
}
 
StaticCondensationFieldSplitPreconditioner::StaticCondensationFieldSplitPreconditioner(
    const InputParameters & params)
  : FieldSplitPreconditionerTempl<MooseStaticCondensationPreconditioner>(params)
{
  std::shared_ptr<Split> top_split = _nl.getSplit(_decomposition_split);
  top_split->setup(_nl, prefix());
}
 
const libMesh::DofMapBase &
StaticCondensationFieldSplitPreconditioner::dofMap() const
{
  return scDofMap();
}
 
const libMesh::System &
StaticCondensationFieldSplitPreconditioner::system() const
{
  return scDofMap().reduced_system();
}
 
std::string
StaticCondensationFieldSplitPreconditioner::prefix() const
{
  return MooseStaticCondensationPreconditioner::prefix();
}
 
void
StaticCondensationFieldSplitPreconditioner::setupDM()
{
  PetscBool ismoose;
  DM dm = LIBMESH_PETSC_NULLPTR;
 
  // Initialize the part of the DM package that's packaged with Moose; in the PETSc source tree this
  // call would be in DMInitializePackage()
  LibmeshPetscCall(DMMooseRegisterAll());
  // Create and set up the DM that will consume the split options and deal with block matrices.
  auto & petsc_solver =
      cast_ref<libMesh::PetscLinearSolver<Number> &>(scSysMat().reduced_system_solver());
  auto ksp = petsc_solver.ksp();
  // if there exists a DMMoose object, then we do not need to recreate one
  LibmeshPetscCall(KSPGetDM(ksp, &dm));
  if (dm)
  {
    LibmeshPetscCall(PetscObjectTypeCompare((PetscObject)dm, DMMOOSE, &ismoose));
    if (ismoose)
      return;
  }
  createMooseDM(&dm);
  LibmeshPetscCall(KSPSetDM(ksp, dm));
// We set the operators ourselves. We do not want the DM to generate the operators
#if PETSC_VERSION_LESS_THAN(3, 25, 0)
  LibmeshPetscCall(KSPSetDMActive(ksp, PETSC_FALSE));
#else
  LibmeshPetscCall(KSPSetDMActive(ksp, KSP_DMACTIVE_OPERATOR, PETSC_FALSE));
#endif
  LibmeshPetscCall(DMDestroy(&dm));
}
 
KSP
StaticCondensationFieldSplitPreconditioner::getKSP()
{
  auto & petsc_solver =
      cast_ref<libMesh::PetscLinearSolver<Number> &>(scSysMat().reduced_system_solver());
  return petsc_solver.ksp();
}