From af5a9aaf84ffb7927049c236e67e0aa7ebc40b19 Mon Sep 17 00:00:00 2001 From: Denis Davydov Date: Mon, 28 Mar 2016 22:54:13 +0200 Subject: [PATCH] extend mpi/step-40_direct_solver to complex-valued PETSc --- tests/mpi/step-40_direct_solver.cc | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) 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); -- 2.39.5