/* paraTestJob - A good test job to run on Parasol.  Can be configured to take a long time or crash. */
#include "common.h"
#include "linefile.h"
#include "hash.h"
#include "options.h"

void usage()
/* Explain usage and exit. */
{
errAbort(
  "paraTestJob - A good test job to run on Parasol.  Can be configured to take a long time or crash\n"
  "usage:\n"
  "   paraTestJob count\n"
  "Run a relatively time consuming algorithm count times.  This\n"
  "algorithm takes about 1/10 per second each time\n"
  "options:\n"
  "   -crash - Try to write to NULL when done.\n"
  "   -err   - Return -1 error code when done\n"
  "   -output=file - make some output in file as well\n"
  "   -heavy=n - make output heavy: n extra 'lumberjack lines\n"
  "   -input = file - make it read in a file too\n"
  "   -sleep=n Sleep for N seconds\n"
  );
}

double cumErr = 0.0;

void noteResults(double a, double b)
/* Don't really do anything. */
{
double diff = a - b;
if (diff < 0) diff = -diff;
cumErr += diff;
}

void compute()
/* Do a relatively time consuming calculation.
 * This should last about 1/10 second. */
{
int i;
double res;

for (i=0; i<1060000; ++i)
    {
    res = sqrt(i);
    noteResults(i, res*res);
    }
}

void paraTestJob(char *countString)
/* paraTestJob - A good test job to run on Parasol.  Can be configured to take a long time or crash. */
{
int i, count = atoi(countString);
char *outName = optionVal("output", NULL);
int heavy = optionInt("heavy", 0);
FILE *f = NULL;
if (outName != NULL)
    f = mustOpen(outName, "w");
for (i=0; i<count; ++i)
    {
    compute();
    if (f != NULL)
        {
	int j;
	fprintf(f, "Computation number %d of %d.  Error so far is %f\n",
		i+1, count, cumErr);
	for (j=0; j<heavy; ++j)
	    {
	    fprintf(f, "  I'm a lumberjack and I'm ok.  I work all night and I sleep all day\n");
	    }
	}
    }
printf("Cumulative error %f\n", cumErr);
if (optionExists("sleep"))
    sleep(optionInt("sleep", 1));
if (optionExists("err"))
    errAbort("Aborting with error");
if (optionExists("crash"))
    {
    char *s = NULL;
    s[0] = 0;
    }
}

int main(int argc, char *argv[])
/* Process command line. */
{
optionHash(&argc, argv);
if (argc != 2)
    usage();
paraTestJob(argv[1]);
return 0;
}