source: trunk/platforms/caba-tsar-virtual_dspin_network-v0/soft/main.c @ 3

#include "system.h"
#include "stdio.h"
#include "stdlib.h"
#include "matrice.h"

#include "../segmentation.h"

#define NPROCS 16
#define SIZE 500 //8000
#define SORT_TYPE 1

volatile int lock=0;
volatile int compteur=NPROCS;
volatile int nprocs=NPROCS;

volatile int write_boost[16*120];

  unsigned int SortArr0[16*(SIZE+200)];

void SORT(int *base, int n, int size, int (*compar) (),int type);
void QSORT(int *base, int n, int size, int (*compar) ());
void simple_sort(int *base, int n, int size, int (*compar) ());
void insertion_sort(int *base, int n, int size, int (*compar) ());
void bubble_sort(int *base, int n, int size, int (*compar) ());
int compare(int *n1, int *n2)
{
  return (*((unsigned int *) n1) - *((unsigned int *) n2));
}

int main()
{
  register int p; 

  p=procnum();

  puts("Hello from processor ");
  putchar(p+'0');
  putchar('\n');


  int i;
  int j;
  int temp;

  if(p+1 <= nprocs){
//      for(i=1 ; i <= 8 ; i++ ){
        puts("Memory copy \n");
        memcpy(SortArr0+p*(SIZE+200), gQSortNum0 + ( (4*(p+1)*i) % 32) ,SIZE);

        puts("Sort... \n");
        SORT((int *) (SortArr0+p*(SIZE+200)), (int) SIZE, sizeof(unsigned int), compare,SORT_TYPE);
        for (j = 1; j < SIZE; j++)
        {
                if (SortArr0[p*(SIZE+200)+j] < SortArr0[p*(SIZE+200) + j - 1])
                {
                        puts("ucbqsort: failed\n");
                        while(1);
                }

        }
        puts("ucbqsort: success\n");
//      }
  }
  

  if(p+1 <= nprocs){
        lock_lock(&lock);
        compteur--;
        if(!compteur){
                putchar('\n');
                puti(temp);
                putchar('\n');
        }
        lock_unlock(&lock);
  }

  while(1);
}

static int (*qcmp) ();
static int qsz;
static int thresh;
static int mthresh;
static void qst(int *, int *);
void QSORT(base, n, size, compar)
     int *base;
     int n;
     int size;
     int (*compar) ();
{
  register int c, *i, *j, *lo, *hi;
  int *min, *max;
  if (n <= 1)
    return;
  qsz = size;
  qcmp = compar;
  thresh = 4;
  mthresh = 6;
  max = base + n;
  if (n >= 4)
    {
      qst(base, max);
      hi = base + thresh;
    } 
  else
    {
      hi = max;
    }

  puts("Quick sort done... \n");

  for (j = lo = base; (lo ++) < hi;)
    if ((*qcmp) (j, lo) > 0)
      j = lo;
  if (j != base)
    {
      for (i = base, hi = base + 1; i < hi;)
        {
          c = *j;
          *j++ = *i;
          *i++ = c;
        }

    }
  for (min = base; (hi = min ++) < max;)
    {
      while ((*qcmp) (hi --, min) > 0)
        ;
      if ((hi ++) != min)
        {
          for (lo = min + 1; --lo >= min;)
            {
              c = *lo;
              for (i = j = lo; (j --) >= hi; i = j)
                *i = *j;
              *i = c;
            }

        }
    }

}


void simple_sort(base, n, size, compar)
     int *base;
     int n;
     int size;
     int (*compar) ();
{
  register int c, *i, *j, *lo, p, *k;
  int *max;
  k=write_boost;
  if (n <= 1)
    return;
  qsz = size;
  qcmp = compar;
  max = base + n;
  puts("Selection Sort\n");
  p=procnum();
//  k=4;
  for(j = base; j < max ; j++)
  {
    putchar('-'); // added for debug
    
    lo=j;
    c=(*lo);

    for(i = j; i < max ; i++)
    {
      *(k+(p * 70)+20)=p; 

      asm("     nop             \n"); 

      if( c-(*i) > 0)
      {

        lo=i;
        c=(*lo);

      }

      asm("     nop             \n"); 
    }

    *lo = *j;
    *j = c;

    putchar('+');

  }

  
}

void insertion_sort(base, n, size, compar)
     int *base;
     int n;
     int size;
     int (*compar) ();
{
  register int c, *i, *j,p;
  int *max;
  if (n <= 1)
    return;
  qsz = size;
  qcmp = compar;
  max = base + n;
  puts("Insertion Sort\n");
  p=procnum();
  for(j = base; j < max-1 ; j++){
    for(i = j; i < max-1 ; i++){
      if( (*i)-(*(i+1)) > 0){
//      if( (*qcmp) (i,i+1) > 0){
        c = *i;
        *i = *(i+1);
        *(i+1)=c;
      }
    }
  }
}


void bubble_sort(base, n, size, compar)
     int *base;
     int n;
     int size;
     int (*compar) ();
{
  register int c, *i;
  register int b_swap=1;
  int *max;
  if (n <= 1)
    return;
  qsz = size;
  qcmp = compar;
  max = base + n;
  puts("Bubble Sort\n");
  while(b_swap){
    b_swap=0;
    for(i = base; i < max-1 ; i++){
      if( (*i)-(*(i+1)) > 0){
//      if( (*qcmp) (i,i+1) > 0){
        b_swap=1;
        c = *i;
        *i = *(i+1);
        *(i+1)=c;
      }
    }
  }
}

static
void qst(base, max)
     int *base, *max;
{
  register int c, *i, *j, *jj;
  register int ii;
  int *mid, *tmp;
  int lo, hi;
  lo = max - base;
  do
    {
      mid = i = base + 1 * ((lo / qsz) >> 1);
      if (lo >= mthresh)
        {
          j = ((*qcmp) ((jj = base), i) > 0 ? jj : i);
          if ((*qcmp) (j, (tmp = max - 1)) > 0)
            {
              j = (j == jj ? i : jj);
              if ((*qcmp) (j, tmp) < 0)
                j = tmp;
            }

          if (j != i)
            {
              ii = 1;
              do
                {
                  c = *i;
                  *i++ = *j;
                  *j++ = c;
                }

              while (--ii);
            }

        }
      for (i = base, j = max - 1;;)
        {
          while (i < mid && (*qcmp) (i, mid) <= 0)
            i ++;
          while (j > mid)
            {
              if ((*qcmp) (mid, j) <= 0)
                {
                  j --;
                  continue;
                }

              tmp = i + 1;
              if (i == mid)
                {
                  mid = jj = j;
                } 
              else
                {
                  jj = j;
                  j --;
                }

              goto swap;
            }

          if (i == mid)
            {
              break;
            } 
          else
            {
              jj = mid;
              tmp = mid = i;
              j --;
            }

        swap:
          ii = 1;
          do
            {
              c = *i;
              *i++ = *jj;
              *jj++ = c;
            }

          while (--ii);
          i = tmp;
        }

      i = (j = mid) + 1;
      if ((lo = j - base) <= (hi = max - i))
        {
          if (lo >= thresh)
            qst(base, j);
          base = i;
          lo = hi;
        } 
      else
        {
          if (hi >= thresh)
            qst(i, max);
          max = j;
        }

    }
  while (lo >= thresh);
}

void SORT(int *base, int n, int size, int (*compar) (),int type){
  switch(type){
  case 0:
    QSORT(base, n, size, compar);
    break;
  case 1:
    simple_sort(base, n, size, compar);
    break;
  case 2:
    insertion_sort(base, n, size, compar);
    break;
  case 3:
    bubble_sort(base, n, size, compar);
    break;
  default:
    break;
  }

}