Edit File - File Manager

Edit File: benchmark.php

<?php

error_reporting(E_ALL);

/**
 * @package JAMA
 */

require_once '../Matrix.php';
require_once 'Stats.php';

/**
 * Example of use of Matrix Class, featuring magic squares.
 */
class Benchmark {
	public $stat;

/**
	 * Simple function to replicate PHP 5 behaviour
	 */
	function microtime_float() {
		list($usec, $sec) = explode(" ", microtime());

return ((float)$usec + (float)$sec);
	}	//	function microtime_float()

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;">&sigma;:<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 displayStats()

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 runEig()

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 runLU()

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 runQR()

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 runCholesky()

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 runSVD()

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;
	}	//	function run()

public function __construct() {
		$this->stat = new Base();
	}	//	function Benchmark()

}  // class Benchmark		(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>
				<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>
			</li>
		</ul>
		<?php
		break;
}

Back to File Manager