Ferrite-FEM / Ferrite-FEM/Ferrite.jl
Reordering of nodal Dirichlet conditions fails
Open
bug
- Dominant language
- Julia
- Stars
- 453
- Forks
- 115
- Avg merge
- 1d 17h
- Merged PRs (30d)
- 11
Description
MWE on master
```julia
using Ferrite
grid = generate_grid(Line, (1,));
ip = Lagrange{RefLine,1}()
dh = DofHandler(grid)
add!(dh, :u, ip)
close!(dh);
ch = ConstraintHandler(dh);
∂Ω = [1] # Triggers error in renumbering
# ∂Ω = getfacetset(grid, "left") # Works
dbc = Dirichlet(:u, ∂Ω, (x, t) -> 0)
add!(ch, dbc);
close!(ch)
perm = [2,1]
Ferrite.renumber!(dh,ch,perm)
```
(Edit by @KnutAM moving ip above dh)
Contributor guide
Assessment
This issue has not been assessed yet.