On Wed, Mar 3, 2010 at 4:48 PM, Raul Miller <[email protected]> wrote:
> Specifically, I am still pondering how to implement deconvolution
> for the general case.
Deconvolution, let's see.
]a=: 8-?17$~6
_3 2 2 4 _8 5
]b=: 8-?17$~6
_5 0 _3 1 8 6
]c=: +//.a*/b
15 _10 _1 _29 12 _37 56 21 _35 _8 30
b 4 :('r=.$0'; 'while. x<:&#y'; 'do. r=.r,t=.y%&{.x';
'y=.}.y-/@:,:t*x'; 'end.'; 'r') c
_3 2 2 4 _8 5
The multi-dimensional case is similar and is exercise to the reader.
Ambrus
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
