From: Denis Davydov Date: Mon, 28 Mar 2016 20:54:13 +0000 (+0200) Subject: extend mpi/step-40_direct_solver to complex-valued PETSc X-Git-Tag: v8.5.0-rc1~1153^2~32 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=af5a9aaf84ffb7927049c236e67e0aa7ebc40b19;p=dealii.git extend mpi/step-40_direct_solver to complex-valued PETSc --- diff --git a/tests/mpi/step-40_direct_solver.cc b/tests/mpi/step-40_direct_solver.cc index 82717cb7ca..26194dd7c1 100644 --- a/tests/mpi/step-40_direct_solver.cc +++ b/tests/mpi/step-40_direct_solver.cc @@ -138,7 +138,7 @@ namespace Step40 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); @@ -182,8 +182,8 @@ namespace Step40 const unsigned int dofs_per_cell = fe.dofs_per_cell; const unsigned int n_q_points = quadrature_formula.size(); - FullMatrix cell_matrix (dofs_per_cell, dofs_per_cell); - Vector cell_rhs (dofs_per_cell); + FullMatrix cell_matrix (dofs_per_cell, dofs_per_cell); + Vector cell_rhs (dofs_per_cell); std::vector local_dof_indices (dofs_per_cell); @@ -193,8 +193,8 @@ namespace Step40 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);