/*
You need to make sure

NMAX is defined to 1000

and 

LARGE is defined

*/

/* you should invoke cleanupCache (or something like it) 
   before the start of each time measurement start=ts_get_ll() in 
   cart2pol.c */ 

void cleanUpCache(void) {
  int i,j, N;

  // clean up cache //
  for(i=0;i<NMAX;i++) {
    for(j=0;j<NMAX;j++) {
      refSetRCART((void*)array1, i, j, N, 0);
      refSetICART((void*)array1, i, j, N, 0);
    }
  }
}

void initTestCases_0(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=10.0;
    test[i].theta=asin(0.1);
    test[i].N=10;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_1(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=10.0;
    test[i].theta=asin(0.1);
    test[i].N=700;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_2(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=10.0;
    test[i].theta=asin(0.1);
    test[i].N=1000;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}




void initTestCases_3(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.0123);
    test[i].N=10;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_4(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.0123);
    test[i].N=700;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_5(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.0123);
    test[i].N=1000;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}




void initTestCases_6(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.2);
    test[i].N=10;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_7(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.2);
    test[i].N=700;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}

void initTestCases_8(void) {
  int i,j, N;

  printf("Initialization can take some time on the XUP.... \n");

  for(i=0;i<2;i++) {
    test[i].R=100.0;
    test[i].theta=asin(0.2);
    test[i].N=1000;
    test[i].array=(void*)array0;
  }

  N=test[0].N;

  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) {
      double dx=((double)i)/N;
      double dy=((double)j)/N;
      double r=sqrt(dx*dx+dy*dy);
 
      refSetRCART((void*)array0, i, j, N, doubleToSHORT16(SHORT16ToDouble(0x7fff)*cos(r*3*3.14159)));
      refSetICART((void*)array0, i, j, N, -doubleToSHORT16(SHORT16ToDouble(0x7fff)*sin(r*3*3.14159)));
    }
  }
}