Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
amdis
amdis-core
Commits
31bc0cd2
Commit
31bc0cd2
authored
Jul 27, 2020
by
Praetorius, Simon
Browse files
correct boundary condition treecontainer access
parent
a7a4301c
Changes
1
Hide whitespace changes
Inline
Side-by-side
amdis/ProblemStat.inc.hpp
View file @
31bc0cd2
...
...
@@ -308,10 +308,10 @@ addDirichletBC(Predicate const& predicate, RowTreePath row, ColTreePath col, Val
"Function passed to addDirichletBC for `predicate` does not model the Functor<bool(WorldVector)> concept"
);
auto
localView
=
globalBasis_
->
localView
();
auto
i
=
child
(
localView
.
tree
(),
makeTreePath
(
row
)
)
;
auto
j
=
child
(
localView
.
tree
(),
makeTreePath
(
col
)
)
;
auto
rowBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
makeTreePath
(
row
)
);
auto
colBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
makeTreePath
(
col
)
);
auto
i
=
makeTreePath
(
row
);
auto
j
=
makeTreePath
(
col
);
auto
rowBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
i
);
auto
colBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
j
);
auto
valueGridFct
=
makeGridFunction
(
values
,
this
->
gridView
());
...
...
@@ -328,10 +328,10 @@ void ProblemStat<Traits>::
addDirichletBC
(
BoundaryType
id
,
RowTreePath
row
,
ColTreePath
col
,
Values
const
&
values
)
{
auto
localView
=
globalBasis_
->
localView
();
auto
i
=
child
(
localView
.
tree
(),
makeTreePath
(
row
)
)
;
auto
j
=
child
(
localView
.
tree
(),
makeTreePath
(
col
)
)
;
auto
rowBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
makeTreePath
(
row
)
);
auto
colBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
makeTreePath
(
col
)
);
auto
i
=
makeTreePath
(
row
);
auto
j
=
makeTreePath
(
col
);
auto
rowBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
i
);
auto
colBasis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
j
);
auto
valueGridFct
=
makeGridFunction
(
values
,
this
->
gridView
());
...
...
@@ -349,7 +349,7 @@ addPeriodicBC(BoundaryType id, WorldMatrix const& matrix, WorldVector const& vec
auto
basis
=
Dune
::
Functions
::
subspaceBasis
(
*
globalBasis_
,
treepath
());
auto
bc
=
makePeriodicBC
<
SystemMatrix
,
SolutionVector
,
SystemVector
>
(
std
::
move
(
basis
),
{
*
boundaryManager_
,
id
},
{
matrix
,
vector
});
boundaryConditions_
[
localView
.
tree
()][
localView
.
tree
()].
push_back
(
makeUniquePtr
(
std
::
move
(
bc
)));
boundaryConditions_
[
treepath
()][
treepath
()].
push_back
(
makeUniquePtr
(
std
::
move
(
bc
)));
}
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment