ASPECT
simulator.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 2024 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 
22 #ifndef _aspect_simulator_h
23 #define _aspect_simulator_h
24 
25 #include <deal.II/base/timer.h>
26 #include <deal.II/base/parameter_handler.h>
27 #include <deal.II/base/conditional_ostream.h>
28 #include <deal.II/base/symmetric_tensor.h>
29 
30 DEAL_II_DISABLE_EXTRA_DIAGNOSTICS
31 
32 #include <deal.II/lac/affine_constraints.h>
33 
34 #include <deal.II/distributed/tria.h>
35 
36 #include <deal.II/dofs/dof_handler.h>
37 #include <deal.II/dofs/dof_tools.h>
38 
39 #include <deal.II/fe/fe_system.h>
40 #include <deal.II/fe/mapping.h>
41 #include <deal.II/base/tensor_function.h>
42 
43 DEAL_II_ENABLE_EXTRA_DIAGNOSTICS
44 
45 #include <aspect/global.h>
69 #include <aspect/advection_field.h>
70 
71 #include <boost/iostreams/tee.hpp>
72 #include <boost/iostreams/stream.hpp>
73 
74 #include <memory>
75 #include <thread>
76 
77 
78 namespace aspect
79 {
80  template <int dim>
81  class MeltHandler;
82 
83  template <int dim>
84  class NewtonHandler;
85 
86  template <int dim>
88 
89  namespace StokesSolver
90  {
91  template <int dim>
92  class Direct;
93  }
94 
95  template <int dim, int velocity_degree>
97 
98  namespace MeshDeformation
99  {
100  template <int dim>
102  }
103 
104  template <int dim>
106 
107  namespace internal
108  {
109  namespace Assembly
110  {
111  namespace Scratch
112  {
113  template <int dim> struct StokesPreconditioner;
114  template <int dim> struct StokesSystem;
115  template <int dim> struct AdvectionSystem;
116  }
117 
118  namespace CopyData
119  {
120  template <int dim> struct StokesPreconditioner;
121  template <int dim> struct StokesSystem;
122  template <int dim> struct AdvectionSystem;
123  }
124  }
125  }
126 
127  namespace Assemblers
128  {
129  template <int dim> class Interface;
130  template <int dim> class Manager;
131  }
132 
134  {
138  double residual;
139  double residual_old;
142  std::pair<double,double> stokes_residuals;
143  };
144 
148  template <int dim>
150  {
152  :
153  scalar_moment_of_inertia(numbers::signaling_nan<double>()),
154  scalar_angular_momentum(numbers::signaling_nan<double>()),
155  scalar_rotation(numbers::signaling_nan<double>()),
156  tensor_moment_of_inertia(numbers::signaling_nan<SymmetricTensor<2,dim>>()),
157  tensor_angular_momentum(numbers::signaling_nan<Tensor<1,dim>>()),
158  tensor_rotation(numbers::signaling_nan<Tensor<1,dim>>())
159  {};
160 
168 
172  SymmetricTensor<2,dim> tensor_moment_of_inertia;
173  Tensor<1,dim> tensor_angular_momentum;
174  Tensor<1,dim> tensor_rotation;
175  };
176 
180  DeclExceptionMsg(ExcNonlinearSolverNoConvergence,
181  "Nonlinear solver failed to converge in the prescribed number of steps. "
182  "Consider changing `Max nonlinear iterations` or `Nonlinear solver failure "
183  "strategy`.");
184 
192  template <int dim>
193  class Simulator
194  {
195  public:
209  Simulator (const MPI_Comm mpi_communicator,
210  ParameterHandler &prm);
211 
216  ~Simulator ();
217 
229  static
230  void declare_parameters (ParameterHandler &prm);
231 
240  void run ();
241 
253  void
254  write_plugin_graph (std::ostream &output_stream) const;
255 
260 
265 
267 
268  private:
269 
289 
307  {
308  IntermediaryConstructorAction (const std::function<void ()> &action);
309  };
310 
325  void setup_dofs ();
326 
335  void setup_introspection ();
336 
348  void set_initial_temperature_and_compositional_fields ();
349 
365  void compute_initial_pressure_field ();
366 
372  void compute_initial_velocity_boundary_constraints (AffineConstraints<double> &constraints);
373 
379  void compute_current_velocity_boundary_constraints (AffineConstraints<double> &constraints);
380 
393  void compute_current_constraints ();
394 
408  double compute_pressure_scaling_factor () const;
409 
421  void start_timestep ();
422 
430  void solve_timestep ();
431 
443  void solve_single_advection_single_stokes ();
444 
457  void solve_no_advection_iterated_stokes ();
458 
469  void solve_no_advection_single_stokes ();
470 
482  void solve_first_timestep_only_single_stokes ();
483 
496  void solve_iterated_advection_and_stokes ();
497 
510  void solve_single_advection_iterated_stokes ();
511 
524  void solve_no_advection_iterated_defect_correction_stokes ();
525 
538  void solve_single_advection_iterated_defect_correction_stokes ();
539 
553  void solve_iterated_advection_and_defect_correction_stokes ();
554 
573  void solve_iterated_advection_and_newton_stokes (bool use_newton_iterations);
574 
594  void solve_single_advection_and_iterated_newton_stokes (bool use_newton_iterations);
595 
607  void solve_single_advection_no_stokes ();
608 
620  void solve_no_advection_no_stokes ();
621 
630  void build_stokes_preconditioner ();
631 
639  void build_advection_preconditioner (const AdvectionField &advection_field,
641  const double diagonal_strengthening);
642 
649  void assemble_stokes_system ();
650 
665  double assemble_and_solve_temperature (const double &initial_residual = 0,
666  double *residual = nullptr);
667 
683  std::vector<double> assemble_and_solve_composition (const std::vector<double> &initial_residual = {},
684  std::vector<double> *residual = nullptr);
685 
700  double assemble_and_solve_stokes (const double &initial_nonlinear_residual = 0,
701  double *nonlinear_residual = nullptr);
702 
715  void do_one_defect_correction_Stokes_step(DefectCorrectionResiduals &dcr,
716  const bool use_picard);
717 
725  void assemble_advection_system (const AdvectionField &advection_field);
726 
737  double solve_advection (const AdvectionField &advection_field);
738 
743  void interpolate_particle_properties (const std::vector<AdvectionField> &advection_fields);
744 
821  std::pair<double,double>
822  solve_stokes (LinearAlgebra::BlockVector &solution_vector);
823 
836  void postprocess ();
837 
853  void refine_mesh (const unsigned int max_grid_level);
854 
870  unsigned int determine_last_good_snapshot() const;
871 
883  void create_snapshot();
884 
898  void resume_from_snapshot();
899 
906  template <class Archive>
907  void serialize (Archive &ar, const unsigned int version);
924  Table<2,DoFTools::Coupling>
925  setup_system_matrix_coupling () const;
926 
934  void setup_system_matrix (const std::vector<IndexSet> &system_partitioning);
935 
947  void setup_system_preconditioner (const std::vector<IndexSet> &system_partitioning);
948 
970  std::unique_ptr<Assemblers::Manager<dim>> assemblers;
971 
980  void set_assemblers ();
981 
992  void set_advection_assemblers ();
993 
1003  void set_stokes_assemblers ();
1004 
1011  void assemble_stokes_preconditioner ();
1012 
1020  void
1021  local_assemble_stokes_preconditioner (const typename DoFHandler<dim>::active_cell_iterator &cell,
1024 
1032  void
1033  copy_local_to_global_stokes_preconditioner (const internal::Assembly::CopyData::StokesPreconditioner<dim> &data);
1034 
1042  void
1043  local_assemble_stokes_system (const typename DoFHandler<dim>::active_cell_iterator &cell,
1046 
1054  void
1055  copy_local_to_global_stokes_system (const internal::Assembly::CopyData::StokesSystem<dim> &data);
1056 
1064  void
1065  local_assemble_advection_face_terms(const AdvectionField &advection_field,
1066  const typename DoFHandler<dim>::active_cell_iterator &cell,
1076  void
1077  local_assemble_advection_system (const AdvectionField &advection_field,
1078  const Vector<double> &viscosity_per_cell,
1079  const typename DoFHandler<dim>::active_cell_iterator &cell,
1082 
1090  void
1091  copy_local_to_global_advection_system (const AdvectionField &advection_field,
1093 
1118  void make_pressure_rhs_compatible(LinearAlgebra::BlockVector &vector);
1119 
1129  template <typename T>
1130  void get_artificial_viscosity (Vector<T> &viscosity_per_cell,
1131  const AdvectionField &advection_field,
1132  const bool skip_interior_cells = false) const;
1133 
1166  double normalize_pressure(LinearAlgebra::BlockVector &vector) const;
1167 
1207  void denormalize_pressure(const double pressure_adjustment,
1208  LinearAlgebra::BlockVector &vector) const;
1209 
1218  void apply_limiter_to_dg_solutions (const AdvectionField &advection_field);
1219 
1236  void compute_unique_advection_support_points (const std::vector<AdvectionField> &advection_fields,
1237  std::vector<Point<dim>> &unique_support_points,
1238  std::vector<std::vector<unsigned int>> &support_point_index_by_field) const;
1239 
1262  void compute_reactions ();
1263 
1273  void update_solution_vectors_with_reaction_results (const unsigned int block_index,
1274  const LinearAlgebra::BlockVector &distributed_vector,
1275  const LinearAlgebra::BlockVector &distributed_reaction_vector);
1276 
1287  void initialize_current_linearization_point ();
1288 
1308  void interpolate_material_output_into_advection_field (const std::vector<AdvectionField> &adv_field);
1309 
1310 
1318  void interpolate_onto_velocity_system(const TensorFunction<1,dim> &func,
1319  LinearAlgebra::Vector &vec) const;
1320 
1321 
1336  void setup_nullspace_constraints(AffineConstraints<double> &constraints);
1337 
1338 
1351  void remove_nullspace(LinearAlgebra::BlockVector &solution,
1352  LinearAlgebra::BlockVector &distributed_stokes_solution) const;
1353 
1369  compute_net_angular_momentum(const bool use_constant_density,
1370  const LinearAlgebra::BlockVector &solution,
1371  const bool limit_to_top_faces = false) const;
1372 
1389  void remove_net_angular_momentum(const bool use_constant_density,
1390  LinearAlgebra::BlockVector &relevant_dst,
1391  LinearAlgebra::BlockVector &tmp_distributed_stokes,
1392  const bool limit_to_top_faces = false) const;
1393 
1401  void replace_outflow_boundary_ids(const unsigned int boundary_id_offset,
1402  const bool is_composition,
1403  const unsigned int composition_index);
1404 
1412  void restore_outflow_boundary_ids(const unsigned int boundary_id_offset);
1413 
1427  void remove_net_linear_momentum(const bool use_constant_density,
1428  LinearAlgebra::BlockVector &relevant_dst,
1429  LinearAlgebra::BlockVector &tmp_distributed_stokes) const;
1430 
1438  double get_maximal_velocity (const LinearAlgebra::BlockVector &solution) const;
1439 
1452  double get_entropy_variation (const double average_field,
1453  const AdvectionField &advection_field) const;
1454 
1463  std::pair<double,double>
1464  get_extrapolated_advection_field_range (const AdvectionField &advection_field) const;
1465 
1474  void exchange_refinement_flags();
1475 
1476 
1485  void maybe_write_timing_output () const;
1486 
1494  bool maybe_write_checkpoint (const time_t last_checkpoint_time,
1495  const bool force_writing_checkpoint);
1496 
1513  bool maybe_do_initial_refinement (const unsigned int max_refinement_level);
1514 
1523  void maybe_refine_mesh (const double new_time_step,
1524  unsigned int &max_refinement_level);
1525 
1533  void advance_time (const double step_size);
1534 
1542  double
1543  compute_viscosity(internal::Assembly::Scratch::AdvectionSystem<dim> &scratch,
1544  const double global_u_infty,
1545  const double global_field_variation,
1546  const double average_field,
1547  const double global_entropy_variation,
1548  const double cell_diameter,
1549  const AdvectionField &advection_field) const;
1550 
1559  void
1560  compute_advection_system_residual(internal::Assembly::Scratch::AdvectionSystem<dim> &scratch,
1561  const double average_field,
1562  const AdvectionField &advection_field,
1563  double &max_residual,
1564  double &max_velocity,
1565  double &max_density,
1566  double &max_specific_heat,
1567  double &conductivity) const;
1568 
1583  bool
1584  stokes_matrix_depends_on_solution () const;
1585 
1595  bool
1596  stokes_A_block_is_symmetric () const;
1597 
1610  void
1611  check_consistency_of_formulation ();
1612 
1618  void
1619  select_default_solver_and_averaging ();
1620 
1631  void
1632  check_consistency_of_boundary_conditions () const;
1633 
1637  double
1638  compute_initial_newton_residual ();
1639 
1648  double
1649  compute_Eisenstat_Walker_linear_tolerance(const bool EisenstatWalkerChoiceOne,
1650  const double maximum_linear_stokes_solver_tolerance,
1651  const double linear_stokes_solver_tolerance,
1652  const double stokes_residual,
1653  const double newton_residual,
1654  const double newton_residual_old);
1655 
1665  void output_statistics();
1666 
1678  double
1679  compute_initial_stokes_residual();
1680 
1691 
1697  std::unique_ptr<MeltHandler<dim>> melt_handler;
1698 
1704  std::unique_ptr<NewtonHandler<dim>> newton_handler;
1705 
1707 
1709 
1717  std::unique_ptr<VolumeOfFluidHandler<dim>> volume_of_fluid_handler;
1718 
1720 
1721 
1723 
1728  std::ofstream log_file_stream;
1729 
1730  using TeeDevice = boost::iostreams::tee_device<std::ostream, std::ofstream>;
1731  using TeeStream = boost::iostreams::stream<TeeDevice>;
1732 
1735 
1740  ConditionalOStream pcout;
1741 
1751  TableHandler statistics;
1752 
1769 
1770  mutable TimerOutput computing_timer;
1771 
1776  Timer wall_timer;
1777 
1783 
1789  unsigned int last_checkpoint_id;
1790 
1799 
1808  const std::unique_ptr<InitialTopographyModel::Interface<dim>> initial_topography_model;
1809  const std::unique_ptr<GeometryModel::Interface<dim>> geometry_model;
1811  const std::unique_ptr<MaterialModel::Interface<dim>> material_model;
1812  const std::unique_ptr<GravityModel::Interface<dim>> gravity_model;
1813 
1817  const std::unique_ptr<PrescribedStokesSolution::Interface<dim>> prescribed_stokes_solution;
1818 
1829  std::shared_ptr<InitialTemperature::Manager<dim>> initial_temperature_manager;
1830  std::shared_ptr<InitialComposition::Manager<dim>> initial_composition_manager;
1831 
1832  const std::unique_ptr<AdiabaticConditions::Interface<dim>> adiabatic_conditions;
1833 #ifdef ASPECT_WITH_WORLD_BUILDER
1834 
1845  std::shared_ptr<WorldBuilder::World> world_builder;
1846 #endif
1849  const std::unique_ptr<BoundaryHeatFlux::Interface<dim>> boundary_heat_flux;
1850 
1854  std::vector<Particle::Manager<dim>> particle_managers;
1855 
1863  double time;
1864  double time_step;
1866  unsigned int timestep_number;
1867  unsigned int pre_refinement_step;
1868  unsigned int nonlinear_iteration;
1896  parallel::distributed::Triangulation<dim> triangulation;
1899 
1902 
1911  std::unique_ptr<Mapping<dim>> mapping;
1912 
1913  const FESystem<dim> finite_element;
1914 
1915  DoFHandler<dim> dof_handler;
1916 
1918 
1930  AffineConstraints<double> constraints;
1931  AffineConstraints<double> current_constraints;
1932 
1938 
1945 
1953 
1973 
1979 
1993 
1998 
2000 
2001  // only used if is_compressible()
2003 
2004  // only used if operator split is enabled
2006 
2007 
2008 
2009  std::unique_ptr<LinearAlgebra::PreconditionAMG> Amg_preconditioner;
2010  std::unique_ptr<LinearAlgebra::PreconditionBase> Mp_preconditioner;
2011 
2017 
2022  private:
2023 
2029  std::unique_ptr<MeshDeformation::MeshDeformationHandler<dim>> mesh_deformation;
2030 
2034  std::unique_ptr<StokesMatrixFreeHandler<dim>> stokes_matrix_free;
2035 
2039  std::unique_ptr<StokesSolver::Direct<dim>> stokes_direct;
2040 
2041 
2042  friend class boost::serialization::access;
2043  friend class SimulatorAccess<dim>;
2045  friend class VolumeOfFluidHandler<dim>;
2046  friend class StokesMatrixFreeHandler<dim>;
2047  template <int dimension, int velocity_degree> friend class StokesMatrixFreeHandlerLocalSmoothingImplementation;
2048  friend struct Parameters<dim>;
2049  };
2050 }
2051 
2052 
2053 #endif
The NullspaceRemoval struct.
Definition: parameters.h:133
unsigned int nonlinear_iteration
Definition: simulator.h:1868
BoundaryVelocity::Manager< dim > boundary_velocity_manager
Definition: simulator.h:1847
const std::unique_ptr< AdiabaticConditions::Interface< dim > > adiabatic_conditions
Definition: simulator.h:1832
std::shared_ptr< InitialComposition::Manager< dim > > initial_composition_manager
Definition: simulator.h:1830
void write_plugin_graph(std::ostream &output_stream)
parallel::distributed::Triangulation< dim > triangulation
Definition: simulator.h:1896
TimerOutput computing_timer
Definition: simulator.h:1770
BoundaryTemperature::Manager< dim > boundary_temperature_manager
Definition: simulator.h:1814
boost::iostreams::tee_device< std::ostream, std::ofstream > TeeDevice
Definition: simulator.h:1730
std::unique_ptr< StokesSolver::Direct< dim > > stokes_direct
Definition: simulator.h:2039
::TrilinosWrappers::MPI::BlockVector BlockVector
Definition: global.h:268
double pressure_scaling
Definition: simulator.h:1944
Tensor< 1, dim > tensor_rotation
Definition: simulator.h:174
double old_time_step
Definition: simulator.h:1865
Tensor< 1, dim > tensor_angular_momentum
Definition: simulator.h:173
AffineConstraints< double > constraints
Definition: simulator.h:1930
bool assemble_newton_stokes_system
Definition: simulator.h:2015
LinearAlgebra::BlockVector old_solution
Definition: simulator.h:1995
const IntermediaryConstructorAction post_geometry_model_creation_action
Definition: simulator.h:1810
std::unique_ptr< VolumeOfFluidHandler< dim > > volume_of_fluid_handler
Definition: simulator.h:1717
bool rebuild_stokes_matrix
Definition: simulator.h:2013
Parameters< dim > parameters
Definition: simulator.h:1690
SymmetricTensor< 2, dim > tensor_moment_of_inertia
Definition: simulator.h:172
AffineConstraints< double > current_constraints
Definition: simulator.h:1931
std::unique_ptr< MeshDeformation::MeshDeformationHandler< dim > > mesh_deformation
Definition: simulator.h:2029
LinearAlgebra::BlockVector system_rhs
Definition: simulator.h:1997
std::pair< double, double > stokes_residuals
Definition: simulator.h:142
std::vector< Particle::Manager< dim > > particle_managers
Definition: simulator.h:1854
std::thread output_statistics_thread
Definition: simulator.h:1798
LateralAveraging< dim > lateral_averaging
Definition: simulator.h:1887
DeclExceptionMsg(ExcNonlinearSolverNoConvergence, "Nonlinear solver failed to converge in the prescribed number of steps. " "Consider changing `Max nonlinear iterations` or `Nonlinear solver failure " "strategy`.")
TeeStream iostream_tee_stream
Definition: simulator.h:1734
::TrilinosWrappers::MPI::Vector Vector
Definition: global.h:262
MeshRefinement::Manager< dim > mesh_refinement_manager
Definition: simulator.h:1900
std::size_t statistics_last_write_size
Definition: simulator.h:1767
const std::unique_ptr< MaterialModel::Interface< dim > > material_model
Definition: simulator.h:1811
unsigned int nonlinear_solver_failures
Definition: simulator.h:1869
double last_pressure_normalization_adjustment
Definition: simulator.h:1937
::TrilinosWrappers::BlockSparseMatrix BlockSparseMatrix
Definition: global.h:279
double total_walltime_until_last_snapshot
Definition: simulator.h:1782
void declare_parameters(ParameterHandler &prm)
const std::unique_ptr< GeometryModel::Interface< dim > > geometry_model
Definition: simulator.h:1809
LinearAlgebra::BlockVector operator_split_reaction_vector
Definition: simulator.h:2005
DoFHandler< dim > dof_handler
Definition: simulator.h:1915
std::unique_ptr< StokesMatrixFreeHandler< dim > > stokes_matrix_free
Definition: simulator.h:2034
unsigned int last_checkpoint_id
Definition: simulator.h:1789
std::unique_ptr< Assemblers::Manager< dim > > assemblers
Definition: simulator.h:970
bool rebuild_sparsity_and_matrices
Definition: simulator.h:2012
const IntermediaryConstructorAction post_signal_creation
Definition: simulator.h:1708
double global_Omega_diameter
Definition: simulator.h:1897
unsigned int timestep_number
Definition: simulator.h:1866
bool assemble_newton_stokes_matrix
Definition: simulator.h:2014
ConditionalOStream pcout
Definition: simulator.h:1740
unsigned int pre_refinement_step
Definition: simulator.h:1867
TimeStepping::Manager< dim > time_stepping_manager
Definition: simulator.h:1878
std::size_t statistics_last_hash
Definition: simulator.h:1768
TableHandler statistics
Definition: simulator.h:1751
Introspection< dim > introspection
Definition: simulator.h:1719
MPI_Comm mpi_communicator
Definition: simulator.h:1722
typename Parameters< dim >::NullspaceRemoval NullspaceRemoval
Definition: simulator.h:264
LinearAlgebra::BlockSparseMatrix system_preconditioner_matrix
Definition: simulator.h:1992
bool do_pressure_rhs_compatibility_modification
Definition: simulator.h:1952
LinearAlgebra::BlockVector current_linearization_point
Definition: simulator.h:1999
BoundaryConvectiveHeating::Manager< dim > boundary_convective_heating_manager
Definition: simulator.h:1815
LinearAlgebra::BlockVector old_old_solution
Definition: simulator.h:1996
BoundaryTraction::Manager< dim > boundary_traction_manager
Definition: simulator.h:1848
LinearAlgebra::BlockVector pressure_shape_function_integrals
Definition: simulator.h:2002
std::ofstream log_file_stream
Definition: simulator.h:1728
std::unique_ptr< LinearAlgebra::PreconditionAMG > Amg_preconditioner
Definition: simulator.h:2009
std::shared_ptr< InitialTemperature::Manager< dim > > initial_temperature_manager
Definition: simulator.h:1829
TeeDevice iostream_tee_device
Definition: simulator.h:1733
boost::iostreams::stream< TeeDevice > TeeStream
Definition: simulator.h:1731
LinearAlgebra::BlockSparseMatrix system_matrix
Definition: simulator.h:1972
std::unique_ptr< MeltHandler< dim > > melt_handler
Definition: simulator.h:1697
LinearAlgebra::BlockVector inverse_lumped_mass_matrix
Definition: simulator.h:1978
bool rebuild_stokes_preconditioner
Definition: simulator.h:2016
std::unique_ptr< Mapping< dim > > mapping
Definition: simulator.h:1911
const std::unique_ptr< PrescribedStokesSolution::Interface< dim > > prescribed_stokes_solution
Definition: simulator.h:1817
HeatingModel::Manager< dim > heating_model_manager
Definition: simulator.h:1901
SimulatorSignals< dim > signals
Definition: simulator.h:1706
const std::unique_ptr< BoundaryHeatFlux::Interface< dim > > boundary_heat_flux
Definition: simulator.h:1849
const FESystem< dim > finite_element
Definition: simulator.h:1913
bool simulator_is_past_initialization
Definition: simulator.h:288
const std::unique_ptr< InitialTopographyModel::Interface< dim > > initial_topography_model
Definition: simulator.h:1808
LinearAlgebra::BlockVector solution
Definition: simulator.h:1994
std::unique_ptr< NewtonHandler< dim > > newton_handler
Definition: simulator.h:1704
typename Parameters< dim >::NonlinearSolver NonlinearSolver
Definition: simulator.h:259
Postprocess::Manager< dim > postprocess_manager
Definition: simulator.h:1917
const std::unique_ptr< GravityModel::Interface< dim > > gravity_model
Definition: simulator.h:1812
double newton_residual_for_derivative_scaling_factor
Definition: simulator.h:141
double global_volume
Definition: simulator.h:1898
BoundaryComposition::Manager< dim > boundary_composition_manager
Definition: simulator.h:1816
std::unique_ptr< LinearAlgebra::PreconditionBase > Mp_preconditioner
Definition: simulator.h:2010
::TrilinosWrappers::PreconditionILU PreconditionILU
Definition: global.h:302