const double old_time_step);
- Epetra_SerialComm trilinos_communicator;
+ Epetra_SerialComm trilinos_communicator;
Triangulation<dim> triangulation;
ConstraintMatrix stokes_constraints;
std::vector<Epetra_Map> stokes_partitioner;
- BlockSparsityPattern stokes_sparsity_pattern;
TrilinosWrappers::BlockSparseMatrix stokes_matrix;
- BlockSparsityPattern stokes_preconditioner_sparsity_pattern;
TrilinosWrappers::BlockSparseMatrix stokes_preconditioner_matrix;
TrilinosWrappers::BlockVector stokes_solution;
ConstraintMatrix temperature_constraints;
Epetra_Map temperature_partitioner;
- SparsityPattern temperature_sparsity_pattern;
TrilinosWrappers::SparseMatrix temperature_mass_matrix;
TrilinosWrappers::SparseMatrix temperature_stiffness_matrix;
TrilinosWrappers::SparseMatrix temperature_matrix;
DoFTools::make_sparsity_pattern (stokes_dof_handler, coupling, csp);
stokes_constraints.condense (csp);
+
+ BlockSparsityPattern stokes_sparsity_pattern;
stokes_sparsity_pattern.copy_from (csp);
stokes_matrix.reinit (stokes_partitioner, stokes_sparsity_pattern);
DoFTools::make_sparsity_pattern (stokes_dof_handler, coupling, csp);
stokes_constraints.condense (csp);
+
+ BlockSparsityPattern stokes_preconditioner_sparsity_pattern;
stokes_preconditioner_sparsity_pattern.copy_from (csp);
stokes_preconditioner_matrix.reinit (stokes_partitioner,
CompressedSetSparsityPattern csp (n_T, n_T);
DoFTools::make_sparsity_pattern (temperature_dof_handler, csp);
temperature_constraints.condense (csp);
+
+ SparsityPattern temperature_sparsity_pattern;
temperature_sparsity_pattern.copy_from (csp);
temperature_matrix.reinit (temperature_partitioner,