problem.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#pragma once
 
#include <prismspf/core/field_attributes.h>
#include <prismspf/core/refinement_manager.h>
#include <prismspf/core/simulation_timer.h>
#include <prismspf/core/types.h>
 
#include <prismspf/nucleation/nucleation_manager.h>
#include <prismspf/nucleation/nucleus_refinement_function.h>
 
#include <prismspf/solvers/solve_context.h>
 
#include <prismspf/utilities/integrator.h>
 
#include <prismspf/config.h>
 
PRISMS_PF_BEGIN_NAMESPACE

template <unsigned int dim, unsigned int degree, typename number>
class Problem
{
public:
  Problem(const std::vector<FieldAttributes>         &field_attributes,
          const std::vector<SolveBlock>              &solve_blocks,
          const UserInputParameters<dim>             &_user_inputs,
          PhaseFieldTools<dim>                       &_pf_tools,
          const PDEOperatorBase<dim, degree, number> &_pde_operator);

  void
  solve();
 
private:
  void
  solve_increment(SimulationTimer &sim_timer);

  void
  init_system();

  void
  reinit_system();

  std::vector<FieldAttributes> field_attributes;

  std::vector<SolveBlock> solve_blocks;

  const UserInputParameters<dim> *user_inputs_ptr;

  PhaseFieldTools<dim> *pf_tools;

  TriangulationManager<dim> triangulation_manager;

  DoFManager<dim, degree> dof_manager;

  ConstraintManager<dim, degree, number> constraint_manager;

  SolveContext<dim, degree, number> solve_context;

  std::vector<std::shared_ptr<SolverBase<dim, degree, number>>> solvers;

  SolutionIndexer<dim, number> solution_indexer;

  RefinementManager<dim, degree, number> grid_refiner;
};
 
PRISMS_PF_END_NAMESPACE