Chaste
Commit::1fd4e48e3990e67db148bc1bc4cf6991a0049d0c
AbstractFeAssemblerInterface.hpp
1
/*
2
3
Copyright (c) 2005-2024, University of Oxford.
4
All rights reserved.
5
6
University of Oxford means the Chancellor, Masters and Scholars of the
7
University of Oxford, having an administrative office at Wellington
8
Square, Oxford OX1 2JD, UK.
9
10
This file is part of Chaste.
11
12
Redistribution and use in source and binary forms, with or without
13
modification, are permitted provided that the following conditions are met:
14
* Redistributions of source code must retain the above copyright notice,
15
this list of conditions and the following disclaimer.
16
* Redistributions in binary form must reproduce the above copyright notice,
17
this list of conditions and the following disclaimer in the documentation
18
and/or other materials provided with the distribution.
19
* Neither the name of the University of Oxford nor the names of its
20
contributors may be used to endorse or promote products derived from this
21
software without specific prior written permission.
22
23
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
24
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
25
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
26
ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
27
LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
28
CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
29
GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30
HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
31
LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
32
OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
33
34
*/
35
36
#ifndef ABSTRACTFEASSEMBLERINTERFACE_HPP_
37
#define ABSTRACTFEASSEMBLERINTERFACE_HPP_
38
39
#include <cassert>
40
#include "
UblasCustomFunctions.hpp
"
41
#include "
PetscTools.hpp
"
42
50
template
<
bool
CAN_ASSEMBLE_VECTOR,
bool
CAN_ASSEMBLE_MATRIX>
51
class
AbstractFeAssemblerInterface
:
private
boost::noncopyable
52
{
53
protected
:
55
Vec
mVectorToAssemble
;
56
58
Mat
mMatrixToAssemble
;
59
65
bool
mAssembleMatrix
;
66
68
bool
mAssembleVector
;
69
71
bool
mZeroMatrixBeforeAssembly
;
72
74
bool
mZeroVectorBeforeAssembly
;
75
77
PetscInt
mOwnershipRangeLo
;
78
80
PetscInt
mOwnershipRangeHi
;
81
89
virtual
void
DoAssemble
()=0;
90
91
public
:
92
96
AbstractFeAssemblerInterface
();
97
105
void
SetMatrixToAssemble
(
Mat
& rMatToAssemble,
bool
zeroMatrixBeforeAssembly=
true
);
106
114
void
SetVectorToAssemble
(
Vec
& rVecToAssemble,
bool
zeroVectorBeforeAssembly);
115
119
void
Assemble
()
120
{
121
mAssembleMatrix
= CAN_ASSEMBLE_MATRIX;
122
mAssembleVector
= CAN_ASSEMBLE_VECTOR;
123
DoAssemble
();
124
}
125
129
void
AssembleMatrix
()
130
{
131
assert(CAN_ASSEMBLE_MATRIX);
132
mAssembleMatrix
=
true
;
133
mAssembleVector
=
false
;
134
DoAssemble
();
135
}
136
140
void
AssembleVector
()
141
{
142
assert(CAN_ASSEMBLE_VECTOR);
143
mAssembleMatrix
=
false
;
144
mAssembleVector
=
true
;
145
DoAssemble
();
146
}
147
151
virtual
~AbstractFeAssemblerInterface
()
152
{
153
}
154
};
155
156
template
<
bool
CAN_ASSEMBLE_VECTOR,
bool
CAN_ASSEMBLE_MATRIX>
157
AbstractFeAssemblerInterface<CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX>::AbstractFeAssemblerInterface
()
158
: mVectorToAssemble(nullptr),
159
mMatrixToAssemble(nullptr),
160
mZeroMatrixBeforeAssembly(true),
161
mZeroVectorBeforeAssembly(true)
162
{
163
assert(CAN_ASSEMBLE_VECTOR || CAN_ASSEMBLE_MATRIX);
164
}
165
166
template
<
bool
CAN_ASSEMBLE_VECTOR,
bool
CAN_ASSEMBLE_MATRIX>
167
void
AbstractFeAssemblerInterface<CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX>::SetMatrixToAssemble
(
Mat
& rMatToAssemble,
bool
zeroMatrixBeforeAssembly)
168
{
169
assert(rMatToAssemble);
170
MatGetOwnershipRange(rMatToAssemble, &mOwnershipRangeLo, &mOwnershipRangeHi);
171
172
mMatrixToAssemble = rMatToAssemble;
173
mZeroMatrixBeforeAssembly = zeroMatrixBeforeAssembly;
174
}
175
176
template
<
bool
CAN_ASSEMBLE_VECTOR,
bool
CAN_ASSEMBLE_MATRIX>
177
void
AbstractFeAssemblerInterface<CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX>::SetVectorToAssemble
(
Vec
& rVecToAssemble,
bool
zeroVectorBeforeAssembly)
178
{
179
assert(rVecToAssemble);
180
VecGetOwnershipRange(rVecToAssemble, &mOwnershipRangeLo, &mOwnershipRangeHi);
181
182
mVectorToAssemble = rVecToAssemble;
183
mZeroVectorBeforeAssembly = zeroVectorBeforeAssembly;
184
}
185
186
#endif
// ABSTRACTFEASSEMBLERINTERFACE_HPP_
PetscTools.hpp
UblasCustomFunctions.hpp
AbstractFeAssemblerInterface
Definition
AbstractFeAssemblerInterface.hpp:52
AbstractFeAssemblerInterface::mAssembleVector
bool mAssembleVector
Definition
AbstractFeAssemblerInterface.hpp:68
AbstractFeAssemblerInterface::mMatrixToAssemble
Mat mMatrixToAssemble
Definition
AbstractFeAssemblerInterface.hpp:58
AbstractFeAssemblerInterface::AssembleMatrix
void AssembleMatrix()
Definition
AbstractFeAssemblerInterface.hpp:129
AbstractFeAssemblerInterface::Assemble
void Assemble()
Definition
AbstractFeAssemblerInterface.hpp:119
AbstractFeAssemblerInterface::mZeroVectorBeforeAssembly
bool mZeroVectorBeforeAssembly
Definition
AbstractFeAssemblerInterface.hpp:74
AbstractFeAssemblerInterface::SetMatrixToAssemble
void SetMatrixToAssemble(Mat &rMatToAssemble, bool zeroMatrixBeforeAssembly=true)
Definition
AbstractFeAssemblerInterface.hpp:167
AbstractFeAssemblerInterface::SetVectorToAssemble
void SetVectorToAssemble(Vec &rVecToAssemble, bool zeroVectorBeforeAssembly)
Definition
AbstractFeAssemblerInterface.hpp:177
AbstractFeAssemblerInterface::AbstractFeAssemblerInterface
AbstractFeAssemblerInterface()
Definition
AbstractFeAssemblerInterface.hpp:157
AbstractFeAssemblerInterface::mOwnershipRangeLo
PetscInt mOwnershipRangeLo
Definition
AbstractFeAssemblerInterface.hpp:77
AbstractFeAssemblerInterface::~AbstractFeAssemblerInterface
virtual ~AbstractFeAssemblerInterface()
Definition
AbstractFeAssemblerInterface.hpp:151
AbstractFeAssemblerInterface::DoAssemble
virtual void DoAssemble()=0
AbstractFeAssemblerInterface::mAssembleMatrix
bool mAssembleMatrix
Definition
AbstractFeAssemblerInterface.hpp:65
AbstractFeAssemblerInterface::mVectorToAssemble
Vec mVectorToAssemble
Definition
AbstractFeAssemblerInterface.hpp:55
AbstractFeAssemblerInterface::mZeroMatrixBeforeAssembly
bool mZeroMatrixBeforeAssembly
Definition
AbstractFeAssemblerInterface.hpp:71
AbstractFeAssemblerInterface::mOwnershipRangeHi
PetscInt mOwnershipRangeHi
Definition
AbstractFeAssemblerInterface.hpp:80
AbstractFeAssemblerInterface::AssembleVector
void AssembleVector()
Definition
AbstractFeAssemblerInterface.hpp:140
Mat
PetscInt
Vec
pde
src
solver
AbstractFeAssemblerInterface.hpp
Generated by
1.9.8