#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*/ okarea2(int, int, int, int);
scalar_t__ pow(float, int);
float chsq2d(int idev, float *x, int ilowb, int ihighb, int ilowa, int ihigha,
             float avg, int (*read2d)(void *, int, float *, int, int),
             void *userdata) {
  int i, j;
  float svd = 0.0;
  if (!okarea2(ilowb, ihighb, ilowa, ihigha))
    return svd;
  for (i = ilowb; i <= ihighb; i++) {
    read2d(userdata, idev, x, ihigha, i);
    for (j = ilowa; j <= ihigha; j++)
      svd += pow(x[j - 1] - avg, 2);
  }
  return svd;
}


x86 assembly (clang, -O0):
<chsq2d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%edi
mov    -0x18(%rbp),%esi
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%ecx
call   1030 <okarea2@plt>
cmp    $0x0,%eax
jne    117b <chsq2d+0x5b>
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1212 <chsq2d+0xf2>
mov    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1208 <chsq2d+0xe8>
mov    0x10(%rbp),%rax
mov    0x18(%rbp),%rdi
mov    -0x8(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x20(%rbp),%ecx
mov    -0x28(%rbp),%r8d
call   *%rax
mov    -0x1c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     11f5 <chsq2d+0xd5>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x24(%rbp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
cvtsi2ss %rax,%xmm0
addss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ab <chsq2d+0x8b>
jmp    11fa <chsq2d+0xda>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1181 <chsq2d+0x61>
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<chsq2d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
movss  %xmm0,0x8(%rsp)
mov    %r9d,%r13d
mov    %r8d,%r15d
mov    %ecx,%r14d
mov    %edx,%ebx
mov    %rsi,%rbp
mov    %edi,0x14(%rsp)
mov    %edx,%edi
mov    %ecx,%esi
mov    %r8d,%edx
mov    %r9d,%ecx
call   1030 <okarea2@plt>
xorps  %xmm0,%xmm0
test   %eax,%eax
je     12a8 <chsq2d+0x188>
cmp    %r14d,%ebx
jg     12a8 <chsq2d+0x188>
movslq %r15d,%rcx
lea    0x1(%r13),%r12d
mov    %r12d,%eax
sub    %r15d,%eax
mov    %eax,0xc(%rsp)
mov    %rcx,%rax
mov    %rcx,0x18(%rsp)
lea    0x1(%rcx),%rax
mov    %rax,0x20(%rsp)
xorps  %xmm0,%xmm0
mov    %r14d,0x10(%rsp)
jmp    11b5 <chsq2d+0x95>
cs nopw 0x0(%rax,%rax,1)
nop
movss  0x4(%rsp),%xmm0
lea    0x1(%rbx),%eax
cmp    0x10(%rsp),%ebx
mov    %eax,%ebx
je     12a8 <chsq2d+0x188>
movss  %xmm0,0x4(%rsp)
mov    0x68(%rsp),%rdi
mov    0x14(%rsp),%esi
mov    %rbp,%rdx
mov    %r13d,%ecx
mov    %ebx,%r8d
call   *0x60(%rsp)
cmp    %r15d,%r13d
jl     11a0 <chsq2d+0x80>
testb  $0x1,0xc(%rsp)
jne    11f0 <chsq2d+0xd0>
mov    0x18(%rsp),%r14
movss  0x4(%rsp),%xmm0
jmp    1225 <chsq2d+0x105>
nopw   0x0(%rax,%rax,1)
mov    0x18(%rsp),%rax
movss  -0x4(%rbp,%rax,4),%xmm0
subss  0x8(%rsp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  0x4(%rsp),%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,%xmm0
mov    0x20(%rsp),%r14
cmp    %r15d,%r13d
movss  0x8(%rsp),%xmm1
je     11a6 <chsq2d+0x86>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  %xmm0,0x4(%rsp)
movss  -0x4(%rbp,%r14,4),%xmm0
subss  %xmm1,%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
addss  0x4(%rsp),%xmm0
movss  %xmm0,0x4(%rsp)
movss  0x0(%rbp,%r14,4),%xmm0
subss  0x8(%rsp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
movss  0x8(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
addss  0x4(%rsp),%xmm0
add    $0x2,%r14
cmp    %r14d,%r12d
jne    1240 <chsq2d+0x120>
jmp    11a6 <chsq2d+0x86>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @chsq2d(i32 %0, float* %1, i32 %2, i32 %3, i32 %4, i32 %5, float %6, i32 (i8*, i32, float*, i32, i32)* nocapture %7, i8* %8) {
%10 = tail call i32 @okarea2(i32 %2, i32 %3, i32 %4, i32 %5)
%11 = icmp eq i32 %10, 0
%12 = icmp sgt i32 %2, %3
%13 = or i1 %11, %12
br i1 %13, label %63, label %14
14: ; preds = %9
%15 = icmp slt i32 %5, %4
%16 = sext i32 %4 to i64
%17 = add i32 %5, 1
%18 = sub i32 %17, %4
%19 = and i32 %18, 1
%20 = icmp eq i32 %19, 0
%21 = add nsw i64 %16, -1
%22 = getelementptr inbounds float, float* %1, i64 %21
%23 = add nsw i64 %16, 1
%24 = icmp eq i32 %5, %4
br label %25
25: ; preds = %14, %59
%26 = phi float [ 0.000000e+00, %14 ], [ %60, %59 ]
%27 = phi i32 [ %2, %14 ], [ %61, %59 ]
%28 = tail call i32 %7(i8* %8, i32 %0, float* %1, i32 %5, i32 %27)
br i1 %15, label %59, label %29
29: ; preds = %25
br i1 %20, label %36, label %30
30: ; preds = %29
%31 = load float, float* %22
%32 = fsub float %31, %6
%33 = tail call i64 @pow(float %32, i32 2)
%34 = sitofp i64 %33 to float
%35 = fadd float %26, %34
br label %36
36: ; preds = %30, %29
%37 = phi float [ undef, %29 ], [ %35, %30 ]
%38 = phi i64 [ %16, %29 ], [ %23, %30 ]
%39 = phi float [ %26, %29 ], [ %35, %30 ]
br i1 %24, label %59, label %40
40: ; preds = %36, %40
%41 = phi i64 [ %56, %40 ], [ %38, %36 ]
%42 = phi float [ %55, %40 ], [ %39, %36 ]
%43 = add nsw i64 %41, -1
%44 = getelementptr inbounds float, float* %1, i64 %43
%45 = load float, float* %44
%46 = fsub float %45, %6
%47 = tail call i64 @pow(float %46, i32 2)
%48 = sitofp i64 %47 to float
%49 = fadd float %42, %48
%50 = getelementptr inbounds float, float* %1, i64 %41
%51 = load float, float* %50
%52 = fsub float %51, %6
%53 = tail call i64 @pow(float %52, i32 2)
%54 = sitofp i64 %53 to float
%55 = fadd float %49, %54
%56 = add nsw i64 %41, 2
%57 = trunc i64 %56 to i32
%58 = icmp eq i32 %17, %57
br i1 %58, label %59, label %40
59: ; preds = %36, %40, %25
%60 = phi float [ %26, %25 ], [ %37, %36 ], [ %55, %40 ]
%61 = add i32 %27, 1
%62 = icmp eq i32 %27, %3
br i1 %62, label %63, label %25
63: ; preds = %59, %9
%64 = phi float [ 0.000000e+00, %9 ], [ %60, %59 ]
ret float %64
}