#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 */
float atan2(float, float);
int /*<<< orphan*/ memcpy(float *, float const *, int);
float sin(float const);
float sqrt(float);
void interpolate_quaternion(const float floor[4], const float ceil[4],
                            const float interpolation, float output[4]) {
  float c, s1, s2, s3, omega;
  float ceil_temp[4];
  c = floor[0] * ceil[0];
  c += floor[1] * ceil[1];
  c += floor[2] * ceil[2];
  c += floor[3] * ceil[3];
  memcpy(ceil_temp, ceil, 4 * sizeof(float));
  if (c < 0.0f) {
    ceil_temp[0] *= -1.0f;
    ceil_temp[1] *= -1.0f;
    ceil_temp[2] *= -1.0f;
    ceil_temp[3] *= -1.0f;
    c *= -1.0f;
  }
  if (c > 0.9999f) {
    output[0] = (1 - interpolation) * floor[0] + interpolation * ceil_temp[0];
    output[1] = (1 - interpolation) * floor[1] + interpolation * ceil_temp[1];
    output[2] = (1 - interpolation) * floor[2] + interpolation * ceil_temp[2];
    output[3] = (1 - interpolation) * floor[3] + interpolation * ceil_temp[3];
  } else {
    s1 = sqrt(1.0f - c * c);
    omega = atan2(s1, c);
    s2 = sin((1 - interpolation) * omega);
    s3 = sin(interpolation * omega);
    output[0] = s2 * floor[0] + s3 * ceil_temp[0];
    output[0] /= s1;
    output[1] = s2 * floor[1] + s3 * ceil_temp[1];
    output[1] /= s1;
    output[2] = s2 * floor[2] + s3 * ceil_temp[2];
    output[2] /= s1;
    output[3] = s2 * floor[3] + s3 * ceil_temp[3];
    output[3] /= s1;
  }
}


x86 assembly (clang, -O0):
<interpolate_quaternion>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0xc(%rax),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
lea    -0x50(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    $0x10,%edx
call   1050 <memcpy@plt>
xorps  %xmm0,%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    1253 <interpolate_quaternion+0x113>
movss  0xdff(%rip),%xmm0
mulss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  0xded(%rip),%xmm0
mulss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  0xddb(%rip),%xmm0
mulss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  0xdc9(%rip),%xmm0
mulss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  0xdb7(%rip),%xmm0
mulss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
movss  0xda4(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    1330 <interpolate_quaternion+0x1f0>
movss  0xd97(%rip),%xmm0
subss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x14(%rbp),%xmm1
mulss  -0x50(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  0xd68(%rip),%xmm0
subss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm2
movss  -0x14(%rbp),%xmm1
mulss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  0xd37(%rip),%xmm0
subss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
movss  -0x14(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  0xd06(%rip),%xmm0
subss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm2
movss  -0x14(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0xc(%rax)
jmp    14a1 <interpolate_quaternion+0x361>
movss  -0x24(%rbp),%xmm0
movss  -0x24(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm1,%xmm0
movss  0xcb5(%rip),%xmm1
addss  %xmm1,%xmm0
call   1060 <sqrt@plt>
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
movss  -0x24(%rbp),%xmm1
call   1030 <atan2@plt>
movss  %xmm0,-0x34(%rbp)
movss  0xc8b(%rip),%xmm0
subss  -0x14(%rbp),%xmm0
mulss  -0x34(%rbp),%xmm0
call   1040 <sin@plt>
movss  %xmm0,-0x2c(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x34(%rbp),%xmm0
call   1040 <sin@plt>
movss  %xmm0,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  -0x50(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x28(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x28(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movss  0x4(%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0x4(%rax)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  -0x28(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movss  0x8(%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0x8(%rax)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0xc(%rax)
movss  -0x28(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movss  0xc(%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0xc(%rax)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interpolate_quaternion>:
push   %r14
push   %rbx
sub    $0x48,%rsp
mov    %rdx,%r14
movss  %xmm0,0xc(%rsp)
mov    %rdi,%rbx
movss  (%rdi),%xmm0
mulss  (%rsi),%xmm0
movss  0x4(%rdi),%xmm1
mulss  0x4(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  0x8(%rdi),%xmm0
mulss  0x8(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0xc(%rdi),%xmm1
mulss  0xc(%rsi),%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,0x30(%rsp)
lea    0x10(%rsp),%rdi
mov    $0x10,%edx
call   1050 <memcpy@plt>
movaps 0x30(%rsp),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    11bd <interpolate_quaternion+0x7d>
movaps 0xe53(%rip),%xmm0
movaps 0x10(%rsp),%xmm2
xorps  %xmm0,%xmm2
movaps %xmm2,0x10(%rsp)
xorps  %xmm0,%xmm1
ucomiss 0xe4c(%rip),%xmm1
jbe    1241 <interpolate_quaternion+0x101>
movss  0xe46(%rip),%xmm3
movss  0xc(%rsp),%xmm2
subss  %xmm2,%xmm3
movss  0x10(%rsp),%xmm0
mulss  %xmm2,%xmm0
movss  (%rbx),%xmm1
mulss  %xmm3,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r14)
movss  0x14(%rsp),%xmm0
mulss  %xmm2,%xmm0
movss  0x4(%rbx),%xmm1
mulss  %xmm3,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x4(%r14)
movss  0x18(%rsp),%xmm0
mulss  %xmm2,%xmm0
movss  0x8(%rbx),%xmm1
mulss  %xmm3,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%r14)
mulss  0x1c(%rsp),%xmm2
mulss  0xc(%rbx),%xmm3
addss  %xmm2,%xmm3
jmp    1333 <interpolate_quaternion+0x1f3>
movaps %xmm1,%xmm2
mulss  %xmm1,%xmm2
movss  0xdc4(%rip),%xmm0
subss  %xmm2,%xmm0
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm0
jb     1262 <interpolate_quaternion+0x122>
sqrtss %xmm0,%xmm0
jmp    1271 <interpolate_quaternion+0x131>
movaps %xmm1,0x30(%rsp)
call   1060 <sqrt@plt>
movaps 0x30(%rsp),%xmm1
movss  %xmm0,0x30(%rsp)
call   1030 <atan2@plt>
movaps %xmm0,%xmm1
movss  %xmm0,0x2c(%rsp)
movss  0xd87(%rip),%xmm0
subss  0xc(%rsp),%xmm0
mulss  %xmm1,%xmm0
call   1040 <sin@plt>
movss  %xmm0,0x28(%rsp)
movss  0x2c(%rsp),%xmm0
mulss  0xc(%rsp),%xmm0
call   1040 <sin@plt>
movss  0x28(%rsp),%xmm3
movss  0x10(%rsp),%xmm1
mulss  %xmm0,%xmm1
movss  (%rbx),%xmm2
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
movss  0x30(%rsp),%xmm4
divss  %xmm4,%xmm2
movss  %xmm2,(%r14)
movss  0x14(%rsp),%xmm1
mulss  %xmm0,%xmm1
movss  0x4(%rbx),%xmm2
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
divss  %xmm4,%xmm2
movss  %xmm2,0x4(%r14)
movss  0x18(%rsp),%xmm1
mulss  %xmm0,%xmm1
movss  0x8(%rbx),%xmm2
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
divss  %xmm4,%xmm2
movss  %xmm2,0x8(%r14)
mulss  0x1c(%rsp),%xmm0
mulss  0xc(%rbx),%xmm3
addss  %xmm0,%xmm3
divss  %xmm4,%xmm3
movss  %xmm3,0xc(%r14)
add    $0x48,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @interpolate_quaternion(float* nocapture %0, float* %1, float %2, float* nocapture writeonly %3) {
%5 = alloca <4 x float>
%6 = bitcast <4 x float>* %5 to i8*
%7 = load float, float* %0
%8 = load float, float* %1
%9 = fmul float %7, %8
%10 = getelementptr inbounds float, float* %0, i64 1
%11 = load float, float* %10
%12 = getelementptr inbounds float, float* %1, i64 1
%13 = load float, float* %12
%14 = tail call float @llvm.fmuladd.f32(float %11, float %13, float %9)
%15 = getelementptr inbounds float, float* %0, i64 2
%16 = load float, float* %15
%17 = getelementptr inbounds float, float* %1, i64 2
%18 = load float, float* %17
%19 = tail call float @llvm.fmuladd.f32(float %16, float %18, float %14)
%20 = getelementptr inbounds float, float* %0, i64 3
%21 = load float, float* %20
%22 = getelementptr inbounds float, float* %1, i64 3
%23 = load float, float* %22
%24 = tail call float @llvm.fmuladd.f32(float %21, float %23, float %19)
%25 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 0
%26 = call i32 @memcpy(float* %25, float* %1, i32 16)
%27 = fcmp olt float %24, 0.000000e+00
br i1 %27, label %28, label %32
28: ; preds = %4
%29 = load <4 x float>, <4 x float>* %5
%30 = fneg <4 x float> %29
store <4 x float> %30, <4 x float>* %5
%31 = fneg float %24
br label %32
32: ; preds = %28, %4
%33 = phi float [ %31, %28 ], [ %24, %4 ]
%34 = fcmp ogt float %33, 0x3FEFFF2E40000000
br i1 %34, label %35, label %58
35: ; preds = %32
%36 = fsub float 1.000000e+00, %2
%37 = load float, float* %0
%38 = load float, float* %25
%39 = fmul float %38, %2
%40 = call float @llvm.fmuladd.f32(float %36, float %37, float %39)
store float %40, float* %3
%41 = load float, float* %10
%42 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 1
%43 = load float, float* %42
%44 = fmul float %43, %2
%45 = call float @llvm.fmuladd.f32(float %36, float %41, float %44)
%46 = getelementptr inbounds float, float* %3, i64 1
store float %45, float* %46
%47 = load float, float* %15
%48 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 2
%49 = load float, float* %48
%50 = fmul float %49, %2
%51 = call float @llvm.fmuladd.f32(float %36, float %47, float %50)
%52 = getelementptr inbounds float, float* %3, i64 2
store float %51, float* %52
%53 = load float, float* %20
%54 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 3
%55 = load float, float* %54
%56 = fmul float %55, %2
%57 = call float @llvm.fmuladd.f32(float %36, float %53, float %56)
br label %93
58: ; preds = %32
%59 = fneg float %33
%60 = call float @llvm.fmuladd.f32(float %59, float %33, float 1.000000e+00)
%61 = call float @sqrt(float %60)
%62 = call float @atan2(float %61, float %33)
%63 = fsub float 1.000000e+00, %2
%64 = fmul float %63, %62
%65 = call float @sin(float %64)
%66 = fmul float %62, %2
%67 = call float @sin(float %66)
%68 = load float, float* %0
%69 = load float, float* %25
%70 = fmul float %67, %69
%71 = call float @llvm.fmuladd.f32(float %65, float %68, float %70)
%72 = fdiv float %71, %61
store float %72, float* %3
%73 = load float, float* %10
%74 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 1
%75 = load float, float* %74
%76 = fmul float %67, %75
%77 = call float @llvm.fmuladd.f32(float %65, float %73, float %76)
%78 = getelementptr inbounds float, float* %3, i64 1
%79 = fdiv float %77, %61
store float %79, float* %78
%80 = load float, float* %15
%81 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 2
%82 = load float, float* %81
%83 = fmul float %67, %82
%84 = call float @llvm.fmuladd.f32(float %65, float %80, float %83)
%85 = getelementptr inbounds float, float* %3, i64 2
%86 = fdiv float %84, %61
store float %86, float* %85
%87 = load float, float* %20
%88 = getelementptr inbounds <4 x float>, <4 x float>* %5, i64 0, i64 3
%89 = load float, float* %88
%90 = fmul float %67, %89
%91 = call float @llvm.fmuladd.f32(float %65, float %87, float %90)
%92 = fdiv float %91, %61
br label %93
93: ; preds = %58, %35
%94 = phi float [ %57, %35 ], [ %92, %58 ]
%95 = getelementptr inbounds float, float* %3, i64 3
store float %94, float* %95
ret void
}