14 lines
521 B
Text
14 lines
521 B
Text
#import std
|
|
#import nat
|
|
|
|
break = ~&r**+ zipt*+ ~&lh*~+ ~&lzyCPrX|\+ -*^|\~&tK33 :^/~& 0!*t
|
|
|
|
band = pad0+ ~&rSS+ zipt^*D(~&r,^lrrSPT/~<K33tx zipt^/~&r ~&lSNyCK33+ zipp0)^/~&rx ~&B->NlNSPC ~&bt
|
|
|
|
deconv = # takes a natural number n to the n-dimensional deconvolution function
|
|
|
|
~&?\math..div! iota; ~&!*; @h|\; (~&al^?\~&ar break@alh2faltPrXPRX)^^(
|
|
~&B->NlC~&bt*++ gang@t+ ~~*,
|
|
lapack..dgelsd^^(
|
|
(~&||0.!**+ ~&B^?a\~&Y@a ^lriFhNSS2iDrlYSK7LS2SL2rQ/~&alt band@alh2faltPrDPMX)^|\~&+ gang,
|
|
@t =>~&l ~&L+@r))
|