void Column_Transformation(double *A, double x, int col1, int col2, int nrows,
                           int ncols) {
  int i;
  double *pA1, *pA2;
  if (col1 != col2) {
    pA1 = A + col1;
    pA2 = A + col2;
    for (i = 0; i < nrows; pA1 += ncols, pA2 += ncols, i++)
      *pA2 += x * *pA1;
  } else {
    pA1 = A + col1;
    for (i = 0; i < nrows; pA1 += ncols, i++)
      *pA1 += x * *pA1;
  }
}


x86 assembly (clang, -O0):
<Column_Transformation>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
je     11bd <Column_Transformation+0xbd>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11b8 <Column_Transformation+0xb8>
movsd  -0x10(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1153 <Column_Transformation+0x53>
jmp    122c <Column_Transformation+0x12c>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1227 <Column_Transformation+0x127>
movsd  -0x10(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11d7 <Column_Transformation+0xd7>
jmp    122c <Column_Transformation+0x12c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Column_Transformation>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %edx,%esi
jne    11f1 <Column_Transformation+0xf1>
test   %ecx,%ecx
jle    12fb <Column_Transformation+0x1fb>
movslq %esi,%rax
lea    (%rdi,%rax,8),%rsi
movslq %r8d,%rax
lea    -0x1(%rcx),%edi
mov    %ecx,%edx
and    $0x3,%edx
cmp    $0x3,%edi
jb     11c2 <Column_Transformation+0xc2>
and    $0xfffffffc,%ecx
mov    %rax,%r8
shl    $0x5,%r8
lea    0x0(,%rax,8),%rbx
nopw   0x0(%rax,%rax,1)
mov    %rsi,%rdi
movsd  (%rsi),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi)
add    %rbx,%rsi
movsd  (%rdi,%rbx,1),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rdi,%rbx,1)
movsd  (%rbx,%rsi,1),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
lea    (%rsi,%rbx,1),%rbp
movsd  %xmm2,(%rbx,%rsi,1)
movsd  (%rbx,%rbp,1),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
lea    (%rbx,%rbp,1),%rsi
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rbp,1)
add    %rbx,%rsi
add    $0xfffffffc,%ecx
jne    1150 <Column_Transformation+0x50>
add    %r8,%rdi
mov    %rdi,%rsi
test   %edx,%edx
je     12fb <Column_Transformation+0x1fb>
shl    $0x3,%rax
xchg   %ax,%ax
movsd  (%rsi),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi)
add    %rax,%rsi
add    $0xffffffff,%edx
jne    11d0 <Column_Transformation+0xd0>
jmp    12fb <Column_Transformation+0x1fb>
test   %ecx,%ecx
jle    12fb <Column_Transformation+0x1fb>
movslq %edx,%r9
movslq %esi,%r10
movslq %r8d,%r8
lea    -0x1(%rcx),%eax
mov    %ecx,%edx
and    $0x3,%edx
cmp    $0x3,%eax
jb     12c9 <Column_Transformation+0x1c9>
and    $0xfffffffc,%ecx
lea    0x0(,%r9,8),%r14
mov    %r8,%r11
shl    $0x5,%r11
lea    0x0(,%r8,8),%rbp
lea    0x0(,%rbp,2),%r12
add    %rbp,%r12
lea    (%r12,%r9,8),%rax
mov    %r8,%r13
shl    $0x4,%r13
lea    0x0(,%r9,8),%rbx
add    %r13,%rbx
lea    0x0(%rbp,%r9,8),%rsi
lea    0x0(,%r10,8),%r15
lea    (%r12,%r10,8),%r12
lea    0x0(%r13,%r10,8),%r13
lea    0x0(%rbp,%r10,8),%rbp
nopl   0x0(%rax)
movsd  (%rdi,%r15,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%r14,1),%xmm1
movsd  %xmm1,(%rdi,%r14,1)
movsd  (%rdi,%rbp,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%rsi,1),%xmm1
movsd  %xmm1,(%rdi,%rsi,1)
movsd  (%rdi,%r13,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%rbx,1),%xmm1
movsd  %xmm1,(%rdi,%rbx,1)
movsd  (%rdi,%r12,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%rax,1),%xmm1
movsd  %xmm1,(%rdi,%rax,1)
add    %r11,%rdi
add    $0xfffffffc,%ecx
jne    1270 <Column_Transformation+0x170>
lea    (%rdi,%r9,8),%rax
lea    (%rdi,%r10,8),%rcx
test   %edx,%edx
je     12fb <Column_Transformation+0x1fb>
shl    $0x3,%r8
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax,%rsi,1),%xmm1
movsd  %xmm1,(%rax,%rsi,1)
add    %r8,%rsi
add    $0xffffffff,%edx
jne    12e0 <Column_Transformation+0x1e0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @Column_Transformation(double* nocapture %0, double %1, i32 %2, i32 %3, i32 %4, i32 %5) {
%7 = icmp eq i32 %2, %3
%8 = icmp sgt i32 %4, 0
br i1 %7, label %47, label %9
9: ; preds = %6
br i1 %8, label %10, label %100
10: ; preds = %9
%11 = sext i32 %3 to i64
%12 = getelementptr inbounds double, double* %0, i64 %11
%13 = sext i32 %2 to i64
%14 = getelementptr inbounds double, double* %0, i64 %13
%15 = sext i32 %5 to i64
%16 = add i32 %4, -1
%17 = and i32 %4, 3
%18 = icmp ult i32 %16, 3
br i1 %18, label %85, label %19
19: ; preds = %10
%20 = and i32 %4, -4
br label %21
21: ; preds = %21, %19
%22 = phi double* [ %12, %19 ], [ %44, %21 ]
%23 = phi double* [ %14, %19 ], [ %43, %21 ]
%24 = phi i32 [ 0, %19 ], [ %45, %21 ]
%25 = load double, double* %23
%26 = load double, double* %22
%27 = tail call double @llvm.fmuladd.f64(double %1, double %25, double %26)
store double %27, double* %22
%28 = getelementptr inbounds double, double* %23, i64 %15
%29 = getelementptr inbounds double, double* %22, i64 %15
%30 = load double, double* %28
%31 = load double, double* %29
%32 = tail call double @llvm.fmuladd.f64(double %1, double %30, double %31)
store double %32, double* %29
%33 = getelementptr inbounds double, double* %28, i64 %15
%34 = getelementptr inbounds double, double* %29, i64 %15
%35 = load double, double* %33
%36 = load double, double* %34
%37 = tail call double @llvm.fmuladd.f64(double %1, double %35, double %36)
store double %37, double* %34
%38 = getelementptr inbounds double, double* %33, i64 %15
%39 = getelementptr inbounds double, double* %34, i64 %15
%40 = load double, double* %38
%41 = load double, double* %39
%42 = tail call double @llvm.fmuladd.f64(double %1, double %40, double %41)
store double %42, double* %39
%43 = getelementptr inbounds double, double* %38, i64 %15
%44 = getelementptr inbounds double, double* %39, i64 %15
%45 = add i32 %24, 4
%46 = icmp eq i32 %45, %20
br i1 %46, label %85, label %21
47: ; preds = %6
br i1 %8, label %48, label %100
48: ; preds = %47
%49 = sext i32 %2 to i64
%50 = getelementptr inbounds double, double* %0, i64 %49
%51 = sext i32 %5 to i64
%52 = add i32 %4, -1
%53 = and i32 %4, 3
%54 = icmp ult i32 %52, 3
br i1 %54, label %74, label %55
55: ; preds = %48
%56 = and i32 %4, -4
br label %57
57: ; preds = %57, %55
%58 = phi double* [ %50, %55 ], [ %71, %57 ]
%59 = phi i32 [ 0, %55 ], [ %72, %57 ]
%60 = load double, double* %58
%61 = tail call double @llvm.fmuladd.f64(double %1, double %60, double %60)
store double %61, double* %58
%62 = getelementptr inbounds double, double* %58, i64 %51
%63 = load double, double* %62
%64 = tail call double @llvm.fmuladd.f64(double %1, double %63, double %63)
store double %64, double* %62
%65 = getelementptr inbounds double, double* %62, i64 %51
%66 = load double, double* %65
%67 = tail call double @llvm.fmuladd.f64(double %1, double %66, double %66)
store double %67, double* %65
%68 = getelementptr inbounds double, double* %65, i64 %51
%69 = load double, double* %68
%70 = tail call double @llvm.fmuladd.f64(double %1, double %69, double %69)
store double %70, double* %68
%71 = getelementptr inbounds double, double* %68, i64 %51
%72 = add i32 %59, 4
%73 = icmp eq i32 %72, %56
br i1 %73, label %74, label %57
74: ; preds = %57, %48
%75 = phi double* [ %50, %48 ], [ %71, %57 ]
%76 = icmp eq i32 %53, 0
br i1 %76, label %100, label %77
77: ; preds = %74, %77
%78 = phi double* [ %82, %77 ], [ %75, %74 ]
%79 = phi i32 [ %83, %77 ], [ 0, %74 ]
%80 = load double, double* %78
%81 = tail call double @llvm.fmuladd.f64(double %1, double %80, double %80)
store double %81, double* %78
%82 = getelementptr inbounds double, double* %78, i64 %51
%83 = add i32 %79, 1
%84 = icmp eq i32 %83, %53
br i1 %84, label %100, label %77
85: ; preds = %21, %10
%86 = phi double* [ %12, %10 ], [ %44, %21 ]
%87 = phi double* [ %14, %10 ], [ %43, %21 ]
%88 = icmp eq i32 %17, 0
br i1 %88, label %100, label %89
89: ; preds = %85, %89
%90 = phi double* [ %97, %89 ], [ %86, %85 ]
%91 = phi double* [ %96, %89 ], [ %87, %85 ]
%92 = phi i32 [ %98, %89 ], [ 0, %85 ]
%93 = load double, double* %91
%94 = load double, double* %90
%95 = tail call double @llvm.fmuladd.f64(double %1, double %93, double %94)
store double %95, double* %90
%96 = getelementptr inbounds double, double* %91, i64 %15
%97 = getelementptr inbounds double, double* %90, i64 %15
%98 = add i32 %92, 1
%99 = icmp eq i32 %98, %17
br i1 %99, label %100, label %89
100: ; preds = %85, %89, %74, %77, %9, %47
ret void
}