source: CIVL/mods/dev.civl.com/examples/omp/dataracebench-1.3.2/micro-benchmarks-fortran/DRB038-truedepseconddimension-var-yes.f95@ cb4d4f4

1.23 2.0 main test-branch
Last change on this file since cb4d4f4 was aad342c, checked in by Stephen Siegel <siegel@…>, 3 years ago

Performing huge refactor to incorporate ABC, GMC, and SARL into CIVL repo and use Java modules.

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

  • Property mode set to 100755
File size: 1.6 KB
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!Only the outmost loop can be parallelized in this program.
9!Data race pair: b[i][j]@51:13 vs. b[i][j-1]@51:22
10
11program DRB038_truedepseconddimension_var_yes
12 use omp_lib
13 implicit none
14
15 integer i, j, n, m, len, argCount, allocStatus, rdErr, ix
16 character(len=80), dimension(:), allocatable :: args
17 real, dimension(:,:), allocatable :: b
18
19 len = 1000
20
21 argCount = command_argument_count()
22 if (argCount == 0) then
23 write (*,'(a)') "No command line arguments provided."
24 end if
25
26 allocate(args(argCount), stat=allocStatus)
27 if (allocStatus > 0) then
28 write (*,'(a)') "Allocation error, program terminated."
29 stop
30 end if
31
32 do ix = 1, argCount
33 call get_command_argument(ix,args(ix))
34 end do
35
36 if (argCount >= 1) then
37 read (args(1), '(i10)', iostat=rdErr) len
38 if (rdErr /= 0 ) then
39 write (*,'(a)') "Error, invalid integer value."
40 end if
41 end if
42
43 n = len
44 m = len
45
46 allocate (b(len,len))
47
48 do i = 1, n
49 !$omp parallel do
50 do j = 2, m
51 b(i,j) = b(i,j-1)
52 end do
53 !$omp end parallel do
54 end do
55! print 100, b(5,5)
56! 100 format ('b(5,5) =', F20.6)
57
58 deallocate(args,b)
59end program
Note: See TracBrowser for help on using the repository browser.