ndarray-blas-level2/test/gbmv.js

BLAS Level 2 operations for ndarrays

'use strict';

var RandMatGen = require('./util/rand-matrix-gen.js');
var ndarray = require('ndarray');
var assertCloseTo = require('./util/close-to');
var constants = require('./util/constants');

var gemv = require('../gemv');
var gbmv = require('../gbmv');

describe('GBMV (general banded matrix-vector product)', function () {
  var n = 15;
  var seed;
  var matGen = new RandMatGen(seed, Float64Array);
  var x = ndarray(new Float64Array(n), [n]);
  var x0 = ndarray(new Float64Array(n), [n]);
  var xn = ndarray(new Float64Array(n), [n]);
  var B = ndarray(new Float64Array(n * n), [n, n]);

  it('gbmv', function () {
    for (var t = 0; t < constants.NUM_TESTS; ++t) {
      seed = matGen.setRandomSeed(36);
      matGen.makeBandedMatrix(n, n, 3, 1, B);
      gbmv(B, 3, 1, x0, x);
      gemv(1, B, x0, 0, xn);
      assertCloseTo(x, xn, constants.TEST_TOLERANCE, 'Failure seed value: "' + seed + '".');
    }
  });
});