The Computer Language
Benchmarks Game

spectral-norm Dart program

source code

// The Computer Language Benchmarks Game
// http://benchmarksgame.alioth.debian.org/
//
// contributed by Jos Hirth
// based on the JavaScript version by Ian Osgood with modifications by Isaac Gouy

import 'dart:math' as Math;
import 'dart:typed_data';

double A(int i, int j) {
  int div = ((i + j) * (i + j + 1) >> 1) + i + 1;
  return 1.0 / div;
}

void Au(Float64List u, Float64List w) {
  int len = u.length;
  for (int i = 0; i < len; ++i) {
    double t = 0.0;
    for (int j = 0; j < len; ++j) {
      t += A(i, j) * u[j];
    }
    w[i] = t;
  }
}

void Atu(Float64List w, Float64List v) {
  int len = w.length;
  for (int i = 0; i < len; ++i) {
    double t = 0.0;
    for (int j = 0; j < len; ++j) {
      t += A(j, i) * w[j];
    }
    v[i] = t;
  }
}

void AtAu(Float64List u, Float64List v, Float64List w) {
  Au(u, w);
  Atu(w, v);
}

double spectralNorm(n) {
  var u = new Float64List(n)..fillRange(0, n, 1.0),
      v = new Float64List(n),
      w = new Float64List(n),
      vv = 0.0,
      vBv = 0.0;

  for (int i = 0; i < 10; ++i) {
    AtAu(u, v, w);
    AtAu(v, u, w);
  }
  for (int i = 0; i < n; ++i) {
    vBv += u[i] * v[i];
    vv  += v[i] * v[i];
  }
  return Math.sqrt(vBv / vv);
}

void main(args) {
  int n = args.length > 0 ? int.parse(args[0]) : 100;
  print(spectralNorm(n).toStringAsFixed(9));
}
    

notes, command-line, and program output

NOTES:
64-bit Ubuntu quad core
Dart VM version: 1.22.0 (Mon Feb 13 17:43:03 2017) on "linux_x64"


Wed, 15 Feb 2017 20:21:48 GMT

MAKE:
mv spectralnorm.dart spectralnorm.dart
mv: 'spectralnorm.dart' and 'spectralnorm.dart' are the same file
/home/dunham/benchmarksgame/nanobench/makefiles/u64q.programs.Makefile:687: recipe for target 'spectralnorm.dart_run' failed
make: [spectralnorm.dart_run] Error 1 (ignored)
0.01s to complete and log all make actions

COMMAND LINE:
/usr/local/src/dart-sdk/bin/dart  spectralnorm.dart 5500

PROGRAM OUTPUT:
1.274224153