mshr fails on simple off file


46
views
0
4 weeks ago by

I am trying to create a very simple mesh based on .off format exported from FREECAD with the following code:

from dolfin import *
from mshr import *

domain = Surface3D("k1.off")
mesh = generate_mesh(domain,2)

plot(mesh, interactive=True)

k1.off only holds a circle.

This leads to  the following output:
Generating mesh with CGAL 3D mesh generator
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[0]PETSC ERROR: to get more information on the crash.
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 59.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

What the issue here? Do I something wrong? Or is this a bug in fenics.

Does anybody have the sample problem?

Thanks, Markus

File attached: k1.off (9.93 KB)

Community: FEniCS Project
Please login to add an answer/comment or follow this question.

Similar posts:
Search »