source: CIVL/examples/omp/dataracebench-1.3.2/micro-benchmarks-fortran/DRB065-pireduction-orig-no.f95

main
Last change on this file was ea777aa, checked in by Alex Wilton <awilton@…>, 3 years ago

Moved examples, include, build_default.properties, common.xml, and README out from dev.civl.com into the root of the repo.

git-svn-id: svn://vsl.cis.udel.edu/civl/trunk@5704 fb995dde-84ed-4084-dfe6-e5aef3e2452c

  • Property mode set to 100755
File size: 983 bytes
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
10program 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
32end program
Note: See TracBrowser for help on using the repository browser.