void scale_beta(double *beta2, double *beta1, double *mean, double *sd, int p) {
  int j;
  double t, s = 0;
  for (j = p - 1; j >= 0; --j) {
    t = beta1[j] * mean[j];
    if (sd[j] != 0)
      t /= sd[j];
    beta2[j] = beta1[j] * sd[j];
    s += t;
  }
  beta2[0] += s;
}


x86 assembly (clang, -O0):
<scale_beta>:
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)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x28(%rbp)
jl     11d6 <scale_beta+0xd6>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1177 <scale_beta+0x77>
jp     1177 <scale_beta+0x77>
jmp    1192 <scale_beta+0x92>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    1129 <scale_beta+0x29>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<scale_beta>:
test   %r8d,%r8d
jle    1171 <scale_beta+0x71>
mov    %r8d,%eax
add    $0x1,%rax
xorpd  %xmm1,%xmm1
movsd  0xee8(%rip),%xmm2
xorpd  %xmm0,%xmm0
nopl   0x0(%rax)
lea    -0x2(%rax),%r8d
movsd  (%rsi,%r8,8),%xmm3
movsd  (%rdx,%r8,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  (%rcx,%r8,8),%xmm5
movapd %xmm5,%xmm6
cmpneqsd %xmm1,%xmm6
movapd %xmm6,%xmm7
andpd  %xmm5,%xmm7
andnpd %xmm2,%xmm6
orpd   %xmm7,%xmm6
divsd  %xmm6,%xmm4
addsd  %xmm4,%xmm0
mulsd  %xmm5,%xmm3
movsd  %xmm3,(%rdi,%r8,8)
add    $0xffffffffffffffff,%rax
cmp    $0x1,%rax
ja     1120 <scale_beta+0x20>
jmp    1175 <scale_beta+0x75>
xorpd  %xmm0,%xmm0
addsd  (%rdi),%xmm0
movsd  %xmm0,(%rdi)
ret

LLVM IR (-O0):
define void @scale_beta(double* %0, double* %1, double* %2, double* %3, i32 %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
store double* %0, double** %6
store double* %1, double** %7
store double* %2, double** %8
store double* %3, double** %9
store i32 %4, i32* %10
store double 0.000000e+00, double* %13
%14 = load i32, i32* %10
%15 = sub nsw i32 %14, 1
store i32 %15, i32* %11
br label %16
16: ; preds = %64, %5
%17 = load i32, i32* %11
%18 = icmp sge i32 %17, 0
br i1 %18, label %19, label %67
19: ; preds = %16
%20 = load double*, double** %7
%21 = load i32, i32* %11
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = load double*, double** %8
%26 = load i32, i32* %11
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = fmul double %24, %29
store double %30, double* %12
%31 = load double*, double** %9
%32 = load i32, i32* %11
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
%36 = fcmp une double %35, 0.000000e+00
br i1 %36, label %37, label %45
37: ; preds = %19
%38 = load double*, double** %9
%39 = load i32, i32* %11
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = load double, double* %12
%44 = fdiv double %43, %42
store double %44, double* %12
br label %45
45: ; preds = %37, %19
%46 = load double*, double** %7
%47 = load i32, i32* %11
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
%50 = load double, double* %49
%51 = load double*, double** %9
%52 = load i32, i32* %11
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
%55 = load double, double* %54
%56 = fmul double %50, %55
%57 = load double*, double** %6
%58 = load i32, i32* %11
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %56, double* %60
%61 = load double, double* %12
%62 = load double, double* %13
%63 = fadd double %62, %61
store double %63, double* %13
br label %64
64: ; preds = %45
%65 = load i32, i32* %11
%66 = add nsw i32 %65, -1
store i32 %66, i32* %11
br label %16
67: ; preds = %16
%68 = load double, double* %13
%69 = load double*, double** %6
%70 = getelementptr inbounds double, double* %69, i64 0
%71 = load double, double* %70
%72 = fadd double %71, %68
store double %72, double* %70
ret void
}

LLVM IR (-O2):
define void @scale_beta(double* nocapture %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %28
7: ; preds = %5
%8 = zext i32 %4 to i64
br label %9
9: ; preds = %7, %9
%10 = phi i64 [ %8, %7 ], [ %12, %9 ]
%11 = phi double [ 0.000000e+00, %7 ], [ %26, %9 ]
%12 = add nsw i64 %10, -1
%13 = and i64 %12, 4294967295
%14 = getelementptr inbounds double, double* %1, i64 %13
%15 = load double, double* %14
%16 = getelementptr inbounds double, double* %2, i64 %13
%17 = load double, double* %16
%18 = fmul double %15, %17
%19 = getelementptr inbounds double, double* %3, i64 %13
%20 = load double, double* %19
%21 = fcmp une double %20, 0.000000e+00
%22 = select i1 %21, double %20, double 1.000000e+00
%23 = fdiv double %18, %22
%24 = fmul double %15, %20
%25 = getelementptr inbounds double, double* %0, i64 %13
store double %24, double* %25
%26 = fadd double %11, %23
%27 = icmp ugt i64 %10, 1
br i1 %27, label %9, label %28
28: ; preds = %9, %5
%29 = phi double [ 0.000000e+00, %5 ], [ %26, %9 ]
%30 = load double, double* %0
%31 = fadd double %29, %30
store double %31, double* %0
ret void
}