#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*/ copyshrinkrange(double *, double *, int, int, int, int);
int /*<<< orphan*/ cov(double *, double *, int, int);
int /*<<< orphan*/ cov2cor(double *, double *, int);
double fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, long);
int /*<<< orphan*/ free(double *);
double *malloc(int);
int /*<<< orphan*/ stderr;
int thin(double *x, int n, int p, int *active, double cormax, int windowsize,
         int stepsize) {
  int j, i, k;
  double *S = NULL, *P = NULL, *xwin = NULL;
  int ws = (windowsize < p) ? windowsize : p, ws2 = ws * ws;
  stepsize = (stepsize > windowsize) ? windowsize : stepsize;
  j = 0;
  while (j < p && ws > 1) {
    if (!(xwin = malloc(sizeof(double) * ws * n))) {
      fprintf(stderr, "can't malloc %ld bytes\n",
              (long)sizeof(double) * ws * n);
      return 0;
    };
    copyshrinkrange(x, xwin, n, p, j, j + ws);
    if (!(S = malloc(sizeof(double) * ws2))) {
      fprintf(stderr, "can't malloc %ld bytes\n", (long)sizeof(double) * ws2);
      return 0;
    };
    if (!(P = malloc(sizeof(double) * ws2))) {
      fprintf(stderr, "can't malloc %ld bytes\n", (long)sizeof(double) * ws2);
      return 0;
    };
    cov(xwin, S, n, ws);
    cov2cor(S, P, ws);
    if (xwin) {
      free(xwin);
      xwin = NULL;
    };
    if (S) {
      free(S);
      S = NULL;
    };
    for (i = 1; i < ws; i++) {
      for (k = 0; k < i; k++) {
        if (fabs(P[i * ws + k]) > cormax && active[i + j] && active[k + j]) {
          active[i + j] = 0;
        }
      }
    }
    if (P) {
      free(P);
      P = NULL;
    };
    stepsize = (stepsize < p - j) ? stepsize : p - j;
    j += stepsize;
    ws = (windowsize < p - j) ? windowsize : p - j;
    ws2 = ws * ws;
  }
  return 1;
}


x86 assembly (clang, -O0):
<thin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
mov    %r9d,-0x30(%rbp)
movq   $0x0,-0x48(%rbp)
movq   $0x0,-0x50(%rbp)
movq   $0x0,-0x58(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11b5 <thin+0x55>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x64(%rbp)
jmp    11bb <thin+0x5b>
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
imul   -0x5c(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jle    11e2 <thin+0x82>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x68(%rbp)
jmp    11e8 <thin+0x88>
mov    -0x30(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%ecx
xor    %eax,%eax
cmp    -0x18(%rbp),%ecx
mov    %al,-0x69(%rbp)
jge    1210 <thin+0xb0>
cmpl   $0x1,-0x5c(%rbp)
setg   %al
mov    %al,-0x69(%rbp)
mov    -0x69(%rbp),%al
test   $0x1,%al
jne    1220 <thin+0xc0>
jmp    14c0 <thin+0x360>
movslq -0x5c(%rbp),%rax
shl    $0x3,%rax
movslq -0x14(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x58(%rbp)
cmp    $0x0,%rax
jne    1276 <thin+0x116>
mov    0x2d94(%rip),%rax
mov    (%rax),%edi
movslq -0x5c(%rbp),%rdx
shl    $0x3,%rdx
movslq -0x14(%rbp),%rax
imul   %rax,%rdx
lea    0xdab(%rip),%rsi
call   1070 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    14c7 <thin+0x367>
mov    -0x10(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x34(%rbp),%r8d
mov    -0x34(%rbp),%r9d
add    -0x5c(%rbp),%r9d
call   1040 <copyshrinkrange@plt>
movslq -0x60(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x48(%rbp)
cmp    $0x0,%rax
jne    12db <thin+0x17b>
mov    0x2d27(%rip),%rax
mov    (%rax),%edi
movslq -0x60(%rbp),%rdx
shl    $0x3,%rdx
lea    0xd46(%rip),%rsi
call   1070 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    14c7 <thin+0x367>
movslq -0x60(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x50(%rbp)
cmp    $0x0,%rax
jne    1321 <thin+0x1c1>
mov    0x2ce1(%rip),%rax
mov    (%rax),%edi
movslq -0x60(%rbp),%rdx
shl    $0x3,%rdx
lea    0xd00(%rip),%rsi
call   1070 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    14c7 <thin+0x367>
mov    -0x58(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x5c(%rbp),%ecx
call   1050 <cov@plt>
mov    -0x48(%rbp),%rdi
mov    -0x50(%rbp),%rsi
mov    -0x5c(%rbp),%edx
call   1060 <cov2cor@plt>
cmpq   $0x0,-0x58(%rbp)
je     1360 <thin+0x200>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x58(%rbp)
cmpq   $0x0,-0x48(%rbp)
je     137c <thin+0x21c>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x48(%rbp)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jge    1436 <thin+0x2d6>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1423 <thin+0x2c3>
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x5c(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xc41(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    1410 <thin+0x2b0>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1410 <thin+0x2b0>
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1410 <thin+0x2b0>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    1415 <thin+0x2b5>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1396 <thin+0x236>
jmp    1428 <thin+0x2c8>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1383 <thin+0x223>
cmpq   $0x0,-0x50(%rbp)
je     1452 <thin+0x2f2>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x50(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    146e <thin+0x30e>
mov    -0x30(%rbp),%eax
mov    %eax,-0x70(%rbp)
jmp    1477 <thin+0x317>
mov    -0x18(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    14a2 <thin+0x342>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x74(%rbp)
jmp    14ab <thin+0x34b>
mov    -0x18(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
imul   -0x5c(%rbp),%eax
mov    %eax,-0x60(%rbp)
jmp    11f5 <thin+0x95>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<thin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%rbx
mov    %esi,0xc(%rsp)
mov    %rdi,0x30(%rsp)
cmp    %edx,%r8d
mov    %edx,%ecx
mov    %r8d,0x10(%rsp)
cmovl  %r8d,%ecx
mov    $0x1,%eax
mov    %edx,0x14(%rsp)
test   %edx,%edx
jle    1364 <thin+0x204>
cmp    $0x2,%ecx
jl     1364 <thin+0x204>
mov    %r9d,%ebp
movsd  %xmm0,0x20(%rsp)
mov    0x10(%rsp),%eax
cmp    %eax,%r9d
cmovg  %eax,%ebp
mov    %ecx,%r14d
imul   %ecx,%r14d
movslq 0xc(%rsp),%rax
mov    %rax,0x28(%rsp)
xor    %eax,%eax
mov    %rax,0x18(%rsp)
nopl   (%rax)
mov    %ecx,%r15d
mov    %r15,%r13
imul   0x28(%rsp),%r13
shl    $0x3,%r13
mov    %r13d,%edi
call   1080 <malloc@plt>
test   %rax,%rax
je     134a <thin+0x1ea>
mov    %rax,%r12
mov    0x18(%rsp),%r8
lea    (%r15,%r8,1),%r9d
mov    0x30(%rsp),%rdi
mov    %rax,%rsi
mov    0xc(%rsp),%edx
mov    0x14(%rsp),%ecx
call   1040 <copyshrinkrange@plt>
mov    %r14d,%r13d
shl    $0x3,%r13
mov    %r13d,%edi
call   1080 <malloc@plt>
test   %rax,%rax
je     134a <thin+0x1ea>
mov    %rax,%r14
mov    %r13d,%edi
call   1080 <malloc@plt>
test   %rax,%rax
je     134a <thin+0x1ea>
mov    %ebp,%r13d
mov    %r12,%rdi
mov    %r14,%rsi
mov    0xc(%rsp),%edx
mov    %r15d,%ecx
mov    %rax,%rbp
call   1050 <cov@plt>
mov    %r14,%rdi
mov    %rbp,%rsi
mov    %r15d,%edx
call   1060 <cov2cor@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
movslq 0x18(%rsp),%r8
lea    (%rbx,%r8,4),%rcx
lea    0x0(,%r15,8),%rdx
add    %rbp,%rdx
lea    0x0(,%r15,8),%r9
mov    $0x1,%esi
movsd  0x20(%rsp),%xmm1
movapd 0xd5a(%rip),%xmm2
jmp    12bc <thin+0x15c>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rsi
add    %r9,%rdx
cmp    %r15,%rsi
je     1300 <thin+0x1a0>
lea    (%rsi,%r8,1),%rbp
xor    %eax,%eax
jmp    12d9 <thin+0x179>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rax
cmp    %rax,%rsi
je     12b0 <thin+0x150>
movsd  (%rdx,%rax,8),%xmm0
andpd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    12d0 <thin+0x170>
cmpl   $0x0,(%rbx,%rbp,4)
je     12d0 <thin+0x170>
cmpl   $0x0,(%rcx,%rax,4)
je     12d0 <thin+0x170>
movl   $0x0,(%rbx,%rbp,4)
jmp    12d0 <thin+0x170>
nopl   (%rax)
call   1030 <free@plt>
mov    0x14(%rsp),%edx
mov    %edx,%eax
mov    0x18(%rsp),%rsi
sub    %esi,%eax
mov    %r13d,%ebp
cmp    %eax,%r13d
cmovge %eax,%ebp
add    %ebp,%esi
mov    %edx,%ecx
sub    %esi,%ecx
mov    0x10(%rsp),%eax
cmp    %eax,%ecx
cmovg  %eax,%ecx
mov    %rsi,0x18(%rsp)
cmp    %esi,%edx
jle    1343 <thin+0x1e3>
mov    %ecx,%r14d
imul   %ecx,%r14d
cmp    $0x1,%ecx
jg     11d0 <thin+0x70>
mov    $0x1,%eax
jmp    1364 <thin+0x204>
mov    0x2c8f(%rip),%rax
mov    (%rax),%edi
lea    0xcb6(%rip),%rsi
mov    %r13,%rdx
call   1070 <fprintf@plt>
xor    %eax,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [24 x i8] c"can't malloc %ld bytes\0A\00"
define i32 @thin(double* %0, i32 %1, i32 %2, i32* %3, double %4, i32 %5, i32 %6) {
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32*
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca double*
%20 = alloca double*
%21 = alloca double*
%22 = alloca i32
%23 = alloca i32
store double* %0, double** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32* %3, i32** %12
store double %4, double* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store double* null, double** %19
store double* null, double** %20
store double* null, double** %21
%24 = load i32, i32* %14
%25 = load i32, i32* %11
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %29
27: ; preds = %7
%28 = load i32, i32* %14
br label %31
29: ; preds = %7
%30 = load i32, i32* %11
br label %31
31: ; preds = %29, %27
%32 = phi i32 [ %28, %27 ], [ %30, %29 ]
store i32 %32, i32* %22
%33 = load i32, i32* %22
%34 = load i32, i32* %22
%35 = mul nsw i32 %33, %34
store i32 %35, i32* %23
%36 = load i32, i32* %15
%37 = load i32, i32* %14
%38 = icmp sgt i32 %36, %37
br i1 %38, label %39, label %41
39: ; preds = %31
%40 = load i32, i32* %14
br label %43
41: ; preds = %31
%42 = load i32, i32* %15
br label %43
43: ; preds = %41, %39
%44 = phi i32 [ %40, %39 ], [ %42, %41 ]
store i32 %44, i32* %15
store i32 0, i32* %16
br label %45
45: ; preds = %219, %43
%46 = load i32, i32* %16
%47 = load i32, i32* %11
%48 = icmp slt i32 %46, %47
br i1 %48, label %49, label %52
49: ; preds = %45
%50 = load i32, i32* %22
%51 = icmp sgt i32 %50, 1
br label %52
52: ; preds = %49, %45
%53 = phi i1 [ false, %45 ], [ %51, %49 ]
br i1 %53, label %54, label %224
54: ; preds = %52
%55 = load i32, i32* %22
%56 = sext i32 %55 to i64
%57 = mul i64 8, %56
%58 = load i32, i32* %10
%59 = sext i32 %58 to i64
%60 = mul i64 %57, %59
%61 = trunc i64 %60 to i32
%62 = call double* @malloc(i32 %61)
store double* %62, double** %21
%63 = icmp ne double* %62, null
br i1 %63, label %73, label %64
64: ; preds = %54
%65 = load i32, i32* @stderr
%66 = load i32, i32* %22
%67 = sext i32 %66 to i64
%68 = mul nsw i64 8, %67
%69 = load i32, i32* %10
%70 = sext i32 %69 to i64
%71 = mul nsw i64 %68, %70
%72 = call i32 @fprintf(i32 %65, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %71)
store i32 0, i32* %8
br label %225
73: ; preds = %54
%74 = load double*, double** %9
%75 = load double*, double** %21
%76 = load i32, i32* %10
%77 = load i32, i32* %11
%78 = load i32, i32* %16
%79 = load i32, i32* %16
%80 = load i32, i32* %22
%81 = add nsw i32 %79, %80
%82 = call i32 @copyshrinkrange(double* %74, double* %75, i32 %76, i32 %77, i32 %78, i32 %81)
%83 = load i32, i32* %23
%84 = sext i32 %83 to i64
%85 = mul i64 8, %84
%86 = trunc i64 %85 to i32
%87 = call double* @malloc(i32 %86)
store double* %87, double** %19
%88 = icmp ne double* %87, null
br i1 %88, label %95, label %89
89: ; preds = %73
%90 = load i32, i32* @stderr
%91 = load i32, i32* %23
%92 = sext i32 %91 to i64
%93 = mul nsw i64 8, %92
%94 = call i32 @fprintf(i32 %90, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %93)
store i32 0, i32* %8
br label %225
95: ; preds = %73
%96 = load i32, i32* %23
%97 = sext i32 %96 to i64
%98 = mul i64 8, %97
%99 = trunc i64 %98 to i32
%100 = call double* @malloc(i32 %99)
store double* %100, double** %20
%101 = icmp ne double* %100, null
br i1 %101, label %108, label %102
102: ; preds = %95
%103 = load i32, i32* @stderr
%104 = load i32, i32* %23
%105 = sext i32 %104 to i64
%106 = mul nsw i64 8, %105
%107 = call i32 @fprintf(i32 %103, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %106)
store i32 0, i32* %8
br label %225
108: ; preds = %95
%109 = load double*, double** %21
%110 = load double*, double** %19
%111 = load i32, i32* %10
%112 = load i32, i32* %22
%113 = call i32 @cov(double* %109, double* %110, i32 %111, i32 %112)
%114 = load double*, double** %19
%115 = load double*, double** %20
%116 = load i32, i32* %22
%117 = call i32 @cov2cor(double* %114, double* %115, i32 %116)
%118 = load double*, double** %21
%119 = icmp ne double* %118, null
br i1 %119, label %120, label %123
120: ; preds = %108
%121 = load double*, double** %21
%122 = call i32 @free(double* %121)
store double* null, double** %21
br label %123
123: ; preds = %120, %108
%124 = load double*, double** %19
%125 = icmp ne double* %124, null
br i1 %125, label %126, label %129
126: ; preds = %123
%127 = load double*, double** %19
%128 = call i32 @free(double* %127)
store double* null, double** %19
br label %129
129: ; preds = %126, %123
store i32 1, i32* %17
br label %130
130: ; preds = %182, %129
%131 = load i32, i32* %17
%132 = load i32, i32* %22
%133 = icmp slt i32 %131, %132
br i1 %133, label %134, label %185
134: ; preds = %130
store i32 0, i32* %18
br label %135
135: ; preds = %178, %134
%136 = load i32, i32* %18
%137 = load i32, i32* %17
%138 = icmp slt i32 %136, %137
br i1 %138, label %139, label %181
139: ; preds = %135
%140 = load double*, double** %20
%141 = load i32, i32* %17
%142 = load i32, i32* %22
%143 = mul nsw i32 %141, %142
%144 = load i32, i32* %18
%145 = add nsw i32 %143, %144
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds double, double* %140, i64 %146
%148 = load double, double* %147
%149 = call double @llvm.fabs.f64(double %148)
%150 = load double, double* %13
%151 = fcmp ogt double %149, %150
br i1 %151, label %152, label %177
152: ; preds = %139
%153 = load i32*, i32** %12
%154 = load i32, i32* %17
%155 = load i32, i32* %16
%156 = add nsw i32 %154, %155
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds i32, i32* %153, i64 %157
%159 = load i32, i32* %158
%160 = icmp ne i32 %159, 0
br i1 %160, label %161, label %177
161: ; preds = %152
%162 = load i32*, i32** %12
%163 = load i32, i32* %18
%164 = load i32, i32* %16
%165 = add nsw i32 %163, %164
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds i32, i32* %162, i64 %166
%168 = load i32, i32* %167
%169 = icmp ne i32 %168, 0
br i1 %169, label %170, label %177
170: ; preds = %161
%171 = load i32*, i32** %12
%172 = load i32, i32* %17
%173 = load i32, i32* %16
%174 = add nsw i32 %172, %173
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds i32, i32* %171, i64 %175
store i32 0, i32* %176
br label %177
177: ; preds = %170, %161, %152, %139
br label %178
178: ; preds = %177
%179 = load i32, i32* %18
%180 = add nsw i32 %179, 1
store i32 %180, i32* %18
br label %135
181: ; preds = %135
br label %182
182: ; preds = %181
%183 = load i32, i32* %17
%184 = add nsw i32 %183, 1
store i32 %184, i32* %17
br label %130
185: ; preds = %130
%186 = load double*, double** %20
%187 = icmp ne double* %186, null
br i1 %187, label %188, label %191
188: ; preds = %185
%189 = load double*, double** %20
%190 = call i32 @free(double* %189)
store double* null, double** %20
br label %191
191: ; preds = %188, %185
%192 = load i32, i32* %15
%193 = load i32, i32* %11
%194 = load i32, i32* %16
%195 = sub nsw i32 %193, %194
%196 = icmp slt i32 %192, %195
br i1 %196, label %197, label %199
197: ; preds = %191
%198 = load i32, i32* %15
br label %203
199: ; preds = %191
%200 = load i32, i32* %11
%201 = load i32, i32* %16
%202 = sub nsw i32 %200, %201
br label %203
203: ; preds = %199, %197
%204 = phi i32 [ %198, %197 ], [ %202, %199 ]
store i32 %204, i32* %15
%205 = load i32, i32* %15
%206 = load i32, i32* %16
%207 = add nsw i32 %206, %205
store i32 %207, i32* %16
%208 = load i32, i32* %14
%209 = load i32, i32* %11
%210 = load i32, i32* %16
%211 = sub nsw i32 %209, %210
%212 = icmp slt i32 %208, %211
br i1 %212, label %213, label %215
213: ; preds = %203
%214 = load i32, i32* %14
br label %219
215: ; preds = %203
%216 = load i32, i32* %11
%217 = load i32, i32* %16
%218 = sub nsw i32 %216, %217
br label %219
219: ; preds = %215, %213
%220 = phi i32 [ %214, %213 ], [ %218, %215 ]
store i32 %220, i32* %22
%221 = load i32, i32* %22
%222 = load i32, i32* %22
%223 = mul nsw i32 %221, %222
store i32 %223, i32* %23
br label %45
224: ; preds = %52
store i32 1, i32* %8
br label %225
225: ; preds = %224, %102, %89, %64
%226 = load i32, i32* %8
ret i32 %226
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [24 x i8] c"can't malloc %ld bytes\0A\00"
define i32 @thin(double* %0, i32 %1, i32 %2, i32* nocapture %3, double %4, i32 %5, i32 %6) {
%8 = icmp slt i32 %5, %2
%9 = select i1 %8, i32 %5, i32 %2
%10 = icmp sgt i32 %2, 0
%11 = icmp sgt i32 %9, 1
%12 = and i1 %10, %11
br i1 %12, label %13, label %91
13: ; preds = %7
%14 = icmp sgt i32 %6, %5
%15 = select i1 %14, i32 %5, i32 %6
%16 = mul nsw i32 %9, %9
%17 = sext i32 %1 to i64
br label %18
18: ; preds = %13, %74
%19 = phi i32 [ %16, %13 ], [ %83, %74 ]
%20 = phi i32 [ %9, %13 ], [ %82, %74 ]
%21 = phi i32 [ 0, %13 ], [ %79, %74 ]
%22 = phi i32 [ %15, %13 ], [ %78, %74 ]
%23 = zext i32 %20 to i64
%24 = shl nuw nsw i64 %23, 3
%25 = mul i64 %24, %17
%26 = trunc i64 %25 to i32
%27 = tail call double* @malloc(i32 %26)
%28 = icmp eq double* %27, null
br i1 %28, label %87, label %29
29: ; preds = %18
%30 = add nsw i32 %20, %21
%31 = tail call i32 @copyshrinkrange(double* %0, double* %27, i32 %1, i32 %2, i32 %21, i32 %30)
%32 = zext i32 %19 to i64
%33 = shl nuw nsw i64 %32, 3
%34 = trunc i64 %33 to i32
%35 = tail call double* @malloc(i32 %34)
%36 = icmp eq double* %35, null
br i1 %36, label %87, label %37
37: ; preds = %29
%38 = tail call double* @malloc(i32 %34)
%39 = icmp eq double* %38, null
br i1 %39, label %87, label %40
40: ; preds = %37
%41 = tail call i32 @cov(double* %27, double* %35, i32 %1, i32 %20)
%42 = tail call i32 @cov2cor(double* %35, double* %38, i32 %20)
%43 = tail call i32 @free(double* %27)
%44 = tail call i32 @free(double* %35)
%45 = sext i32 %21 to i64
%46 = zext i32 %20 to i64
br label %47
47: ; preds = %40, %71
%48 = phi i64 [ 1, %40 ], [ %72, %71 ]
%49 = mul nsw i64 %48, %46
%50 = add nsw i64 %48, %45
%51 = getelementptr inbounds i32, i32* %3, i64 %50
br label %52
52: ; preds = %47, %68
%53 = phi i64 [ 0, %47 ], [ %69, %68 ]
%54 = add nuw nsw i64 %53, %49
%55 = getelementptr inbounds double, double* %38, i64 %54
%56 = load double, double* %55
%57 = tail call double @llvm.fabs.f64(double %56)
%58 = fcmp ogt double %57, %4
br i1 %58, label %59, label %68
59: ; preds = %52
%60 = load i32, i32* %51
%61 = icmp eq i32 %60, 0
br i1 %61, label %68, label %62
62: ; preds = %59
%63 = add nsw i64 %53, %45
%64 = getelementptr inbounds i32, i32* %3, i64 %63
%65 = load i32, i32* %64
%66 = icmp eq i32 %65, 0
br i1 %66, label %68, label %67
67: ; preds = %62
store i32 0, i32* %51
br label %68
68: ; preds = %52, %59, %62, %67
%69 = add nuw nsw i64 %53, 1
%70 = icmp eq i64 %69, %48
br i1 %70, label %71, label %52
71: ; preds = %68
%72 = add nuw nsw i64 %48, 1
%73 = icmp eq i64 %72, %23
br i1 %73, label %74, label %47
74: ; preds = %71
%75 = tail call i32 @free(double* %38)
%76 = sub nsw i32 %2, %21
%77 = icmp slt i32 %22, %76
%78 = select i1 %77, i32 %22, i32 %76
%79 = add nsw i32 %78, %21
%80 = sub nsw i32 %2, %79
%81 = icmp sgt i32 %80, %5
%82 = select i1 %81, i32 %5, i32 %80
%83 = mul nsw i32 %82, %82
%84 = icmp slt i32 %79, %2
%85 = icmp sgt i32 %82, 1
%86 = select i1 %84, i1 %85, i1 false
br i1 %86, label %18, label %91
87: ; preds = %37, %29, %18
%88 = phi i64 [ %25, %18 ], [ %33, %29 ], [ %33, %37 ]
%89 = load i32, i32* @stderr
%90 = tail call i32 @fprintf(i32 %89, i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i64 %88)
br label %91
91: ; preds = %74, %87, %7
%92 = phi i32 [ 1, %7 ], [ 0, %87 ], [ 1, %74 ]
ret i32 %92
}