search

gdmcbain / fenics-tuto-in-skfem

reproducing the examples from the FEniCS tutorial in scikit-fem
14 stars 1 forks source link

pinned dmsh for cylinder in crossflow #24

Closed gdmcbain closed 3 years ago

gdmcbain commented 3 years ago

While learning adaptmesh #23, let's see if we can get a pinned version of dmsh working for the cylinder in cross-flow.

gdmcbain commented 3 years ago

To begin with, pin dmsh at 

dmsh==0.2.2
gdmcbain commented 3 years ago

Then with something ilke


geo = dmsh.Difference(
    dmsh.Rectangle(0.0, length, 0.0, height),
    dmsh.Circle([.2, .2], radius)
)

points, triangles = dmsh.generate(geo, 0.025, tol=1e-9)

we hit

AttributeError: 'MeshTri' object has no attribute 'node_coords'

which does not refer to skfem.MeshTri but rather a class of the same name from meshplex. We read in its CHANGELOG for 0.14.0 that

node_coords is now points

so let's pin

meshplex<0.14
gdmcbain commented 3 years ago

The issue with dmsh is duplicaye nodes nschloe/dmsh#66.