static char rcsid[] = "$Id: for009.c,v 1.12 2002/04/14 13:58:19 a-hasega Exp $";
/* 
 * $RWC_Release: Omni-1.6 $
 * $RWC_Copyright:
 *  Omni Compiler Software Version 1.5-1.6
 *  Copyright (C) 2002 PC Cluster Consortium
 *  
 *  This software is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License version
 *  2.1 published by the Free Software Foundation.
 *  
 *  Omni Compiler Software Version 1.0-1.4
 *  Copyright (C) 1999, 2000, 2001.
 *   Tsukuba Research Center, Real World Computing Partnership, Japan.
 *  
 *  Please check the Copyright and License information in the files named
 *  COPYRIGHT and LICENSE under the top  directory of the Omni Compiler
 *  Software release kit.
 *  
 *  
 *  $
 */
/* for 009:
 * check no implicit barrier at the end of for loop with nowait
 */

#include "omniConfig.h"
#include <omp.h>
#include "omni.h"


int	thds;
int	*buf;


void
barrier (int n)
{
  static int cnt[2] = { 0, 0 };
  static int cur = 0;

  int	c;

  #pragma omp critical
  {
    c = cur;
    if (cnt[c]+1 == n) {
      cur = (cur + 1) % 2;
      cnt[c] = 0;
    } else {
      cnt[c] += 1;
    }
  }

  while (cnt[c] != 0) {
    #pragma omp flush
  }
}


void
clear ()
{
  int lp;
  
  for (lp=0; lp<thds; lp++) {
    buf[lp] = -1;
  }
}


main ()
{
  int	lp, lp0;

  int	errors = 0;


  thds = omp_get_max_threads ();
  if (thds == 1) {
    printf ("should be run this program on multi threads.\n");
    exit (0);
  }
  buf = (int *) malloc (sizeof (int) * thds);
  if (buf == NULL) {
    printf ("can not allocate memory.\n");
    exit (1);
  }
  omp_set_dynamic (0);


  clear ();
  #pragma omp parallel
  {
    int	id = omp_get_thread_num ();
    int	i;

    #pragma omp for schedule (static,1) nowait
    for (lp=0; lp<thds; lp++) {
      barrier (thds);
      if (id == 0) {
	for (i=1; i<thds; i++) {
	  while (buf[i] != i) {
	    #pragma omp flush
	  }
	}
	buf[id] = id;
	#pragma omp flush
      }
    }

    if (id != 0) {
      #pragma omp flush
      if (buf[0] != -1) {
        #pragma omp critical
	{
	  ERROR (errors);
	}
      }
      buf[id] = id;
      #pragma omp flush
    }
  }

  clear ();
  #pragma omp parallel
  {
    int	id = omp_get_thread_num ();
    int	i;

    #pragma omp for schedule (dynamic,1) nowait
    for (lp=0; lp<thds; lp++) {
      barrier (thds);
      if (id == 0) {
	for (i=1; i<thds; i++) {
	  while (buf[i] != i) {
	    #pragma omp flush
	  }
	}
	buf[id] = id;
	#pragma omp flush
      }
    }

    if (id != 0) {
      #pragma omp flush
      if (buf[0] != -1) {
        #pragma omp critical
	{
	  ERROR (errors);
	}
      }
      buf[id] = id;
      #pragma omp flush
    }
  }


  lp0 = -1;
  clear ();
  #pragma omp parallel
  {
    int	id = omp_get_thread_num ();
    int	i;

    #pragma omp for schedule (guided,1) nowait
    for (lp=0; lp<thds*4; lp++) {
      if (lp == 0) {
	lp0 = id;
        #pragma omp flush
	for (i=0; i<thds; i++) {
	  if (id == i) {
	    continue;
	  }
	  while (buf[i] != i) {
	    #pragma omp flush
	  }
	}
	buf[id] = id;
	#pragma omp flush
      }
    }

    while (lp0 == -1) {
      #pragma omp flush
    }
    if (id != lp0) {
      #pragma omp flush
      if (buf[lp0] != -1) {
        #pragma omp critical
	{
	  ERROR (errors);
	}
      }
      buf[id] = id;
      #pragma omp flush
    }
  }


  if (errors == 0) {
    printf ("for 009 : SUCCESS\n");
    return 0;
  } else {
    printf ("for 009 : FAILED\n");
    return 1;
  }
}