Add package extension for BlockArrays.jl. #567
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Codecov ReportBase: 92.87% // Head: 92.45% // Decreases project coverage by
Additional details and impacted files@@ Coverage Diff @@
## master #567 +/- ##
==========================================
- Coverage 92.87% 92.45% -0.43%
==========================================
Files 27 28 +1
Lines 4127 4202 +75
==========================================
+ Hits 3833 3885 +52
- Misses 294 317 +23
Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here. ☔ View full report at Codecov. |
|
Nice! |
|
No change, the local matrix is always blocked even when not using a globally blocked matrix. |
fredrikekre
force-pushed
the
fe/block-glue
branch
from
4dfccd9
to
e889c67
Compare
fredrikekre
force-pushed
the
fe/block-glue
branch
from
e889c67
to
334e9f8
Compare
This patch adds support for global block matrices. A global block matrix can be obtained by `create_sparsity_pattern(BlockMatrix, dh, ch)` or by `create_sparsity_pattern(A::BlockMatrix, dh, ch)`. The former infers the block sizes by looking at the fields of the DofHandler, and the latters uses the blocksizes of the (partially) initilized BlockMatrix (i.e. by `undef_blocks`). For global assembly, `start_assemble(::BlockMatrix, ...)` returns a `BlockAssembler` which pre allocates a vector for pre-computing the (block number, local index) pairs for each DoF. BlockArrays only support local condensation of constraints currently and you have to use `apply_assemble!` instead of `assemble!` + `apply!`.
fredrikekre
force-pushed
the
fe/block-glue
branch
from
334e9f8
to
53fb5f7
Compare
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This patch adds support for global block matrices. A global block matrix
can be obtained by
create_sparsity_pattern(BlockMatrix, dh, ch)or bycreate_sparsity_pattern(A::BlockMatrix, dh, ch). The former infers theblock sizes by looking at the fields of the DofHandler, and the latters
uses the blocksizes of the (partially) initilized BlockMatrix (i.e. by
undef_blocks).For global assembly,
start_assemble(::BlockMatrix, ...)returns aBlockAssemblerwhich pre allocates a vector for pre-computing the(block number, local index) pairs for each DoF. BlockArrays only support
local condensation of constraints currently and you have to use
apply_assemble!instead ofassemble!+apply!.