I'd like to Fold over higher dimensional objects than lists. Semantically, the format would be something like FoldMulti[f, array, list]. As a precondition, Dimensions[array] == Length[list]. The motivating example is from the recent discussions about inverse of CoefficientList. Let poly be so...

I'd like to know how, given symbols f and g, and a list {a,b,c,...}, to produce the expression: f[f[f[g,a],b],c] I can get the result if I do: Apply[f,{{{g,a},b},c},{0,2}] but then how do I manipulate the list {a,b,c} into {{{g,a},b},c}? Is there a better way of doing it?