#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 */
double exp(double);
void initialize(int n, double *dx, double *dy, double **u, double **f) {
  int i, j;
  double xxexp, yy;
  *dx = 1.0 / (n - 1);
  *dy = 1.0 / (n - 1);
#pragma omp parallel for default(shared) private(i, j, xxexp, yy)
  for (i = 0; i < n; i++) {
    xxexp = exp((*dx) * (i - 1));
    for (j = 0; j < n; j++) {
      yy = (*dy) * (j - 1);
      (u[i][j]) = 0.0;
      (f[i][j]) = xxexp * yy;
    }
  }
}


x86 assembly (clang, -O0):
<initialize>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xec3(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xea5(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1220 <initialize+0x110>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
call   1030 <exp@plt>
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    120d <initialize+0xfd>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11a1 <initialize+0x91>
jmp    1212 <initialize+0x102>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    116e <initialize+0x5e>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<initialize>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,0x20(%rsp)
mov    %rcx,0x18(%rsp)
mov    %edi,%r13d
lea    -0x1(%r13),%eax
cvtsi2sd %eax,%xmm0
movsd  0xec5(%rip),%xmm1
divsd  %xmm0,%xmm1
mov    %rsi,0x10(%rsp)
movsd  %xmm1,(%rsi)
movsd  %xmm1,(%rdx)
test   %edi,%edi
jle    133f <initialize+0x22f>
mov    %rdx,%rbx
mov    %r13d,%ebp
lea    0x8(%rdx),%r14
mov    %ebp,%eax
and    $0xfffffffc,%eax
mov    %rax,(%rsp)
mov    %rbp,%rax
neg    %rax
mov    %rax,0x8(%rsp)
lea    0x0(,%rbp,8),%r15
and    $0xffffffffffffffe0,%r15
xor    %r12d,%r12d
jmp    119d <initialize+0x8d>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r12
cmp    %rbp,%r12
je     133f <initialize+0x22f>
lea    -0x1(%r12),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x10(%rsp),%rax
mulsd  (%rax),%xmm0
call   1030 <exp@plt>
mov    0x18(%rsp),%rax
mov    (%rax,%r12,8),%rax
mov    0x20(%rsp),%rcx
mov    (%rcx,%r12,8),%rcx
cmp    $0x4,%r13d
jae    1260 <initialize+0x150>
xor    %esi,%esi
mov    %rsi,%rdx
test   $0x1,%bpl
je     1204 <initialize+0xf4>
lea    -0x1(%rsi),%edx
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
mulsd  (%rbx),%xmm1
movq   $0x0,(%rax,%rsi,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx,%rsi,8)
mov    %rsi,%rdx
or     $0x1,%rdx
not    %rsi
cmp    0x8(%rsp),%rsi
je     1190 <initialize+0x80>
xchg   %ax,%ax
lea    -0x1(%rdx),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  (%rbx),%xmm1
movq   $0x0,(%rax,%rdx,8)
mulsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
movsd  %xmm1,(%rcx,%rdx,8)
mulsd  (%rbx),%xmm2
movq   $0x0,0x8(%rax,%rdx,8)
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rcx,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rbp
jne    1210 <initialize+0x100>
jmp    1190 <initialize+0x80>
nopl   0x0(%rax,%rax,1)
lea    (%rax,%rbp,8),%rsi
lea    (%rcx,%rbp,8),%rdi
cmp    %rdi,%rax
setb   %r10b
cmp    %rsi,%rcx
setb   %r11b
cmp    %r14,%rax
setb   %dl
cmp    %rbx,%rsi
seta   %sil
cmp    %r14,%rcx
setb   %r9b
cmp    %rbx,%rdi
seta   %r8b
test   %r11b,%r10b
jne    11d3 <initialize+0xc3>
and    %sil,%dl
jne    11d3 <initialize+0xc3>
mov    $0x0,%esi
and    %r8b,%r9b
pcmpeqd %xmm6,%xmm6
xorps  %xmm7,%xmm7
movdqa 0xd65(%rip),%xmm8
jne    11d5 <initialize+0xc5>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
movsd  (%rbx),%xmm2
unpcklpd %xmm2,%xmm2
xor    %edx,%edx
movdqa 0xd35(%rip),%xmm3
nopl   0x0(%rax,%rax,1)
movdqa %xmm3,%xmm4
paddd  %xmm6,%xmm4
movdqa %xmm3,%xmm5
psubd  %xmm6,%xmm5
cvtdq2pd %xmm4,%xmm4
cvtdq2pd %xmm5,%xmm5
mulpd  %xmm2,%xmm4
mulpd  %xmm2,%xmm5
movups %xmm7,(%rax,%rdx,1)
movups %xmm7,0x10(%rax,%rdx,1)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,(%rcx,%rdx,1)
movupd %xmm5,0x10(%rcx,%rdx,1)
paddd  %xmm8,%xmm3
add    $0x20,%rdx
cmp    %rdx,%r15
jne    12e0 <initialize+0x1d0>
mov    (%rsp),%rdx
mov    %rdx,%rsi
cmp    %rbp,%rdx
je     1190 <initialize+0x80>
jmp    11d5 <initialize+0xc5>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @initialize(i32 %0, double* %1, double* %2, double** %3, double** %4) {
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca double**
%10 = alloca double**
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
store i32 %0, i32* %6
store double* %1, double** %7
store double* %2, double** %8
store double** %3, double*** %9
store double** %4, double*** %10
%15 = load i32, i32* %6
%16 = sub nsw i32 %15, 1
%17 = sitofp i32 %16 to double
%18 = fdiv double 1.000000e+00, %17
%19 = load double*, double** %7
store double %18, double* %19
%20 = load i32, i32* %6
%21 = sub nsw i32 %20, 1
%22 = sitofp i32 %21 to double
%23 = fdiv double 1.000000e+00, %22
%24 = load double*, double** %8
store double %23, double* %24
store i32 0, i32* %11
br label %25
25: ; preds = %71, %5
%26 = load i32, i32* %11
%27 = load i32, i32* %6
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %74
29: ; preds = %25
%30 = load double*, double** %7
%31 = load double, double* %30
%32 = load i32, i32* %11
%33 = sub nsw i32 %32, 1
%34 = sitofp i32 %33 to double
%35 = fmul double %31, %34
%36 = call double @exp(double %35)
store double %36, double* %13
store i32 0, i32* %12
br label %37
37: ; preds = %67, %29
%38 = load i32, i32* %12
%39 = load i32, i32* %6
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %70
41: ; preds = %37
%42 = load double*, double** %8
%43 = load double, double* %42
%44 = load i32, i32* %12
%45 = sub nsw i32 %44, 1
%46 = sitofp i32 %45 to double
%47 = fmul double %43, %46
store double %47, double* %14
%48 = load double**, double*** %9
%49 = load i32, i32* %11
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double*, double** %48, i64 %50
%52 = load double*, double** %51
%53 = load i32, i32* %12
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
store double 0.000000e+00, double* %55
%56 = load double, double* %13
%57 = load double, double* %14
%58 = fmul double %56, %57
%59 = load double**, double*** %10
%60 = load i32, i32* %11
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double*, double** %59, i64 %61
%63 = load double*, double** %62
%64 = load i32, i32* %12
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
store double %58, double* %66
br label %67
67: ; preds = %41
%68 = load i32, i32* %12
%69 = add nsw i32 %68, 1
store i32 %69, i32* %12
br label %37
70: ; preds = %37
br label %71
71: ; preds = %70
%72 = load i32, i32* %11
%73 = add nsw i32 %72, 1
store i32 %73, i32* %11
br label %25
74: ; preds = %25
ret void
}

LLVM IR (-O2):
define void @initialize(i32 %0, double* nocapture %1, double* nocapture %2, double** nocapture %3, double** nocapture %4) {
%6 = add nsw i32 %0, -1
%7 = sitofp i32 %6 to double
%8 = fdiv double 1.000000e+00, %7
store double %8, double* %1
store double %8, double* %2
%9 = icmp sgt i32 %0, 0
br i1 %9, label %10, label %118
10: ; preds = %5
%11 = zext i32 %0 to i64
%12 = getelementptr double, double* %2, i64 1
%13 = icmp ult i32 %0, 4
%14 = and i64 %11, 4294967292
%15 = icmp eq i64 %14, %11
%16 = and i64 %11, 1
%17 = icmp eq i64 %16, 0
%18 = sub nsw i64 0, %11
br label %19
19: ; preds = %115, %10
%20 = phi i64 [ 0, %10 ], [ %116, %115 ]
%21 = load double, double* %1
%22 = trunc i64 %20 to i32
%23 = add i32 %22, -1
%24 = sitofp i32 %23 to double
%25 = fmul double %21, %24
%26 = tail call double @exp(double %25)
%27 = getelementptr inbounds double*, double** %3, i64 %20
%28 = load double*, double** %27
%29 = getelementptr inbounds double*, double** %4, i64 %20
%30 = load double*, double** %29
br i1 %13, label %78, label %31
31: ; preds = %19
%32 = getelementptr double, double* %28, i64 %11
%33 = getelementptr double, double* %30, i64 %11
%34 = icmp ult double* %28, %33
%35 = icmp ult double* %30, %32
%36 = and i1 %34, %35
%37 = icmp ult double* %28, %12
%38 = icmp ugt double* %32, %2
%39 = and i1 %37, %38
%40 = or i1 %36, %39
%41 = icmp ult double* %30, %12
%42 = icmp ugt double* %33, %2
%43 = and i1 %41, %42
%44 = or i1 %40, %43
br i1 %44, label %78, label %45
45: ; preds = %31
%46 = insertelement <2 x double> poison, double %26, i64 0
%47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
%48 = insertelement <2 x double> poison, double %26, i64 0
%49 = shufflevector <2 x double> %48, <2 x double> poison, <2 x i32> zeroinitializer
%50 = load double, double* %2
%51 = insertelement <2 x double> poison, double %50, i64 0
%52 = shufflevector <2 x double> %51, <2 x double> poison, <2 x i32> zeroinitializer
%53 = insertelement <2 x double> poison, double %50, i64 0
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> zeroinitializer
br label %55
55: ; preds = %55, %45
%56 = phi i64 [ 0, %45 ], [ %74, %55 ]
%57 = phi <2 x i32> [ <i32 0, i32 1>, %45 ], [ %75, %55 ]
%58 = add <2 x i32> %57, <i32 -1, i32 -1>
%59 = add <2 x i32> %57, <i32 1, i32 1>
%60 = sitofp <2 x i32> %58 to <2 x double>
%61 = sitofp <2 x i32> %59 to <2 x double>
%62 = fmul <2 x double> %52, %60
%63 = fmul <2 x double> %54, %61
%64 = getelementptr inbounds double, double* %28, i64 %56
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %65
%66 = getelementptr inbounds double, double* %64, i64 2
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %67
%68 = fmul <2 x double> %47, %62
%69 = fmul <2 x double> %49, %63
%70 = getelementptr inbounds double, double* %30, i64 %56
%71 = bitcast double* %70 to <2 x double>*
store <2 x double> %68, <2 x double>* %71
%72 = getelementptr inbounds double, double* %70, i64 2
%73 = bitcast double* %72 to <2 x double>*
store <2 x double> %69, <2 x double>* %73
%74 = add nuw i64 %56, 4
%75 = add <2 x i32> %57, <i32 4, i32 4>
%76 = icmp eq i64 %74, %14
br i1 %76, label %77, label %55
77: ; preds = %55
br i1 %15, label %115, label %78
78: ; preds = %31, %19, %77
%79 = phi i64 [ 0, %31 ], [ 0, %19 ], [ %14, %77 ]
%80 = xor i64 %79, -1
br i1 %17, label %91, label %81
81: ; preds = %78
%82 = load double, double* %2
%83 = trunc i64 %79 to i32
%84 = add i32 %83, -1
%85 = sitofp i32 %84 to double
%86 = fmul double %82, %85
%87 = getelementptr inbounds double, double* %28, i64 %79
store double 0.000000e+00, double* %87
%88 = fmul double %26, %86
%89 = getelementptr inbounds double, double* %30, i64 %79
store double %88, double* %89
%90 = or i64 %79, 1
br label %91
91: ; preds = %81, %78
%92 = phi i64 [ %79, %78 ], [ %90, %81 ]
%93 = icmp eq i64 %80, %18
br i1 %93, label %115, label %94
94: ; preds = %91, %94
%95 = phi i64 [ %113, %94 ], [ %92, %91 ]
%96 = load double, double* %2
%97 = trunc i64 %95 to i32
%98 = add i32 %97, -1
%99 = sitofp i32 %98 to double
%100 = fmul double %96, %99
%101 = getelementptr inbounds double, double* %28, i64 %95
store double 0.000000e+00, double* %101
%102 = fmul double %26, %100
%103 = getelementptr inbounds double, double* %30, i64 %95
store double %102, double* %103
%104 = add nuw nsw i64 %95, 1
%105 = load double, double* %2
%106 = trunc i64 %104 to i32
%107 = add i32 %106, -1
%108 = sitofp i32 %107 to double
%109 = fmul double %105, %108
%110 = getelementptr inbounds double, double* %28, i64 %104
store double 0.000000e+00, double* %110
%111 = fmul double %26, %109
%112 = getelementptr inbounds double, double* %30, i64 %104
store double %111, double* %112
%113 = add nuw nsw i64 %95, 2
%114 = icmp eq i64 %113, %11
br i1 %114, label %115, label %94
115: ; preds = %91, %94, %77
%116 = add nuw nsw i64 %20, 1
%117 = icmp eq i64 %116, %11
br i1 %117, label %118, label %19
118: ; preds = %115, %5
ret void
}