void ladd_matrix(matrix, depth, scalar) double *matrix;
long depth;
double scalar;
{
  long i, j;
  for (i = 0; i < depth; i++) {
    for (j = 0; j < depth; j++) {
      matrix[i * depth + j] += scalar;
    }
  }
}


x86 assembly (clang, -O0):
<ladd_matrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1184 <ladd_matrix+0x84>
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    116e <ladd_matrix+0x6e>
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
imul   -0x10(%rbp),%rcx
add    -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    112f <ladd_matrix+0x2f>
jmp    1173 <ladd_matrix+0x73>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1119 <ladd_matrix+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ladd_matrix>:
push   %r15
push   %r14
push   %rbx
test   %rsi,%rsi
jle    123c <ladd_matrix+0x13c>
lea    -0x4(%rsi),%r9
mov    %r9,%r10
shr    $0x2,%r10
add    $0x1,%r10
mov    %rsi,%r11
and    $0xfffffffffffffffc,%r11
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %r10,%r8
and    $0xfffffffffffffffe,%r8
lea    0x30(%rdi),%rax
lea    0x0(,%rsi,8),%r14
xor    %r15d,%r15d
mov    %rdi,%rbx
jmp    1163 <ladd_matrix+0x63>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r15
add    %r14,%rax
add    %r14,%rbx
cmp    %rsi,%r15
je     123c <ladd_matrix+0x13c>
cmp    $0x4,%rsi
jae    1170 <ladd_matrix+0x70>
xor    %ecx,%ecx
jmp    1220 <ladd_matrix+0x120>
cmp    $0x4,%r9
jae    1180 <ladd_matrix+0x80>
xor    %ecx,%ecx
jmp    11d8 <ladd_matrix+0xd8>
nopw   0x0(%rax,%rax,1)
mov    %r8,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movupd -0x30(%rax,%rcx,8),%xmm2
movupd -0x20(%rax,%rcx,8),%xmm3
movupd -0x10(%rax,%rcx,8),%xmm4
movupd (%rax,%rcx,8),%xmm5
addpd  %xmm1,%xmm2
addpd  %xmm1,%xmm3
movupd %xmm2,-0x30(%rax,%rcx,8)
movupd %xmm3,-0x20(%rax,%rcx,8)
addpd  %xmm1,%xmm4
addpd  %xmm1,%xmm5
movupd %xmm4,-0x10(%rax,%rcx,8)
movupd %xmm5,(%rax,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1190 <ladd_matrix+0x90>
test   $0x1,%r10b
je     1206 <ladd_matrix+0x106>
mov    %r15,%rdx
imul   %rsi,%rdx
add    %rdx,%rcx
movupd (%rdi,%rcx,8),%xmm2
movupd 0x10(%rdi,%rcx,8),%xmm3
addpd  %xmm1,%xmm2
addpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rcx,8)
movupd %xmm3,0x10(%rdi,%rcx,8)
mov    %r11,%rcx
cmp    %rsi,%r11
je     1150 <ladd_matrix+0x50>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rbx,%rcx,8),%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rsi
jne    1220 <ladd_matrix+0x120>
jmp    1150 <ladd_matrix+0x50>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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