The Computer Language
Benchmarks Game

spectral-norm PHP program

source code

<? 
/* 
   The Computer Language Benchmarks Game
   http://benchmarksgame.alioth.debian.org/

   contributed by Nikita Popov
 */


function A($i, $j){
   return 1.0 / ( ( ( ($i+$j) * ($i+$j+1) ) >> 1 ) + $i + 1 );
}

function Av($n, $v, $start, $end){
   global $_tpl;
   $Av = $_tpl;
   for ($i = $start; $i < $end; ++$i) {
      $sum = 0.0;
      foreach($v as $j=>$v_j) {
         $sum += A($i,$j) * $v_j;
      }
      $Av[$i] = $sum;
   }
   return $Av;
}

function Atv($n, $v, $start, $end){
   global $_tpl;
   $Atv = $_tpl;
   for($i = $start; $i < $end; ++$i) {
      $sum = 0.0;
      foreach($v as $j=>$v_j) {
         $sum += A($j,$i) * $v_j;
      }
      $Atv[$i] = $sum;
   }
   return $Atv;
}

function AtAv($n, $v, $start, $end, $sync){

   $tmp = Av($n, $v, $start, $end);
   if ($sync) $tmp = sync($tmp);

   $tmp = Atv($n, $tmp, $start, $end);
   if ($sync) $tmp = sync($tmp);

   return $tmp;
}

function sync($tmp) {

   global $parent,$chunk_data_size,$total_data_size,$pipe,$pipes;

   if (!$parent) {
      $data = pack('d*', ...$tmp);
      safe_write($pipe, $data);
      return array_values(unpack('d*', safe_read($pipe, $total_data_size)));
   } else {
      $tmps = array();
      foreach($pipes as $pipe) {
         $tmps[] = unpack('d*', safe_read($pipe, $chunk_data_size));
      }
      $tmps[] = $tmp;
      $tmp = array_merge(...$tmps);

      $data = pack('d*', ...$tmp);
      foreach($pipes as $pipe) {
         safe_write($pipe, $data);
      }
      return $tmp;
   }
}

function safe_write($fd, $data) {
   $len = strlen($data);
   do {
      $w = fwrite($fd, $data);
      $len -= $w;
   } while($len && ($data = substr($data, $w)) !== FALSE);
}
function safe_read($fd, $len) {
   $data = '';
   while ($len > 0) {
      $d = fread($fd, $len);
      $len -= strlen($d);
      $data .= $d;
   }
   return $data;
}
function pipe() {
   return stream_socket_pair(STREAM_PF_UNIX, STREAM_SOCK_STREAM, 0);
}


$n = (int) (($argc == 2) ? $argv[1] : 1);

$procs = 1;
if (file_exists('/proc/cpuinfo')) {
   $procs = preg_match_all('/^processor\s/m', file_get_contents('/proc/cpuinfo'), $discard);
}

if ($n < $procs) {
   $procs = 1;
}

$chunk_size = (int) ($n / $procs);
$double_size = strlen(pack('d', 0.0));
$chunk_data_size = $double_size * $chunk_size;
$total_data_size = $double_size * $n;

$pipes = array();
$parent = FALSE;
for($i = 0; $i < $procs; ++$i) {
   $range_begin = $i * $chunk_size;
   if ($i < ($procs - 1)) {
      $pipe = pipe();
      $pipes[] = $pipe[0];
      $pipe = $pipe[1];
      $range_end = $range_begin + $chunk_size;
      $pid = pcntl_fork();
      if ($pid === -1) {
         die('could not fork');
      } else if ($pid) {
         continue;
      }
      break;
   } else {
      $range_end = $n;
      $parent = TRUE;
   }
}

$u = array_fill(0, $n, 1.0);
$_tpl = array_fill($range_begin, $range_end - $range_begin, 0.0);
$sync = $procs > 0;

for ($i=0; $i<10; $i++){
   $v = AtAv($n,$u,$range_begin,$range_end,$sync);
   $u = AtAv($n,$v,$range_begin,$range_end,$sync);
}

if (!$parent) {
   exit(0);
}

$childs = $procs - 1;
while ($childs--) {
   pcntl_wait($s);
}

$vBv = 0.0;
$vv = 0.0;
$i = 0;
foreach($v as $val) {
   $vBv += $u[$i]*$val;
   $vv += $val*$val;
   ++$i;
}
printf("%0.9f\n", sqrt($vBv/$vv));
    

notes, command-line, and program output

NOTES:
64-bit Ubuntu quad core
PHP 7.1.4 (cli) (built: Apr 16 2017 16:17:54) ( NTS )
Copyright (c) 1997-2017 The PHP Group
Zend Engine v3.1.0, Copyright (c) 1998-2017 Zend Technologies


Mon, 17 Apr 2017 01:34:20 GMT

COMMAND LINE:
/usr/local/src/php-7.1.4/bin/php -n  spectralnorm.php 5500

PROGRAM OUTPUT:
1.274224153