Centre Blaise Pascal » Bench4GPU

LENGTH splitter(LENGTH iterations,unsigned int process,

		unsigned int seed_w,unsigned int seed_z)

{

  LENGTH inside[1024],insides=0;

#pragma omp parallel for

  for (int i=0 ; i<process; i++) {

    inside[i]=MainLoopGlobal(iterations/process,seed_w+process,seed_z+process);

    printf("\tFound %lld for process %i\n",(long long)inside[i],i);

  }

  printf("\n");

  for (int i=0 ; i<process; i++) {

    insides+=inside[i];

  }

  return(insides);

} 

  LENGTH iterations=ITERATIONS,insides=0;

  insides=splitter(iterations,process,seed_w,seed_z);

//

// Estimation of Pi using Monte Carlo exploration process

// gcc -std=c99 -O3 -o Pi Pi.c -lm 

//

#include <math.h>

#include <stdio.h>

#include <stdlib.h>

#include <omp.h>

#include <limits.h>

// Marsaglia RNG very simple implementation

#define znew  ((z=36969*(z&65535)+(z>>16))<<16)

#define wnew  ((w=18000*(w&65535)+(w>>16))&65535)

#define MWC   (znew+wnew)

#define SHR3  (jsr=(jsr=(jsr=jsr^(jsr<<17))^(jsr>>13))^(jsr<<5))

#define CONG  (jcong=69069*jcong+1234567)

#define KISS  ((MWC^CONG)+SHR3)

#define MWCfp MWC * 2.328306435454494e-10f

#define KISSfp KISS * 2.328306435454494e-10f

#define ITERATIONS 1000000000

#define PROCESS 4

#ifdef LONG

#define LENGTH long long

#else

#define LENGTH int

#endif

#pragma omp declare target

LENGTH splitter(int,int,int,LENGTH);

LENGTH MainLoopGlobal(LENGTH iterations,unsigned int seed_w,unsigned int seed_z)

{

   unsigned int z=seed_z;

   unsigned int w=seed_w;

   LENGTH total=0;

   for (LENGTH i=0;i<iterations;i++) {

      float x=MWCfp ;

      float y=MWCfp ;

      // Matching test

      int inside=((x*x+y*y) < 1.0f) ? 1:0;

      total+=inside;

   }

   return(total);

}

LENGTH splitter(int process,int seed_w,int seed_z,LENGTH iterations) {

  LENGTH inside[1024],insides=0;

  int i;

#pragma omp target device(0)

#pragma omp teams num_teams(60) thread_limit(4)

#pragma omp distribute

  for (int i=0 ; i<process; i++) {

    inside[i]=MainLoopGlobal(iterations/process,seed_w+process,seed_z+process);

    //    printf("\tFound %lld for process %i\n",(long long)inside[i],i);

  }

  //printf("\n");

  for (int i=0 ; i<process; i++) {

    insides+=inside[i];

  }

  return(insides);

}

int main(int argc, char *argv[]) {

  unsigned int seed_w=10,seed_z=10,process=PROCESS;

  LENGTH iterations=ITERATIONS;

  LENGTH inside[1024],insides=0;

  if (argc > 1) {

    iterations=(LENGTH)atoll(argv[1]);

    process=atoi(argv[2]);

  }

  else {

    printf("\n\tPi : Estimate Pi with Monte Carlo exploration\n\n");

    printf("\t\t#1 : number of iterations (default 1 billion)\n");

    printf("\t\t#2 : number of process (default 4)\n\n");

  }

  printf ("\n\tInformation about architecture:\n\n");

  printf ("\tSizeof int = %lld bytes.\n", (long long)sizeof(int));

  printf ("\tSizeof long = %lld bytes.\n", (long long)sizeof(long));

  printf ("\tSizeof long long = %lld bytes.\n\n", (long long)sizeof(long long));

  printf ("\tMax int = %u\n", INT_MAX);

  printf ("\tMax long = %ld\n", LONG_MAX);

  printf ("\tMax long long = %lld\n\n", LLONG_MAX);

  insides=splitter(process,seed_w,seed_z,iterations);

  float pi=4.*(float)insides/(float)iterations;

  printf("\tPi=%f with error %f and %lld iterations\n\n",pi,

         fabs(pi-4*atan(1))/pi,(long long)iterations);

}

#!/bin/bash

EXE=Pi_OpenMP_LONG

ITERATIONS=10000000000

MYTIME=/usr/bin/time

REPEAT=10

PROCESS=480

[ ! $1 == '' ] && EXE=$1

[ ! $2 == '' ] && ITERATIONS=$2

[ ! $3 == '' ] && PROCESS=$3

LOGFILE=${EXE}_${HOSTNAME}_${ITERATIONS}.log

> $LOGFILE

p=1

while [ $p -le $PROCESS ]

do

    export OMP_NUM_THREADS=$p

    echo -e "Process $p" >> $LOGFILE

    echo -ne "Start $EXE with $ITERATIONS and $p : "

    i=1

    while [ $i -le $REPEAT ]

    do 

        echo -ne "$i "

        $MYTIME ./$EXE $ITERATIONS $p >> $LOGFILE 2>&1 

        i=$(($i+1))

    done

    echo

    p=$(($p+1))

done

EXECUTABLE=Pi_OpenMP_LONG Pi_OpenMP_INT

SOURCE=Pi_OpenMP.c

CC=icpc

CFLAGS=-Wall -O3 -openmp -g

LIBRARY=-lm

all: $(EXECUTABLE)

Pi_OpenMP_LONG: $(SOURCE)

	$(CC) $(CFLAGS) $(DIRECTIVES) -DLONG -o $@ $< $(LIBRARY)

Pi_OpenMP_INT: $(SOURCE)

	$(CC) $(CFLAGS) $(DIRECTIVES) -DINTEGER -o $@ $< $(LIBRARY)

.PHONY: clean check mrproper

mrproper: 

	rm -rf $(EXECUTABLE)

	find . -name "*~" -exec rm {} \;

clean:

	find . -name "*~" -exec rm {} \;

check: $(EXECUTABLE)

	@echo "To be Defined"