Re: [Rd] akima 0.4-5, interpp() bug = COMMON block problem

2006-02-01 Thread Albrecht Gebhardt
After reading my own question once again, I think I can answer it myself: The code had too _few_ COMMON blocks. It seems the code relied on that, that several local variables would carry their values after being initialised in the first subroutine call until the second call. This is not the case.

[Rd] akima 0.4-5, interpp() bug = COMMON block problem

2006-02-01 Thread Albrecht Gebhardt
Hi, I'm currently hunting a bug in the akima library, especially in the code behind the interpp.old function (bi-variate linear interpolation). It is based on a triangulation algorithm, interpolation at a given point needs to know the triangle which contains this point, then the interpolation is a