mexfunction toward octave : computing trouble

[Jérôme LEFEVRE]

Hi,

I try to transpose a Mex function toward octave. The function compute
gradient at node on a unstructured mesh.
I build the linked function with mkoctfile and test it with a simple 2
triangle mesh, but successive result are never homogeneous.
I suspect unsuitable format or allocation gap, but my experience is poor in
c++.
Spurious loop (with no homogenous successive result) in grad.cc are : line
165 to 176 and line 180 to 188

Anyone could help me ? What is wrong ?

Many thanks

grad.cc : gradient function for unstructured mesh
test_grad.m : a simple gradient test with 2 element

test_grad.m
Description: Binary data

grad.cc
Description: Binary data