Mercurial > hg > CbC > CbC_gcc

comparison gcc/testsuite/gfortran.dg/pr65450.f90 @ 111:04ced10e8804

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
gcc 7
author kono
date Fri, 27 Oct 2017 22:46:09 +0900
parents
children 84e7813d76e9
comparison
equal deleted inserted replaced
68:561a7518be6b 111:04ced10e8804
1 ! PR tree-optimization/65450
2 ! { dg-do run }
3 ! { dg-additional-options "-mtune=amdfam10" { target x86_64-*-* i?86-*-* } }
4
5 program pr65450
6 integer :: n, m, o, i, k
7 double precision :: u(500,60,3), h(500,60,3)
8 double precision :: v(500,60)
9 u = 0
10 h = 0
11 o = 1
12 m = 2
13 n = 3
14 do k = 1, 50
15 v = foo (u(:,:,m))
16 u(2:499,1:60,n) = u(2:499,1:60,o)+16.d0
17 h(1:500,2:59,n) = h(1:500,2:59,o)-4.d0*v(1:500,2:59)-32.0d0
18 i = o
19 o = m
20 m = n
21 n = i
22 end do
23 if (abs (v(17, 23) + h(17, 23, 2) + 768.0d0) > 0.5d0) call abort
24 contains
25 function foo(a)
26 double precision :: a(:,:)
27 double precision :: foo(size(a,dim=1),size(a,dim=2))
28 integer :: i, j
29 i = size(a,dim=1)
30 j = size(a,dim=2)
31 foo(2:i-1,1:j) = a(3:i,1:j)-a(1:i-2,1:j)
32 foo(1,1:j) = 2*(a(2,1:j)-a(1,1:j))
33 foo(i,1:j) = 2*(a(i,1:j)-a(i-1,1:j))
34 end function foo
35 end program pr65450