Ferrite-FEM / Ferrite-FEM/Ferrite.jl

CellIterator iterates over cells with no DOFs

Open
#976 10 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
453
Forks
115
Avg merge
1d 17h
Merged PRs (30d)
11

Description

I've tried out the fix for #966. However, it seems like the CellIterator still iterates over every cell in the Grid instead of only the cells belonging to the DofHandler. When Iterating over a cell that doesn't belong to the DofHandler, the CellCache will be reinitialized with an empty dof vector.

```julia
using Ferrite

grid = generate_grid(Hexahedron, (3,3,3))

dh = DofHandler(grid)
sdh = SubDofHandler(dh, [1,2,3])
add!(sdh, :u, Lagrange{RefHexahedron,1}()^3)
close!(dh)

for cell in CellIterator(sdh) # works
@assert celldofs(cell) |> !isempty
end

for cell in CellIterator(dh, [1,2,3]) # works
@assert celldofs(cell) |> !isempty
end

for cell in CellIterator(dh) # errors
@assert celldofs(cell) |> !isempty
end
```

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.