🚀 go-pugleaf

In article <ben.714157437@galaxy>, ben@ph.tn.tudelft.nl (Ben Verwer)
writes:
>The basic question is how to program for images with variable sizes. The 
>solution presented by the author here (malloc a linear array of
>xsize*ysize) is
>imo the best. Those numerical recipies in C readers who propose the famous
>>array of pointers solution are thinking rather 2d. What if you go to 3d:
>an 
>image of pointers? Which is of course implemented as an array of pointers
>to 
>pointers. Good luck!
>As for the author who proposed the function call in the innerloop: you
>must 
>drink lots of coffee at work :)
>
Why not do what I do and have the best of both worlds? Allocate one long
array of nx*ny doubles, then an array of nx double*s, each of which points
to the appropriate section of the previous array. Operations that naturally
run across the entire array can be coded like that, and operations that
make use of the 2Dness of the array can be coded like that. You also have
something you can pass to Fortran. Of course I wrap it up in a C++ class to
hide the ugly details of creating and destroying the objects each time I
use them.

Maynard Handley