#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ assert(int);
void v_p_2() {
  int i, j, k;
  double ina[1024], inb[1024], out[1024], acc;
  for (j = 0; j < 32; j++)
    for (i = 0; i < 32; i++) {
      ina[j * 32 + i] = (i == j) ? 2.0 : 0.0;
      inb[j * 32 + i] = (double)(i + j);
    }
#pragma acc parallel num_gangs(32) num_workers(32)                             \
    vector_length(32) private(acc) copyin(ina, inb) copyout(out)
  {
#pragma acc loop gang worker
    for (k = 0; k < 32; k++)
      for (j = 0; j < 32; j++) {
        acc = 0;
#pragma acc loop vector reduction(+ : acc)
        for (i = 0; i < 32; i++)
          acc += ina[k * 32 + i] * inb[i * 32 + j];
        out[k * 32 + j] = acc;
      }
  }
  for (j = 0; j < 32; j++)
    for (i = 0; i < 32; i++)
      assert(out[j * 32 + i] == (i + j) * 2);
}


x86 assembly (clang, -O0):
<v_p_2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x6030,%rsp
movl   $0x0,-0x8(%rbp)
cmpl   $0x20,-0x8(%rbp)
jge    11d1 <v_p_2+0xc1>
movl   $0x0,-0x4(%rbp)
cmpl   $0x20,-0x4(%rbp)
jge    11be <v_p_2+0xae>
mov    -0x4(%rbp),%eax
mov    -0x8(%rbp),%ecx
movsd  0xeb5(%rip),%xmm0
xorps  %xmm1,%xmm1
movsd  %xmm1,-0x6028(%rbp)
ff
cmp    %ecx,%eax
movsd  %xmm0,-0x6020(%rbp)
ff
je     1176 <v_p_2+0x66>
movsd  -0x6028(%rbp),%xmm0
ff
movsd  %xmm0,-0x6020(%rbp)
ff
movsd  -0x6020(%rbp),%xmm0
ff
mov    -0x8(%rbp),%eax
shl    $0x5,%eax
add    -0x4(%rbp),%eax
cltq
movsd  %xmm0,-0x2010(%rbp,%rax,8)
ff ff
mov    -0x4(%rbp),%eax
add    -0x8(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x8(%rbp),%eax
shl    $0x5,%eax
add    -0x4(%rbp),%eax
cltq
movsd  %xmm0,-0x4010(%rbp,%rax,8)
ff ff
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1133 <v_p_2+0x23>
jmp    11c3 <v_p_2+0xb3>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    1122 <v_p_2+0x12>
movl   $0x0,-0xc(%rbp)
cmpl   $0x20,-0xc(%rbp)
jge    129a <v_p_2+0x18a>
movl   $0x0,-0x8(%rbp)
cmpl   $0x20,-0x8(%rbp)
jge    1287 <v_p_2+0x177>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x6018(%rbp)
ff
movl   $0x0,-0x4(%rbp)
cmpl   $0x20,-0x4(%rbp)
jge    125d <v_p_2+0x14d>
mov    -0xc(%rbp),%eax
shl    $0x5,%eax
add    -0x4(%rbp),%eax
cltq
movsd  -0x2010(%rbp,%rax,8),%xmm0
ff ff
mov    -0x4(%rbp),%eax
shl    $0x5,%eax
add    -0x8(%rbp),%eax
cltq
movsd  -0x4010(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x6018(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x6018(%rbp)
ff
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1205 <v_p_2+0xf5>
movsd  -0x6018(%rbp),%xmm0
ff
mov    -0xc(%rbp),%eax
shl    $0x5,%eax
add    -0x8(%rbp),%eax
cltq
movsd  %xmm0,-0x6010(%rbp,%rax,8)
ff ff
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    11e9 <v_p_2+0xd9>
jmp    128c <v_p_2+0x17c>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    11d8 <v_p_2+0xc8>
movl   $0x0,-0x8(%rbp)
cmpl   $0x20,-0x8(%rbp)
jge    1314 <v_p_2+0x204>
movl   $0x0,-0x4(%rbp)
cmpl   $0x20,-0x4(%rbp)
jge    1301 <v_p_2+0x1f1>
mov    -0x8(%rbp),%eax
shl    $0x5,%eax
add    -0x4(%rbp),%eax
cltq
movsd  -0x6010(%rbp,%rax,8),%xmm0
ff ff
mov    -0x4(%rbp),%eax
add    -0x8(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
sete   %al
setnp  %cl
and    %cl,%al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    12b2 <v_p_2+0x1a2>
jmp    1306 <v_p_2+0x1f6>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    12a1 <v_p_2+0x191>
add    $0x6030,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<v_p_2>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x6008,%rsp
lea    0x2008(%rsp),%rax
lea    0x4008(%rsp),%rcx
xor    %edx,%edx
movsd  0xec8(%rip),%xmm0
jmp    1157 <v_p_2+0x47>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
add    $0x100,%rax
add    $0x100,%rcx
cmp    $0x20,%rdx
je     11b3 <v_p_2+0xa3>
lea    -0x1(%rdx),%rsi
xor    %edi,%edi
jmp    1181 <v_p_2+0x71>
nop
lea    (%rdx,%rdi,1),%ebx
add    $0x1,%ebx
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
movsd  %xmm1,(%rcx,%rdi,8)
movsd  %xmm2,(%rax,%rdi,8)
add    $0x2,%rdi
cmp    $0x20,%rdi
je     1140 <v_p_2+0x30>
movapd %xmm0,%xmm1
cmp    %rdi,%rdx
je     118e <v_p_2+0x7e>
xorpd  %xmm1,%xmm1
movsd  %xmm1,-0x8(%rcx,%rdi,8)
lea    (%rdx,%rdi,1),%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
movsd  %xmm1,-0x8(%rax,%rdi,8)
movapd %xmm0,%xmm1
cmp    %rdi,%rsi
je     1160 <v_p_2+0x50>
xorpd  %xmm1,%xmm1
jmp    1160 <v_p_2+0x50>
lea    0x4008(%rsp),%rax
lea    0x2100(%rsp),%r8
xor    %r9d,%r9d
cs nopw 0x0(%rax,%rax,1)
mov    %r9,%rsi
shl    $0x5,%rsi
mov    %r8,%rcx
xor    %ebx,%ebx
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
mov    %rcx,%rdi
xor    %edx,%edx
nopl   0x0(%rax)
movsd  -0x8(%rax,%rdx,8),%xmm1
mulsd  -0x100(%rdi),%xmm1
ff
addsd  %xmm0,%xmm1
movsd  (%rax,%rdx,8),%xmm0
mulsd  (%rdi),%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%rdx
add    $0x200,%rdi
cmp    $0x20,%rdx
jne    11f0 <v_p_2+0xe0>
lea    (%rbx,%rsi,1),%rdx
movsd  %xmm0,(%rsp,%rdx,8)
add    $0x1,%rbx
add    $0x8,%rcx
cmp    $0x20,%rbx
jne    11e0 <v_p_2+0xd0>
add    $0x1,%r9
add    $0x100,%rax
cmp    $0x20,%r9
jne    11d0 <v_p_2+0xc0>
mov    %rsp,%r14
xor    %r12d,%r12d
xor    %r15d,%r15d
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%r12,%rbx,1),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
cmpeqsd (%r14,%rbx,4),%xmm0
movq   %xmm0,%rdi
and    $0x1,%edi
call   1030 <assert@plt>
add    $0x2,%rbx
cmp    $0x40,%rbx
jne    1260 <v_p_2+0x150>
add    $0x1,%r15
add    $0x2,%r12
add    $0x100,%r14
cmp    $0x20,%r15
jne    1250 <v_p_2+0x140>
add    $0x6008,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @v_p_2() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca [1024 x double]
%5 = alloca [1024 x double]
%6 = alloca [1024 x double]
%7 = alloca double
store i32 0, i32* %2
br label %8
8: ; preds = %41, %0
%9 = load i32, i32* %2
%10 = icmp slt i32 %9, 32
br i1 %10, label %11, label %44
11: ; preds = %8
store i32 0, i32* %1
br label %12
12: ; preds = %37, %11
%13 = load i32, i32* %1
%14 = icmp slt i32 %13, 32
br i1 %14, label %15, label %40
15: ; preds = %12
%16 = load i32, i32* %1
%17 = load i32, i32* %2
%18 = icmp eq i32 %16, %17
%19 = zext i1 %18 to i64
%20 = select i1 %18, double 2.000000e+00, double 0.000000e+00
%21 = load i32, i32* %2
%22 = mul nsw i32 %21, 32
%23 = load i32, i32* %1
%24 = add nsw i32 %22, %23
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds [1024 x double], [1024 x double]* %4, i64 0, i64 %25
store double %20, double* %26
%27 = load i32, i32* %1
%28 = load i32, i32* %2
%29 = add nsw i32 %27, %28
%30 = sitofp i32 %29 to double
%31 = load i32, i32* %2
%32 = mul nsw i32 %31, 32
%33 = load i32, i32* %1
%34 = add nsw i32 %32, %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [1024 x double], [1024 x double]* %5, i64 0, i64 %35
store double %30, double* %36
br label %37
37: ; preds = %15
%38 = load i32, i32* %1
%39 = add nsw i32 %38, 1
store i32 %39, i32* %1
br label %12
40: ; preds = %12
br label %41
41: ; preds = %40
%42 = load i32, i32* %2
%43 = add nsw i32 %42, 1
store i32 %43, i32* %2
br label %8
44: ; preds = %8
store i32 0, i32* %3
br label %45
45: ; preds = %88, %44
%46 = load i32, i32* %3
%47 = icmp slt i32 %46, 32
br i1 %47, label %48, label %91
48: ; preds = %45
store i32 0, i32* %2
br label %49
49: ; preds = %84, %48
%50 = load i32, i32* %2
%51 = icmp slt i32 %50, 32
br i1 %51, label %52, label %87
52: ; preds = %49
store double 0.000000e+00, double* %7
store i32 0, i32* %1
br label %53
53: ; preds = %73, %52
%54 = load i32, i32* %1
%55 = icmp slt i32 %54, 32
br i1 %55, label %56, label %76
56: ; preds = %53
%57 = load i32, i32* %3
%58 = mul nsw i32 %57, 32
%59 = load i32, i32* %1
%60 = add nsw i32 %58, %59
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds [1024 x double], [1024 x double]* %4, i64 0, i64 %61
%63 = load double, double* %62
%64 = load i32, i32* %1
%65 = mul nsw i32 %64, 32
%66 = load i32, i32* %2
%67 = add nsw i32 %65, %66
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds [1024 x double], [1024 x double]* %5, i64 0, i64 %68
%70 = load double, double* %69
%71 = load double, double* %7
%72 = call double @llvm.fmuladd.f64(double %63, double %70, double %71)
store double %72, double* %7
br label %73
73: ; preds = %56
%74 = load i32, i32* %1
%75 = add nsw i32 %74, 1
store i32 %75, i32* %1
br label %53
76: ; preds = %53
%77 = load double, double* %7
%78 = load i32, i32* %3
%79 = mul nsw i32 %78, 32
%80 = load i32, i32* %2
%81 = add nsw i32 %79, %80
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds [1024 x double], [1024 x double]* %6, i64 0, i64 %82
store double %77, double* %83
br label %84
84: ; preds = %76
%85 = load i32, i32* %2
%86 = add nsw i32 %85, 1
store i32 %86, i32* %2
br label %49
87: ; preds = %49
br label %88
88: ; preds = %87
%89 = load i32, i32* %3
%90 = add nsw i32 %89, 1
store i32 %90, i32* %3
br label %45
91: ; preds = %45
store i32 0, i32* %2
br label %92
92: ; preds = %119, %91
%93 = load i32, i32* %2
%94 = icmp slt i32 %93, 32
br i1 %94, label %95, label %122
95: ; preds = %92
store i32 0, i32* %1
br label %96
96: ; preds = %115, %95
%97 = load i32, i32* %1
%98 = icmp slt i32 %97, 32
br i1 %98, label %99, label %118
99: ; preds = %96
%100 = load i32, i32* %2
%101 = mul nsw i32 %100, 32
%102 = load i32, i32* %1
%103 = add nsw i32 %101, %102
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds [1024 x double], [1024 x double]* %6, i64 0, i64 %104
%106 = load double, double* %105
%107 = load i32, i32* %1
%108 = load i32, i32* %2
%109 = add nsw i32 %107, %108
%110 = mul nsw i32 %109, 2
%111 = sitofp i32 %110 to double
%112 = fcmp oeq double %106, %111
%113 = zext i1 %112 to i32
%114 = call i32 @assert(i32 %113)
br label %115
115: ; preds = %99
%116 = load i32, i32* %1
%117 = add nsw i32 %116, 1
store i32 %117, i32* %1
br label %96
118: ; preds = %96
br label %119
119: ; preds = %118
%120 = load i32, i32* %2
%121 = add nsw i32 %120, 1
store i32 %121, i32* %2
br label %92
122: ; preds = %92
ret void
}

LLVM IR (-O2):
define void @v_p_2() {
%1 = alloca [1024 x double]
%2 = alloca [1024 x double]
%3 = alloca [1024 x double]
%4 = bitcast [1024 x double]* %1 to i8*
%5 = bitcast [1024 x double]* %2 to i8*
%6 = bitcast [1024 x double]* %3 to i8*
br label %7
7: ; preds = %0, %31
%8 = phi i64 [ 0, %0 ], [ %32, %31 ]
%9 = shl nsw i64 %8, 5
br label %10
10: ; preds = %10, %7
%11 = phi i64 [ 0, %7 ], [ %29, %10 ]
%12 = icmp eq i64 %11, %8
%13 = select i1 %12, double 2.000000e+00, double 0.000000e+00
%14 = add nuw nsw i64 %11, %9
%15 = getelementptr inbounds [1024 x double], [1024 x double]* %1, i64 0, i64 %14
store double %13, double* %15
%16 = add nuw nsw i64 %11, %8
%17 = trunc i64 %16 to i32
%18 = sitofp i32 %17 to double
%19 = getelementptr inbounds [1024 x double], [1024 x double]* %2, i64 0, i64 %14
store double %18, double* %19
%20 = or i64 %11, 1
%21 = icmp eq i64 %20, %8
%22 = select i1 %21, double 2.000000e+00, double 0.000000e+00
%23 = add nuw nsw i64 %20, %9
%24 = getelementptr inbounds [1024 x double], [1024 x double]* %1, i64 0, i64 %23
store double %22, double* %24
%25 = add nuw nsw i64 %20, %8
%26 = trunc i64 %25 to i32
%27 = sitofp i32 %26 to double
%28 = getelementptr inbounds [1024 x double], [1024 x double]* %2, i64 0, i64 %23
store double %27, double* %28
%29 = add nuw nsw i64 %11, 2
%30 = icmp eq i64 %29, 32
br i1 %30, label %31, label %10
31: ; preds = %10
%32 = add nuw nsw i64 %8, 1
%33 = icmp eq i64 %32, 32
br i1 %33, label %34, label %7
34: ; preds = %31, %66
%35 = phi i64 [ %67, %66 ], [ 0, %31 ]
%36 = shl nsw i64 %35, 5
br label %37
37: ; preds = %34, %61
%38 = phi i64 [ 0, %34 ], [ %64, %61 ]
br label %39
39: ; preds = %39, %37
%40 = phi i64 [ 0, %37 ], [ %59, %39 ]
%41 = phi double [ 0.000000e+00, %37 ], [ %58, %39 ]
%42 = add nuw nsw i64 %40, %36
%43 = getelementptr inbounds [1024 x double], [1024 x double]* %1, i64 0, i64 %42
%44 = load double, double* %43
%45 = shl nuw nsw i64 %40, 5
%46 = add nuw nsw i64 %45, %38
%47 = getelementptr inbounds [1024 x double], [1024 x double]* %2, i64 0, i64 %46
%48 = load double, double* %47
%49 = tail call double @llvm.fmuladd.f64(double %44, double %48, double %41)
%50 = or i64 %40, 1
%51 = add nuw nsw i64 %50, %36
%52 = getelementptr inbounds [1024 x double], [1024 x double]* %1, i64 0, i64 %51
%53 = load double, double* %52
%54 = shl nuw nsw i64 %50, 5
%55 = add nuw nsw i64 %54, %38
%56 = getelementptr inbounds [1024 x double], [1024 x double]* %2, i64 0, i64 %55
%57 = load double, double* %56
%58 = tail call double @llvm.fmuladd.f64(double %53, double %57, double %49)
%59 = add nuw nsw i64 %40, 2
%60 = icmp eq i64 %59, 32
br i1 %60, label %61, label %39
61: ; preds = %39
%62 = add nuw nsw i64 %38, %36
%63 = getelementptr inbounds [1024 x double], [1024 x double]* %3, i64 0, i64 %62
store double %58, double* %63
%64 = add nuw nsw i64 %38, 1
%65 = icmp eq i64 %64, 32
br i1 %65, label %66, label %37
66: ; preds = %61
%67 = add nuw nsw i64 %35, 1
%68 = icmp eq i64 %67, 32
br i1 %68, label %69, label %34
69: ; preds = %66, %86
%70 = phi i64 [ %87, %86 ], [ 0, %66 ]
%71 = shl nsw i64 %70, 5
br label %72
72: ; preds = %69, %72
%73 = phi i64 [ 0, %69 ], [ %84, %72 ]
%74 = add nuw nsw i64 %73, %71
%75 = getelementptr inbounds [1024 x double], [1024 x double]* %3, i64 0, i64 %74
%76 = load double, double* %75
%77 = add nuw nsw i64 %73, %70
%78 = trunc i64 %77 to i32
%79 = shl i32 %78, 1
%80 = sitofp i32 %79 to double
%81 = fcmp oeq double %76, %80
%82 = zext i1 %81 to i32
%83 = tail call i32 @assert(i32 %82)
%84 = add nuw nsw i64 %73, 1
%85 = icmp eq i64 %84, 32
br i1 %85, label %86, label %72
86: ; preds = %72
%87 = add nuw nsw i64 %70, 1
%88 = icmp eq i64 %87, 32
br i1 %88, label %89, label %69
89: ; preds = %86
ret void
}