#!/bin/bash

ecode=0

FAILED=""

if test -z "${OPTIONS}";then

OPTIONS="sha1 sha256 sha384 sha256-new sha1-one"

fi

echo "Running for: $OPTIONS"

TIME=ns

function plot()

{

	hash=$1

	R --no-save <<_EOF_

#plot

require(lattice)

z=read.csv("out-$hash.txt", header=TRUE);

options(scipen=999)

m <- matrix(ncol=3,nrow=256)

png(filename = "$hash-density.png",width=3072,height=4096,units="px",bg="white");

d <- density(unlist(z[1,]));

plot(d, ylab="Delta", xlab="${TIME}", col="red", main=sprintf("estimated PDF for %d", 1));

for (row in c(4,72,253,256)){

  lines(density(unlist(z[row,])), col = sample(colours(), 5));

}

dev.off();

colnames(m) <- c("Delta","TimeMedian","TimeAverage")

for (row in 1:256){

  v=tail(z[row,], length(z[row,])-1)

  m[row,][1] = unlist(z[row,][1]);

  m[row,][2] = median(unlist(v));

  m[row,][3] = mean(unlist(v));

}

png(filename = "$hash-timings-med.png",width=1024,height=1024,units="px",bg="white");

plot(m[,1],m[,2],xlab="Delta",ylab="Median timings (${TIME})");

dev.off();

png(filename = "$hash-timings-ave.png",width=1024,height=1024,units="px",bg="white");

plot(m[,1],m[,3],xlab="Delta",ylab="Average timings (${TIME})");

dev.off();

estatus=0;

for (row in 2:256){

  v=ks.test(unlist(z[1,]),unlist(z[row,]));

  p = 0.05/255

  if (v[2] < p) {

    print (row);

    print (v[2]);

    estatus=1;

  }

}

q(status=estatus);

_EOF_

if test $? != 0;then

	echo "*** Failed test for $hash"

	FAILED="${FAILED} $hash"

	ecode=1

fi

}

for hash in ${OPTIONS};do

	if ! test -f out-$hash.txt;then

		./mini-record-timing $hash

	fi

	plot $hash

done

echo "Failed: ${FAILED}"

exit $ecode