mtimes

import mir.ndslice;
import mir.math;

auto a = mininitRcslice!double(3, 5);
auto b = mininitRcslice!double(5, 4);

a[] =
[[-5,  1,  7, 7, -4],
 [-1, -5,  6, 3, -3],
 [-5, -2, -3, 6,  0]];

b[] =
[[-5, -3,  3,  1],
 [ 4,  3,  6,  4],
 [-4, -2, -2,  2],
 [-1,  9,  4,  8],
 [ 9,  8,  3, -2]];

assert(mtimes!(double, double)(a, b) ==
    [[-42,  35,  -7, 77],
     [-69, -21, -42, 21],
     [ 23,  69,   3, 29]]);

import mir.ndslice;
import mir.math;

auto a = mininitRcslice!cdouble(3, 5);
auto b = mininitRcslice!cdouble(5, 4);

a[] =
[[-5 + 0i,  1,  7, 7, -4],
 [-1 + 0i, -5,  6, 3, -3],
 [-5 + 0i, -2, -3, 6,  0]];

b[] =
[[-5 + 0i, -3,  3,  1],
 [ 4 + 0i,  3,  6,  4],
 [-4 + 0i, -2, -2,  2],
 [-1 + 0i,  9,  4,  8],
 [ 9 + 0i, 8,  3, -2]];

assert(mtimes!(cdouble, cdouble)(a, b) ==
    [[-42,  35,  -7, 77],
     [-69, -21, -42, 21],
     [ 23,  69,   3, 29]]);