Reference documentation for deal.II version 9.1.1
\(\newcommand{\dealcoloneq}{\mathrel{\vcenter{:}}=}\)
petsc_full_matrix.cc
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 2004 - 2018 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE.md at
12 // the top level directory of deal.II.
13 //
14 // ---------------------------------------------------------------------
15 
16 #include <deal.II/lac/petsc_full_matrix.h>
17 
18 #ifdef DEAL_II_WITH_PETSC
19 
20 # include <deal.II/lac/exceptions.h>
21 # include <deal.II/lac/petsc_compatibility.h>
22 
23 DEAL_II_NAMESPACE_OPEN
24 
25 namespace PETScWrappers
26 {
28  {
29  // empty constructor generate an empty matrix
30  do_reinit(0, 0);
31  }
32 
34  {
35  do_reinit(m, n);
36  }
37 
38  void
40  {
41  // get rid of old matrix and generate a
42  // new one
43  const PetscErrorCode ierr = destroy_matrix(matrix);
44  AssertThrow(ierr == 0, ExcPETScError(ierr));
45 
46  do_reinit(m, n);
47  }
48 
49  void
51  {
52  // use the call sequence indicating only a maximal number of
53  // elements per row for all rows globally
54  const PetscErrorCode ierr =
55  MatCreateSeqDense(PETSC_COMM_SELF, m, n, nullptr, &matrix);
56  AssertThrow(ierr == 0, ExcPETScError(ierr));
57  }
58 
59 
60  const MPI_Comm &
62  {
63  static const MPI_Comm communicator = MPI_COMM_SELF;
64  return communicator;
65  }
66 } // namespace PETScWrappers
67 
68 
69 DEAL_II_NAMESPACE_CLOSE
70 
71 #endif // DEAL_II_WITH_PETSC
PetscErrorCode destroy_matrix(Mat &matrix)
virtual const MPI_Comm & get_mpi_communicator() const override
types::global_dof_index size_type
#define AssertThrow(cond, exc)
Definition: exceptions.h:1519
void do_reinit(const size_type m, const size_type n)
void reinit(const size_type m, const size_type n)