void update_test_point_dense_SAG(double *w, double *g, long d, double stepSize,
                                 double lambda) {
  for (long j = 0; j < d; j++) {
    w[j] -= stepSize * (g[j] + lambda * w[j]);
  }
}


x86 assembly (clang, -O0):
<update_test_point_dense_SAG>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    119e <update_test_point_dense_SAG+0x9e>
movsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1122 <update_test_point_dense_SAG+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<update_test_point_dense_SAG>:
test   %rdx,%rdx
jle    123d <update_test_point_dense_SAG+0x13d>
xorpd  0xeef(%rip),%xmm0
cmp    $0x4,%rdx
jb     1131 <update_test_point_dense_SAG+0x31>
lea    (%rsi,%rdx,8),%rax
cmp    %rdi,%rax
jbe    11bf <update_test_point_dense_SAG+0xbf>
lea    (%rdi,%rdx,8),%rax
cmp    %rsi,%rax
jbe    11bf <update_test_point_dense_SAG+0xbf>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%dl
je     1161 <update_test_point_dense_SAG+0x61>
movsd  (%rdi,%rax,8),%xmm2
movapd %xmm1,%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rsi,%rax,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
or     $0x1,%rax
add    %rdx,%rcx
je     123d <update_test_point_dense_SAG+0x13d>
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm2
movapd %xmm1,%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rsi,%rax,8),%xmm3
movsd  0x8(%rdi,%rax,8),%xmm4
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
movapd %xmm1,%xmm2
mulsd  %xmm4,%xmm2
addsd  0x8(%rsi,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,0x8(%rdi,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdx
jne    1170 <update_test_point_dense_SAG+0x70>
jmp    123d <update_test_point_dense_SAG+0x13d>
mov    %rdx,%rax
and    $0xfffffffffffffffc,%rax
movapd %xmm1,%xmm8
unpcklpd %xmm1,%xmm8
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd (%rsi,%rcx,8),%xmm4
movupd 0x10(%rsi,%rcx,8),%xmm5
movupd (%rdi,%rcx,8),%xmm6
movupd 0x10(%rdi,%rcx,8),%xmm7
movapd %xmm8,%xmm2
mulpd  %xmm6,%xmm2
addpd  %xmm4,%xmm2
movapd %xmm8,%xmm4
mulpd  %xmm7,%xmm4
addpd  %xmm5,%xmm4
mulpd  %xmm3,%xmm2
addpd  %xmm6,%xmm2
mulpd  %xmm3,%xmm4
addpd  %xmm7,%xmm4
movupd %xmm2,(%rdi,%rcx,8)
movupd %xmm4,0x10(%rdi,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    11e0 <update_test_point_dense_SAG+0xe0>
cmp    %rdx,%rax
jne    1133 <update_test_point_dense_SAG+0x33>
ret

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

LLVM IR (-O2):
define void @update_test_point_dense_SAG(double* nocapture %0, double* nocapture %1, i64 %2, double %3, double %4) {
%6 = icmp sgt i64 %2, 0
br i1 %6, label %7, label %67
7: ; preds = %5
%8 = fneg double %3
%9 = icmp ult i64 %2, 4
br i1 %9, label %50, label %10
10: ; preds = %7
%11 = getelementptr double, double* %0, i64 %2
%12 = getelementptr double, double* %1, i64 %2
%13 = icmp ugt double* %12, %0
%14 = icmp ugt double* %11, %1
%15 = and i1 %13, %14
br i1 %15, label %50, label %16
16: ; preds = %10
%17 = and i64 %2, -4
%18 = insertelement <2 x double> poison, double %4, i64 0
%19 = shufflevector <2 x double> %18, <2 x double> poison, <2 x i32> zeroinitializer
%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 = insertelement <2 x double> poison, double %8, i64 0
%23 = shufflevector <2 x double> %22, <2 x double> poison, <2 x i32> zeroinitializer
%24 = insertelement <2 x double> poison, double %8, i64 0
%25 = shufflevector <2 x double> %24, <2 x double> poison, <2 x i32> zeroinitializer
br label %26
26: ; preds = %26, %16
%27 = phi i64 [ 0, %16 ], [ %46, %26 ]
%28 = getelementptr inbounds double, double* %1, i64 %27
%29 = bitcast double* %28 to <2 x double>*
%30 = load <2 x double>, <2 x double>* %29
%31 = getelementptr inbounds double, double* %28, i64 2
%32 = bitcast double* %31 to <2 x double>*
%33 = load <2 x double>, <2 x double>* %32
%34 = getelementptr inbounds double, double* %0, i64 %27
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = getelementptr inbounds double, double* %34, i64 2
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %19, <2 x double> %36, <2 x double> %30)
%41 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %21, <2 x double> %39, <2 x double> %33)
%42 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %23, <2 x double> %40, <2 x double> %36)
%43 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %25, <2 x double> %41, <2 x double> %39)
%44 = bitcast double* %34 to <2 x double>*
store <2 x double> %42, <2 x double>* %44
%45 = bitcast double* %37 to <2 x double>*
store <2 x double> %43, <2 x double>* %45
%46 = add nuw i64 %27, 4
%47 = icmp eq i64 %46, %17
br i1 %47, label %48, label %26
48: ; preds = %26
%49 = icmp eq i64 %17, %2
br i1 %49, label %67, label %50
50: ; preds = %10, %7, %48
%51 = phi i64 [ 0, %10 ], [ 0, %7 ], [ %17, %48 ]
%52 = xor i64 %51, -1
%53 = and i64 %2, 1
%54 = icmp eq i64 %53, 0
br i1 %54, label %63, label %55
55: ; preds = %50
%56 = getelementptr inbounds double, double* %1, i64 %51
%57 = load double, double* %56
%58 = getelementptr inbounds double, double* %0, i64 %51
%59 = load double, double* %58
%60 = tail call double @llvm.fmuladd.f64(double %4, double %59, double %57)
%61 = tail call double @llvm.fmuladd.f64(double %8, double %60, double %59)
store double %61, double* %58
%62 = or i64 %51, 1
br label %63
63: ; preds = %55, %50
%64 = phi i64 [ %51, %50 ], [ %62, %55 ]
%65 = sub i64 0, %2
%66 = icmp eq i64 %52, %65
br i1 %66, label %67, label %68
67: ; preds = %63, %68, %48, %5
ret void
68: ; preds = %63, %68
%69 = phi i64 [ %83, %68 ], [ %64, %63 ]
%70 = getelementptr inbounds double, double* %1, i64 %69
%71 = load double, double* %70
%72 = getelementptr inbounds double, double* %0, i64 %69
%73 = load double, double* %72
%74 = tail call double @llvm.fmuladd.f64(double %4, double %73, double %71)
%75 = tail call double @llvm.fmuladd.f64(double %8, double %74, double %73)
store double %75, double* %72
%76 = add nuw nsw i64 %69, 1
%77 = getelementptr inbounds double, double* %1, i64 %76
%78 = load double, double* %77
%79 = getelementptr inbounds double, double* %0, i64 %76
%80 = load double, double* %79
%81 = tail call double @llvm.fmuladd.f64(double %4, double %80, double %78)
%82 = tail call double @llvm.fmuladd.f64(double %8, double %81, double %80)
store double %82, double* %79
%83 = add nuw nsw i64 %69, 2
%84 = icmp eq i64 %83, %2
br i1 %84, label %67, label %68
}