1 files changed, 47 insertions, 0 deletions
diff --git a/libgomp/testsuite/libgomp.c/appendix-a/a.16.1.c b/libgomp/testsuite/libgomp.c/appendix-a/a.16.1.c
new file mode 100644
index 00000000000..28a994a28fa
--- /dev/null
+++ b/libgomp/testsuite/libgomp.c/appendix-a/a.16.1.c
@@ -0,0 +1,47 @@
+/* { dg-do run } */
+
+#include <stdio.h>
+
+float
+work1 (int i)
+{
+  return 1.0 * i;
+}
+
+float
+work2 (int i)
+{
+  return 2.0 * i;
+}
+
+void
+a16 (float *x, float *y, int *index, int n)
+{
+  int i;
+#pragma omp parallel for shared(x, y, index, n)
+  for (i = 0; i < n; i++)
+    {
+#pragma omp atomic
+      x[index[i]] += work1 (i);
+      y[i] += work2 (i);
+    }
+}
+int
+main ()
+{
+  float x[1000];
+  float y[10000];
+  int index[10000];
+  int i;
+  for (i = 0; i < 10000; i++)
+    {
+      index[i] = i % 1000;
+      y[i] = 0.0;
+    }
+  for (i = 0; i < 1000; i++)
+    x[i] = 0.0;
+  a16 (x, y, index, 10000);
+  for (i = 0; i < 10; i++)
+    printf ("x[%d] = %f, y[%d] = %f\n", i, x[i], i, y[i]);
+  return 0;
+}