#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 cos(float);
double sin(float);
int /*<<< orphan*/ zapit(float *, int);
int gen_ucsb2_stf(float *slip, float *t0, float *stf, int nt, float *dt,
                  float *z0) {
  int it, nstf;
  float tau, tau1, tau2, tau1x2, arg1, arg2;
  float sum, t, alpha;
  float pi = 3.141592654;
  zapit(stf, nt);
  tau = 2.0 * (*t0);
  tau1 = 0.5 * 0.13 * tau;
  tau2 = tau - tau1;
  tau1x2 = 2.0 * tau1;
  nstf = (int)((tau) / (*dt) + 0.5);
  if (nstf > nt)
    nstf = nt;
  if (nstf == 0)
    return (0);
  for (it = 0; it < nstf; it++) {
    t = it * (*dt);
    alpha = 0.0;
    if (t < tau1) {
      arg1 = pi * t / tau1;
      arg2 = 0.5 * arg1;
      alpha = 0.7 - 0.7 * cos(arg1) + 0.6 * sin(arg2);
    } else if (t < tau1x2) {
      arg1 = pi * t / tau1;
      arg2 = pi * (t - tau1) / tau2;
      alpha = 1.0 - 0.7 * cos(arg1) + 0.3 * cos(arg2);
    } else if (t < tau) {
      arg1 = pi * (t - tau1) / tau2;
      alpha = 0.3 + 0.3 * cos(arg1);
    }
    stf[it] = alpha;
  }
  sum = 0.0;
  for (it = 0; it < nstf; it++)
    sum = sum + (*dt) * stf[it];
  if (sum <= 0.0)
    return (0);
  sum = (*slip) / sum;
  for (it = 0; it < nstf; it++)
    stf[it] = stf[it] * sum;
  return (nstf);
}


x86 assembly (clang, -O0):
<gen_ucsb2_stf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movss  0xee6(%rip),%xmm0
movss  %xmm0,-0x68(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x24(%rbp),%esi
call   1040 <zapit@plt>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe89(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe73(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x44(%rbp),%xmm0
subss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x48(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe3e(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x44(%rbp),%xmm0
mov    -0x30(%rbp),%rax
divss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe10(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jle    120d <gen_ucsb2_stf+0xdd>
mov    -0x24(%rbp),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x40(%rbp)
jne    1223 <gen_ucsb2_stf+0xf3>
movl   $0x0,-0x4(%rbp)
jmp    14c8 <gen_ucsb2_stf+0x398>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    13fe <gen_ucsb2_stf+0x2ce>
cvtsi2ssl -0x3c(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x60(%rbp),%xmm1
movss  -0x48(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    12e9 <gen_ucsb2_stf+0x1b9>
movss  -0x68(%rbp),%xmm0
mulss  -0x60(%rbp),%xmm0
divss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd78(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x54(%rbp),%xmm0
call   1050 <cos@plt>
movsd  0xd79(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd85(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movss  -0x58(%rbp),%xmm0
call   1030 <sin@plt>
movsd  -0x70(%rbp),%xmm1
movsd  0xd5d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
jmp    13de <gen_ucsb2_stf+0x2ae>
movss  -0x60(%rbp),%xmm1
movss  -0x50(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1381 <gen_ucsb2_stf+0x251>
movss  -0x68(%rbp),%xmm0
mulss  -0x60(%rbp),%xmm0
divss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x68(%rbp),%xmm0
movss  -0x60(%rbp),%xmm1
subss  -0x48(%rbp),%xmm1
mulss  %xmm1,%xmm0
divss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x54(%rbp),%xmm0
call   1050 <cos@plt>
movsd  0xce1(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xcdd(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movss  -0x58(%rbp),%xmm0
call   1050 <cos@plt>
movsd  -0x78(%rbp),%xmm1
movsd  0xcad(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
jmp    13d9 <gen_ucsb2_stf+0x2a9>
movss  -0x60(%rbp),%xmm1
movss  -0x44(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    13d4 <gen_ucsb2_stf+0x2a4>
movss  -0x68(%rbp),%xmm0
movss  -0x60(%rbp),%xmm1
subss  -0x48(%rbp),%xmm1
mulss  %xmm1,%xmm0
divss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
call   1050 <cos@plt>
movsd  0xc55(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
jmp    13d9 <gen_ucsb2_stf+0x2a9>
jmp    13de <gen_ucsb2_stf+0x2ae>
movss  -0x64(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    122a <gen_ucsb2_stf+0xfa>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    144e <gen_ucsb2_stf+0x31e>
movss  -0x5c(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    140d <gen_ucsb2_stf+0x2dd>
movss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     1470 <gen_ucsb2_stf+0x340>
movl   $0x0,-0x4(%rbp)
jmp    14c8 <gen_ucsb2_stf+0x398>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
divss  -0x5c(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    14c2 <gen_ucsb2_stf+0x392>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x5c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1489 <gen_ucsb2_stf+0x359>
mov    -0x40(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gen_ucsb2_stf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,%r15
mov    %ecx,%ebp
mov    %rdx,%rbx
mov    %rsi,%r12
mov    %rdi,%r14
mov    %rdx,%rdi
mov    %ecx,%esi
call   1040 <zapit@plt>
movss  (%r12),%xmm1
addss  %xmm1,%xmm1
movss  (%r15),%xmm0
movss  %xmm1,0x4(%rsp)
divss  %xmm0,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  0xe8d(%rip),%xmm1
cvttsd2si %xmm1,%r12d
cmp    %ebp,%r12d
cmovg  %ebp,%r12d
xor    %eax,%eax
test   %r12d,%r12d
je     14aa <gen_ucsb2_stf+0x37a>
jle    14aa <gen_ucsb2_stf+0x37a>
mov    %r14,0x20(%rsp)
movss  0x4(%rsp),%xmm2
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
mulsd  0xe4e(%rip),%xmm1
cvtsd2ss %xmm1,%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x8(%rsp)
movaps %xmm1,%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,0x1c(%rsp)
mov    %r12d,%r13d
lea    -0x1(%r13),%rbp
xor    %r14d,%r14d
movss  %xmm1,0xc(%rsp)
nopl   (%rax)
xorps  %xmm2,%xmm2
cvtsi2ss %r14d,%xmm2
mulss  %xmm0,%xmm2
ucomiss %xmm2,%xmm1
jbe    1250 <gen_ucsb2_stf+0x120>
mulss  0xe3f(%rip),%xmm2
divss  %xmm1,%xmm2
movaps %xmm2,%xmm0
mulss  0xe34(%rip),%xmm0
movss  %xmm0,(%rsp)
movaps %xmm2,%xmm0
call   1050 <cos@plt>
mulsd  0xdfb(%rip),%xmm0
addsd  0xe03(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movss  (%rsp),%xmm0
call   1030 <sin@plt>
mulsd  0xdf3(%rip),%xmm0
addsd  0x10(%rsp),%xmm0
jmp    1300 <gen_ucsb2_stf+0x1d0>
nopl   0x0(%rax,%rax,1)
movss  0x1c(%rsp),%xmm0
ucomiss %xmm2,%xmm0
jbe    12c0 <gen_ucsb2_stf+0x190>
movaps %xmm2,%xmm0
movaps %xmm2,%xmm3
movss  0xdcf(%rip),%xmm2
mulss  %xmm2,%xmm0
divss  %xmm1,%xmm0
subss  %xmm1,%xmm3
mulss  %xmm2,%xmm3
divss  0x8(%rsp),%xmm3
movss  %xmm3,(%rsp)
call   1050 <cos@plt>
mulsd  0xd87(%rip),%xmm0
addsd  0xd87(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movss  (%rsp),%xmm0
call   1050 <cos@plt>
mulsd  0xd5f(%rip),%xmm0
addsd  0x10(%rsp),%xmm0
jmp    1300 <gen_ucsb2_stf+0x1d0>
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
movss  0x4(%rsp),%xmm1
ucomiss %xmm2,%xmm1
jbe    1304 <gen_ucsb2_stf+0x1d4>
subss  0xc(%rsp),%xmm2
mulss  0xd5c(%rip),%xmm2
divss  0x8(%rsp),%xmm2
movaps %xmm2,%xmm0
call   1050 <cos@plt>
movsd  0xd1e(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
nopw   0x0(%rax,%rax,1)
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rbx,%r14,4)
cmp    %r14,%rbp
je     1323 <gen_ucsb2_stf+0x1f3>
movss  (%r15),%xmm0
add    $0x1,%r14
movss  0xc(%rsp),%xmm1
jmp    11e0 <gen_ucsb2_stf+0xb0>
test   %r12d,%r12d
jle    1345 <gen_ucsb2_stf+0x215>
movss  (%r15),%xmm0
mov    %r13d,%eax
and    $0x3,%eax
cmp    $0x3,%rbp
jae    134c <gen_ucsb2_stf+0x21c>
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
mov    0x20(%rsp),%rsi
jmp    13a0 <gen_ucsb2_stf+0x270>
xor    %eax,%eax
jmp    14aa <gen_ucsb2_stf+0x37a>
mov    %r13d,%edx
and    $0xfffffffc,%edx
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
mov    0x20(%rsp),%rsi
nopl   0x0(%rax)
movss  (%rbx,%rcx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  0x4(%rbx,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
addss  %xmm2,%xmm1
movss  0x8(%rbx,%rcx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  0xc(%rbx,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
addss  %xmm2,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1360 <gen_ucsb2_stf+0x230>
test   %rax,%rax
je     13c6 <gen_ucsb2_stf+0x296>
lea    (%rbx,%rcx,4),%rcx
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movss  (%rcx,%rdx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm2,%xmm1
add    $0x1,%rdx
cmp    %rdx,%rax
jne    13b0 <gen_ucsb2_stf+0x280>
xor    %eax,%eax
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jae    14aa <gen_ucsb2_stf+0x37a>
test   %r12d,%r12d
jle    14a7 <gen_ucsb2_stf+0x377>
movss  (%rsi),%xmm0
divss  %xmm1,%xmm0
cmp    $0x8,%r12d
jae    13f2 <gen_ucsb2_stf+0x2c2>
xor    %eax,%eax
jmp    1490 <gen_ucsb2_stf+0x360>
mov    %r13d,%eax
and    $0xfffffff8,%eax
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%rax),%rcx
mov    %rcx,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     145e <gen_ucsb2_stf+0x32e>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopl   0x0(%rax)
movups (%rbx,%rcx,4),%xmm2
movups 0x10(%rbx,%rcx,4),%xmm3
movups 0x20(%rbx,%rcx,4),%xmm4
movups 0x30(%rbx,%rcx,4),%xmm5
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,(%rbx,%rcx,4)
movups %xmm3,0x10(%rbx,%rcx,4)
mulps  %xmm1,%xmm4
mulps  %xmm1,%xmm5
movups %xmm4,0x20(%rbx,%rcx,4)
movups %xmm5,0x30(%rbx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1420 <gen_ucsb2_stf+0x2f0>
jmp    1460 <gen_ucsb2_stf+0x330>
xor    %ecx,%ecx
test   $0x1,%dl
je     147d <gen_ucsb2_stf+0x34d>
movups (%rbx,%rcx,4),%xmm2
movups 0x10(%rbx,%rcx,4),%xmm3
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,(%rbx,%rcx,4)
movups %xmm3,0x10(%rbx,%rcx,4)
cmp    %r13,%rax
je     14a7 <gen_ucsb2_stf+0x377>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rbx,%rax,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rbx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r13
jne    1490 <gen_ucsb2_stf+0x360>
mov    %r12d,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @gen_ucsb2_stf(float* %0, float* %1, float* %2, i32 %3, float* %4, float* %5) {
%7 = alloca i32
%8 = alloca float*
%9 = alloca float*
%10 = alloca float*
%11 = alloca i32
%12 = alloca float*
%13 = alloca float*
%14 = alloca i32
%15 = alloca i32
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
%22 = alloca float
%23 = alloca float
%24 = alloca float
%25 = alloca float
store float* %0, float** %8
store float* %1, float** %9
store float* %2, float** %10
store i32 %3, i32* %11
store float* %4, float** %12
store float* %5, float** %13
store float 0x400921FB60000000, float* %25
%26 = load float*, float** %10
%27 = load i32, i32* %11
%28 = call i32 @zapit(float* %26, i32 %27)
%29 = load float*, float** %9
%30 = load float, float* %29
%31 = fpext float %30 to double
%32 = fmul double 2.000000e+00, %31
%33 = fptrunc double %32 to float
store float %33, float* %16
%34 = load float, float* %16
%35 = fpext float %34 to double
%36 = fmul double 6.500000e-02, %35
%37 = fptrunc double %36 to float
store float %37, float* %17
%38 = load float, float* %16
%39 = load float, float* %17
%40 = fsub float %38, %39
store float %40, float* %18
%41 = load float, float* %17
%42 = fpext float %41 to double
%43 = fmul double 2.000000e+00, %42
%44 = fptrunc double %43 to float
store float %44, float* %19
%45 = load float, float* %16
%46 = load float*, float** %12
%47 = load float, float* %46
%48 = fdiv float %45, %47
%49 = fpext float %48 to double
%50 = fadd double %49, 5.000000e-01
%51 = fptosi double %50 to i32
store i32 %51, i32* %15
%52 = load i32, i32* %15
%53 = load i32, i32* %11
%54 = icmp sgt i32 %52, %53
br i1 %54, label %55, label %57
55: ; preds = %6
%56 = load i32, i32* %11
store i32 %56, i32* %15
br label %57
57: ; preds = %55, %6
%58 = load i32, i32* %15
%59 = icmp eq i32 %58, 0
br i1 %59, label %60, label %61
60: ; preds = %57
store i32 0, i32* %7
br label %192
61: ; preds = %57
store i32 0, i32* %14
br label %62
62: ; preds = %140, %61
%63 = load i32, i32* %14
%64 = load i32, i32* %15
%65 = icmp slt i32 %63, %64
br i1 %65, label %66, label %143
66: ; preds = %62
%67 = load i32, i32* %14
%68 = sitofp i32 %67 to float
%69 = load float*, float** %12
%70 = load float, float* %69
%71 = fmul float %68, %70
store float %71, float* %23
store float 0.000000e+00, float* %24
%72 = load float, float* %23
%73 = load float, float* %17
%74 = fcmp olt float %72, %73
br i1 %74, label %75, label %92
75: ; preds = %66
%76 = load float, float* %25
%77 = load float, float* %23
%78 = fmul float %76, %77
%79 = load float, float* %17
%80 = fdiv float %78, %79
store float %80, float* %20
%81 = load float, float* %20
%82 = fpext float %81 to double
%83 = fmul double 5.000000e-01, %82
%84 = fptrunc double %83 to float
store float %84, float* %21
%85 = load float, float* %20
%86 = call double @cos(float %85)
%87 = call double @llvm.fmuladd.f64(double 0xBFE6666666666666, double %86, double 0x3FE6666666666666)
%88 = load float, float* %21
%89 = call double @sin(float %88)
%90 = call double @llvm.fmuladd.f64(double 6.000000e-01, double %89, double %87)
%91 = fptrunc double %90 to float
store float %91, float* %24
br label %134
92: ; preds = %66
%93 = load float, float* %23
%94 = load float, float* %19
%95 = fcmp olt float %93, %94
br i1 %95, label %96, label %116
96: ; preds = %92
%97 = load float, float* %25
%98 = load float, float* %23
%99 = fmul float %97, %98
%100 = load float, float* %17
%101 = fdiv float %99, %100
store float %101, float* %20
%102 = load float, float* %25
%103 = load float, float* %23
%104 = load float, float* %17
%105 = fsub float %103, %104
%106 = fmul float %102, %105
%107 = load float, float* %18
%108 = fdiv float %106, %107
store float %108, float* %21
%109 = load float, float* %20
%110 = call double @cos(float %109)
%111 = call double @llvm.fmuladd.f64(double 0xBFE6666666666666, double %110, double 1.000000e+00)
%112 = load float, float* %21
%113 = call double @cos(float %112)
%114 = call double @llvm.fmuladd.f64(double 3.000000e-01, double %113, double %111)
%115 = fptrunc double %114 to float
store float %115, float* %24
br label %133
116: ; preds = %92
%117 = load float, float* %23
%118 = load float, float* %16
%119 = fcmp olt float %117, %118
br i1 %119, label %120, label %132
120: ; preds = %116
%121 = load float, float* %25
%122 = load float, float* %23
%123 = load float, float* %17
%124 = fsub float %122, %123
%125 = fmul float %121, %124
%126 = load float, float* %18
%127 = fdiv float %125, %126
store float %127, float* %20
%128 = load float, float* %20
%129 = call double @cos(float %128)
%130 = call double @llvm.fmuladd.f64(double 3.000000e-01, double %129, double 3.000000e-01)
%131 = fptrunc double %130 to float
store float %131, float* %24
br label %132
132: ; preds = %120, %116
br label %133
133: ; preds = %132, %96
br label %134
134: ; preds = %133, %75
%135 = load float, float* %24
%136 = load float*, float** %10
%137 = load i32, i32* %14
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds float, float* %136, i64 %138
store float %135, float* %139
br label %140
140: ; preds = %134
%141 = load i32, i32* %14
%142 = add nsw i32 %141, 1
store i32 %142, i32* %14
br label %62
143: ; preds = %62
store float 0.000000e+00, float* %22
store i32 0, i32* %14
br label %144
144: ; preds = %158, %143
%145 = load i32, i32* %14
%146 = load i32, i32* %15
%147 = icmp slt i32 %145, %146
br i1 %147, label %148, label %161
148: ; preds = %144
%149 = load float, float* %22
%150 = load float*, float** %12
%151 = load float, float* %150
%152 = load float*, float** %10
%153 = load i32, i32* %14
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds float, float* %152, i64 %154
%156 = load float, float* %155
%157 = call float @llvm.fmuladd.f32(float %151, float %156, float %149)
store float %157, float* %22
br label %158
158: ; preds = %148
%159 = load i32, i32* %14
%160 = add nsw i32 %159, 1
store i32 %160, i32* %14
br label %144
161: ; preds = %144
%162 = load float, float* %22
%163 = fpext float %162 to double
%164 = fcmp ole double %163, 0.000000e+00
br i1 %164, label %165, label %166
165: ; preds = %161
store i32 0, i32* %7
br label %192
166: ; preds = %161
%167 = load float*, float** %8
%168 = load float, float* %167
%169 = load float, float* %22
%170 = fdiv float %168, %169
store float %170, float* %22
store i32 0, i32* %14
br label %171
171: ; preds = %187, %166
%172 = load i32, i32* %14
%173 = load i32, i32* %15
%174 = icmp slt i32 %172, %173
br i1 %174, label %175, label %190
175: ; preds = %171
%176 = load float*, float** %10
%177 = load i32, i32* %14
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds float, float* %176, i64 %178
%180 = load float, float* %179
%181 = load float, float* %22
%182 = fmul float %180, %181
%183 = load float*, float** %10
%184 = load i32, i32* %14
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds float, float* %183, i64 %185
store float %182, float* %186
br label %187
187: ; preds = %175
%188 = load i32, i32* %14
%189 = add nsw i32 %188, 1
store i32 %189, i32* %14
br label %171
190: ; preds = %171
%191 = load i32, i32* %15
store i32 %191, i32* %7
br label %192
192: ; preds = %190, %165, %60
%193 = load i32, i32* %7
ret i32 %193
}

LLVM IR (-O2):
define i32 @gen_ucsb2_stf(float* nocapture %0, float* nocapture %1, float* %2, i32 %3, float* nocapture %4, float* nocapture readnone %5) {
%7 = tail call i32 @zapit(float* %2, i32 %3)
%8 = load float, float* %1
%9 = fmul float %8, 2.000000e+00
%10 = fpext float %9 to double
%11 = fmul double %10, 6.500000e-02
%12 = fptrunc double %11 to float
%13 = fsub float %9, %12
%14 = fmul float %12, 2.000000e+00
%15 = load float, float* %4
%16 = fdiv float %9, %15
%17 = fpext float %16 to double
%18 = fadd double %17, 5.000000e-01
%19 = fptosi double %18 to i32
%20 = icmp sgt i32 %19, %3
%21 = select i1 %20, i32 %3, i32 %19
%22 = icmp eq i32 %21, 0
br i1 %22, label %191, label %23
23: ; preds = %6
%24 = icmp sgt i32 %21, 0
br i1 %24, label %25, label %191
25: ; preds = %23
%26 = zext i32 %21 to i64
br label %35
27: ; preds = %73
br i1 %24, label %28, label %191
28: ; preds = %27
%29 = load float, float* %4
%30 = add nsw i64 %26, -1
%31 = and i64 %26, 3
%32 = icmp ult i64 %30, 3
br i1 %32, label %102, label %33
33: ; preds = %28
%34 = and i64 %26, 4294967292
br label %80
35: ; preds = %78, %25
%36 = phi float [ %15, %25 ], [ %79, %78 ]
%37 = phi i64 [ 0, %25 ], [ %76, %78 ]
%38 = trunc i64 %37 to i32
%39 = sitofp i32 %38 to float
%40 = fmul float %36, %39
%41 = fcmp olt float %40, %12
br i1 %41, label %42, label %51
42: ; preds = %35
%43 = fmul float %40, 0x400921FB60000000
%44 = fdiv float %43, %12
%45 = fmul float %44, 5.000000e-01
%46 = tail call double @cos(float %44)
%47 = tail call double @llvm.fmuladd.f64(double %46, double 0xBFE6666666666666, double 0x3FE6666666666666)
%48 = tail call double @sin(float %45)
%49 = tail call double @llvm.fmuladd.f64(double %48, double 6.000000e-01, double %47)
%50 = fptrunc double %49 to float
br label %73
51: ; preds = %35
%52 = fcmp olt float %40, %14
br i1 %52, label %53, label %64
53: ; preds = %51
%54 = fmul float %40, 0x400921FB60000000
%55 = fdiv float %54, %12
%56 = fsub float %40, %12
%57 = fmul float %56, 0x400921FB60000000
%58 = fdiv float %57, %13
%59 = tail call double @cos(float %55)
%60 = tail call double @llvm.fmuladd.f64(double %59, double 0xBFE6666666666666, double 1.000000e+00)
%61 = tail call double @cos(float %58)
%62 = tail call double @llvm.fmuladd.f64(double %61, double 3.000000e-01, double %60)
%63 = fptrunc double %62 to float
br label %73
64: ; preds = %51
%65 = fcmp olt float %40, %9
br i1 %65, label %66, label %73
66: ; preds = %64
%67 = fsub float %40, %12
%68 = fmul float %67, 0x400921FB60000000
%69 = fdiv float %68, %13
%70 = tail call double @cos(float %69)
%71 = tail call double @llvm.fmuladd.f64(double %70, double 3.000000e-01, double 3.000000e-01)
%72 = fptrunc double %71 to float
br label %73
73: ; preds = %53, %66, %64, %42
%74 = phi float [ %50, %42 ], [ %63, %53 ], [ %72, %66 ], [ 0.000000e+00, %64 ]
%75 = getelementptr inbounds float, float* %2, i64 %37
store float %74, float* %75
%76 = add nuw nsw i64 %37, 1
%77 = icmp eq i64 %76, %26
br i1 %77, label %27, label %78
78: ; preds = %73
%79 = load float, float* %4
br label %35
80: ; preds = %80, %33
%81 = phi i64 [ 0, %33 ], [ %99, %80 ]
%82 = phi float [ 0.000000e+00, %33 ], [ %98, %80 ]
%83 = phi i64 [ 0, %33 ], [ %100, %80 ]
%84 = getelementptr inbounds float, float* %2, i64 %81
%85 = load float, float* %84
%86 = tail call float @llvm.fmuladd.f32(float %29, float %85, float %82)
%87 = or i64 %81, 1
%88 = getelementptr inbounds float, float* %2, i64 %87
%89 = load float, float* %88
%90 = tail call float @llvm.fmuladd.f32(float %29, float %89, float %86)
%91 = or i64 %81, 2
%92 = getelementptr inbounds float, float* %2, i64 %91
%93 = load float, float* %92
%94 = tail call float @llvm.fmuladd.f32(float %29, float %93, float %90)
%95 = or i64 %81, 3
%96 = getelementptr inbounds float, float* %2, i64 %95
%97 = load float, float* %96
%98 = tail call float @llvm.fmuladd.f32(float %29, float %97, float %94)
%99 = add nuw nsw i64 %81, 4
%100 = add i64 %83, 4
%101 = icmp eq i64 %100, %34
br i1 %101, label %102, label %80
102: ; preds = %80, %28
%103 = phi float [ undef, %28 ], [ %98, %80 ]
%104 = phi i64 [ 0, %28 ], [ %99, %80 ]
%105 = phi float [ 0.000000e+00, %28 ], [ %98, %80 ]
%106 = icmp eq i64 %31, 0
br i1 %106, label %117, label %107
107: ; preds = %102, %107
%108 = phi i64 [ %114, %107 ], [ %104, %102 ]
%109 = phi float [ %113, %107 ], [ %105, %102 ]
%110 = phi i64 [ %115, %107 ], [ 0, %102 ]
%111 = getelementptr inbounds float, float* %2, i64 %108
%112 = load float, float* %111
%113 = tail call float @llvm.fmuladd.f32(float %29, float %112, float %109)
%114 = add nuw nsw i64 %108, 1
%115 = add i64 %110, 1
%116 = icmp eq i64 %115, %31
br i1 %116, label %117, label %107
117: ; preds = %107, %102
%118 = phi float [ %103, %102 ], [ %113, %107 ]
%119 = fcmp ugt float %118, 0.000000e+00
br i1 %119, label %120, label %191
120: ; preds = %117
%121 = load float, float* %0
%122 = fdiv float %121, %118
br i1 %24, label %123, label %191
123: ; preds = %120
%124 = zext i32 %21 to i64
%125 = icmp ult i32 %21, 8
br i1 %125, label %182, label %126
126: ; preds = %123
%127 = and i64 %26, 4294967288
%128 = insertelement <4 x float> poison, float %122, i64 0
%129 = shufflevector <4 x float> %128, <4 x float> poison, <4 x i32> zeroinitializer
%130 = insertelement <4 x float> poison, float %122, i64 0
%131 = shufflevector <4 x float> %130, <4 x float> poison, <4 x i32> zeroinitializer
%132 = add nsw i64 %127, -8
%133 = lshr exact i64 %132, 3
%134 = add nuw nsw i64 %133, 1
%135 = and i64 %134, 1
%136 = icmp eq i64 %132, 0
br i1 %136, label %166, label %137
137: ; preds = %126
%138 = and i64 %134, 4611686018427387902
br label %139
139: ; preds = %139, %137
%140 = phi i64 [ 0, %137 ], [ %163, %139 ]
%141 = phi i64 [ 0, %137 ], [ %164, %139 ]
%142 = getelementptr inbounds float, float* %2, i64 %140
%143 = bitcast float* %142 to <4 x float>*
%144 = load <4 x float>, <4 x float>* %143
%145 = getelementptr inbounds float, float* %142, i64 4
%146 = bitcast float* %145 to <4 x float>*
%147 = load <4 x float>, <4 x float>* %146
%148 = fmul <4 x float> %129, %144
%149 = fmul <4 x float> %131, %147
%150 = bitcast float* %142 to <4 x float>*
store <4 x float> %148, <4 x float>* %150
%151 = bitcast float* %145 to <4 x float>*
store <4 x float> %149, <4 x float>* %151
%152 = or i64 %140, 8
%153 = getelementptr inbounds float, float* %2, i64 %152
%154 = bitcast float* %153 to <4 x float>*
%155 = load <4 x float>, <4 x float>* %154
%156 = getelementptr inbounds float, float* %153, i64 4
%157 = bitcast float* %156 to <4 x float>*
%158 = load <4 x float>, <4 x float>* %157
%159 = fmul <4 x float> %129, %155
%160 = fmul <4 x float> %131, %158
%161 = bitcast float* %153 to <4 x float>*
store <4 x float> %159, <4 x float>* %161
%162 = bitcast float* %156 to <4 x float>*
store <4 x float> %160, <4 x float>* %162
%163 = add nuw i64 %140, 16
%164 = add i64 %141, 2
%165 = icmp eq i64 %164, %138
br i1 %165, label %166, label %139
166: ; preds = %139, %126
%167 = phi i64 [ 0, %126 ], [ %163, %139 ]
%168 = icmp eq i64 %135, 0
br i1 %168, label %180, label %169
169: ; preds = %166
%170 = getelementptr inbounds float, float* %2, i64 %167
%171 = bitcast float* %170 to <4 x float>*
%172 = load <4 x float>, <4 x float>* %171
%173 = getelementptr inbounds float, float* %170, i64 4
%174 = bitcast float* %173 to <4 x float>*
%175 = load <4 x float>, <4 x float>* %174
%176 = fmul <4 x float> %129, %172
%177 = fmul <4 x float> %131, %175
%178 = bitcast float* %170 to <4 x float>*
store <4 x float> %176, <4 x float>* %178
%179 = bitcast float* %173 to <4 x float>*
store <4 x float> %177, <4 x float>* %179
br label %180
180: ; preds = %166, %169
%181 = icmp eq i64 %127, %26
br i1 %181, label %191, label %182
182: ; preds = %123, %180
%183 = phi i64 [ 0, %123 ], [ %127, %180 ]
br label %184
184: ; preds = %182, %184
%185 = phi i64 [ %189, %184 ], [ %183, %182 ]
%186 = getelementptr inbounds float, float* %2, i64 %185
%187 = load float, float* %186
%188 = fmul float %122, %187
store float %188, float* %186
%189 = add nuw nsw i64 %185, 1
%190 = icmp eq i64 %189, %124
br i1 %190, label %191, label %184
191: ; preds = %184, %180, %23, %27, %120, %117, %6
%192 = phi i32 [ 0, %6 ], [ 0, %117 ], [ %21, %120 ], [ 0, %27 ], [ 0, %23 ], [ %21, %180 ], [ %21, %184 ]
ret i32 %192
}