float tfunk(float p[]) {
  int j;
  float q, r, sumd = 0.0, sumr = 0.0;
  static float wid[4 + 1] = {0.0, 1.0, 3.0, 10.0, 30.0};
  for (j = 1; j <= 4; j++) {
    q = p[j] * wid[j];
    r = (float)(q >= 0 ? (int)(q + 0.5) : -(int)(0.5 - q));
    sumr += q * q;
    sumd += (q - r) * (q - r);
  }
  return 1 + sumr * (1 + (sumd > 0.3 * 0.3 ? 2.0 : 2.0 * sumd / (0.3 * 0.3)));
}


x86 assembly (clang, -O0):
<tfunk>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x1,-0xc(%rbp)
cmpl   $0x4,-0xc(%rbp)
jg     11f9 <tfunk+0xf9>
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0xc(%rbp),%rcx
lea    0x2eef(%rip),%rax
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     117d <tfunk+0x7d>
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeab(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x20(%rbp)
jmp    119d <tfunk+0x9d>
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe8a(%rip),%xmm0
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
movss  -0x1c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x10(%rbp),%xmm0
subss  -0x14(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
subss  -0x14(%rbp),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    111f <tfunk+0x1f>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xde8(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1234 <tfunk+0x134>
movsd  0xdde(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    125a <tfunk+0x15a>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdc3(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0xdaf(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  0xda4(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd94(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<tfunk>:
movss  0x4(%rdi),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
jae    1128 <tfunk+0x28>
movsd  0xee4(%rip),%xmm2
subsd  %xmm1,%xmm2
cvttsd2si %xmm2,%eax
neg    %eax
jmp    1134 <tfunk+0x34>
addsd  0xed0(%rip),%xmm1
cvttsd2si %xmm1,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movaps %xmm0,%xmm2
subss  %xmm1,%xmm2
mulss  %xmm2,%xmm2
xorps  %xmm3,%xmm3
movss  0x8(%rdi),%xmm1
mulss  0xeca(%rip),%xmm1
addss  %xmm3,%xmm2
ucomiss %xmm3,%xmm1
xorps  %xmm3,%xmm3
cvtss2sd %xmm1,%xmm3
jae    117a <tfunk+0x7a>
movsd  0xe92(%rip),%xmm4
subsd  %xmm3,%xmm4
cvttsd2si %xmm4,%eax
neg    %eax
jmp    1186 <tfunk+0x86>
addsd  0xe7e(%rip),%xmm3
cvttsd2si %xmm3,%eax
xorps  %xmm4,%xmm4
cvtsi2ss %eax,%xmm4
movaps %xmm1,%xmm3
subss  %xmm4,%xmm3
mulss  %xmm3,%xmm3
addss  %xmm2,%xmm3
movss  0xc(%rdi),%xmm2
mulss  0xe7b(%rip),%xmm2
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm2
xorps  %xmm4,%xmm4
cvtss2sd %xmm2,%xmm4
jae    11cc <tfunk+0xcc>
movsd  0xe40(%rip),%xmm5
subsd  %xmm4,%xmm5
cvttsd2si %xmm5,%eax
neg    %eax
jmp    11d8 <tfunk+0xd8>
addsd  0xe2c(%rip),%xmm4
cvttsd2si %xmm4,%eax
xorps  %xmm5,%xmm5
cvtsi2ss %eax,%xmm5
movaps %xmm2,%xmm4
subss  %xmm5,%xmm4
mulss  %xmm4,%xmm4
addss  %xmm3,%xmm4
movss  0x10(%rdi),%xmm3
mulss  0xe2d(%rip),%xmm3
xorps  %xmm5,%xmm5
ucomiss %xmm5,%xmm3
xorps  %xmm5,%xmm5
cvtss2sd %xmm3,%xmm5
jae    121e <tfunk+0x11e>
movsd  0xdee(%rip),%xmm6
subsd  %xmm5,%xmm6
cvttsd2si %xmm6,%eax
neg    %eax
jmp    122a <tfunk+0x12a>
addsd  0xdda(%rip),%xmm5
cvttsd2si %xmm5,%eax
xorps  %xmm5,%xmm5
cvtsi2ss %eax,%xmm5
movaps %xmm3,%xmm6
subss  %xmm5,%xmm6
mulss  %xmm6,%xmm6
addss  %xmm4,%xmm6
xorps  %xmm4,%xmm4
cvtss2sd %xmm6,%xmm4
ucomisd 0xdc1(%rip),%xmm4
jbe    125b <tfunk+0x15b>
movsd  0xdaf(%rip),%xmm4
jmp    126f <tfunk+0x16f>
addsd  %xmm4,%xmm4
divsd  0xda9(%rip),%xmm4
addsd  0xda9(%rip),%xmm4
mulss  %xmm0,%xmm0
xorps  %xmm5,%xmm5
addss  %xmm0,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm5,%xmm1
mulss  %xmm2,%xmm2
addss  %xmm1,%xmm2
mulss  %xmm3,%xmm3
addss  %xmm2,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
mulsd  %xmm4,%xmm0
addsd  0xd73(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
ret

LLVM IR (-O0):
@tfunk.wid = internal global [5 x float] [float 0.000000e+00, float 1.000000e+00, float 3.000000e+00, float 1.000000e+01, float 3.000000e+01]
define float @tfunk(float* %0) {
%2 = alloca float*
%3 = alloca i32
%4 = alloca float
%5 = alloca float
%6 = alloca float
%7 = alloca float
store float* %0, float** %2
store float 0.000000e+00, float* %6
store float 0.000000e+00, float* %7
store i32 1, i32* %3
br label %8
8: ; preds = %50, %1
%9 = load i32, i32* %3
%10 = icmp sle i32 %9, 4
br i1 %10, label %11, label %53
11: ; preds = %8
%12 = load float*, float** %2
%13 = load i32, i32* %3
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds float, float* %12, i64 %14
%16 = load float, float* %15
%17 = load i32, i32* %3
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds [5 x float], [5 x float]* @tfunk.wid, i64 0, i64 %18
%20 = load float, float* %19
%21 = fmul float %16, %20
store float %21, float* %4
%22 = load float, float* %4
%23 = fcmp oge float %22, 0.000000e+00
br i1 %23, label %24, label %29
24: ; preds = %11
%25 = load float, float* %4
%26 = fpext float %25 to double
%27 = fadd double %26, 5.000000e-01
%28 = fptosi double %27 to i32
br label %35
29: ; preds = %11
%30 = load float, float* %4
%31 = fpext float %30 to double
%32 = fsub double 5.000000e-01, %31
%33 = fptosi double %32 to i32
%34 = sub nsw i32 0, %33
br label %35
35: ; preds = %29, %24
%36 = phi i32 [ %28, %24 ], [ %34, %29 ]
%37 = sitofp i32 %36 to float
store float %37, float* %5
%38 = load float, float* %4
%39 = load float, float* %4
%40 = load float, float* %7
%41 = call float @llvm.fmuladd.f32(float %38, float %39, float %40)
store float %41, float* %7
%42 = load float, float* %4
%43 = load float, float* %5
%44 = fsub float %42, %43
%45 = load float, float* %4
%46 = load float, float* %5
%47 = fsub float %45, %46
%48 = load float, float* %6
%49 = call float @llvm.fmuladd.f32(float %44, float %47, float %48)
store float %49, float* %6
br label %50
50: ; preds = %35
%51 = load i32, i32* %3
%52 = add nsw i32 %51, 1
store i32 %52, i32* %3
br label %8
53: ; preds = %8
%54 = load float, float* %7
%55 = fpext float %54 to double
%56 = load float, float* %6
%57 = fpext float %56 to double
%58 = fcmp ogt double %57, 0x3FB70A3D70A3D70A
br i1 %58, label %59, label %60
59: ; preds = %53
br label %65
60: ; preds = %53
%61 = load float, float* %6
%62 = fpext float %61 to double
%63 = fmul double 2.000000e+00, %62
%64 = fdiv double %63, 0x3FB70A3D70A3D70A
br label %65
65: ; preds = %60, %59
%66 = phi double [ 2.000000e+00, %59 ], [ %64, %60 ]
%67 = fadd double 1.000000e+00, %66
%68 = call double @llvm.fmuladd.f64(double %55, double %67, double 1.000000e+00)
%69 = fptrunc double %68 to float
ret float %69
}

LLVM IR (-O2):
define float @tfunk(float* nocapture %0) {
%2 = getelementptr inbounds float, float* %0, i64 1
%3 = load float, float* %2
%4 = fcmp ult float %3, 0.000000e+00
%5 = fpext float %3 to double
br i1 %4, label %9, label %6
6: ; preds = %1
%7 = fadd double %5, 5.000000e-01
%8 = fptosi double %7 to i32
br label %13
9: ; preds = %1
%10 = fsub double 5.000000e-01, %5
%11 = fptosi double %10 to i32
%12 = sub nsw i32 0, %11
br label %13
13: ; preds = %9, %6
%14 = phi i32 [ %8, %6 ], [ %12, %9 ]
%15 = sitofp i32 %14 to float
%16 = fsub float %3, %15
%17 = tail call float @llvm.fmuladd.f32(float %16, float %16, float 0.000000e+00)
%18 = getelementptr inbounds float, float* %0, i64 2
%19 = load float, float* %18
%20 = fmul float %19, 3.000000e+00
%21 = fcmp ult float %20, 0.000000e+00
%22 = fpext float %20 to double
br i1 %21, label %26, label %23
23: ; preds = %13
%24 = fadd double %22, 5.000000e-01
%25 = fptosi double %24 to i32
br label %30
26: ; preds = %13
%27 = fsub double 5.000000e-01, %22
%28 = fptosi double %27 to i32
%29 = sub nsw i32 0, %28
br label %30
30: ; preds = %26, %23
%31 = phi i32 [ %25, %23 ], [ %29, %26 ]
%32 = sitofp i32 %31 to float
%33 = fsub float %20, %32
%34 = tail call float @llvm.fmuladd.f32(float %33, float %33, float %17)
%35 = getelementptr inbounds float, float* %0, i64 3
%36 = load float, float* %35
%37 = fmul float %36, 1.000000e+01
%38 = fcmp ult float %37, 0.000000e+00
%39 = fpext float %37 to double
br i1 %38, label %43, label %40
40: ; preds = %30
%41 = fadd double %39, 5.000000e-01
%42 = fptosi double %41 to i32
br label %47
43: ; preds = %30
%44 = fsub double 5.000000e-01, %39
%45 = fptosi double %44 to i32
%46 = sub nsw i32 0, %45
br label %47
47: ; preds = %43, %40
%48 = phi i32 [ %42, %40 ], [ %46, %43 ]
%49 = sitofp i32 %48 to float
%50 = fsub float %37, %49
%51 = tail call float @llvm.fmuladd.f32(float %50, float %50, float %34)
%52 = getelementptr inbounds float, float* %0, i64 4
%53 = load float, float* %52
%54 = fmul float %53, 3.000000e+01
%55 = fcmp ult float %54, 0.000000e+00
%56 = fpext float %54 to double
br i1 %55, label %60, label %57
57: ; preds = %47
%58 = fadd double %56, 5.000000e-01
%59 = fptosi double %58 to i32
br label %64
60: ; preds = %47
%61 = fsub double 5.000000e-01, %56
%62 = fptosi double %61 to i32
%63 = sub nsw i32 0, %62
br label %64
64: ; preds = %60, %57
%65 = phi i32 [ %59, %57 ], [ %63, %60 ]
%66 = sitofp i32 %65 to float
%67 = fsub float %54, %66
%68 = tail call float @llvm.fmuladd.f32(float %67, float %67, float %51)
%69 = fpext float %68 to double
%70 = fcmp ogt double %69, 0x3FB70A3D70A3D70A
br i1 %70, label %75, label %71
71: ; preds = %64
%72 = fmul double %69, 2.000000e+00
%73 = fdiv double %72, 0x3FB70A3D70A3D70A
%74 = fadd double %73, 1.000000e+00
br label %75
75: ; preds = %64, %71
%76 = phi double [ %74, %71 ], [ 3.000000e+00, %64 ]
%77 = tail call float @llvm.fmuladd.f32(float %3, float %3, float 0.000000e+00)
%78 = tail call float @llvm.fmuladd.f32(float %20, float %20, float %77)
%79 = tail call float @llvm.fmuladd.f32(float %37, float %37, float %78)
%80 = tail call float @llvm.fmuladd.f32(float %54, float %54, float %79)
%81 = fpext float %80 to double
%82 = tail call double @llvm.fmuladd.f64(double %81, double %76, double 1.000000e+00)
%83 = fptrunc double %82 to float
ret float %83
}