#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 ceil(float);
int /*<<< orphan*/ cumulativeSum(float *, float *, int);
int /*<<< orphan*/ diffArray(float *, float *, int);
int fabs(float);
int /*<<< orphan*/ grd_trapezoid(float, float, float, float, float **, float **,
                                 float **, float **, int *, int *);
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *, ...);
int /*<<< orphan*/ scalefloats(float *, int, int);
float sqrt(float);
int sumfloats(float *, int);
void grd_triangle(float GMAX, float SMAX, float Ktgt, float Ts, float **kx,
                  float **gx, float **slew, float **time, int *N, int *N_ramp) {
  float T_ramp;
  int Ngx;
  float scaling;
  int n;
  if (fabs(Ktgt / 4257) > (GMAX * GMAX) / SMAX) {
    if (1)
      printf("grd_triangle: using a trapezoid lobe\n");
    grd_trapezoid(GMAX, SMAX, Ktgt, Ts, kx, gx, slew, time, N, N_ramp);
    return;
  }
  T_ramp = sqrt(fabs(Ktgt / 4257) / SMAX);
  *N_ramp = ceil(T_ramp / Ts) + 3;
  Ngx = 2 * *N_ramp;
  *gx = (float *)malloc(Ngx * sizeof(float));
  for (n = 0; n < *N_ramp; n++) {
    (*gx)[n] = n * Ts * SMAX;
    (*gx)[Ngx - 1 - n] = (*gx)[n];
  }
  scaling = Ktgt / (sumfloats((*gx), Ngx) * 4257 * Ts);
  if (fabs(scaling) > 1 && 1)
    printf("grd_triangle: abs(scaling)=%.2f> 1 !\n", scaling);
  scalefloats(*gx, Ngx, scaling);
  if (kx) {
    *kx = (float *)malloc(Ngx * sizeof(float));
    cumulativeSum(*gx, *kx, Ngx);
    scalefloats(*kx, Ngx, 4257 * Ts);
  }
  if (slew) {
    *slew = (float *)malloc(Ngx * sizeof(float));
    diffArray(*gx, *slew, Ngx);
    scalefloats(*slew, Ngx, 1 / Ts);
  }
  if (time) {
    *time = (float *)malloc(Ngx * sizeof(float));
    for (n = 0; n < Ngx; n++)
      (*time)[n] = n * Ts;
  }
  *N = Ngx;
  return;
}


x86 assembly (clang, -O0):
<grd_triangle>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
movss  %xmm2,-0xc(%rbp)
movss  %xmm3,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
movss  -0xc(%rbp),%xmm0
movss  0xe1f(%rip),%xmm1
divss  %xmm1,%xmm0
call   1060 <fabs@plt>
cvtsi2ss %eax,%xmm0
movss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
divss  -0x8(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    124a <grd_triangle+0xaa>
lea    0xdfb(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
movss  -0x4(%rbp),%xmm0
movss  -0x8(%rbp),%xmm1
movss  -0xc(%rbp),%xmm2
movss  -0x10(%rbp),%xmm3
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%r8
mov    -0x40(%rbp),%r9
call   10a0 <grd_trapezoid@plt>
jmp    14b2 <grd_triangle+0x312>
movss  -0xc(%rbp),%xmm0
movss  0xda9(%rip),%xmm1
divss  %xmm1,%xmm0
call   1060 <fabs@plt>
cvtsi2ss %eax,%xmm0
divss  -0x8(%rbp),%xmm0
call   10b0 <sqrt@plt>
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
divss  -0x10(%rbp),%xmm0
call   1040 <ceil@plt>
mov    %eax,%ecx
add    $0x3,%ecx
mov    -0x40(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x40(%rbp),%rax
mov    (%rax),%eax
shl    $0x1,%eax
mov    %eax,-0x48(%rbp)
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    -0x40(%rbp),%rcx
cmp    (%rcx),%eax
jge    131d <grd_triangle+0x17d>
cvtsi2ssl -0x50(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
sub    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    12b9 <grd_triangle+0x119>
movss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%rdi
mov    -0x48(%rbp),%esi
call   1070 <sumfloats@plt>
movss  -0x54(%rbp),%xmm0
imul   $0x10a1,%eax,%eax
cvtsi2ss %eax,%xmm1
mulss  -0x10(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x4c(%rbp),%xmm0
call   1060 <fabs@plt>
cmp    $0x1,%eax
jle    137d <grd_triangle+0x1dd>
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcb8(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x20(%rbp),%rax
mov    (%rax),%rdi
mov    -0x48(%rbp),%esi
cvttss2si -0x4c(%rbp),%edx
call   1030 <scalefloats@plt>
cmpq   $0x0,-0x18(%rbp)
je     13eb <grd_triangle+0x24b>
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,%rcx
mov    -0x18(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x20(%rbp),%rax
mov    (%rax),%rdi
mov    -0x18(%rbp),%rax
mov    (%rax),%rsi
mov    -0x48(%rbp),%edx
call   10c0 <cumulativeSum@plt>
mov    -0x18(%rbp),%rax
mov    (%rax),%rdi
mov    -0x48(%rbp),%esi
movss  0xc23(%rip),%xmm0
mulss  -0x10(%rbp),%xmm0
cvttss2si %xmm0,%edx
call   1030 <scalefloats@plt>
cmpq   $0x0,-0x28(%rbp)
je     1445 <grd_triangle+0x2a5>
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x20(%rbp),%rax
mov    (%rax),%rdi
mov    -0x28(%rbp),%rax
mov    (%rax),%rsi
mov    -0x48(%rbp),%edx
call   1080 <diffArray@plt>
mov    -0x28(%rbp),%rax
mov    (%rax),%rdi
mov    -0x48(%rbp),%esi
movss  0xbcd(%rip),%xmm0
divss  -0x10(%rbp),%xmm0
cvttss2si %xmm0,%edx
call   1030 <scalefloats@plt>
cmpq   $0x0,-0x30(%rbp)
je     14a9 <grd_triangle+0x309>
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,%rcx
mov    -0x30(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    14a4 <grd_triangle+0x304>
cvtsi2ssl -0x50(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    (%rax),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1470 <grd_triangle+0x2d0>
jmp    14a9 <grd_triangle+0x309>
mov    -0x48(%rbp),%ecx
mov    -0x38(%rbp),%rax
mov    %ecx,(%rax)
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<grd_triangle>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,%r15
mov    %r8,%r13
mov    %rcx,%rbp
mov    %rdx,%rbx
mov    %rsi,%r14
mov    %rdi,%r12
movdqa %xmm3,0x10(%rsp)
movss  %xmm1,0x8(%rsp)
movss  %xmm0,0xc(%rsp)
movss  %xmm2,0x4(%rsp)
movaps %xmm2,%xmm0
divss  0xe0d(%rip),%xmm0
movss  %xmm0,0x2c(%rsp)
call   1070 <fabs@plt>
movss  0x8(%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  0xc(%rsp),%xmm1
mulss  %xmm1,%xmm1
divss  %xmm2,%xmm1
ucomiss %xmm1,%xmm0
jbe    1266 <grd_triangle+0xb6>
lea    0xe61(%rip),%rdi
call   1030 <puts@plt>
movss  0xc(%rsp),%xmm0
movss  0x8(%rsp),%xmm1
movss  0x4(%rsp),%xmm2
movaps 0x10(%rsp),%xmm3
mov    %r12,%rdi
mov    %r14,%rsi
mov    %rbx,%rdx
mov    %rbp,%rcx
mov    %r13,%r8
mov    %r15,%r9
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    10b0 <grd_trapezoid@plt>
movss  0x2c(%rsp),%xmm0
call   1070 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  0x8(%rsp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
mov    %rbp,0x40(%rsp)
mov    %r13,0x38(%rsp)
mov    %rbx,0x30(%rsp)
jb     129b <grd_triangle+0xeb>
sqrtss %xmm0,%xmm0
jmp    12a0 <grd_triangle+0xf0>
call   10c0 <sqrt@plt>
movaps 0x10(%rsp),%xmm1
divss  %xmm1,%xmm0
call   1050 <ceil@plt>
mov    %eax,%r13d
lea    0x3(%r13),%eax
mov    %eax,(%r15)
lea    0x6(%r13,%r13,1),%ebp
lea    0x18(,%r13,8),%edi
mov    %edi,%ebx
call   10a0 <malloc@plt>
mov    %rax,(%r14)
mov    (%r15),%ecx
mov    %rbp,%r15
test   %ecx,%ecx
movss  0x8(%rsp),%xmm2
jle    1373 <grd_triangle+0x1c3>
cmp    $0x1,%ecx
jne    12ec <grd_triangle+0x13c>
xor    %edx,%edx
jmp    134b <grd_triangle+0x19b>
mov    %ecx,%esi
and    $0xfffffffe,%esi
lea    -0x2(%r15),%edi
xor    %edx,%edx
movaps 0x10(%rsp),%xmm1
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm2,%xmm0
movss  %xmm0,(%rax,%rdx,4)
lea    0x1(%rdi),%ebp
movslq %ebp,%rbp
movss  %xmm0,(%rax,%rbp,4)
lea    0x1(%rdx),%ebp
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm2,%xmm0
movss  %xmm0,0x4(%rax,%rdx,4)
movslq %edi,%rdi
movss  %xmm0,(%rax,%rdi,4)
add    $0x2,%rdx
add    $0xfffffffe,%edi
cmp    %rsi,%rdx
jne    1300 <grd_triangle+0x150>
test   $0x1,%cl
je     1373 <grd_triangle+0x1c3>
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
mulss  0x10(%rsp),%xmm0
mulss  %xmm2,%xmm0
movss  %xmm0,(%rax,%rdx,4)
not    %edx
add    %r15d,%edx
movslq %edx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    %rax,%rdi
mov    %r15d,%esi
call   1080 <sumfloats@plt>
imul   $0x10a1,%eax,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
mulss  0x10(%rsp),%xmm1
movss  0x4(%rsp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0x4(%rsp)
call   1070 <fabs@plt>
cmp    $0x2,%eax
jl     13c3 <grd_triangle+0x213>
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xca4(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    (%r14),%rdi
cvttss2si 0x4(%rsp),%edx
mov    %r15d,%esi
call   1040 <scalefloats@plt>
test   %r12,%r12
mov    %ebx,%ebp
je     1412 <grd_triangle+0x262>
mov    %ebp,%edi
call   10a0 <malloc@plt>
mov    %rax,(%r12)
mov    (%r14),%rdi
mov    %rax,%rsi
mov    %r15d,%edx
call   10d0 <cumulativeSum@plt>
mov    (%r12),%rdi
movss  0xc00(%rip),%xmm0
mulss  0x10(%rsp),%xmm0
cvttss2si %xmm0,%edx
mov    %r15d,%esi
call   1040 <scalefloats@plt>
mov    0x30(%rsp),%rbx
test   %rbx,%rbx
je     1451 <grd_triangle+0x2a1>
mov    %ebp,%edi
call   10a0 <malloc@plt>
mov    %rax,(%rbx)
mov    (%r14),%rdi
mov    %rax,%rsi
mov    %r15d,%edx
call   1090 <diffArray@plt>
mov    (%rbx),%rdi
movss  0xbc5(%rip),%xmm0
divss  0x10(%rsp),%xmm0
cvttss2si %xmm0,%edx
mov    %r15d,%esi
call   1040 <scalefloats@plt>
mov    0x40(%rsp),%rbx
test   %rbx,%rbx
je     1589 <grd_triangle+0x3d9>
mov    %ebp,%edi
call   10a0 <malloc@plt>
movaps 0x10(%rsp),%xmm2
mov    %rax,(%rbx)
cmp    $0xfffffffe,%r13d
jl     1589 <grd_triangle+0x3d9>
cmp    $0x2,%r15d
mov    $0x1,%ecx
cmovge %r15d,%ecx
cmp    $0x8,%ecx
jae    1491 <grd_triangle+0x2e1>
xor    %edx,%edx
jmp    1570 <grd_triangle+0x3c0>
mov    %ecx,%edx
and    $0x7ffffff8,%edx
movaps %xmm2,%xmm0
shufps $0x0,%xmm2,%xmm0
lea    -0x8(%rdx),%rsi
mov    %rsi,%rbp
shr    $0x3,%rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     15a0 <grd_triangle+0x3f0>
mov    %rbp,%rdi
and    $0xfffffffffffffffe,%rdi
movdqa 0xb49(%rip),%xmm1
xor    %esi,%esi
movdqa 0xb4e(%rip),%xmm8
movdqa 0xb56(%rip),%xmm3
movdqa 0xb5e(%rip),%xmm4
movdqa 0xb66(%rip),%xmm5
nopw   0x0(%rax,%rax,1)
movdqa %xmm1,%xmm6
paddd  %xmm8,%xmm6
cvtdq2ps %xmm1,%xmm7
cvtdq2ps %xmm6,%xmm6
mulps  %xmm0,%xmm7
mulps  %xmm0,%xmm6
movups %xmm7,(%rax,%rsi,4)
movups %xmm6,0x10(%rax,%rsi,4)
movdqa %xmm1,%xmm6
paddd  %xmm3,%xmm6
movdqa %xmm1,%xmm7
paddd  %xmm4,%xmm7
cvtdq2ps %xmm6,%xmm6
cvtdq2ps %xmm7,%xmm7
mulps  %xmm0,%xmm6
mulps  %xmm0,%xmm7
movups %xmm6,0x20(%rax,%rsi,4)
movups %xmm7,0x30(%rax,%rsi,4)
add    $0x10,%rsi
paddd  %xmm5,%xmm1
add    $0xfffffffffffffffe,%rdi
jne    14f0 <grd_triangle+0x340>
test   $0x1,%bpl
je     1565 <grd_triangle+0x3b5>
cvtdq2ps %xmm1,%xmm3
paddd  0xacd(%rip),%xmm1
cvtdq2ps %xmm1,%xmm1
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm1
movups %xmm3,(%rax,%rsi,4)
movups %xmm1,0x10(%rax,%rsi,4)
cmp    %rcx,%rdx
je     1589 <grd_triangle+0x3d9>
nopw   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
mulss  %xmm2,%xmm0
movss  %xmm0,(%rax,%rdx,4)
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1570 <grd_triangle+0x3c0>
mov    0x38(%rsp),%rax
mov    %r15d,(%rax)
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movdqa 0xa68(%rip),%xmm1
xor    %esi,%esi
test   $0x1,%bpl
jne    1548 <grd_triangle+0x398>
jmp    1565 <grd_triangle+0x3b5>

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

LLVM IR (-O2):
@.str.1 = [38 x i8] c"grd_triangle: abs(scaling)=%.2f> 1 !\0A\00"
@str = [37 x i8] c"grd_triangle: using a trapezoid lobe\00"
define void @grd_triangle(float %0, float %1, float %2, float %3, float** %4, float** %5, float** %6, float** %7, i32* %8, i32* %9) {
%11 = fdiv float %2, 4.257000e+03
%12 = tail call i32 @fabs(float %11)
%13 = sitofp i32 %12 to float
%14 = fmul float %0, %0
%15 = fdiv float %14, %1
%16 = fcmp olt float %15, %13
br i1 %16, label %17, label %20
17: ; preds = %10
%18 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([37 x i8], [37 x i8]* @str, i64 0, i64 0))
%19 = tail call i32 @grd_trapezoid(float %0, float %1, float %2, float %3, float** %4, float** %5, float** %6, float** %7, i32* %8, i32* %9)
br label %195
20: ; preds = %10
%21 = tail call i32 @fabs(float %11)
%22 = sitofp i32 %21 to float
%23 = fdiv float %22, %1
%24 = tail call float @sqrt(float %23)
%25 = fdiv float %24, %3
%26 = tail call i32 @ceil(float %25)
%27 = add nsw i32 %26, 3
store i32 %27, i32* %9
%28 = shl i32 %27, 1
%29 = shl i32 %27, 3
%30 = tail call i64 @malloc(i32 %29)
%31 = inttoptr i64 %30 to float*
store float* %31, float** %5
%32 = load i32, i32* %9
%33 = icmp sgt i32 %32, 0
br i1 %33, label %34, label %81
34: ; preds = %20
%35 = zext i32 %32 to i64
%36 = and i64 %35, 1
%37 = icmp eq i32 %32, 1
br i1 %37, label %67, label %38
38: ; preds = %34
%39 = and i64 %35, 4294967294
br label %40
40: ; preds = %40, %38
%41 = phi i64 [ 0, %38 ], [ %64, %40 ]
%42 = phi i64 [ 0, %38 ], [ %65, %40 ]
%43 = trunc i64 %41 to i32
%44 = sitofp i32 %43 to float
%45 = fmul float %44, %3
%46 = fmul float %45, %1
%47 = getelementptr inbounds float, float* %31, i64 %41
store float %46, float* %47
%48 = trunc i64 %41 to i32
%49 = xor i32 %48, -1
%50 = add i32 %28, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %31, i64 %51
store float %46, float* %52
%53 = or i64 %41, 1
%54 = trunc i64 %53 to i32
%55 = sitofp i32 %54 to float
%56 = fmul float %55, %3
%57 = fmul float %56, %1
%58 = getelementptr inbounds float, float* %31, i64 %53
store float %57, float* %58
%59 = trunc i64 %53 to i32
%60 = xor i32 %59, -1
%61 = add i32 %28, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %31, i64 %62
store float %57, float* %63
%64 = add nuw nsw i64 %41, 2
%65 = add i64 %42, 2
%66 = icmp eq i64 %65, %39
br i1 %66, label %67, label %40
67: ; preds = %40, %34
%68 = phi i64 [ 0, %34 ], [ %64, %40 ]
%69 = icmp eq i64 %36, 0
br i1 %69, label %81, label %70
70: ; preds = %67
%71 = trunc i64 %68 to i32
%72 = sitofp i32 %71 to float
%73 = fmul float %72, %3
%74 = fmul float %73, %1
%75 = getelementptr inbounds float, float* %31, i64 %68
store float %74, float* %75
%76 = trunc i64 %68 to i32
%77 = xor i32 %76, -1
%78 = add i32 %28, %77
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %31, i64 %79
store float %74, float* %80
br label %81
81: ; preds = %70, %67, %20
%82 = tail call i32 @sumfloats(float* %31, i32 %28)
%83 = mul nsw i32 %82, 4257
%84 = sitofp i32 %83 to float
%85 = fmul float %84, %3
%86 = fdiv float %2, %85
%87 = tail call i32 @fabs(float %86)
%88 = icmp sgt i32 %87, 1
br i1 %88, label %89, label %92
89: ; preds = %81
%90 = fpext float %86 to double
%91 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([38 x i8], [38 x i8]* @.str.1, i64 0, i64 0), double %90)
br label %92
92: ; preds = %89, %81
%93 = load float*, float** %5
%94 = fptosi float %86 to i32
%95 = tail call i32 @scalefloats(float* %93, i32 %28, i32 %94)
%96 = icmp eq float** %4, null
br i1 %96, label %106, label %97
97: ; preds = %92
%98 = tail call i64 @malloc(i32 %29)
%99 = inttoptr i64 %98 to float*
store float* %99, float** %4
%100 = load float*, float** %5
%101 = tail call i32 @cumulativeSum(float* %100, float* %99, i32 %28)
%102 = load float*, float** %4
%103 = fmul float %3, 4.257000e+03
%104 = fptosi float %103 to i32
%105 = tail call i32 @scalefloats(float* %102, i32 %28, i32 %104)
br label %106
106: ; preds = %97, %92
%107 = icmp eq float** %6, null
br i1 %107, label %117, label %108
108: ; preds = %106
%109 = tail call i64 @malloc(i32 %29)
%110 = inttoptr i64 %109 to float*
store float* %110, float** %6
%111 = load float*, float** %5
%112 = tail call i32 @diffArray(float* %111, float* %110, i32 %28)
%113 = load float*, float** %6
%114 = fdiv float 1.000000e+00, %3
%115 = fptosi float %114 to i32
%116 = tail call i32 @scalefloats(float* %113, i32 %28, i32 %115)
br label %117
117: ; preds = %108, %106
%118 = icmp eq float** %7, null
br i1 %118, label %194, label %119
119: ; preds = %117
%120 = tail call i64 @malloc(i32 %29)
%121 = inttoptr i64 %120 to float*
store float* %121, float** %7
%122 = icmp sgt i32 %26, -3
br i1 %122, label %123, label %194
123: ; preds = %119
%124 = call i32 @llvm.smax.i32(i32 %28, i32 1)
%125 = zext i32 %124 to i64
%126 = icmp ult i32 %124, 8
br i1 %126, label %184, label %127
127: ; preds = %123
%128 = and i64 %125, 2147483640
%129 = insertelement <4 x float> poison, float %3, i64 0
%130 = shufflevector <4 x float> %129, <4 x float> poison, <4 x i32> zeroinitializer
%131 = insertelement <4 x float> poison, float %3, i64 0
%132 = shufflevector <4 x float> %131, <4 x float> poison, <4 x i32> zeroinitializer
%133 = add nsw i64 %128, -8
%134 = lshr exact i64 %133, 3
%135 = add nuw nsw i64 %134, 1
%136 = and i64 %135, 1
%137 = icmp eq i64 %133, 0
br i1 %137, label %168, label %138
138: ; preds = %127
%139 = and i64 %135, 4611686018427387902
br label %140
140: ; preds = %140, %138
%141 = phi i64 [ 0, %138 ], [ %164, %140 ]
%142 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %138 ], [ %165, %140 ]
%143 = phi i64 [ 0, %138 ], [ %166, %140 ]
%144 = add <4 x i32> %142, <i32 4, i32 4, i32 4, i32 4>
%145 = sitofp <4 x i32> %142 to <4 x float>
%146 = sitofp <4 x i32> %144 to <4 x float>
%147 = fmul <4 x float> %130, %145
%148 = fmul <4 x float> %132, %146
%149 = getelementptr inbounds float, float* %121, i64 %141
%150 = bitcast float* %149 to <4 x float>*
store <4 x float> %147, <4 x float>* %150
%151 = getelementptr inbounds float, float* %149, i64 4
%152 = bitcast float* %151 to <4 x float>*
store <4 x float> %148, <4 x float>* %152
%153 = or i64 %141, 8
%154 = add <4 x i32> %142, <i32 8, i32 8, i32 8, i32 8>
%155 = add <4 x i32> %142, <i32 12, i32 12, i32 12, i32 12>
%156 = sitofp <4 x i32> %154 to <4 x float>
%157 = sitofp <4 x i32> %155 to <4 x float>
%158 = fmul <4 x float> %130, %156
%159 = fmul <4 x float> %132, %157
%160 = getelementptr inbounds float, float* %121, i64 %153
%161 = bitcast float* %160 to <4 x float>*
store <4 x float> %158, <4 x float>* %161
%162 = getelementptr inbounds float, float* %160, i64 4
%163 = bitcast float* %162 to <4 x float>*
store <4 x float> %159, <4 x float>* %163
%164 = add nuw i64 %141, 16
%165 = add <4 x i32> %142, <i32 16, i32 16, i32 16, i32 16>
%166 = add i64 %143, 2
%167 = icmp eq i64 %166, %139
br i1 %167, label %168, label %140
168: ; preds = %140, %127
%169 = phi i64 [ 0, %127 ], [ %164, %140 ]
%170 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %127 ], [ %165, %140 ]
%171 = icmp eq i64 %136, 0
br i1 %171, label %182, label %172
172: ; preds = %168
%173 = add <4 x i32> %170, <i32 4, i32 4, i32 4, i32 4>
%174 = sitofp <4 x i32> %170 to <4 x float>
%175 = sitofp <4 x i32> %173 to <4 x float>
%176 = fmul <4 x float> %130, %174
%177 = fmul <4 x float> %132, %175
%178 = getelementptr inbounds float, float* %121, i64 %169
%179 = bitcast float* %178 to <4 x float>*
store <4 x float> %176, <4 x float>* %179
%180 = getelementptr inbounds float, float* %178, i64 4
%181 = bitcast float* %180 to <4 x float>*
store <4 x float> %177, <4 x float>* %181
br label %182
182: ; preds = %168, %172
%183 = icmp eq i64 %128, %125
br i1 %183, label %194, label %184
184: ; preds = %123, %182
%185 = phi i64 [ 0, %123 ], [ %128, %182 ]
br label %186
186: ; preds = %184, %186
%187 = phi i64 [ %192, %186 ], [ %185, %184 ]
%188 = trunc i64 %187 to i32
%189 = sitofp i32 %188 to float
%190 = fmul float %189, %3
%191 = getelementptr inbounds float, float* %121, i64 %187
store float %190, float* %191
%192 = add nuw nsw i64 %187, 1
%193 = icmp eq i64 %192, %125
br i1 %193, label %194, label %186
194: ; preds = %186, %182, %119, %117
store i32 %28, i32* %8
br label %195
195: ; preds = %194, %17
ret void
}