Tom Gustafsson

Yes, it's again something to do with 1D plotting. Here is a workaround: ```python ax = plot(basis, phi, nrefs=0) ```

I believe that the difference is because ElementLineP2 uses nodal basis and ElementLinePp uses hierarchical basis.

In most applications you can use either one of those but they don't give exactly the same result because the basis functions are different. There are infinite different bases you...

We can leave this open because of the bug in skfem.visuals.matplotlib.plot for MeshLine1.

I cannot remember what was the remaining work and cannot deduce it by reading this, closing. The workaround is to specify `nrefs=0`, hierarchical bases cannot use `nrefs>0`.

Thanks for the report. Based on this, all named boundaries should be removed from the mesh during adaptive refine: https://github.com/kinnala/scikit-fem/blob/c6b36757a2ba3ebb048abfc9d49a9a9205155a35/skfem/mesh/mesh_tri_1.py#L389 What is check_mesh doing?

Anyhow, don’t expect this to be fixed any time soon. It is a somewhat challenging indexing problem to solve in general.

I don't quite follow what happens in this issue. Can you provide the code of `check_mesh`?

I'm unable to reproduce this with my own minimal example, perhaps it was fixed by another issue, or then I misunderstood the issue: ```python from skfem import * import numpy...