#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 sqrt(float);
float sd_it(float sd, int t, int v, float m, float mnext) {
  float ret = 0;
  if (t > 0)
    ret = (1 - 1 / (double)t) * sd * sd + (t + 1) * (mnext - m) * (mnext - m);
  return sqrt(ret);
}


x86 assembly (clang, -O0):
<sd_it>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movss  %xmm0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm1,-0x10(%rbp)
movss  %xmm2,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
cmpl   $0x0,-0x8(%rbp)
jle    11ad <sd_it+0x9d>
cvtsi2sdl -0x8(%rbp),%xmm0
movsd  0xeb4(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  0xea8(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm1
movss  -0x14(%rbp),%xmm3
subss  -0x10(%rbp),%xmm3
mulss  %xmm3,%xmm1
movss  -0x14(%rbp),%xmm3
subss  -0x10(%rbp),%xmm3
mulss  %xmm3,%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
call   1030 <sqrt@plt>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sd_it>:
xorps  %xmm4,%xmm4
xorps  %xmm3,%xmm3
test   %edi,%edi
jle    1169 <sd_it+0x59>
xorps  %xmm3,%xmm3
cvtsi2sd %edi,%xmm3
movsd  0xed7(%rip),%xmm5
movapd %xmm5,%xmm6
divsd  %xmm3,%xmm6
subsd  %xmm6,%xmm5
cvtss2sd %xmm0,%xmm0
mulsd  %xmm0,%xmm5
add    $0x1,%edi
xorps  %xmm3,%xmm3
cvtsi2ss %edi,%xmm3
subss  %xmm1,%xmm2
mulss  %xmm2,%xmm3
mulss  %xmm2,%xmm3
xorps  %xmm1,%xmm1
cvtss2sd %xmm3,%xmm1
mulsd  %xmm0,%xmm5
addsd  %xmm1,%xmm5
xorps  %xmm3,%xmm3
cvtsd2ss %xmm5,%xmm3
ucomiss %xmm4,%xmm3
jb     1176 <sd_it+0x66>
xorps  %xmm0,%xmm0
sqrtss %xmm3,%xmm0
ret
movaps %xmm3,%xmm0
jmp    1030 <sqrt@plt>

LLVM IR (-O0):
define float @sd_it(float %0, i32 %1, i32 %2, float %3, float %4) {
%6 = alloca float
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
%10 = alloca float
%11 = alloca float
store float %0, float* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store float %3, float* %9
store float %4, float* %10
store float 0.000000e+00, float* %11
%12 = load i32, i32* %7
%13 = icmp sgt i32 %12, 0
br i1 %13, label %14, label %38
14: ; preds = %5
%15 = load i32, i32* %7
%16 = sitofp i32 %15 to double
%17 = fdiv double 1.000000e+00, %16
%18 = fsub double 1.000000e+00, %17
%19 = load float, float* %6
%20 = fpext float %19 to double
%21 = fmul double %18, %20
%22 = load float, float* %6
%23 = fpext float %22 to double
%24 = load i32, i32* %7
%25 = add nsw i32 %24, 1
%26 = sitofp i32 %25 to float
%27 = load float, float* %10
%28 = load float, float* %9
%29 = fsub float %27, %28
%30 = fmul float %26, %29
%31 = load float, float* %10
%32 = load float, float* %9
%33 = fsub float %31, %32
%34 = fmul float %30, %33
%35 = fpext float %34 to double
%36 = call double @llvm.fmuladd.f64(double %21, double %23, double %35)
%37 = fptrunc double %36 to float
store float %37, float* %11
br label %38
38: ; preds = %14, %5
%39 = load float, float* %11
%40 = call float @sqrt(float %39)
ret float %40
}

LLVM IR (-O2):
define float @sd_it(float %0, i32 %1, i32 %2, float %3, float %4) {
%6 = icmp sgt i32 %1, 0
br i1 %6, label %7, label %21
7: ; preds = %5
%8 = sitofp i32 %1 to double
%9 = fdiv double 1.000000e+00, %8
%10 = fsub double 1.000000e+00, %9
%11 = fpext float %0 to double
%12 = fmul double %10, %11
%13 = add nuw nsw i32 %1, 1
%14 = sitofp i32 %13 to float
%15 = fsub float %4, %3
%16 = fmul float %15, %14
%17 = fmul float %15, %16
%18 = fpext float %17 to double
%19 = tail call double @llvm.fmuladd.f64(double %12, double %11, double %18)
%20 = fptrunc double %19 to float
br label %21
21: ; preds = %7, %5
%22 = phi float [ %20, %7 ], [ 0.000000e+00, %5 ]
%23 = tail call float @sqrt(float %22)
ret float %23
}