LU decomposition: Difference between revisions

m
→‎{{header|Raku}}: oops, fix variable name
m (→‎{{header|Raku}}: bit more idiomatic)
m (→‎{{header|Raku}}: oops, fix variable name)
Line 3,987:
my @U = matrix-zero $n;
for ^$n X ^$n -> ($i,$j) {
if $j ≥ $i { @U[$i;$j] = @AP[$i;$j] - [+] map { @U[$_;$j] × @L[$i;$_] }, ^$i }
else { @L[$i;$j] = (@AP[$i;$j] - [+] map { @U[$_;$j] × @L[$i;$_] }, ^$j) / @U[$j;$j] }
}
@P, @Aʼ, @L, @U;
2,392

edits