1.23
2.0
acw/focus-triggers
main
test-branch
| Line | |
|---|
| 1 | !!!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!!!
|
|---|
| 2 | !!! Copyright (c) 2017-20, Lawrence Livermore National Security, LLC
|
|---|
| 3 | !!! and DataRaceBench project contributors. See the DataRaceBench/COPYRIGHT file for details.
|
|---|
| 4 | !!!
|
|---|
| 5 | !!! SPDX-License-Identifier: (BSD-3-Clause)
|
|---|
| 6 | !!!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!!!
|
|---|
| 7 |
|
|---|
| 8 | !Classic PI calculation using reduction. No data race pairs.
|
|---|
| 9 |
|
|---|
| 10 | program DRB065_pireduction_orig_no
|
|---|
| 11 | use omp_lib
|
|---|
| 12 | implicit none
|
|---|
| 13 |
|
|---|
| 14 | real(kind = 16) :: x, interval_width, pi
|
|---|
| 15 | integer(kind = 8) :: i, num_steps
|
|---|
| 16 |
|
|---|
| 17 | pi = 0.0
|
|---|
| 18 | num_steps = 2000000000
|
|---|
| 19 | interval_width = 1.0/num_steps
|
|---|
| 20 |
|
|---|
| 21 | !$omp parallel do reduction(+:pi) private(x)
|
|---|
| 22 | do i = 1, num_steps
|
|---|
| 23 | x = (i + 0.5) * interval_width
|
|---|
| 24 | pi = pi + 1.0 / (x*x + 1.0);
|
|---|
| 25 | end do
|
|---|
| 26 | !$omp end parallel do
|
|---|
| 27 |
|
|---|
| 28 | pi = pi * 4.0 * interval_width;
|
|---|
| 29 | print 100, pi
|
|---|
| 30 | 100 format ("PI =",F24.20)
|
|---|
| 31 |
|
|---|
| 32 | end program
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
