#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 calc_max(double *, unsigned long);
double calc_min(double *, unsigned long);
int ceil(int);
int /*<<< orphan*/ do_vect(double **, unsigned long);
int /*<<< orphan*/ free(double *);
unsigned long pindex(int, int, int, int, unsigned int, unsigned int,
                     unsigned int);
void update_parameters(int slicemode, double **ppro, int *rangea, int *rangeb,
                       int *rangec, unsigned *extx, unsigned *exty,
                       unsigned *extz, double *origx, double *origy,
                       double *origz, int maxx, int minx, int maxy, int miny,
                       int maxz, int minz, double width, int *k, double *cutoff,
                       int *stride, double *pphi, unsigned long *nvox) {
  double maxdensity, mindensity;
  int i, j, l;
  unsigned long ind1, ind2;
  *extx = maxx - minx + 1;
  *exty = maxy - miny + 1;
  *extz = maxz - minz + 1;
  *origx += minx * width;
  *origy += miny * width;
  *origz += minz * width;
  switch (slicemode) {
  case 1:
    *rangea = *extx;
    *rangeb = *exty;
    *rangec = *extz;
    break;
  case 2:
    *rangea = *extz;
    *rangeb = *extx;
    *rangec = *exty;
    break;
  case 3:
    *rangea = *exty;
    *rangeb = *extz;
    *rangec = *extx;
    break;
  }
  *k = (*rangec) / 2;
  *nvox = (*extx) * (*exty) * (*extz);
  mindensity = calc_min(pphi, *nvox);
  maxdensity = calc_max(pphi, *nvox);
  *cutoff = 0.5 * (maxdensity + mindensity);
  *stride = ceil(*rangea / 100.0);
  if (*stride < ceil(*rangeb / 100.0))
    *stride = ceil(*rangeb / 100.0);
  if (*stride < 1)
    *stride = 1;
  free(*ppro);
  do_vect(ppro, (unsigned long)(*rangea * *rangeb));
  for (i = 0; i < *rangea; i++)
    for (j = 0; j < *rangeb; j++) {
      ind2 = i + j * *rangea;
      for (l = 0; l < *rangec; l++) {
        ind1 = pindex(slicemode, i, j, l, *extx, *exty, *extz);
        *(*ppro + ind2) += *(pphi + ind1);
      }
      *(*ppro + ind2) /= ((double)*rangec);
    }
}


x86 assembly (clang, -O0):
<update_parameters>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x88(%rbp),%rax
mov    0x80(%rbp),%rax
mov    0x78(%rbp),%rax
mov    0x70(%rbp),%rax
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%eax
mov    0x58(%rbp),%eax
mov    0x50(%rbp),%eax
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
mov    0x38(%rbp),%ecx
mov    0x40(%rbp),%eax
sub    %eax,%ecx
add    $0x1,%ecx
mov    -0x30(%rbp),%rax
mov    %ecx,(%rax)
mov    0x48(%rbp),%ecx
mov    0x50(%rbp),%eax
sub    %eax,%ecx
add    $0x1,%ecx
mov    0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    0x58(%rbp),%ecx
mov    0x60(%rbp),%eax
sub    %eax,%ecx
add    $0x1,%ecx
mov    0x18(%rbp),%rax
mov    %ecx,(%rax)
cvtsi2sdl 0x40(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl 0x50(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl 0x60(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
mov    %eax,-0x6c(%rbp)
sub    $0x1,%eax
je     128a <update_parameters+0x12a>
jmp    1268 <update_parameters+0x108>
mov    -0x6c(%rbp),%eax
sub    $0x2,%eax
je     12b3 <update_parameters+0x153>
jmp    1279 <update_parameters+0x119>
mov    -0x6c(%rbp),%eax
sub    $0x3,%eax
je     12dc <update_parameters+0x17c>
jmp    1300 <update_parameters+0x1a0>
mov    -0x30(%rbp),%rax
mov    (%rax),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    0x10(%rbp),%rax
mov    (%rax),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    0x18(%rbp),%rax
mov    (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
jmp    1300 <update_parameters+0x1a0>
mov    0x18(%rbp),%rax
mov    (%rax),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%rax
mov    (%rax),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    0x10(%rbp),%rax
mov    (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
jmp    1300 <update_parameters+0x1a0>
mov    0x10(%rbp),%rax
mov    (%rax),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    0x18(%rbp),%rax
mov    (%rax),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%rax
mov    (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    0x68(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
mov    0x10(%rbp),%rcx
imul   (%rcx),%eax
mov    0x18(%rbp),%rcx
imul   (%rcx),%eax
mov    %eax,%eax
mov    %eax,%ecx
mov    0x88(%rbp),%rax
mov    %rcx,(%rax)
mov    0x80(%rbp),%rdi
mov    0x88(%rbp),%rax
mov    (%rax),%rsi
call   1060 <calc_min@plt>
movsd  %xmm0,-0x48(%rbp)
mov    0x80(%rbp),%rdi
mov    0x88(%rbp),%rax
mov    (%rax),%rsi
call   1070 <calc_max@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm1
addsd  -0x48(%rbp),%xmm1
movsd  0xc88(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    0x70(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movsd  0xc64(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
mov    %eax,%ecx
mov    0x78(%rbp),%rax
mov    %ecx,(%rax)
mov    0x78(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x70(%rbp)
mov    -0x20(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movsd  0xc36(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
mov    %eax,%ecx
mov    -0x70(%rbp),%eax
cmp    %ecx,%eax
jge    1409 <update_parameters+0x2a9>
mov    -0x20(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movsd  0xc0c(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
mov    %eax,%ecx
mov    0x78(%rbp),%rax
mov    %ecx,(%rax)
mov    0x78(%rbp),%rax
cmpl   $0x1,(%rax)
jge    1420 <update_parameters+0x2c0>
mov    0x78(%rbp),%rax
movl   $0x1,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%rdi
call   1030 <free@plt>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    -0x20(%rbp),%rcx
imul   (%rcx),%eax
movslq %eax,%rsi
call   1080 <do_vect@plt>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    153d <update_parameters+0x3dd>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    152a <update_parameters+0x3ca>
mov    -0x4c(%rbp),%eax
mov    -0x50(%rbp),%ecx
mov    -0x18(%rbp),%rdx
imul   (%rdx),%ecx
add    %ecx,%eax
cltq
mov    %rax,-0x68(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jge    14fb <update_parameters+0x39b>
mov    -0x4(%rbp),%edi
mov    -0x4c(%rbp),%esi
mov    -0x50(%rbp),%edx
mov    -0x54(%rbp),%ecx
mov    -0x30(%rbp),%rax
mov    (%rax),%r8d
mov    0x10(%rbp),%rax
mov    (%rax),%r9d
mov    0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,(%rsp)
call   1040 <pindex@plt>
mov    %rax,-0x60(%rbp)
mov    0x80(%rbp),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x68(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    148d <update_parameters+0x32d>
mov    -0x28(%rbp),%rax
cvtsi2sdl (%rax),%xmm1
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1462 <update_parameters+0x302>
jmp    152f <update_parameters+0x3cf>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    144c <update_parameters+0x2ec>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<update_parameters>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%r12
mov    %rsi,0x8(%rsp)
mov    0xe8(%rsp),%r14
mov    0x78(%rsp),%r11
mov    0x70(%rsp),%rbx
mov    0x80(%rsp),%r15
mov    0xc0(%rsp),%r10d
mov    0xb0(%rsp),%esi
mov    0xa0(%rsp),%eax
mov    0x98(%rsp),%ecx
sub    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,(%r9)
mov    0xa8(%rsp),%ecx
sub    %esi,%ecx
add    $0x1,%ecx
mov    %ecx,(%rbx)
mov    0xb8(%rsp),%ecx
sub    %r10d,%ecx
add    $0x1,%ecx
cvtsi2sd %eax,%xmm1
mov    %ecx,(%r11)
mulsd  %xmm0,%xmm1
addsd  (%r15),%xmm1
movsd  %xmm1,(%r15)
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mov    0x88(%rsp),%rsi
mulsd  %xmm0,%xmm1
addsd  (%rsi),%xmm1
mov    0xd8(%rsp),%r15
movsd  %xmm1,(%rsi)
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
mov    0x90(%rsp),%rax
mulsd  %xmm0,%xmm1
addsd  (%rax),%xmm1
mov    0xd0(%rsp),%rbp
movsd  %xmm1,(%rax)
mov    0xc8(%rsp),%rax
cmp    $0x1,%edi
mov    %r8,0x20(%rsp)
mov    %edi,0x1c(%rsp)
je     1261 <update_parameters+0x101>
cmp    $0x2,%edi
je     1276 <update_parameters+0x116>
cmp    $0x3,%edi
mov    0x70(%rsp),%rsi
jne    128b <update_parameters+0x12b>
mov    (%rsi),%ecx
mov    %ecx,(%rdx)
mov    (%r11),%ecx
mov    %ecx,(%r12)
mov    (%r9),%ecx
jmp    1286 <update_parameters+0x126>
mov    (%r9),%ecx
mov    %ecx,(%rdx)
mov    0x70(%rsp),%rsi
mov    (%rsi),%ecx
mov    %ecx,(%r12)
mov    (%r11),%ecx
jmp    1286 <update_parameters+0x126>
mov    %ecx,(%rdx)
mov    (%r9),%ecx
mov    %ecx,(%r12)
mov    0x70(%rsp),%rsi
mov    (%rsi),%ecx
mov    %ecx,(%r8)
jmp    128e <update_parameters+0x12e>
mov    (%r8),%ecx
mov    %rdx,%r13
mov    %ecx,%edx
shr    $0x1f,%edx
add    %ecx,%edx
sar    %edx
mov    %edx,(%rax)
mov    (%rsi),%esi
mov    %r9,0x30(%rsp)
imul   (%r9),%esi
imul   (%r11),%esi
mov    %rsi,(%r14)
mov    0xe0(%rsp),%rbx
mov    %rbx,%rdi
call   1060 <calc_min@plt>
movsd  %xmm0,0x10(%rsp)
mov    (%r14),%rsi
mov    %rbx,%rdi
call   1070 <calc_max@plt>
addsd  0x10(%rsp),%xmm0
mulsd  0xd23(%rip),%xmm0
movsd  %xmm0,0x0(%rbp)
mov    %r13,0x10(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x0(%r13),%xmm0
divsd  0xd10(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
mov    %eax,%r14d
mov    %eax,(%r15)
xorps  %xmm0,%xmm0
cvtsi2sdl (%r12),%xmm0
divsd  0xcf0(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
cmp    %eax,%r14d
jge    1345 <update_parameters+0x1e5>
xorps  %xmm0,%xmm0
cvtsi2sdl (%r12),%xmm0
divsd  0xcd1(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <ceil@plt>
mov    %eax,(%r15)
jmp    1348 <update_parameters+0x1e8>
mov    (%r15),%eax
mov    0x8(%rsp),%rcx
test   %eax,%eax
jg     1358 <update_parameters+0x1f8>
movl   $0x1,(%r15)
mov    (%rcx),%rdi
call   1030 <free@plt>
mov    0x10(%rsp),%r14
movslq (%r14),%rax
movslq (%r12),%rsi
imul   %rax,%rsi
mov    0x8(%rsp),%rdi
call   1080 <do_vect@plt>
mov    0x8(%rsp),%rdi
mov    (%r14),%ecx
test   %ecx,%ecx
mov    0x20(%rsp),%rbp
jle    147e <update_parameters+0x31e>
mov    (%r12),%ebx
xor    %esi,%esi
mov    %r12,0x28(%rsp)
mov    0x30(%rsp),%r15
jmp    13bb <update_parameters+0x25b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%esi
cmp    %ecx,%esi
jge    147e <update_parameters+0x31e>
test   %ebx,%ebx
jle    13b0 <update_parameters+0x250>
mov    0x0(%rbp),%eax
xor    %edx,%edx
jmp    13f4 <update_parameters+0x294>
cs nopw 0x0(%rax,%rax,1)
mov    (%rdi),%rcx
movsd  (%rcx,%r13,8),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%r13,8)
add    $0x1,%edx
mov    (%r14),%ecx
cmp    %ebx,%edx
jge    13b0 <update_parameters+0x250>
imul   %edx,%ecx
add    %esi,%ecx
movslq %ecx,%r13
test   %eax,%eax
jle    13d0 <update_parameters+0x270>
xor    %r14d,%r14d
mov    0x1c(%rsp),%ebx
nopw   0x0(%rax,%rax,1)
mov    (%r15),%r8d
mov    0x70(%rsp),%rax
mov    (%rax),%r9d
mov    0x78(%rsp),%rax
mov    (%rax),%eax
mov    %eax,(%rsp)
mov    %ebx,%edi
mov    %r14d,%ecx
mov    %esi,%ebp
mov    %edx,%r12d
call   1040 <pindex@plt>
mov    %r12d,%edx
mov    %ebp,%esi
mov    0x8(%rsp),%rdi
mov    0x20(%rsp),%rbp
mov    0xe0(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
mov    (%rdi),%rcx
addsd  (%rcx,%r13,8),%xmm0
movsd  %xmm0,(%rcx,%r13,8)
add    $0x1,%r14d
mov    0x0(%rbp),%eax
cmp    %eax,%r14d
jl     1410 <update_parameters+0x2b0>
mov    0x28(%rsp),%r12
mov    (%r12),%ebx
mov    0x10(%rsp),%r14
jmp    13d9 <update_parameters+0x279>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @update_parameters(i32 %0, double** %1, i32* %2, i32* %3, i32* %4, i32* %5, i32* %6, i32* %7, double* %8, double* %9, double* %10, i32 %11, i32 %12, i32 %13, i32 %14, i32 %15, i32 %16, double %17, i32* %18, double* %19, i32* %20, double* %21, i64* %22) {
%24 = alloca i32
%25 = alloca double**
%26 = alloca i32*
%27 = alloca i32*
%28 = alloca i32*
%29 = alloca i32*
%30 = alloca i32*
%31 = alloca i32*
%32 = alloca double*
%33 = alloca double*
%34 = alloca double*
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = alloca i32
%41 = alloca double
%42 = alloca i32*
%43 = alloca double*
%44 = alloca i32*
%45 = alloca double*
%46 = alloca i64*
%47 = alloca double
%48 = alloca double
%49 = alloca i32
%50 = alloca i32
%51 = alloca i32
%52 = alloca i64
%53 = alloca i64
store i32 %0, i32* %24
store double** %1, double*** %25
store i32* %2, i32** %26
store i32* %3, i32** %27
store i32* %4, i32** %28
store i32* %5, i32** %29
store i32* %6, i32** %30
store i32* %7, i32** %31
store double* %8, double** %32
store double* %9, double** %33
store double* %10, double** %34
store i32 %11, i32* %35
store i32 %12, i32* %36
store i32 %13, i32* %37
store i32 %14, i32* %38
store i32 %15, i32* %39
store i32 %16, i32* %40
store double %17, double* %41
store i32* %18, i32** %42
store double* %19, double** %43
store i32* %20, i32** %44
store double* %21, double** %45
store i64* %22, i64** %46
%54 = load i32, i32* %35
%55 = load i32, i32* %36
%56 = sub nsw i32 %54, %55
%57 = add nsw i32 %56, 1
%58 = load i32*, i32** %29
store i32 %57, i32* %58
%59 = load i32, i32* %37
%60 = load i32, i32* %38
%61 = sub nsw i32 %59, %60
%62 = add nsw i32 %61, 1
%63 = load i32*, i32** %30
store i32 %62, i32* %63
%64 = load i32, i32* %39
%65 = load i32, i32* %40
%66 = sub nsw i32 %64, %65
%67 = add nsw i32 %66, 1
%68 = load i32*, i32** %31
store i32 %67, i32* %68
%69 = load i32, i32* %36
%70 = sitofp i32 %69 to double
%71 = load double, double* %41
%72 = load double*, double** %32
%73 = load double, double* %72
%74 = call double @llvm.fmuladd.f64(double %70, double %71, double %73)
store double %74, double* %72
%75 = load i32, i32* %38
%76 = sitofp i32 %75 to double
%77 = load double, double* %41
%78 = load double*, double** %33
%79 = load double, double* %78
%80 = call double @llvm.fmuladd.f64(double %76, double %77, double %79)
store double %80, double* %78
%81 = load i32, i32* %40
%82 = sitofp i32 %81 to double
%83 = load double, double* %41
%84 = load double*, double** %34
%85 = load double, double* %84
%86 = call double @llvm.fmuladd.f64(double %82, double %83, double %85)
store double %86, double* %84
%87 = load i32, i32* %24
switch i32 %87, label %118 [
i32 1, label %88
i32 2, label %98
i32 3, label %108
]
88: ; preds = %23
%89 = load i32*, i32** %29
%90 = load i32, i32* %89
%91 = load i32*, i32** %26
store i32 %90, i32* %91
%92 = load i32*, i32** %30
%93 = load i32, i32* %92
%94 = load i32*, i32** %27
store i32 %93, i32* %94
%95 = load i32*, i32** %31
%96 = load i32, i32* %95
%97 = load i32*, i32** %28
store i32 %96, i32* %97
br label %118
98: ; preds = %23
%99 = load i32*, i32** %31
%100 = load i32, i32* %99
%101 = load i32*, i32** %26
store i32 %100, i32* %101
%102 = load i32*, i32** %29
%103 = load i32, i32* %102
%104 = load i32*, i32** %27
store i32 %103, i32* %104
%105 = load i32*, i32** %30
%106 = load i32, i32* %105
%107 = load i32*, i32** %28
store i32 %106, i32* %107
br label %118
108: ; preds = %23
%109 = load i32*, i32** %30
%110 = load i32, i32* %109
%111 = load i32*, i32** %26
store i32 %110, i32* %111
%112 = load i32*, i32** %31
%113 = load i32, i32* %112
%114 = load i32*, i32** %27
store i32 %113, i32* %114
%115 = load i32*, i32** %29
%116 = load i32, i32* %115
%117 = load i32*, i32** %28
store i32 %116, i32* %117
br label %118
118: ; preds = %23, %108, %98, %88
%119 = load i32*, i32** %28
%120 = load i32, i32* %119
%121 = sdiv i32 %120, 2
%122 = load i32*, i32** %42
store i32 %121, i32* %122
%123 = load i32*, i32** %29
%124 = load i32, i32* %123
%125 = load i32*, i32** %30
%126 = load i32, i32* %125
%127 = mul i32 %124, %126
%128 = load i32*, i32** %31
%129 = load i32, i32* %128
%130 = mul i32 %127, %129
%131 = zext i32 %130 to i64
%132 = load i64*, i64** %46
store i64 %131, i64* %132
%133 = load double*, double** %45
%134 = load i64*, i64** %46
%135 = load i64, i64* %134
%136 = call double @calc_min(double* %133, i64 %135)
store double %136, double* %48
%137 = load double*, double** %45
%138 = load i64*, i64** %46
%139 = load i64, i64* %138
%140 = call double @calc_max(double* %137, i64 %139)
store double %140, double* %47
%141 = load double, double* %47
%142 = load double, double* %48
%143 = fadd double %141, %142
%144 = fmul double 5.000000e-01, %143
%145 = load double*, double** %43
store double %144, double* %145
%146 = load i32*, i32** %26
%147 = load i32, i32* %146
%148 = sitofp i32 %147 to double
%149 = fdiv double %148, 1.000000e+02
%150 = fptosi double %149 to i32
%151 = call i32 @ceil(i32 %150)
%152 = load i32*, i32** %44
store i32 %151, i32* %152
%153 = load i32*, i32** %44
%154 = load i32, i32* %153
%155 = load i32*, i32** %27
%156 = load i32, i32* %155
%157 = sitofp i32 %156 to double
%158 = fdiv double %157, 1.000000e+02
%159 = fptosi double %158 to i32
%160 = call i32 @ceil(i32 %159)
%161 = icmp slt i32 %154, %160
br i1 %161, label %162, label %170
162: ; preds = %118
%163 = load i32*, i32** %27
%164 = load i32, i32* %163
%165 = sitofp i32 %164 to double
%166 = fdiv double %165, 1.000000e+02
%167 = fptosi double %166 to i32
%168 = call i32 @ceil(i32 %167)
%169 = load i32*, i32** %44
store i32 %168, i32* %169
br label %170
170: ; preds = %162, %118
%171 = load i32*, i32** %44
%172 = load i32, i32* %171
%173 = icmp slt i32 %172, 1
br i1 %173, label %174, label %176
174: ; preds = %170
%175 = load i32*, i32** %44
store i32 1, i32* %175
br label %176
176: ; preds = %174, %170
%177 = load double**, double*** %25
%178 = load double*, double** %177
%179 = call i32 @free(double* %178)
%180 = load double**, double*** %25
%181 = load i32*, i32** %26
%182 = load i32, i32* %181
%183 = load i32*, i32** %27
%184 = load i32, i32* %183
%185 = mul nsw i32 %182, %184
%186 = sext i32 %185 to i64
%187 = call i32 @do_vect(double** %180, i64 %186)
store i32 0, i32* %49
br label %188
188: ; preds = %251, %176
%189 = load i32, i32* %49
%190 = load i32*, i32** %26
%191 = load i32, i32* %190
%192 = icmp slt i32 %189, %191
br i1 %192, label %193, label %254
193: ; preds = %188
store i32 0, i32* %50
br label %194
194: ; preds = %247, %193
%195 = load i32, i32* %50
%196 = load i32*, i32** %27
%197 = load i32, i32* %196
%198 = icmp slt i32 %195, %197
br i1 %198, label %199, label %250
199: ; preds = %194
%200 = load i32, i32* %49
%201 = load i32, i32* %50
%202 = load i32*, i32** %26
%203 = load i32, i32* %202
%204 = mul nsw i32 %201, %203
%205 = add nsw i32 %200, %204
%206 = sext i32 %205 to i64
store i64 %206, i64* %53
store i32 0, i32* %51
br label %207
207: ; preds = %234, %199
%208 = load i32, i32* %51
%209 = load i32*, i32** %28
%210 = load i32, i32* %209
%211 = icmp slt i32 %208, %210
br i1 %211, label %212, label %237
212: ; preds = %207
%213 = load i32, i32* %24
%214 = load i32, i32* %49
%215 = load i32, i32* %50
%216 = load i32, i32* %51
%217 = load i32*, i32** %29
%218 = load i32, i32* %217
%219 = load i32*, i32** %30
%220 = load i32, i32* %219
%221 = load i32*, i32** %31
%222 = load i32, i32* %221
%223 = call i64 @pindex(i32 %213, i32 %214, i32 %215, i32 %216, i32 %218, i32 %220, i32 %222)
store i64 %223, i64* %52
%224 = load double*, double** %45
%225 = load i64, i64* %52
%226 = getelementptr inbounds double, double* %224, i64 %225
%227 = load double, double* %226
%228 = load double**, double*** %25
%229 = load double*, double** %228
%230 = load i64, i64* %53
%231 = getelementptr inbounds double, double* %229, i64 %230
%232 = load double, double* %231
%233 = fadd double %232, %227
store double %233, double* %231
br label %234
234: ; preds = %212
%235 = load i32, i32* %51
%236 = add nsw i32 %235, 1
store i32 %236, i32* %51
br label %207
237: ; preds = %207
%238 = load i32*, i32** %28
%239 = load i32, i32* %238
%240 = sitofp i32 %239 to double
%241 = load double**, double*** %25
%242 = load double*, double** %241
%243 = load i64, i64* %53
%244 = getelementptr inbounds double, double* %242, i64 %243
%245 = load double, double* %244
%246 = fdiv double %245, %240
store double %246, double* %244
br label %247
247: ; preds = %237
%248 = load i32, i32* %50
%249 = add nsw i32 %248, 1
store i32 %249, i32* %50
br label %194
250: ; preds = %194
br label %251
251: ; preds = %250
%252 = load i32, i32* %49
%253 = add nsw i32 %252, 1
store i32 %253, i32* %49
br label %188
254: ; preds = %188
ret void
}

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