#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 fabs(double);
double pow(double, double);
double duu_dx(double **U, int i, int j, double dx, double gamma) {
  return 1.0 / dx *
             (pow((U[i][j] + U[i + 1][j]) / 2.0, 2.0) -
              pow((U[i - 1][j] + U[i][j]) / 2.0, 2.0)) +
         gamma / dx *
             (fabs(U[i][j] + U[i + 1][j]) / 2.0 * (U[i][j] - U[i + 1][j]) /
                  2.0 -
              fabs(U[i - 1][j] + U[i][j]) / 2.0 * (U[i - 1][j] - U[i][j]) /
                  2.0);
}


x86 assembly (clang, -O0):
<duu_dx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  0xed4(%rip),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xe8b(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xe7f(%rip),%xmm1
call   1030 <pow@plt>
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xe3e(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xe32(%rip),%xmm1
call   1030 <pow@plt>
movsd  -0x30(%rbp),%xmm2
movaps %xmm0,%xmm1
movsd  -0x28(%rbp),%xmm0
subsd  %xmm1,%xmm2
movsd  -0x20(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm3
movaps 0xdec(%rip),%xmm4
pand   %xmm4,%xmm3
movsd  0xdd0(%rip),%xmm4
divsd  %xmm4,%xmm3
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm4
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm4
mulsd  %xmm4,%xmm3
movsd  0xd91(%rip),%xmm4
divsd  %xmm4,%xmm3
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm4
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm4
movaps 0xd67(%rip),%xmm5
pand   %xmm5,%xmm4
movsd  0xd4b(%rip),%xmm5
divsd  %xmm5,%xmm4
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm5
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm5
mulsd  %xmm5,%xmm4
movsd  0xd0c(%rip),%xmm5
divsd  %xmm5,%xmm4
subsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<duu_dx>:
movapd %xmm0,%xmm9
movsd  0xef3(%rip),%xmm4
divsd  %xmm0,%xmm4
movslq %esi,%rax
mov    -0x8(%rdi,%rax,8),%rcx
mov    (%rdi,%rax,8),%rsi
movslq %edx,%rdx
mov    0x8(%rdi,%rax,8),%rax
movsd  (%rsi,%rdx,8),%xmm0
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movsd  (%rax,%rdx,8),%xmm5
movsd  (%rcx,%rdx,8),%xmm6
movapd %xmm5,%xmm7
unpcklpd %xmm6,%xmm7
addpd  %xmm3,%xmm7
movapd 0xebf(%rip),%xmm8
movapd %xmm7,%xmm2
mulpd  %xmm8,%xmm2
mulpd  %xmm2,%xmm2
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm4,%xmm2
divsd  %xmm9,%xmm1
andpd  0xea5(%rip),%xmm7
mulpd  %xmm8,%xmm7
unpcklpd %xmm0,%xmm5
unpcklpd %xmm6,%xmm0
subpd  %xmm5,%xmm0
mulpd  %xmm7,%xmm0
mulpd  %xmm8,%xmm0
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
subsd  %xmm3,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
ret

LLVM IR (-O0):
define double @duu_dx(double** %0, i32 %1, i32 %2, double %3, double %4) {
%6 = alloca double**
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca double
store double** %0, double*** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double %3, double* %9
store double %4, double* %10
%11 = load double, double* %9
%12 = fdiv double 1.000000e+00, %11
%13 = load double**, double*** %6
%14 = load i32, i32* %7
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double*, double** %13, i64 %15
%17 = load double*, double** %16
%18 = load i32, i32* %8
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %17, i64 %19
%21 = load double, double* %20
%22 = load double**, double*** %6
%23 = load i32, i32* %7
%24 = add nsw i32 %23, 1
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double*, double** %22, i64 %25
%27 = load double*, double** %26
%28 = load i32, i32* %8
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
%32 = fadd double %21, %31
%33 = fdiv double %32, 2.000000e+00
%34 = call double @pow(double %33, double 2.000000e+00)
%35 = load double**, double*** %6
%36 = load i32, i32* %7
%37 = sub nsw i32 %36, 1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double*, double** %35, i64 %38
%40 = load double*, double** %39
%41 = load i32, i32* %8
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
%44 = load double, double* %43
%45 = load double**, double*** %6
%46 = load i32, i32* %7
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double*, double** %45, i64 %47
%49 = load double*, double** %48
%50 = load i32, i32* %8
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = fadd double %44, %53
%55 = fdiv double %54, 2.000000e+00
%56 = call double @pow(double %55, double 2.000000e+00)
%57 = fsub double %34, %56
%58 = load double, double* %10
%59 = load double, double* %9
%60 = fdiv double %58, %59
%61 = load double**, double*** %6
%62 = load i32, i32* %7
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double*, double** %61, i64 %63
%65 = load double*, double** %64
%66 = load i32, i32* %8
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %65, i64 %67
%69 = load double, double* %68
%70 = load double**, double*** %6
%71 = load i32, i32* %7
%72 = add nsw i32 %71, 1
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double*, double** %70, i64 %73
%75 = load double*, double** %74
%76 = load i32, i32* %8
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = fadd double %69, %79
%81 = call double @llvm.fabs.f64(double %80)
%82 = fdiv double %81, 2.000000e+00
%83 = load double**, double*** %6
%84 = load i32, i32* %7
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double*, double** %83, i64 %85
%87 = load double*, double** %86
%88 = load i32, i32* %8
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = load double**, double*** %6
%93 = load i32, i32* %7
%94 = add nsw i32 %93, 1
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double*, double** %92, i64 %95
%97 = load double*, double** %96
%98 = load i32, i32* %8
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %97, i64 %99
%101 = load double, double* %100
%102 = fsub double %91, %101
%103 = fmul double %82, %102
%104 = fdiv double %103, 2.000000e+00
%105 = load double**, double*** %6
%106 = load i32, i32* %7
%107 = sub nsw i32 %106, 1
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double*, double** %105, i64 %108
%110 = load double*, double** %109
%111 = load i32, i32* %8
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds double, double* %110, i64 %112
%114 = load double, double* %113
%115 = load double**, double*** %6
%116 = load i32, i32* %7
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double*, double** %115, i64 %117
%119 = load double*, double** %118
%120 = load i32, i32* %8
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %119, i64 %121
%123 = load double, double* %122
%124 = fadd double %114, %123
%125 = call double @llvm.fabs.f64(double %124)
%126 = fdiv double %125, 2.000000e+00
%127 = load double**, double*** %6
%128 = load i32, i32* %7
%129 = sub nsw i32 %128, 1
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds double*, double** %127, i64 %130
%132 = load double*, double** %131
%133 = load i32, i32* %8
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = load double**, double*** %6
%138 = load i32, i32* %7
%139 = sext i32 %138 to i64
%140 = getelementptr inbounds double*, double** %137, i64 %139
%141 = load double*, double** %140
%142 = load i32, i32* %8
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %141, i64 %143
%145 = load double, double* %144
%146 = fsub double %136, %145
%147 = fmul double %126, %146
%148 = fdiv double %147, 2.000000e+00
%149 = fsub double %104, %148
%150 = fmul double %60, %149
%151 = call double @llvm.fmuladd.f64(double %12, double %57, double %150)
ret double %151
}

LLVM IR (-O2):
define double @duu_dx(double** nocapture %0, i32 %1, i32 %2, double %3, double %4) {
%6 = fdiv double 1.000000e+00, %3
%7 = sext i32 %1 to i64
%8 = getelementptr inbounds double*, double** %0, i64 %7
%9 = load double*, double** %8
%10 = sext i32 %2 to i64
%11 = getelementptr inbounds double, double* %9, i64 %10
%12 = load double, double* %11
%13 = add nsw i32 %1, 1
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double*, double** %0, i64 %14
%16 = load double*, double** %15
%17 = getelementptr inbounds double, double* %16, i64 %10
%18 = load double, double* %17
%19 = add nsw i32 %1, -1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double*, double** %0, i64 %20
%22 = load double*, double** %21
%23 = getelementptr inbounds double, double* %22, i64 %10
%24 = load double, double* %23
%25 = insertelement <2 x double> poison, double %12, i64 0
%26 = shufflevector <2 x double> %25, <2 x double> poison, <2 x i32> zeroinitializer
%27 = insertelement <2 x double> poison, double %18, i64 0
%28 = insertelement <2 x double> %27, double %24, i64 1
%29 = fadd <2 x double> %26, %28
%30 = fmul <2 x double> %29, <double 5.000000e-01, double 5.000000e-01>
%31 = fmul <2 x double> %30, %30
%32 = shufflevector <2 x double> %31, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%33 = fsub <2 x double> %31, %32
%34 = extractelement <2 x double> %33, i64 0
%35 = fdiv double %4, %3
%36 = call <2 x double> @llvm.fabs.v2f64(<2 x double> %29)
%37 = fmul <2 x double> %36, <double 5.000000e-01, double 5.000000e-01>
%38 = insertelement <2 x double> %25, double %24, i64 1
%39 = insertelement <2 x double> %27, double %12, i64 1
%40 = fsub <2 x double> %38, %39
%41 = fmul <2 x double> %40, %37
%42 = fmul <2 x double> %41, <double 5.000000e-01, double 5.000000e-01>
%43 = shufflevector <2 x double> %42, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%44 = fsub <2 x double> %42, %43
%45 = extractelement <2 x double> %44, i64 0
%46 = fmul double %35, %45
%47 = tail call double @llvm.fmuladd.f64(double %6, double %34, double %46)
ret double %47
}