int matrix_add_two_rows(double *input, int rows, int cols, int row_location_1,
                        int row_location_2, double scalar) {
  int j;
  for (j = 0; j < cols; j++) {
    input[row_location_1 + j * cols] =
        input[row_location_1 + j * cols] +
        scalar * input[row_location_2 + j * cols];
  }
  return 0;
}


x86 assembly (clang, -O0):
<matrix_add_two_rows>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1190 <matrix_add_two_rows+0x90>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <matrix_add_two_rows+0x21>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrix_add_two_rows>:
push   %rbx
test   %edx,%edx
jle    1193 <matrix_add_two_rows+0x93>
mov    %edx,%r9d
movslq %ecx,%rax
movslq %r8d,%r8
cmp    $0x1,%edx
jne    111b <matrix_add_two_rows+0x1b>
xor    %edx,%edx
jmp    1174 <matrix_add_two_rows+0x74>
mov    %r9d,%r11d
and    $0xfffffffe,%r11d
lea    (%r8,%r9,1),%r10
mov    %r9,%rcx
shl    $0x4,%rcx
lea    (%rax,%r9,1),%rsi
xor    %edx,%edx
mov    %rdi,%rbx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rbx,%r8,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  (%rbx,%r10,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx,%rsi,8),%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
add    $0x2,%rdx
add    %rcx,%rbx
cmp    %rdx,%r11
jne    1140 <matrix_add_two_rows+0x40>
test   $0x1,%r9b
je     1193 <matrix_add_two_rows+0x93>
imul   %r9,%rdx
add    %rdx,%rax
add    %r8,%rdx
mulsd  (%rdi,%rdx,8),%xmm0
addsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
xor    %eax,%eax
pop    %rbx
ret

LLVM IR (-O0):
define i32 @matrix_add_two_rows(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, double %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca i32
store double* %0, double** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store double %5, double* %12
store i32 0, i32* %13
br label %14
14: ; preds = %47, %6
%15 = load i32, i32* %13
%16 = load i32, i32* %9
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %50
18: ; preds = %14
%19 = load double*, double** %7
%20 = load i32, i32* %10
%21 = load i32, i32* %13
%22 = load i32, i32* %9
%23 = mul nsw i32 %21, %22
%24 = add nsw i32 %20, %23
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %19, i64 %25
%27 = load double, double* %26
%28 = load double, double* %12
%29 = load double*, double** %7
%30 = load i32, i32* %11
%31 = load i32, i32* %13
%32 = load i32, i32* %9
%33 = mul nsw i32 %31, %32
%34 = add nsw i32 %30, %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %29, i64 %35
%37 = load double, double* %36
%38 = call double @llvm.fmuladd.f64(double %28, double %37, double %27)
%39 = load double*, double** %7
%40 = load i32, i32* %10
%41 = load i32, i32* %13
%42 = load i32, i32* %9
%43 = mul nsw i32 %41, %42
%44 = add nsw i32 %40, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %39, i64 %45
store double %38, double* %46
br label %47
47: ; preds = %18
%48 = load i32, i32* %13
%49 = add nsw i32 %48, 1
store i32 %49, i32* %13
br label %14
50: ; preds = %14
ret i32 0
}

LLVM IR (-O2):
define i32 @matrix_add_two_rows(double* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, double %5) {
%7 = icmp sgt i32 %2, 0
br i1 %7, label %8, label %51
8: ; preds = %6
%9 = zext i32 %2 to i64
%10 = sext i32 %3 to i64
%11 = sext i32 %4 to i64
%12 = and i64 %9, 1
%13 = icmp eq i32 %2, 1
br i1 %13, label %39, label %14
14: ; preds = %8
%15 = and i64 %9, 4294967294
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 0, %14 ], [ %36, %16 ]
%18 = phi i64 [ 0, %14 ], [ %37, %16 ]
%19 = mul nsw i64 %17, %9
%20 = add nsw i64 %19, %10
%21 = getelementptr inbounds double, double* %0, i64 %20
%22 = load double, double* %21
%23 = add nsw i64 %19, %11
%24 = getelementptr inbounds double, double* %0, i64 %23
%25 = load double, double* %24
%26 = tail call double @llvm.fmuladd.f64(double %5, double %25, double %22)
store double %26, double* %21
%27 = or i64 %17, 1
%28 = mul nsw i64 %27, %9
%29 = add nsw i64 %28, %10
%30 = getelementptr inbounds double, double* %0, i64 %29
%31 = load double, double* %30
%32 = add nsw i64 %28, %11
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = load double, double* %33
%35 = tail call double @llvm.fmuladd.f64(double %5, double %34, double %31)
store double %35, double* %30
%36 = add nuw nsw i64 %17, 2
%37 = add i64 %18, 2
%38 = icmp eq i64 %37, %15
br i1 %38, label %39, label %16
39: ; preds = %16, %8
%40 = phi i64 [ 0, %8 ], [ %36, %16 ]
%41 = icmp eq i64 %12, 0
br i1 %41, label %51, label %42
42: ; preds = %39
%43 = mul nsw i64 %40, %9
%44 = add nsw i64 %43, %10
%45 = getelementptr inbounds double, double* %0, i64 %44
%46 = load double, double* %45
%47 = add nsw i64 %43, %11
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = tail call double @llvm.fmuladd.f64(double %5, double %49, double %46)
store double %50, double* %45
br label %51
51: ; preds = %42, %39, %6
ret i32 0
}