[GRASS-user] v.surf.rst error: Input window change

[Daniel Victoria]

Hamish,

The output did look fine.
I reduced the region a bit and rounded the resolution and now I don't
see the message. But I'm almost sure that my old region settings had
round resolution and coordinates.

As for compiling, I'm using the WinGrass 7 binaries so I can't do much
about it. I don't have a build environment set up in my win machine.

Daniel

On Mon, Jun 11, 2012 at 2:09 AM, Hamish <hamish_b at yahoo.com> wrote:
> Daniel wrote:
>> As an update, if I remove the mask raster option, the interpolation
>> finishes ok but it gives me the message at the end...
>>
>> G_ludcmp() failed! n=0  d=0.00
>> G_ludcmp() failed! n=0  d=0.00
>> G_ludcmp() failed! n=0  d=0.00
>> v.surf.rst complete.
>
> does the output look ok?
> what if you change your region slightly? (e.g. round it outwards to nice
> whole numbers with 'g.region -a res=100' followed by resetting back to
> your original resolution with 'g.region res=$orig')
>
>
> If you compiled grass7 with OpenMP support, try 'make distclean' +
> recompiling without and see if you still get the error. You might be able
> to try the same without recompiling by doing
>   export OMP_NUM_THREADS=1
> on the command line before v.surf.rst (on the same command line shell).
>
> the specific error messages is from lib/rst/interp_float/matrix.c
> + lib/gmath/lu.c. I guess a matrix with 0 elements counts as 'singular'
> and so can't be decomposed. (and so I suspect it is happening before the
> OpenMP pragmas are reached, so that could be a red herring)
>
>
> Hamish