UserMatrix system_matrix (system_sparsity,
parameters.preconditioning);
system_matrix.copy_from (mass_matrix);
- constraints.condense (system_matrix);
+ constraints.condense (static_cast<SparseMatrix<double>&>(system_matrix));
const unsigned int
solver_steps1 = solve (system_matrix, u, tmp_u_bar),
solver_steps2 = solve (system_matrix, v, tmp_v_bar);
parameters.theta *
parameters.theta,
laplace_matrix);
- constraints.condense (system_matrix);
+ constraints.condense (static_cast<SparseMatrix<double>&>(system_matrix));
if (parameters.extrapolate_old_solutions)
// solve with a hopefully good guess
const unsigned int solver_steps1 = solve (system_matrix, v, right_hand_side1);
system_matrix.copy_from (mass_matrix);
- constraints.condense (system_matrix);
+ constraints.condense (static_cast<SparseMatrix<double>&>(system_matrix));
if (true)
{
Vector<double> tmp (right_hand_side2.size());
parameters.theta *
parameters.theta,
laplace_matrix);
- constraints.condense (system_matrix);
+ constraints.condense (static_cast<SparseMatrix<double>&>(system_matrix));
if (parameters.extrapolate_old_solutions)
// solve with a hopefully good guess
const unsigned int solver_steps1 = solve (system_matrix, u, right_hand_side1);
system_matrix.copy_from (mass_matrix);
- constraints.condense (system_matrix);
+ constraints.condense (static_cast<SparseMatrix<double>&>(system_matrix));
if (true)
{
Vector<double> tmp (right_hand_side2.size());