Re: crash in mshm_refine

[Carlo De Falco]

On 6 Jul 2013, at 14:12, Marco Vassallo <address@hidden> wrote:

> > 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

that is better that what I get, at least it doesn't crash …

> which seems to occur here:
> 
>  line 334:         dolfin::refine (*r_mesh, mesh, cell_markers);

sounds strange to me, I believe that error will be issued later in the 
function when you resize the output matrices.

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

you should run Octave in a debugger and place a breakpoint in "Fmshm_refine" 
(the "F" is automatically prepended to the function name by the DEFUN_DLD macro)

If you build Octave from source to run in debug mode you can use the 
"run_octave"
scritpt provided in the source tree with the "-g" command-line option.

> Thanks
> 
> marco

c.