locally_relevant_dofs,
mpi_communicator);
system_rhs.reinit (locally_owned_dofs, mpi_communicator);
- system_rhs = 0;
+ system_rhs = PetscScalar();
constraints.clear ();
constraints.reinit (locally_relevant_dofs);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
- FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
- Vector<double> cell_rhs (dofs_per_cell);
+ FullMatrix<PetscScalar> cell_matrix (dofs_per_cell, dofs_per_cell);
+ Vector<PetscScalar> cell_rhs (dofs_per_cell);
std::vector<types::global_dof_index> local_dof_indices (dofs_per_cell);
for (; cell!=endc; ++cell)
if (cell->is_locally_owned())
{
- cell_matrix = 0;
- cell_rhs = 0;
+ cell_matrix = PetscScalar();
+ cell_rhs = PetscScalar();
fe_values.reinit (cell);