This is apparently a problem in mshr. I'l take a look. Best regards
Benjamin Kehlet On 20 August 2015 at 22:43, bill wilson <[email protected]> wrote: > This used to work and only recently started producing the following error: > > ubuntu@ip-172-31-23-251:~$ python > Python 2.7.6 (default, Jun 22 2015, 17:58:13) > [GCC 4.8.2] on linux2 > Type "help", "copyright", "credits" or "license" for more information. > >>> import mshr > >>> from dolfin import * > >>> domain=mshr.Rectangle(Point(0.,0.),Point(8.,8.)) > >>> domain.set_subdomain(1,mshr.Circle(Point(4.,4.),1)) > >>> t4=mshr.generate_mesh(domain,45) > >>> V0 = FunctionSpace(t4, 'DG', 0) > *** Warning: Mesh is empty, unable to create entities of dimension 2. > *** Warning: Mesh is empty, unable to create entities of dimension 2. > *** Warning: Mesh is empty, unable to create entities of dimension 2. > *** Warning: Mesh is empty, unable to create entities of dimension 2. > Traceback (most recent call last): > File "<stdin>", line 1, in <module> > File > "/usr/lib/python2.7/dist-packages/dolfin/functions/functionspace.py", line > 468, in __init__ > FunctionSpaceBase.__init__(self, mesh, element, constrained_domain) > File > "/usr/lib/python2.7/dist-packages/dolfin/functions/functionspace.py", line > 168, in __init__ > dolfin_dofmap = cpp.DofMap(ufc_dofmap, mesh) > File "/usr/lib/python2.7/dist-packages/dolfin/cpp/fem.py", line 660, in > __init__ > _fem.DofMap_swiginit(self,_fem.new_DofMap(*args)) > 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 > *** > *** [email protected] > *** > *** 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 entity. > *** Reason: Mesh entity index 0 out of range [0, 0] for entity of > dimension 2. > *** Where: This error was encountered inside MeshEntity.cpp. > *** Process: 0 > *** > *** DOLFIN version: 1.6.0 > *** Git changeset: unknown > *** > ------------------------------------------------------------------------- > > > _______________________________________________ > fenics mailing list > [email protected] > http://fenicsproject.org/mailman/listinfo/fenics > >
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
