9 lines
256 B
Text
9 lines
256 B
Text
|
var [const] GSL=Import("zklGSL"); // libGSL (GNU Scientific Library)
|
||
|
|
fcn dconv1D(f,g){
|
||
|
|
fsz,hsz:=f.len(), g.len() - fsz +1;
|
||
|
|
A:=GSL.Matrix(g.len(),hsz);
|
||
|
|
foreach n,fn in ([0..].zip(f)){ foreach rc in (hsz){ A[rc+n,rc]=fn } }
|
||
|
|
h:=A.AxEQb(g);
|
||
|
|
h
|
||
|
|
}
|