main.c

Go to the documentation of this file.

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <unistd.h>
#include <stdint.h>
#include <pthread.h>

#define NUM_ACCEL 4
#define RUN 32
#define OPS_PER_ACCEL RUN/NUM_ACCEL

#include "black_scholes.h"
#include "fixedptc.h"
#include "fixedptc.c"
#include "black_scholes.c"
#include "seed.h"

fixedpt asset_path_test ( int seed );

//pthread_mutex_t mutex1 = PTHREAD_MUTEX_INITIALIZER;
//int dummySeed = 123;

struct thread_data
{
   int  startidx;
   int  maxidx;
};

void *black_scholes(void* threadarg)
{

   int i;
   fixedpt sum = 0;
   struct thread_data* arg = (struct thread_data*) threadarg;
   int startidx = arg->startidx;
   int maxidx = arg->maxidx;
   fixedpt u;
   for (i = startidx; i < maxidx; i++)
   {
      fixedpt seed = seeds[i];
      u = asset_path_test (seed);
      //printf("u = %d\n", u >> (FIXEDPT_BITS - FIXEDPT_WBITS));
      result[i] = u;
      sum += (u >> FIXEDPT_BITS - FIXEDPT_WBITS);
   }

   pthread_exit((void*)sum);
}

int main ( void )
{

   legup_start_counter(0);
   int run = RUN;
   int i;
   fixedpt sum = 0;
   //create the thread variables
   fixedpt result[NUM_ACCEL] = {0};
   pthread_t threads[NUM_ACCEL];
   struct thread_data data[NUM_ACCEL];

   for (i = 0; i < NUM_ACCEL; i++)
   {
      //initialize structs to pass into accels
      data[i].startidx = i * OPS_PER_ACCEL;
      data[i].maxidx = (i + 1) * OPS_PER_ACCEL;
   }

   //launch threads
   for (i = 0; i < NUM_ACCEL; i++)
   {
      pthread_create(&threads[i], NULL, black_scholes, (void *)&data[i]);
   }

   //join the threads
   for (i = 0; i < NUM_ACCEL; i++)
   {
      pthread_join(threads[i], (void**)&result[i]);
   }

   for (i = 0; i < NUM_ACCEL; i++)
   {
      sum += result[i];
   }

   int perf_counter = legup_stop_counter(0);
   printf("perf_counter = %d\n", perf_counter);
   printf("sum = %d\n", sum); // the golden result will be the SUM of the prices
   if (sum == 10752)
   {
      printf("RESULT: PASS\n");
   }
   else
   {
      printf("RESULT: FAIL\n");
   }

   return (int)sum;
}

fixedpt asset_path_test (int seed )
{

   int n = 100;
   fixedpt mu, s0, sigma, t1;
   fixedpt s;
   //int holdSeed = seed;

   s0 = 13107200; // fixedpt_rconst(200.0);
//        printf("%d\n", s0>> (FIXEDPT_BITS - FIXEDPT_WBITS));

   mu = 16384; // fixedpt_rconst(0.25);
//        printf("%d\n", mu>> (FIXEDPT_BITS - FIXEDPT_WBITS));

   sigma = 4391; //fixedpt_rconst(0.067);
//        printf("%d\n", sigma>> (FIXEDPT_BITS - FIXEDPT_WBITS));

   t1 = 131072; // fixedpt_rconst(2.0);
//        printf("%d\n", t1>> (FIXEDPT_BITS - FIXEDPT_WBITS));

   s = asset_path_fixed_simplified ( s0, mu, sigma, t1, n, &seed);

   return s;
}
/******************************************************************************/