double ch_normalize(double *vec, int beg, int end) {
  int i;
  double scale;
  double ch_norm(double *vec, int beg, int end);
  scale = ch_norm(vec, beg, end);
  vec = vec + beg;
  for (i = end - beg + 1; i; i--) {
    *vec = *vec / scale;
    vec++;
  }
  return (scale);
}


x86 assembly (clang, -O0):
<ch_normalize>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x10(%rbp),%edx
call   1030 <ch_norm@plt>
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
je     118e <ch_normalize+0x7e>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
divsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x8(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    1155 <ch_normalize+0x45>
movsd  -0x20(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ch_normalize>:
push   %rbp
push   %r14
push   %rbx
mov    %edx,%ebx
mov    %esi,%ebp
mov    %rdi,%r14
call   1030 <ch_norm@plt>
sub    %ebp,%ebx
mov    %ebx,%ecx
inc    %ecx
je     11e6 <ch_normalize+0xd6>
movslq %ebp,%rdx
lea    (%r14,%rdx,8),%rax
test   %ebx,%ebx
je     11c8 <ch_normalize+0xb8>
add    $0x1,%rbx
mov    %rbx,%r8
and    $0xfffffffffffffffe,%r8
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%r8),%rsi
mov    %rsi,%rbp
shr    %rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     11eb <ch_normalize+0xdb>
mov    %rbp,%rsi
and    $0xfffffffffffffffe,%rsi
lea    (%r14,%rdx,8),%rdx
add    $0x10,%rdx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movupd -0x10(%rdx,%rdi,8),%xmm2
movupd (%rdx,%rdi,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdx,%rdi,8)
divpd  %xmm1,%xmm3
movupd %xmm3,(%rdx,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1180 <ch_normalize+0x70>
test   $0x1,%bpl
je     11bc <ch_normalize+0xac>
movupd (%rax,%rdi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rax,%rdi,8)
cmp    %r8,%rbx
je     11e6 <ch_normalize+0xd6>
lea    (%rax,%r8,8),%rax
sub    %r8d,%ecx
mov    %ecx,%ecx
xor    %edx,%edx
nopl   0x0(%rax)
movsd  (%rax,%rdx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rdx,8)
add    $0x1,%rdx
cmp    %edx,%ecx
jne    11d0 <ch_normalize+0xc0>
pop    %rbx
pop    %r14
pop    %rbp
ret
xor    %edi,%edi
test   $0x1,%bpl
jne    11ae <ch_normalize+0x9e>
jmp    11bc <ch_normalize+0xac>

LLVM IR (-O0):
define double @ch_normalize(double* %0, i32 %1, i32 %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
store double* %0, double** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
%9 = load double*, double** %4
%10 = load i32, i32* %5
%11 = load i32, i32* %6
%12 = call double @ch_norm(double* %9, i32 %10, i32 %11)
store double %12, double* %8
%13 = load double*, double** %4
%14 = load i32, i32* %5
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %13, i64 %15
store double* %16, double** %4
%17 = load i32, i32* %6
%18 = load i32, i32* %5
%19 = sub nsw i32 %17, %18
%20 = add nsw i32 %19, 1
store i32 %20, i32* %7
br label %21
21: ; preds = %32, %3
%22 = load i32, i32* %7
%23 = icmp ne i32 %22, 0
br i1 %23, label %24, label %35
24: ; preds = %21
%25 = load double*, double** %4
%26 = load double, double* %25
%27 = load double, double* %8
%28 = fdiv double %26, %27
%29 = load double*, double** %4
store double %28, double* %29
%30 = load double*, double** %4
%31 = getelementptr inbounds double, double* %30, i32 1
store double* %31, double** %4
br label %32
32: ; preds = %24
%33 = load i32, i32* %7
%34 = add nsw i32 %33, -1
store i32 %34, i32* %7
br label %21
35: ; preds = %21
%36 = load double, double* %8
ret double %36
}

LLVM IR (-O2):
define double @ch_normalize(double* %0, i32 %1, i32 %2) {
%4 = tail call double @ch_norm(double* %0, i32 %1, i32 %2)
%5 = sub i32 1, %1
%6 = add i32 %5, %2
%7 = icmp eq i32 %6, 0
br i1 %7, label %68, label %8
8: ; preds = %3
%9 = sext i32 %1 to i64
%10 = getelementptr inbounds double, double* %0, i64 %9
%11 = sub i32 %2, %1
%12 = zext i32 %11 to i64
%13 = add nuw nsw i64 %12, 1
%14 = icmp eq i32 %11, 0
br i1 %14, label %57, label %15
15: ; preds = %8
%16 = and i64 %13, 8589934590
%17 = getelementptr double, double* %10, i64 %16
%18 = trunc i64 %16 to i32
%19 = sub i32 %6, %18
%20 = insertelement <2 x double> poison, double %4, i64 0
%21 = shufflevector <2 x double> %20, <2 x double> poison, <2 x i32> zeroinitializer
%22 = add nsw i64 %16, -2
%23 = lshr exact i64 %22, 1
%24 = add nuw i64 %23, 1
%25 = and i64 %24, 1
%26 = icmp eq i64 %22, 0
br i1 %26, label %46, label %27
27: ; preds = %15
%28 = and i64 %24, -2
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %43, %29 ]
%31 = phi i64 [ 0, %27 ], [ %44, %29 ]
%32 = getelementptr double, double* %10, i64 %30
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = fdiv <2 x double> %34, %21
%36 = bitcast double* %32 to <2 x double>*
store <2 x double> %35, <2 x double>* %36
%37 = or i64 %30, 2
%38 = getelementptr double, double* %10, i64 %37
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fdiv <2 x double> %40, %21
%42 = bitcast double* %38 to <2 x double>*
store <2 x double> %41, <2 x double>* %42
%43 = add nuw i64 %30, 4
%44 = add nuw i64 %31, 2
%45 = icmp eq i64 %44, %28
br i1 %45, label %46, label %29
46: ; preds = %29, %15
%47 = phi i64 [ 0, %15 ], [ %43, %29 ]
%48 = icmp eq i64 %25, 0
br i1 %48, label %55, label %49
49: ; preds = %46
%50 = getelementptr double, double* %10, i64 %47
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = fdiv <2 x double> %52, %21
%54 = bitcast double* %50 to <2 x double>*
store <2 x double> %53, <2 x double>* %54
br label %55
55: ; preds = %46, %49
%56 = icmp eq i64 %13, %16
br i1 %56, label %68, label %57
57: ; preds = %8, %55
%58 = phi double* [ %10, %8 ], [ %17, %55 ]
%59 = phi i32 [ %6, %8 ], [ %19, %55 ]
br label %60
60: ; preds = %57, %60
%61 = phi double* [ %65, %60 ], [ %58, %57 ]
%62 = phi i32 [ %66, %60 ], [ %59, %57 ]
%63 = load double, double* %61
%64 = fdiv double %63, %4
store double %64, double* %61
%65 = getelementptr inbounds double, double* %61, i64 1
%66 = add nsw i32 %62, -1
%67 = icmp eq i32 %66, 0
br i1 %67, label %68, label %60
68: ; preds = %60, %55, %3
ret double %4
}