+
+ /* Note that I've transposed c__, length, and f relative to the
+ original Fortran code. e.g. length was length(maxfunc,n)
+ in Fortran [ or actually length(maxfunc, maxdims), but by
+ using malloc I can just allocate n ], corresponding to
+ length[n][maxfunc] in C, but I've changed the code to access
+ it as length[maxfunc][n]. That is, the maxfunc direction
+ is the discontiguous one. This makes it easier to resize
+ dynamically (by adding contiguous rows) using realloc, without
+ having to move data around manually. */