Meshing after applying CSGRotation fails

Issue #22 resolved
Simon Funke created an issue

It seems that generating a mesh from a rotated CSG fails.

A simple example is this:

from mshr import *                                                                                                                                                                                                                             
from dolfin import *                                                                                                                                                                                                                           

left_vessel = Rectangle(dolfin.Point(-0.5, 10.), dolfin.Point(0.5, 20.))                                                                                                                                                                       
left_vessel = CSGRotation(left_vessel, 10)                                                                                                                                                                                                     

mesh = generate_mesh(left_vessel, 10) 

Running this script with FEniCS 1.5 crashes with

Traceback (most recent call last):
  File "generate_bifurcation_mesh.py", line 7, in <module>
    mesh = generate_mesh(left_vessel, 10)
  File "/home/simon/.hashdist/bld/profile/vzrgs5vtcmxt/lib/python2.7/site-packages/mshr.py", line 494, in generate_mesh
    _generate(m, geometry, resolution, backend)
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@fenicsproject.org
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error:   Unable to create mesh from CSG geometry.
*** Reason:  Unhandled geometry dimension -732706896.
*** Where:   This error was encountered inside MeshGenerator.cpp.
*** Process: unknown
*** 
*** DOLFIN version: 1.5.0
*** Git changeset:  
*** -------------------------------------------------------------------------

Comments (1)

  1. Log in to comment