This may be a problem with your physical ids. Can you remove all physical ids from the geo file, and try to export a “bare” mesh, without material id/boundary id information?
I currently use vtk for all my needs. L. > On 24 Jun 2020, at 14:53, Victoria W. <wiedo...@gmail.com> wrote: > > Hi David, > > Thank you for getting back to me. I've tried with a .inp mesh produced by > both gmsh and freeCAD, but I get the same error with my .inp gmsh export and > when I use the FreeCAD export I have to use the read_abaqus() read in method > which has been giving me a bad_alloc error message. I even tried writing the > 1D mesh from scratch by looking at the tutorial mesh, but that also produced > a GridIn error. Do you have any other program suggestions for exporting the > mesh? I'm currently installing Salome and looking to try that next, but > please let me know if there's something else I should be using. > > Thank you and all the best, > Victoria > > On Tuesday, June 23, 2020 at 12:36:01 PM UTC-4, David Wells wrote: > Hi Victoria, > > This is definitely a bug with our gmsh reader - the problem is it assigns the > boundary ids of 0 to all nodes (since they are associated with a material id > of 0), which doesn't make sense since none of the given elements are actually > boundary elements. My best guess is that we never tested this code with 1D > meshes with nonzero codimension. > > The fix to the GMSH code is not immediately obvious. As a temporary > workaround, would it be possible for you to convert your mesh to a .inp? That > reader seems to work correctly for the dim = 1 spacedim = 2 combination. > > Best, > David > > On Thu, Jun 18, 2020 at 10:17 AM Victoria W. <wied...@gmail.com> wrote: > Hi David, > > I've been using read_msh() with a version 1 gmsh 1D mesh for the boundary > definition. Attaching the code snipped of my read_domain() function and the > mesh here. Thank you for your help. > > Best, Victoria > > > On Thursday, June 18, 2020 at 9:35:59 AM UTC-4, David Wells wrote: > Hi Victoria, > > It's difficult to tell what is going on without seeing some code but this > looks like a bug in deal.II: that's the correct GridIn function and > everything should be at the boundary but the assertion still fails for > reasons that appear to be wrong. > > Can you reproduce the same error if you use that mesh with step-34, where you > only modify the call to gi.read_ucd() to gi.read_msh()? If so, would you > please post the mesh to this thread? > > Best, > David > > On Wed, Jun 17, 2020 at 1:28 PM Victoria W. <wied...@gmail.com> wrote: > Hello, > > I've been getting the following error message when I try and input a 1D mesh > in 2D space for a boundary element method simulation. I'm working with > step-34 tutorial program to model an electric field on an electrode, but I > can't get past this mesh input step. I have tried all of the mesh types that > gmsh can output and that dealii will except. Has anyone seen this before, > and does anyone have advice for me? Thank you to anyone who can help or give > advice! > > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dea...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/63e62794-13a9-4562-9d60-4568681c077eo%40googlegroups.com. > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dea...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/0c04e288-82dd-4f2b-b873-ebafbe49987co%40googlegroups.com. > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/688f383d-d44c-47de-819c-fced52432518o%40googlegroups.com. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/588CA626-AFF3-4C49-BF93-3DED45D6AAFB%40gmail.com.
