message: Change Add and ParAdd for hypre matrices to use hypre parcsradd
diff: @@ -1563,18 +1563,11 @@ void HypreParMatrix::Destroy()
  HypreParMatrix *Add(double alpha, const HypreParMatrix &A,
  double beta, const HypreParMatrix &B)
  {
  - hypre_ParCSRMatrix *C_hypre =
  - internal::hypre_ParCSRMatrixAdd(const_cast<HypreParMatrix &>(A),
  - const_cast<HypreParMatrix &>(B));
  - MFEM_VERIFY(C_hypre, "error in hypre_ParCSRMatrixAdd");
  -
  - hypre_MatvecCommPkgCreate(C_hypre);
  - HypreParMatrix *C = new HypreParMatrix(C_hypre);
  - *C = 0.0;
  - C->Add(alpha, A);
  - C->Add(beta, B);
  + hypre_ParCSRMatrix *C;
  + hypre_ParcsrAdd(alpha, A, beta, B, &C);
  + hypre_MatvecCommPkgCreate(C);
  - return C;
  + return new HypreParMatrix(C);
  }
  HypreParMatrix * ParMult(const HypreParMatrix *A, const HypreParMatrix *B)
  @@ -1590,7 +1583,8 @@ HypreParMatrix * ParMult(const HypreParMatrix *A, const HypreParMatrix *B)
  HypreParMatrix * ParAdd(const HypreParMatrix *A, const HypreParMatrix *B)
  {
  - hypre_ParCSRMatrix * C = internal::hypre_ParCSRMatrixAdd(*A,*B);
  + hypre_ParCSRMatrix *C;
  + hypre_ParcsrAdd(1.0, *A, 1.0, *B, &C);
  hypre_MatvecCommPkgCreate(C);
  