Smooth triangles for 3d mesh not working?

Hello,

I tried to use the smooth method for 3d mesh of triangles and I got this error:

RuntimeError:

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
***     fenics-support@googlegroups.com
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error:   Unable to find normal.
*** Reason:  Normal vector is not defined in dimension 3 (only defined when the triangle is in R^2.
*** Where:   This error was encountered inside TriangleCell.cpp.
*** Process: 0
***
*** DOLFIN version: 2019.1.0
*** Git changeset:
*** -------------------------------------------------------------------------

The code to reproduce this error:

from fenics import *
mesh_cube = UnitCubeMesh(2,2,2)
bnd_mesh = BoundaryMesh(mesh_cube, 'exterior')
bnd_mesh.smooth(10) # Error

Does anyone if there is a way to use smooth for this kind of mesh ?

1 Like