Ah, I see -- the shape is constrained to be the same between different classes of conditions, just the amplitude varies.
As they point out in the paper, in the general case this is a nonlinear problem, which is outside the computational range of 3dDeconvolve. I once started to write such a program, but abandoned in in December 2007 when no one I talked to was interested in it.
So there are no "future plans" for this capability, unless I can be talked into it, but there were "past plans" -- that never ripened.