Skip to content

Commit a19254f

Browse files
chmerdonchmerdon
andauthored
Two bugfixes for ItemIntegratorDG (#93)
Co-authored-by: chmerdon <christian.merdon@wias-berlin.de>
1 parent 21ee674 commit a19254f

File tree

2 files changed

+9
-2
lines changed

2 files changed

+9
-2
lines changed

CHANGELOG.md

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -4,6 +4,7 @@
44

55
### Fixed
66
- fixed type issues in LinearOperator and ItemIntegrator
7+
- fixed region awareness and output array sizing for boundary entities in ItemIntegratorDG
78

89
## v1.7.0
910

src/common_operators/item_integrator_dg.jl

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -234,6 +234,12 @@ function build_assembler!(O::ItemIntegratorDG{Tv}, FE_args::Array{<:FEVectorBloc
234234
end
235235
end
236236

237+
if itemregions[item] > 0
238+
if !(visit_region[itemregions[item]])
239+
continue
240+
end
241+
end
242+
237243
for qp in 1:nweights
238244
## evaluate arguments at quadrature points
239245
fill!(input_args, 0)
@@ -364,9 +370,9 @@ function evaluate(O::ItemIntegratorDG{Tv, UT}, sol; kwargs...) where {Tv, UT}
364370
elseif AT <: ON_EDGES
365371
nitems = size(grid[EdgeNodes], 2)
366372
elseif AT <: ON_BFACES
367-
nitems = size(grid[BFaceNodes], 2)
373+
nitems = size(grid[FaceNodes], 2)
368374
elseif AT <: ON_BEDGES
369-
nitems = size(grid[BEdgeNodes], 2)
375+
nitems = size(grid[EdgeNodes], 2)
370376
end
371377
if O.parameters[:piecewise]
372378
b = zeros(eltype(sol[1].entries), O.parameters[:resultdim], nitems)

0 commit comments

Comments
 (0)