static char rcsid[] = "$Id: atomic003.c,v 1.6 2002/04/03 10:38:33 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.
 *  
 *  
 *  $
 */
/* atomic 003:
 * 異るatomic間での動作確認
 */

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


#ifdef __OMNI_SCASH__
#define LOOPNUM		(thds * 5)
#else
#define LOOPNUM		(thds * 1000)
#endif
#define INT_BIT		(sizeof(int) * 8)

int		errors = 0;
int		thds;
int		atom_incr, atom_incr2, atom_decr, atom_decr2;
int		atom_plus, atom_mul, atom_minus, atom_div;
int		atom_and, atom_or, atom_xor;
unsigned int	atom_lshift, atom_rshift;


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
  }
}


int
sameas (int v)
{
  return v;
}


void
clear ()
{
  int	i;

  atom_incr = 0;
  atom_incr2 = 0;
  atom_decr = 0;
  atom_decr2 = 0;

  atom_plus = 0;
  atom_mul = 1;
  atom_minus = 0;
  atom_div = 1 << (INT_BIT - 2);

  for (i=atom_and=0; i<INT_BIT; i++)
    atom_and |= (1<<i);
  for (i=atom_xor=0; i<INT_BIT; i++)
    atom_xor |= (1<<i);
  atom_or = 0;
  atom_lshift = 1;
  atom_rshift = 1 << (INT_BIT - 2);
}


int
check ()
{
  int	i, tmp;

  int	err = 0;


  if (atom_incr != LOOPNUM) {
    err ++;
  }
  if (atom_incr2 != LOOPNUM) {
    err ++;
  }
  if (atom_decr != -LOOPNUM) {
    err ++;
  }
  if (atom_decr2 != -LOOPNUM) {
    err ++;
  }

  if (atom_plus != LOOPNUM*2) {
    err ++;
  }
  if (atom_minus != -LOOPNUM*3) {
    err ++;
  }

  if (atom_mul != 1) {
    err ++;
  }
  if (atom_div != (1 << (INT_BIT - 2))) {
    err ++;
  }
  for (i=tmp=0; i<INT_BIT; i++)
    tmp |= (1<<i);
  for (i=0; i<LOOPNUM; i++)
    tmp &= ~(1<<(i%(INT_BIT-1)));
  if (atom_and != tmp) {
    err ++;
  }
  for (i=tmp=0; i<INT_BIT; i++)
    tmp |= (1<<i);
  for (i=0; i<LOOPNUM; i++) {
    tmp ^= (1<<(i%(INT_BIT-1)));
  }
  if (atom_xor != tmp) {
    err ++;
  }
  for(i=tmp=0; i<LOOPNUM; i++)
    tmp |= (1<<(i%(INT_BIT-1)));  
  if (atom_or != tmp) {
    err ++;
  }
  if (atom_lshift != 1) {
    err ++;
  }
  if (atom_rshift != (1 << (INT_BIT - 2))) {
    err ++;
  }

  return err;
}


void
func_atomic (int id)
{
  int i;


  switch (id) {
  case 0:
    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_incr ++;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      ++ atom_incr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_decr --;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      -- atom_decr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_plus += sameas(1) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_minus -= sameas(2) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_mul *= sameas(3) - 1;
      #pragma omp atomic	
      atom_mul /= 1 - sameas(3);
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_div /= 1 + sameas(3);
        #pragma omp atomic	
        atom_div *= 3 + sameas(1);
      }
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_and &= ~(1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_xor ^= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_or |= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_lshift <<= 1;

        #pragma omp atomic	
        atom_lshift >>= 1;
      }
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_rshift >>= 2;

        #pragma omp atomic	
        atom_rshift <<= 2;
      }
    }
    break;

  case 1:
    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_incr ++;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      ++ atom_incr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_decr --;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      -- atom_decr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_plus += sameas(1) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_minus -= sameas(2) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_mul *= sameas(3) - 1;
      #pragma omp atomic	
      atom_mul /= 1 - sameas(3);
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_div /= 1 + sameas(3);
        #pragma omp atomic	
        atom_div *= 3 + sameas(1);
      }
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_and &= ~(1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_xor ^= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_or |= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_lshift <<= 1;

      #pragma omp atomic	
      atom_lshift >>= 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_rshift >>= 1;

      #pragma omp atomic	
      atom_rshift <<= 1;
    }
    break;

  case 2:
    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_incr ++;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      ++ atom_incr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_decr --;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      -- atom_decr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_plus += sameas(1) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_minus -= sameas(2) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_mul *= sameas(3) - 1;
      #pragma omp atomic	
      atom_mul /= 1 - sameas(3);
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_div /= 1 + sameas(3);
        #pragma omp atomic	
        atom_div *= 3 + sameas(1);
      }
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_and &= ~(1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_xor ^= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_or |= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_lshift <<= 1;

        #pragma omp atomic	
        atom_lshift >>= 1;
      }
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      if (omp_get_thread_num () < (INT_BIT - 2)) {
        #pragma omp atomic	
        atom_rshift >>= 1;

        #pragma omp atomic	
        atom_rshift <<= 1;
      }
    }
    break;

  case 3:
    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_incr ++;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      ++ atom_incr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      atom_decr --;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic
      -- atom_decr2;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_plus += sameas(1) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_minus -= sameas(2) + 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_mul *= sameas(3) - 1;
      #pragma omp atomic	
      atom_mul /= 1 - sameas(3);
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_div /= 1 + sameas(3);
      #pragma omp atomic	
      atom_div *= 3 + sameas(1);
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_and &= ~(1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_xor ^= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
      #pragma omp atomic	
      atom_or |= (1<<(i%(INT_BIT-1)));
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_lshift <<= 1;

      #pragma omp atomic	
      atom_lshift >>= 1;
    }

    barrier (thds);
    for (i=0; i<LOOPNUM/thds; i++) {
      #pragma omp atomic	
      atom_rshift >>= 2;

      #pragma omp atomic	
      atom_rshift <<= 2;
    }
    break;

  default:
    break;
  }
}


main ()
{
  thds = omp_get_max_threads ();
  if (thds == 1) {
    printf ("should be run this program on multi threads.\n");
    exit (0);
  }
  if (4 < thds) {
    thds = 4;
    omp_set_num_threads (4);
  }

  omp_set_dynamic (0);


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

    switch (id) {
    case 0:
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_incr ++;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        ++ atom_incr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_decr --;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        -- atom_decr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_plus += sameas(1) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_minus -= sameas(2) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_mul *= sameas(3) - 1;  
        #pragma omp atomic	  
        atom_mul /= 1 - sameas(3);  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_div /= 1 + sameas(3);  
        #pragma omp atomic	  
        atom_div *= 3 + sameas(1);  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_and &= ~(1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_xor ^= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_or |= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_lshift <<= 1;  
  
        #pragma omp atomic	  
        atom_lshift >>= 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_rshift >>= 2;  
  
        #pragma omp atomic	  
        atom_rshift <<= 2;  
      }
      break;

    case 1:
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_incr ++;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        ++ atom_incr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_decr --;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        -- atom_decr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_plus += sameas(1) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_minus -= sameas(2) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_mul *= sameas(3) - 1;  
        #pragma omp atomic	  
        atom_mul /= 1 - sameas(3);  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_div /= 1 + sameas(3);  
        #pragma omp atomic	  
        atom_div *= 3 + sameas(1);  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_and &= ~(1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_xor ^= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_or |= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_lshift <<= 1;  
  
        #pragma omp atomic	  
        atom_lshift >>= 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_rshift >>= 2;  
  
        #pragma omp atomic	  
        atom_rshift <<= 2;  
      }
      break;

    case 2:
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_incr ++;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        ++ atom_incr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_decr --;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        -- atom_decr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_plus += sameas(1) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_minus -= sameas(2) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_mul *= sameas(3) - 1;  
        #pragma omp atomic	  
        atom_mul /= 1 - sameas(3);  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_div /= 1 + sameas(3);  
        #pragma omp atomic	  
        atom_div *= 3 + sameas(1);  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_and &= ~(1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_xor ^= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_or |= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_lshift <<= 1;  
  
        #pragma omp atomic	  
        atom_lshift >>= 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_rshift >>= 2;  
  
        #pragma omp atomic	  
        atom_rshift <<= 2;  
      }
      break;

    case 3:
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_incr ++;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        ++ atom_incr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        atom_decr --;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic  
        -- atom_decr2;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_plus += sameas(1) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_minus -= sameas(2) + 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_mul *= sameas(3) - 1;  
        #pragma omp atomic	  
        atom_mul /= 1 - sameas(3);  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_div /= 1 + sameas(3);  
        #pragma omp atomic	  
        atom_div *= 3 + sameas(1);  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_and &= ~(1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_xor ^= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=LOOPNUM/thds*id; i<LOOPNUM/thds*(id+1); i++) {
        #pragma omp atomic	  
        atom_or |= (1<<(i%(INT_BIT-1)));  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_lshift <<= 1;  
  
        #pragma omp atomic	  
        atom_lshift >>= 1;  
      }  
  
      barrier (thds);
      for (i=0; i<LOOPNUM/thds; i++) {  
        #pragma omp atomic	  
        atom_rshift >>= 2;  
  
        #pragma omp atomic	  
        atom_rshift <<= 2;  
      }
      break;

    default:
      break;
    }
  }
  errors += check ();

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

    func_atomic (id);
  }
  errors += check ();

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