this was indeed due to a bug in the gstat code, I think I found it. I will upload new sources the end of next week; if you want to proceed now, below you'll find the diff -c output. The example you sent me looks nice!
Thanks for notifying me.
Best regards, -- Edzer
Marco Hugentobler wrote:
...I try to interpolate scattered points with a breakline. When i run gstat, i always get errors. When i take out the 'edges:...' command, everything works fine, otherwise, i get the errors.
My command file looks like this:
~/gstat/src$ diff -c gls.c /tmp/gstat-2.4.3/src/gls.c *** gls.c Fri Sep 26 21:14:53 2003 --- /tmp/gstat-2.4.3/src/gls.c Mon Jun 23 20:30:23 2003 *************** *** 223,231 **** }
for (i = 0, global = 1; i < n_vars && global; i++) ! global = (d[i]->sel == d[i]->list ! && d[i]->n_list == d[i]->n_original ! && d[i]->n_list == d[i]->n_sel);
/* * global things: enter whenever (a) first time, (b) local selections or --- 223,229 ---- }
for (i = 0, global = 1; i < n_vars && global; i++)
! global = (d[i]->sel == d[i]->list && d[i]->n_list == d[i]->n_original);
/* * global things: enter whenever (a) first time, (b) local selections or