RE: crash in mshm_refine

[Marco Vassallo]

> From: address@hidden

> Subject: crash in mshm_refine
> Date: Sat, 6 Jul 2013 09:48:11 +0200
> To: address@hidden
> CC: address@hidden
>
> Hi Marco,
>
> try the following:
>
> msh = msh3m_structured_mesh (1:80,1:80,1:80,1,1:6);
> msh = mshm_refine (msh, 1);
>
> on my system this crashes Octave.

on my system this give an error

   error: out of memory or dimension too large for Octave's index type

which seems to occur here:

 line 334:         dolfin::refine (*r_mesh, mesh, cell_markers);

by the way, how can I use a debugger on a DLD function?

Thanks

marco
>
> Even though this problem could actually be too large to handle,
> the crash should not happen, I think you should check the validity
> of your pointers before writing the output mesh.
>
> c.