I've been running it through a debugger, but what I can glean from it is minimal.
I'll give some code:
for(...){
/* initializing matrix (matrix subroutine below) */
X=matrix(1, x, 1, y);
...
/* freeing matrix */
free_matrix(X, 1, x, 1, y);
}
float **matrix(long nrl, long nrh, long ncl...