Ferrite-FEM / Ferrite-FEM/Ferrite.jl

Reordering of nodal Dirichlet conditions fails

Open
#1,116 1 comment 0 reactions 0 assignees View on GitHub
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

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.