void HTS_set_duration(int *duration, double *mean, double *vari, double *remain,
                      int size, double frame_length) {
  int i;
  double temp1, temp2;
  double rho = 0.0;
  if (frame_length != 0.0) {
    temp1 = 0.0;
    temp2 = 0.0;
    for (i = 0; i < size; i++) {
      temp1 += mean[i];
      temp2 += vari[i];
    }
    rho = (frame_length - temp1) / temp2;
  }
  for (i = 0; i < size; i++) {
    temp1 = mean[i] + rho * vari[i] + *remain;
    duration[i] = (int)(temp1 + 0.5);
    if (duration[i] < 1)
      duration[i] = 1;
    *remain = temp1 - (double)duration[i];
  }
}


x86 assembly (clang, -O0):
<HTS_set_duration>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1142 <HTS_set_duration+0x42>
jp     1142 <HTS_set_duration+0x42>
jmp    11b5 <HTS_set_duration+0xb5>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11a1 <HTS_set_duration+0xa1>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1159 <HTS_set_duration+0x59>
movsd  -0x30(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1265 <HTS_set_duration+0x165>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xdfc(%rip),%xmm0
addsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edx
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jge    1239 <HTS_set_duration+0x139>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x1,(%rax,%rcx,4)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
subsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11bc <HTS_set_duration+0xbc>
pop    %rbp
ret

x86 assembly (clang, -O2):
<HTS_set_duration>:
push   %rbx
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
mov    %r8d,%r10d
jne    1114 <HTS_set_duration+0x14>
jnp    11d0 <HTS_set_duration+0xd0>
test   %r8d,%r8d
jle    1136 <HTS_set_duration+0x36>
lea    -0x1(%r10),%rax
mov    %r10d,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jae    1143 <HTS_set_duration+0x43>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
xorpd  %xmm2,%xmm2
jmp    1197 <HTS_set_duration+0x97>
xorpd  %xmm2,%xmm2
xorpd  %xmm1,%xmm1
jmp    11c4 <HTS_set_duration+0xc4>
mov    %r10d,%r11d
and    $0xfffffffc,%r11d
xorpd  %xmm1,%xmm1
xor    %eax,%eax
xorpd  %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addsd  (%rsi,%rax,8),%xmm2
addsd  (%rdx,%rax,8),%xmm1
addsd  0x8(%rsi,%rax,8),%xmm2
addsd  0x8(%rdx,%rax,8),%xmm1
addsd  0x10(%rsi,%rax,8),%xmm2
addsd  0x10(%rdx,%rax,8),%xmm1
addsd  0x18(%rsi,%rax,8),%xmm2
addsd  0x18(%rdx,%rax,8),%xmm1
add    $0x4,%rax
cmp    %rax,%r11
jne    1160 <HTS_set_duration+0x60>
test   %r9,%r9
je     11c4 <HTS_set_duration+0xc4>
lea    (%rdx,%rax,8),%r11
lea    (%rsi,%rax,8),%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
addsd  (%rax,%rbx,8),%xmm2
addsd  (%r11,%rbx,8),%xmm1
add    $0x1,%rbx
cmp    %rbx,%r9
jne    11b0 <HTS_set_duration+0xb0>
subsd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
test   %r8d,%r8d
jle    1233 <HTS_set_duration+0x133>
movsd  (%rcx),%xmm2
xor    %eax,%eax
movsd  0xe1d(%rip),%xmm0
mov    $0x1,%r8d
nopl   0x0(%rax)
movsd  (%rdx,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  (%rsi,%rax,8),%xmm3
addsd  %xmm2,%xmm3
movapd %xmm3,%xmm2
addsd  %xmm0,%xmm2
cvttsd2si %xmm2,%ebx
test   %ebx,%ebx
cmovle %r8d,%ebx
mov    %ebx,(%rdi,%rax,4)
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
subsd  %xmm2,%xmm3
movsd  %xmm3,(%rcx)
add    $0x1,%rax
movapd %xmm3,%xmm2
cmp    %rax,%r10
jne    11f0 <HTS_set_duration+0xf0>
pop    %rbx
ret

LLVM IR (-O0):
define void @HTS_set_duration(i32* %0, double* %1, double* %2, double* %3, i32 %4, double %5) {
%7 = alloca i32*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca double
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
store i32* %0, i32** %7
store double* %1, double** %8
store double* %2, double** %9
store double* %3, double** %10
store i32 %4, i32* %11
store double %5, double* %12
store double 0.000000e+00, double* %16
%17 = load double, double* %12
%18 = fcmp une double %17, 0.000000e+00
br i1 %18, label %19, label %48
19: ; preds = %6
store double 0.000000e+00, double* %14
store double 0.000000e+00, double* %15
store i32 0, i32* %13
br label %20
20: ; preds = %39, %19
%21 = load i32, i32* %13
%22 = load i32, i32* %11
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %42
24: ; preds = %20
%25 = load double*, double** %8
%26 = load i32, i32* %13
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double, double* %14
%31 = fadd double %30, %29
store double %31, double* %14
%32 = load double*, double** %9
%33 = load i32, i32* %13
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double, double* %15
%38 = fadd double %37, %36
store double %38, double* %15
br label %39
39: ; preds = %24
%40 = load i32, i32* %13
%41 = add nsw i32 %40, 1
store i32 %41, i32* %13
br label %20
42: ; preds = %20
%43 = load double, double* %12
%44 = load double, double* %14
%45 = fsub double %43, %44
%46 = load double, double* %15
%47 = fdiv double %45, %46
store double %47, double* %16
br label %48
48: ; preds = %42, %6
store i32 0, i32* %13
br label %49
49: ; preds = %97, %48
%50 = load i32, i32* %13
%51 = load i32, i32* %11
%52 = icmp slt i32 %50, %51
br i1 %52, label %53, label %100
53: ; preds = %49
%54 = load double*, double** %8
%55 = load i32, i32* %13
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
%58 = load double, double* %57
%59 = load double, double* %16
%60 = load double*, double** %9
%61 = load i32, i32* %13
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
%65 = call double @llvm.fmuladd.f64(double %59, double %64, double %58)
%66 = load double*, double** %10
%67 = load double, double* %66
%68 = fadd double %65, %67
store double %68, double* %14
%69 = load double, double* %14
%70 = fadd double %69, 5.000000e-01
%71 = fptosi double %70 to i32
%72 = load i32*, i32** %7
%73 = load i32, i32* %13
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds i32, i32* %72, i64 %74
store i32 %71, i32* %75
%76 = load i32*, i32** %7
%77 = load i32, i32* %13
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds i32, i32* %76, i64 %78
%80 = load i32, i32* %79
%81 = icmp slt i32 %80, 1
br i1 %81, label %82, label %87
82: ; preds = %53
%83 = load i32*, i32** %7
%84 = load i32, i32* %13
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds i32, i32* %83, i64 %85
store i32 1, i32* %86
br label %87
87: ; preds = %82, %53
%88 = load double, double* %14
%89 = load i32*, i32** %7
%90 = load i32, i32* %13
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds i32, i32* %89, i64 %91
%93 = load i32, i32* %92
%94 = sitofp i32 %93 to double
%95 = fsub double %88, %94
%96 = load double*, double** %10
store double %95, double* %96
br label %97
97: ; preds = %87
%98 = load i32, i32* %13
%99 = add nsw i32 %98, 1
store i32 %99, i32* %13
br label %49
100: ; preds = %49
ret void
}

LLVM IR (-O2):
define void @HTS_set_duration(i32* nocapture writeonly %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, i32 %4, double %5) {
%7 = fcmp une double %5, 0.000000e+00
br i1 %7, label %8, label %78
8: ; preds = %6
%9 = icmp sgt i32 %4, 0
br i1 %9, label %10, label %73
10: ; preds = %8
%11 = zext i32 %4 to i64
%12 = add nsw i64 %11, -1
%13 = and i64 %11, 3
%14 = icmp ult i64 %12, 3
br i1 %14, label %52, label %15
15: ; preds = %10
%16 = and i64 %11, 4294967292
br label %17
17: ; preds = %17, %15
%18 = phi i64 [ 0, %15 ], [ %49, %17 ]
%19 = phi double [ 0.000000e+00, %15 ], [ %48, %17 ]
%20 = phi double [ 0.000000e+00, %15 ], [ %45, %17 ]
%21 = phi i64 [ 0, %15 ], [ %50, %17 ]
%22 = getelementptr inbounds double, double* %1, i64 %18
%23 = load double, double* %22
%24 = fadd double %20, %23
%25 = getelementptr inbounds double, double* %2, i64 %18
%26 = load double, double* %25
%27 = fadd double %19, %26
%28 = or i64 %18, 1
%29 = getelementptr inbounds double, double* %1, i64 %28
%30 = load double, double* %29
%31 = fadd double %24, %30
%32 = getelementptr inbounds double, double* %2, i64 %28
%33 = load double, double* %32
%34 = fadd double %27, %33
%35 = or i64 %18, 2
%36 = getelementptr inbounds double, double* %1, i64 %35
%37 = load double, double* %36
%38 = fadd double %31, %37
%39 = getelementptr inbounds double, double* %2, i64 %35
%40 = load double, double* %39
%41 = fadd double %34, %40
%42 = or i64 %18, 3
%43 = getelementptr inbounds double, double* %1, i64 %42
%44 = load double, double* %43
%45 = fadd double %38, %44
%46 = getelementptr inbounds double, double* %2, i64 %42
%47 = load double, double* %46
%48 = fadd double %41, %47
%49 = add nuw nsw i64 %18, 4
%50 = add i64 %21, 4
%51 = icmp eq i64 %50, %16
br i1 %51, label %52, label %17
52: ; preds = %17, %10
%53 = phi double [ undef, %10 ], [ %45, %17 ]
%54 = phi double [ undef, %10 ], [ %48, %17 ]
%55 = phi i64 [ 0, %10 ], [ %49, %17 ]
%56 = phi double [ 0.000000e+00, %10 ], [ %48, %17 ]
%57 = phi double [ 0.000000e+00, %10 ], [ %45, %17 ]
%58 = icmp eq i64 %13, 0
br i1 %58, label %73, label %59
59: ; preds = %52, %59
%60 = phi i64 [ %70, %59 ], [ %55, %52 ]
%61 = phi double [ %69, %59 ], [ %56, %52 ]
%62 = phi double [ %66, %59 ], [ %57, %52 ]
%63 = phi i64 [ %71, %59 ], [ 0, %52 ]
%64 = getelementptr inbounds double, double* %1, i64 %60
%65 = load double, double* %64
%66 = fadd double %62, %65
%67 = getelementptr inbounds double, double* %2, i64 %60
%68 = load double, double* %67
%69 = fadd double %61, %68
%70 = add nuw nsw i64 %60, 1
%71 = add i64 %63, 1
%72 = icmp eq i64 %71, %13
br i1 %72, label %73, label %59
73: ; preds = %52, %59, %8
%74 = phi double [ 0.000000e+00, %8 ], [ %53, %52 ], [ %66, %59 ]
%75 = phi double [ 0.000000e+00, %8 ], [ %54, %52 ], [ %69, %59 ]
%76 = fsub double %5, %74
%77 = fdiv double %76, %75
br label %78
78: ; preds = %73, %6
%79 = phi double [ %77, %73 ], [ 0.000000e+00, %6 ]
%80 = icmp sgt i32 %4, 0
br i1 %80, label %81, label %102
81: ; preds = %78
%82 = zext i32 %4 to i64
%83 = load double, double* %3
br label %84
84: ; preds = %81, %84
%85 = phi double [ %83, %81 ], [ %99, %84 ]
%86 = phi i64 [ 0, %81 ], [ %100, %84 ]
%87 = getelementptr inbounds double, double* %1, i64 %86
%88 = load double, double* %87
%89 = getelementptr inbounds double, double* %2, i64 %86
%90 = load double, double* %89
%91 = tail call double @llvm.fmuladd.f64(double %79, double %90, double %88)
%92 = fadd double %91, %85
%93 = fadd double %92, 5.000000e-01
%94 = fptosi double %93 to i32
%95 = getelementptr inbounds i32, i32* %0, i64 %86
%96 = icmp sgt i32 %94, 1
%97 = select i1 %96, i32 %94, i32 1
store i32 %97, i32* %95
%98 = sitofp i32 %97 to double
%99 = fsub double %92, %98
store double %99, double* %3
%100 = add nuw nsw i64 %86, 1
%101 = icmp eq i64 %100, %82
br i1 %101, label %102, label %84
102: ; preds = %84, %78
ret void
}