File Source for benchmark.php

<?php
/**
* @package JAMA
*/
require_once "../Matrix.php";
require_once "Stats.php";
/**
* Example of use of Matrix Class, featuring magic squares.
*/
class Benchmark {
var $stat;
/**
* Simple function to replicate PHP 5 behaviour
*/
function microtime_float() {
list($usec, $sec) = explode(" ", microtime());
return ((float)$usec + (float)$sec);
}
function Benchmark() {
$this->stat = new Base();
}
function displayStats( $times = null ) {
$this->stat->setData($times);
$stats = $this->stat->calcFull();
echo '<table style="margin-left:32px;">';
echo '<tr><td style="text-align:right;">n:<td style="text-align:right;">' . $stats['count'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">Mean:<td style="text-align:right;">' . $stats['mean'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">Min.:<td style="text-align:right;">' . $stats['min'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">Max.:<td style="text-align:right;">' . $stats['max'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">σ:<td style="text-align:right;">' . $stats['stdev'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">Variance:<td style="text-align:right;">' . $stats['variance'] . ' </td></tr>';
echo '<tr><td style="text-align:right;">Range:<td style="text-align:right;">' . $stats['range'] . ' </td></tr>';
echo '</table>';
return $stats;
}
function runEig($n = 4, $t = 100) {
$times = array();
for( $i = 0; $i < $t; $i++ ) {
$M = Matrix::random($n, $n);
$start_time = $this->microtime_float();
$E = new EigenvalueDecomposition($M);
$stop_time = $this->microtime_float();
$times[] = $stop_time - $start_time;
}
return $times;
}
function runLU($n = 4, $t = 100) {
$times = array();
for( $i = 0; $i < $t; $i++ ) {
$M = Matrix::random($n, $n);
$start_time = $this->microtime_float();
$E = new LUDecomposition($M);
$stop_time = $this->microtime_float();
$times[] = $stop_time - $start_time;
}
return $times;
}
function runQR($n = 4, $t = 100) {
$times = array();
for( $i = 0; $i < $t; $i++ ) {
$M = Matrix::random($n, $n);
$start_time = $this->microtime_float();
$E = new QRDecomposition($M);
$stop_time = $this->microtime_float();
$times[] = $stop_time - $start_time;
}
return $times;
}
function runCholesky($n = 4, $t = 100) {
$times = array();
for( $i = 0; $i < $t; $i++ ) {
$M = Matrix::random($n, $n);
$start_time = $this->microtime_float();
$E = new CholeskyDecomposition($M);
$stop_time = $this->microtime_float();
$times[] = $stop_time - $start_time;
}
return $times;
}
function runSVD($n = 4, $t = 100) {
$times = array();
for( $i = 0; $i < $t; $i++ ) {
$M = Matrix::random($n, $n);
$start_time = $this->microtime_float();
$E = new SingularValueDecomposition($M);
$stop_time = $this->microtime_float();
$times[] = $stop_time - $start_time;
}
return $times;
}
function run() {
$n = 8;
$t = 16;
$sum = 0;
echo "Cholesky decomposition: $t random {$n}x{$n} matrices ";
$r = $this->displayStats($this->runCholesky($n, $t));
$sum += $r['mean'] * $n;
echo '<hr />';
echo "Eigenvalue decomposition: $t random {$n}x{$n} matrices ";
$r = $this->displayStats($this->runEig($n, $t));
$sum += $r['mean'] * $n;
echo '<hr />';
echo "LU decomposition: $t random {$n}x{$n} matrices ";
$r = $this->displayStats($this->runLU($n, $t));
$sum += $r['mean'] * $n;
echo '<hr />';
echo "QR decomposition: $t random {$n}x{$n} matrices ";
$r = $this->displayStats($this->runQR($n, $t));
$sum += $r['mean'] * $n;
echo '<hr />';
echo "Singular Value decomposition: $t random {$n}x{$n} matrices ";
$r = $this->displayStats($this->runSVD($n, $t));
$sum += $r['mean'] * $n;
return $sum;
}
} // end MagicSquareExample
$benchmark = new Benchmark();
switch($_REQUEST['decomposition']) {
case 'cholesky':
$m = array();
for($i = 2; $i <= 8; $i *= 2) {
$t = 32/$i;
echo "Cholesky decomposition: $t random {$i}x{$i} matrices ";
$s = $benchmark->displayStats($benchmark->runCholesky($i, $t));
$m[$i] = $s['mean'];
echo " ";
}
echo '<pre>';
foreach($m as $x => $y)
echo "$x\t" . 1000*$y . "\n";
echo '</pre>';
break;
case 'eigenvalue':
$m = array();
for($i = 2; $i <= 8; $i *= 2) {
$t = 32/$i;
echo "Eigenvalue decomposition: $t random {$i}x{$i} matrices ";
$s = $benchmark->displayStats($benchmark->runEig($i, $t));
$m[$i] = $s['mean'];
echo " ";
}
echo '<pre>';
foreach($m as $x => $y)
echo "$x\t" . 1000*$y . "\n";
echo '</pre>';
break;
case 'lu':
$m = array();
for($i = 2; $i <= 8; $i *= 2) {
$t = 32/$i;
echo "LU decomposition: $t random {$i}x{$i} matrices ";
$s = $benchmark->displayStats($benchmark->runLU($i, $t));
$m[$i] = $s['mean'];
echo " ";
}
echo '<pre>';
foreach($m as $x => $y)
echo "$x\t" . 1000*$y . "\n";
echo '</pre>';
break;
case 'qr':
$m = array();
for($i = 2; $i <= 8; $i *= 2) {
$t = 32/$i;
echo "QR decomposition: $t random {$i}x{$i} matrices ";
$s = $benchmark->displayStats($benchmark->runQR($i, $t));
$m[$i] = $s['mean'];
echo " ";
}
echo '<pre>';
foreach($m as $x => $y)
echo "$x\t" . 1000*$y . "\n";
echo '</pre>';
break;
case 'svd':
$m = array();
for($i = 2; $i <= 8; $i *= 2) {
$t = 32/$i;
echo "Singular value decomposition: $t random {$i}x{$i} matrices ";
$s = $benchmark->displayStats($benchmark->runSVD($i, $t));
$m[$i] = $s['mean'];
echo " ";
}
echo '<pre>';
foreach($m as $x => $y)
echo "$x\t" . 1000*$y . "\n";
echo '</pre>';
break;
case 'all':
$s = $benchmark->run();
print(" Total: {$s}s ");
break;
default:
?>
<ul>
<li><a href="benchmark.php?decomposition=all">Complete Benchmark</a></li>
<ul>
<li><a href="benchmark.php?decomposition=cholesky">Cholesky</a></li>
<li><a href="benchmark.php?decomposition=eigenvalue">Eigenvalue</a></li>
<li><a href="benchmark.php?decomposition=lu">LU</a></li>
<li><a href="benchmark.php?decomposition=qr">QR</a></li>
<li><a href="benchmark.php?decomposition=svd">Singular Value</a></li>
</ul>
</ul>
<?php
break;
}
?>

Source for file benchmark.php