#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 */
short ModDiff(float, float);
void IntegratePhase(float **psi, float **phi, short **flows, long nrow,
                    long ncol) {
  long row, col;
  short **rowflow, **colflow;
  rowflow = flows;
  colflow = &(flows[nrow - 1]);
  phi[0][0] = psi[0][0];
  for (col = 1; col < ncol; col++) {
    phi[0][col] =
        phi[0][col - 1] + (ModDiff(psi[0][col], psi[0][col - 1]) +
                           colflow[0][col - 1] * 6.28318530717958647692);
  }
  for (row = 1; row < nrow; row++) {
    for (col = 0; col < ncol; col++) {
      phi[row][col] =
          phi[row - 1][col] + (ModDiff(psi[row][col], psi[row - 1][col]) -
                               rowflow[row - 1][col] * 6.28318530717958647692);
    }
  }
}


x86 assembly (clang, -O0):
<IntegratePhase>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
movss  %xmm0,(%rax)
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1217 <IntegratePhase+0x107>
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm1
call   1030 <ModDiff@plt>
movsd  -0x50(%rbp),%xmm0
cwtl
cvtsi2sd %eax,%xmm2
mov    -0x48(%rbp),%rax
mov    (%rax),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
movsd  0xe1a(%rip),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1169 <IntegratePhase+0x59>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    1329 <IntegratePhase+0x219>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1313 <IntegratePhase+0x203>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
call   1030 <ModDiff@plt>
movsd  -0x58(%rbp),%xmm0
cwtl
cvtsi2sd %eax,%xmm2
mov    -0x40(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xd23(%rip),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1235 <IntegratePhase+0x125>
jmp    1318 <IntegratePhase+0x208>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    121f <IntegratePhase+0x10f>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<IntegratePhase>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8,%r15
mov    %rcx,%r14
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%rbx
mov    (%rdi),%rax
movss  (%rax),%xmm0
mov    (%rsi),%rcx
movss  %xmm0,(%rcx)
cmp    $0x2,%r8
jl     11b3 <IntegratePhase+0xa3>
lea    -0x2(%r15),%rcx
mov    %rcx,0x8(%rsp)
xor    %ebp,%ebp
nopl   0x0(%rax)
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%rsp)
movss  (%rax,%rbp,4),%xmm1
movss  0x4(%rax,%rbp,4),%xmm0
call   1030 <ModDiff@plt>
cwtl
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x8(%r12,%r14,8),%rax
movswl (%rax,%rbp,2),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  0xe77(%rip),%xmm1
addsd  %xmm0,%xmm1
addsd  (%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
mov    0x0(%r13),%rax
movss  %xmm0,0x4(%rax,%rbp,4)
cmp    %rbp,0x8(%rsp)
je     11b3 <IntegratePhase+0xa3>
mov    (%rbx),%rax
add    $0x1,%rbp
jmp    1150 <IntegratePhase+0x40>
cmp    $0x2,%r14
jl     126c <IntegratePhase+0x15c>
mov    $0x1,%ebp
mov    %r14,0x10(%rsp)
jmp    11e2 <IntegratePhase+0xd2>
nopl   0x0(%rax)
add    $0x1,%rbp
mov    0x10(%rsp),%r14
cmp    %r14,%rbp
je     126c <IntegratePhase+0x15c>
test   %r15,%r15
jle    11d0 <IntegratePhase+0xc0>
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
mov    -0x8(%r13,%rbp,8),%rax
movss  (%rax,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%rsp)
mov    -0x8(%rbx,%rbp,8),%rax
mov    (%rbx,%rbp,8),%rcx
movss  (%rcx,%r14,4),%xmm0
movss  (%rax,%r14,4),%xmm1
call   1030 <ModDiff@plt>
movsd  0xde2(%rip),%xmm2
cwtl
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x8(%r12,%rbp,8),%rax
movswl (%rax,%r14,2),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
addsd  (%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
mov    0x0(%r13,%rbp,8),%rax
movss  %xmm0,(%rax,%r14,4)
add    $0x1,%r14
cmp    %r14,%r15
jne    11f0 <IntegratePhase+0xe0>
jmp    11d0 <IntegratePhase+0xc0>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @IntegratePhase(float** nocapture %0, float** nocapture %1, i16** nocapture %2, i64 %3, i64 %4) {
%6 = add nsw i64 %3, -1
%7 = getelementptr inbounds i16*, i16** %2, i64 %6
%8 = load float*, float** %0
%9 = load float, float* %8
%10 = load float*, float** %1
store float %9, float* %10
%11 = icmp sgt i64 %4, 1
br i1 %11, label %16, label %12
12: ; preds = %16, %5
%13 = icmp sgt i64 %3, 1
br i1 %13, label %14, label %79
14: ; preds = %12
%15 = icmp sgt i64 %4, 0
br label %41
16: ; preds = %5, %39
%17 = phi float* [ %40, %39 ], [ %8, %5 ]
%18 = phi float [ %34, %39 ], [ %9, %5 ]
%19 = phi i64 [ %37, %39 ], [ 1, %5 ]
%20 = add nsw i64 %19, -1
%21 = fpext float %18 to double
%22 = getelementptr inbounds float, float* %17, i64 %19
%23 = load float, float* %22
%24 = getelementptr inbounds float, float* %17, i64 %20
%25 = load float, float* %24
%26 = tail call signext i16 @ModDiff(float %23, float %25)
%27 = sitofp i16 %26 to double
%28 = load i16*, i16** %7
%29 = getelementptr inbounds i16, i16* %28, i64 %20
%30 = load i16, i16* %29
%31 = sitofp i16 %30 to double
%32 = tail call double @llvm.fmuladd.f64(double %31, double 0x401921FB54442D18, double %27)
%33 = fadd double %32, %21
%34 = fptrunc double %33 to float
%35 = load float*, float** %1
%36 = getelementptr inbounds float, float* %35, i64 %19
store float %34, float* %36
%37 = add nuw nsw i64 %19, 1
%38 = icmp eq i64 %37, %4
br i1 %38, label %12, label %39
39: ; preds = %16
%40 = load float*, float** %0
br label %16
41: ; preds = %14, %76
%42 = phi i64 [ 1, %14 ], [ %77, %76 ]
br i1 %15, label %43, label %76
43: ; preds = %41
%44 = add nsw i64 %42, -1
%45 = getelementptr inbounds float*, float** %1, i64 %44
%46 = getelementptr inbounds float*, float** %0, i64 %42
%47 = getelementptr inbounds float*, float** %0, i64 %44
%48 = getelementptr inbounds i16*, i16** %2, i64 %44
%49 = getelementptr inbounds float*, float** %1, i64 %42
br label %50
50: ; preds = %43, %50
%51 = phi i64 [ 0, %43 ], [ %74, %50 ]
%52 = load float*, float** %45
%53 = getelementptr inbounds float, float* %52, i64 %51
%54 = load float, float* %53
%55 = fpext float %54 to double
%56 = load float*, float** %46
%57 = getelementptr inbounds float, float* %56, i64 %51
%58 = load float, float* %57
%59 = load float*, float** %47
%60 = getelementptr inbounds float, float* %59, i64 %51
%61 = load float, float* %60
%62 = tail call signext i16 @ModDiff(float %58, float %61)
%63 = sitofp i16 %62 to double
%64 = load i16*, i16** %48
%65 = getelementptr inbounds i16, i16* %64, i64 %51
%66 = load i16, i16* %65
%67 = sitofp i16 %66 to double
%68 = fneg double %67
%69 = tail call double @llvm.fmuladd.f64(double %68, double 0x401921FB54442D18, double %63)
%70 = fadd double %69, %55
%71 = fptrunc double %70 to float
%72 = load float*, float** %49
%73 = getelementptr inbounds float, float* %72, i64 %51
store float %71, float* %73
%74 = add nuw nsw i64 %51, 1
%75 = icmp eq i64 %74, %4
br i1 %75, label %76, label %50
76: ; preds = %50, %41
%77 = add nuw nsw i64 %42, 1
%78 = icmp eq i64 %77, %3
br i1 %78, label %79, label %41
79: ; preds = %76, %12
ret void
}