I put up some basic code along these lines here:
Data cube. Turned out going the PDL route didn't offer much -- wasn't a math issue as much as a data access issue -- so I followed the suggestions of japhy and gryng and
built a module. It is slow, but it did the trick for me. nop