Processing math: 25%
Reference documentation for deal.II version 9.2.0
\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}
\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}
\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}
\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}
Main Page
Tutorial
Code gallery
Modules
Namespaces
Namespace List
Namespace Members
All
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
z
Functions
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
Variables
a
b
c
d
e
f
i
l
m
n
o
p
r
s
t
u
v
w
z
Typedefs
a
b
c
g
i
l
m
o
p
s
t
v
Enumerations
Enumerator
a
b
c
d
e
f
g
h
i
l
m
n
o
p
r
s
t
u
v
w
Classes
Class List
Class Index
Class Hierarchy
Class Members
All
:
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
~
Functions
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
z
~
Variables
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Typedefs
a
b
c
d
e
f
g
h
i
j
l
m
n
o
p
q
r
s
t
v
w
Enumerations
a
b
c
d
e
f
g
i
m
o
p
r
s
t
u
v
w
z
Enumerator
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
Related Functions
:
a
b
c
d
f
g
h
i
l
m
o
p
r
s
t
u
v
Related Pages
Files
File List
File Members
All
a
b
c
d
e
f
g
h
i
l
m
n
o
p
q
r
s
t
u
v
z
Functions
a
b
c
d
e
f
g
h
i
l
m
n
o
p
s
t
u
v
z
Variables
a
c
d
f
g
i
l
m
p
q
s
t
u
v
Typedefs
Enumerations
Enumerator
m
u
Macros
a
b
d
e
h
i
m
s
t
dealii.org
•
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Modules
Pages
source
lac
petsc_full_matrix.cc
Go to the documentation of this file.
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
{
27
FullMatrix::FullMatrix
()
28
{
29
// empty constructor generate an empty matrix
30
do_reinit
(0, 0);
31
}
32
33
FullMatrix::FullMatrix
(
const
size_type
m,
const
size_type
n)
34
{
35
do_reinit
(
m
,
n
);
36
}
37
38
void
39
FullMatrix::reinit
(
const
size_type
m,
const
size_type
n)
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
50
FullMatrix::do_reinit
(
const
size_type
m,
const
size_type
n)
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
&
61
FullMatrix::get_mpi_communicator
()
const
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
PETScWrappers::destroy_matrix
PetscErrorCode destroy_matrix(Mat &matrix)
Definition:
petsc_compatibility.h:71
PETScWrappers
Definition:
petsc_block_sparse_matrix.h:36
petsc_compatibility.h
PETScWrappers::FullMatrix::get_mpi_communicator
virtual const MPI_Comm & get_mpi_communicator() const override
Definition:
petsc_full_matrix.cc:61
PETScWrappers::MatrixBase::n
size_type n() const
Definition:
petsc_matrix_base.cc:250
MPI_Comm
exceptions.h
PETScWrappers::FullMatrix::do_reinit
void do_reinit(const size_type m, const size_type n)
Definition:
petsc_full_matrix.cc:50
PETScWrappers::FullMatrix::FullMatrix
FullMatrix()
Definition:
petsc_full_matrix.cc:27
DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_OPEN
Definition:
config.h:358
PETScWrappers::MatrixBase::matrix
Mat matrix
Definition:
petsc_matrix_base.h:968
unsigned int
petsc_full_matrix.h
LACExceptions::ExcPETScError
Definition:
exceptions.h:56
PETScWrappers::MatrixBase::m
size_type m() const
Definition:
petsc_matrix_base.cc:237
DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_NAMESPACE_CLOSE
Definition:
config.h:359
AssertThrow
#define AssertThrow(cond, exc)
Definition:
exceptions.h:1531
PETScWrappers::FullMatrix::reinit
void reinit(const size_type m, const size_type n)
Definition:
petsc_full_matrix.cc:39
Generated by
1.8.17