/+dub.sdl:

dependency "lubeck" version="~>0.1"

libs "lapack" "blas"

​

+/

import std.stdio;

import mir.ndslice: magic, repeat, as, slice;

import lubeck: mtimes;

​

void main()

{

    auto n = 5;

    // Magic Square

    auto matrix = n.magic.as!double.slice;

    // [1 1 1 1 1]

    auto vec = 1.repeat(n).as!double.slice;

    // Uses CBLAS for multiplication

    matrix.mtimes(vec).writeln;

    matrix.mtimes(matrix).writeln;

}