Deconvolution/2D+: Difference between revisions

m
→‎{{header|Perl}}: a few simplifications
m (→‎{{header|Perl}}: oops, need 'say' for output)
m (→‎{{header|Perl}}: a few simplifications)
Line 859:
my @dg = shape(@g);
my @hsize;
push @hsize, ($dg[$_] - $df[$_]) + 1 for 0..$#df;
my @toSolve = map { [row(\@g, \@f, \@hsize, $_)] } coords(shape(@g));
rref( \@toSolve );
 
my @h;
my $mag = $hsize[0]*$hsize[1]*$hsize[2];
for (0.. $mag-1) {
push @h, $toSolve[$_][$mag];
}
my @soln;
my $n = 0;
for (coords(@hsize)) {
my($ik,$j,$ki) = split ' ', $_;
$solnh[$ki][$j][$ik] = $htoSolve[$n++][-1];
$n++;
}
@solnh;
}
 
Line 886 ⟶ 880:
my @row;
my @fdim = shape(@f);
@gc = @gc[0..$#fdim]; # clip extraneous values
for (coords(@hsize)) {
my @hc = reverse split ' ', $_;
Line 997 ⟶ 990:
]
);
 
my $d = scalar shape(@g);
print "${d}D array:\n";
print "h =\n";
print "\nff =\n";
 
my @h = deconvolve_N( \@g, \@f );
my @ff = deconvolve_N( \@g, \@h );
 
my $d = scalar shape(@g);
print "${d}D arrayarrays:\n";
print "h =\n";
pretty_print(0,@h);
print "\nff =\n";
pretty_print(0,@ff);</lang>
{{out}}
<pre>3D arrayarrays:
h =
[
2,392

edits