OpenBlasSharp 0.1.0

OpenBlasSharp

The purpose of this project is to provide a .NET wrapper for OpenBLAS.

Installation

The NuGet package is available. This package does not include the native DLL. You must download the compiled binary and put libopenblas.dll at the directory of the executable file.

All the classes are in the OpenBlasSharp namespace.

using OpenBlasSharp;

Examples

Dot product

var len = 3;

var rnd = new Random(42);

var x = Enumerable.Range(0, len).Select(i => rnd.NextDouble()).ToArray();
var y = Enumerable.Range(0, len).Select(i => rnd.NextDouble()).ToArray();

fixed (double* px = x)
fixed (double* py = y)
{
    var result = Blas.Ddot(len, px, 1, py, 1);
}

Matrix multiplication

var m = 3;
var n = 5;
var k = 4;

var rnd = new Random(42);

var a = Enumerable.Range(0, m * k).Select(i => rnd.NextDouble()).ToArray();
var b = Enumerable.Range(0, k * n).Select(i => rnd.NextDouble()).ToArray();
var c = new double[m * n];

fixed (double* pa = a)
fixed (double* pb = b)
fixed (double* pc = c)
{
    Blas.Dgemm(
        Order.ColMajor,
        Transpose.NoTrans,
        Transpose.NoTrans,
        m, n, k,
        1.0,
        pa, m,
        pb, k,
        1.0,
        pc, m);
}

Singular value decomposition

var m = 4;
var n = 3;

var rnd = new Random(42);

var a = Enumerable.Range(0, m * n).Select(i => rnd.NextDouble()).ToArray();
var s = new double[Math.Min(m, n)];
var u = new double[m * m];
var vt = new double[n * n];
var work = new double[Math.Min(m, n) - 1];

fixed (double* pa = a)
fixed (double* ps = s)
fixed (double* pu = u)
fixed (double* pvt = vt)
fixed (double* pwork = work)
{
    Lapack.Dgesvd(
        MatrixLayout.ColMajor,
        'A', 'A',
        m, n,
        pa, m,
        ps,
        pu, m,
        pvt, n,
        pwork);
}

OpenBLAS specific functions

var numThreads = OpenBlas.GetNumThreads();

Development status

• Doc comments for the BLAS functions are missing.
• Low-level LAPACK functions with the _work suffix are not supported.
• LAPACK functions that require function pointers are not supported.

License

OpenBlasSharp is available under the BSD-3-Clause license.