#include "pidomus.h"#include <deal.II/lac/trilinos_precondition.h>#include <deal.II/numerics/solution_transfer.h>#include <deal.II/base/quadrature_lib.h>#include <deal.II/base/logstream.h>#include <deal.II/base/function.h>#include <deal.II/base/utilities.h>#include <deal.II/base/conditional_ostream.h>#include <deal.II/base/work_stream.h>#include <deal.II/fe/mapping_q_eulerian.h>#include <deal.II/lac/full_matrix.h>#include <deal.II/lac/solver_cg.h>#include <deal.II/lac/solver_gmres.h>#include <deal.II/lac/sparse_direct.h>#include <deal.II/lac/constraint_matrix.h>#include <deal.II/lac/block_sparsity_pattern.h>#include <deal.II/lac/trilinos_block_vector.h>#include <deal.II/lac/trilinos_sparse_matrix.h>#include <deal.II/lac/trilinos_solver.h>#include <deal.II/lac/block_linear_operator.h>#include <deal.II/lac/linear_operator.h>#include <deal.II/grid/tria.h>#include <deal.II/grid/grid_generator.h>#include <deal.II/grid/tria_accessor.h>#include <deal.II/grid/tria_iterator.h>#include <deal.II/grid/filtered_iterator.h>#include <deal.II/grid/tria_boundary_lib.h>#include <deal.II/grid/grid_tools.h>#include <deal.II/grid/grid_refinement.h>#include <deal.II/dofs/dof_handler.h>#include <deal.II/dofs/dof_renumbering.h>#include <deal.II/dofs/dof_accessor.h>#include <deal.II/dofs/dof_tools.h>#include <deal.II/numerics/vector_tools.h>#include <deal.II/numerics/matrix_tools.h>#include <deal.II/numerics/error_estimator.h>#include <deal.II/distributed/solution_transfer.h>#include <deal.II/distributed/tria.h>#include <deal.II/distributed/grid_refinement.h>#include <deal2lkit/utilities.h>#include <typeinfo>#include <fstream>#include <iostream>#include <sstream>#include <limits>#include <numeric>#include <locale>#include <string>#include <math.h>#include "lac/lac_initializer.h"