///////////////////////////////////////////////////////////////////////////////
//    Copyright (C) 2005 Intel Corp.
//
//    Subject to the terms and conditions set forth below, Intel hereby
//    grants you a nonexclusive, nontransferable license, to use,
//    reproduce and distribute the example code sequences contained
//    herein, in object code format, solely as part of your computer
//    program(s) and solely in order to allow your computer program(s) to
//    implement the multimedia instruction extensions contained in such
//    sequences solely with respect to the Intel instruction set
//    architecture.  No other license, express, implied, statutory, by
//    estoppel or otherwise, to any other intellectual property rights is
//    granted herein.
//
//    ALL INFORMATION, SAMPLES AND OTHER MATERIALS PROVIDED HEREIN
//    INCLUDING, WITHOUT LIMITATION, THE EXAMPLE CODE SEQUENCES ARE
//    PROVIDED "AS IS" WITH NO WARRANTIES, EXPRESS, IMPLIED, STATUTORY OR
//    OTHERWISE, AND INTEL SPECIFICALLY DISCLAIMS ANY IMPLIED WARRANTY OF
//    MERCHANTABILITY,  NONINFRINGEMENT OR FITNESS FOR ANY PARTICULAR
//    PURPOSE.
//
//    THE MATERIALS PROVIDED HEREIN ARE PROVIDED WITHOUT CHARGE.
//    THEREFORE, IN NO EVENT WILL INTEL BE LIABLE FOR ANY DAMAGES OF ANY
//    KIND, INCLUDING DIRECT OR INDIRECT DAMAGES, LOSS OF DATA, LOST
//    PROFITS, COST OF COVER OR SPECIAL, INCIDENTAL, CONSEQUENTIAL,
//    DAMAGES ARISING FROM THE USE OF THE MATERIALS PROVIDED HEREIN,
//    INCLUDING WITHOUT LIMITATION THE EXAMPLE CODE SEQUENCES, HOWEVER
//    CAUSED AND ON ANY THEORY OF LIABILITY.  THIS LIMITATION WILL APPLY
//    EVEN IF INTEL OR ANY AUTHORIZED AGENT OF INTEL HAS BEEN ADVISED OF
//    THE POSSIBILITY OF SUCH DAMAGE.
//
///////////////////////////////////////////////////////////////////////////////#ifndef THREADPOOL_H
///////////////////////////////////////////////////////////////////////////////
//
//  This implments a brute force method of determining the prime numbers
//  in a given range 1..n.
//
///////////////////////////////////////////////////////////////////////////////
#include <stdio.h>
#include <windows.h>
#include <time.h>
#include <math.h>
static const int gMaxThreads  = 32;
static       int gProgress    = 0,
                 gPrimesFound = 0;
int              globalPrimes[100000];


///////////////////////////////////////////////////////////////////////////////
//
//  GetCommandLineNumThreads
//
//  Grabs the number of threads from the commandline
//
///////////////////////////////////////////////////////////////////////////////
void GetCommandLineArguments(int argc, char **argv, int &start, int &end)
{
    if( argc == 3 )
    {
        start = atoi(argv[1]);
        end   = atoi(argv[2]);
    }
   else
    {
        printf("Usage:- %s <start range> <end range> \n", argv[0]);
        exit(-1);
    }
}

void ShowProgress( int val, int range )
{
    int percentDone = 0;
#pragma omp critical (one)
	{
       gProgress++;

       percentDone = (int)((float)gProgress/(float)range *200.0f + 0.5f);
	}

    if( percentDone % 10 == 0 )
        printf("\b\b\b\b%3d%%", percentDone);
}

bool TestForPrime(int val)
{
    int limit, factor = 3;

    if( val == 1 )
        return true;
    else if( val == 2 )
        return true;

    limit = (long)(sqrtf((float)val)+0.5f);
    while( (factor <= limit) && (val % factor))
        factor ++;

    return (factor > limit);
}
#                         National PARAM Supercomputing Facility

void FindPrimes(int start, int end)
{
    // start is always odd
    int range = end - start + 1; 

#pragma omp parallel for
    for( int i = start; i <= end; i += 2 )
    {
        if( TestForPrime(i) )
#pragma omp critical (two)
            globalPrimes[gPrimesFound++] = i;
        
        ShowProgress(i, range);
    }
}

int main(int argc, char **argv)
{
	int     start, end;
    clock_t before, after;

    GetCommandLineArguments(argc, argv, start, end);

    if( start <= 2 )
        globalPrimes[gPrimesFound++] = 2; // Special case

    if((start % 2) == 0 )
        start = start + 1;

    before = clock();
    FindPrimes(start, end);
    after = clock();


    printf("\n\n%8d primes found between %6d and %6d in %7.2f secs\n",
           gPrimesFound,
           start, end, (float)(after - before)/ CLOCKS_PER_SEC);
}