18 lines
582 B
D
18 lines
582 B
D
import std.stdio, std.range, std.numeric, std.algorithm;
|
|
|
|
T[][] matMul(T)(immutable T[][] A, immutable T[][] B) pure nothrow {
|
|
immutable Bt = B[0].length.iota.map!(i=> B.transversal(i).array)
|
|
.array;
|
|
return A.map!((in a) => Bt.map!(b => a.dotProduct(b)).array).array;
|
|
}
|
|
|
|
void main() {
|
|
immutable a = [[1, 2], [3, 4], [3, 6]];
|
|
immutable b = [[-3, -8, 3,], [-2, 1, 4]];
|
|
|
|
immutable form = "[%([%(%d, %)],\n %)]]";
|
|
writefln("A = \n" ~ form ~ "\n", a);
|
|
writefln("B = \n" ~ form ~ "\n", b);
|
|
writefln("A * B = \n" ~ form, matMul(a, b));
|
|
}
|