void xgerc(int m, int n, double alpha1, int ix0, const double y_data[],
           double A_data[], int ia0, int lda) {
  double temp;
  int i;
  int ijA;
  int ix;
  int j;
  int jA;
  int jy;
  if (alpha1 != 0.0) {
    jA = ia0;
    jy = 0;
    for (j = 0; j < n; j++) {
      if (y_data[jy] != 0.0) {
        temp = y_data[jy] * alpha1;
        ix = ix0;
        i = m + jA;
        for (ijA = jA; ijA < i; ijA++) {
          A_data[ijA - 1] += A_data[ix - 1] * temp;
          ix++;
        }
      }
      jy++;
      jA += lda;
    }
  }
}


x86 assembly (clang, -O0):
<xgerc>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    113e <xgerc+0x3e>
jp     113e <xgerc+0x3e>
jmp    1232 <xgerc+0x132>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x50(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    122d <xgerc+0x12d>
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1183 <xgerc+0x83>
jp     1183 <xgerc+0x83>
jmp    120d <xgerc+0x10d>
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x4(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1208 <xgerc+0x108>
mov    -0x28(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11af <xgerc+0xaf>
jmp    120d <xgerc+0x10d>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    0x10(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1152 <xgerc+0x52>
jmp    1232 <xgerc+0x132>
pop    %rbp
ret

x86 assembly (clang, -O2):
<xgerc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
mov    %edi,-0x64(%rsp)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1123 <xgerc+0x23>
jnp    134b <xgerc+0x24b>
test   %esi,%esi
jle    134b <xgerc+0x24b>
movslq %edx,%rax
movslq %r9d,%rcx
movslq 0x38(%rsp),%rdx
movslq -0x64(%rsp),%rdi
mov    %esi,%esi
mov    %rsi,-0x18(%rsp)
lea    -0x1(%rcx),%rsi
mov    %rsi,-0x48(%rsp)
mov    %rdi,-0x40(%rsp)
lea    (%rcx,%rdi,1),%r13
lea    0x1(%rcx),%r10
mov    %rcx,%rsi
not    %rsi
mov    %rsi,-0x50(%rsp)
mov    %rax,-0x58(%rsp)
lea    (%r8,%rax,8),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x60(%rsp)
mov    %rcx,%rbx
neg    %rbx
lea    (%r8,%rcx,8),%rbp
add    $0x8,%rbp
lea    0x0(,%rdx,8),%rax
mov    %rax,-0x30(%rsp)
xor    %edi,%edi
xorpd  %xmm1,%xmm1
mov    %r10,-0x28(%rsp)
mov    %r13,-0x20(%rsp)
mov    %rcx,-0x38(%rsp)
mov    %rdx,-0x10(%rsp)
jmp    11d5 <xgerc+0xd5>
nopl   0x0(%rax)
add    $0x1,%rdi
mov    -0x10(%rsp),%rdx
add    %rdx,%rcx
add    %rdx,%r13
add    %rdx,%r10
sub    %rdx,%rbx
add    -0x30(%rsp),%rbp
cmp    -0x18(%rsp),%rdi
je     134b <xgerc+0x24b>
cmp    %r10,%r13
mov    %r10,%rax
cmovg  %r13,%rax
mov    %rdi,%r15
imul   %rdx,%r15
mov    -0x20(%rsp),%rdx
add    %r15,%rdx
mov    %rbx,%rsi
mov    -0x28(%rsp),%rbx
lea    (%rbx,%r15,1),%r12
mov    %rsi,%rbx
cmp    %r12,%rdx
mov    %r12,%r11
cmovg  %rdx,%r11
cmovg  %rdx,%r12
mov    -0x8(%rsp),%rdx
movsd  (%rdx,%rdi,8),%xmm2
ucomisd %xmm1,%xmm2
jne    121d <xgerc+0x11d>
jnp    11b0 <xgerc+0xb0>
cmpl   $0x0,-0x64(%rsp)
jle    11b0 <xgerc+0xb0>
mov    -0x38(%rsp),%rdx
add    %r15,%rdx
sub    %rdx,%r11
mulsd  %xmm0,%xmm2
mov    %rcx,%r14
mov    -0x58(%rsp),%r9
cmp    $0x4,%r11
jb     1310 <xgerc+0x210>
mov    -0x48(%rsp),%rdx
add    %r15,%rdx
lea    (%r8,%rdx,8),%r9
mov    -0x50(%rsp),%rdx
sub    %r15,%rdx
add    %rdx,%r12
mov    -0x60(%rsp),%rdx
lea    (%rdx,%r12,8),%rdx
add    $0x8,%rdx
cmp    %rdx,%r9
jae    1289 <xgerc+0x189>
lea    (%r8,%rcx,8),%rdx
lea    (%rdx,%r12,8),%rdx
mov    %rcx,%r14
mov    -0x58(%rsp),%r9
cmp    %rdx,-0x60(%rsp)
jb     1310 <xgerc+0x210>
add    %rbx,%rax
and    $0xfffffffffffffffc,%rax
mov    %r11,%r15
and    $0xfffffffffffffffc,%r15
lea    (%rcx,%r15,1),%r14
mov    -0x58(%rsp),%rdx
lea    (%r15,%rdx,1),%r9
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
xor    %edx,%edx
mov    -0x60(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%rsi,%rdx,8),%xmm4
movupd 0x10(%rsi,%rdx,8),%xmm5
movupd -0x10(%rbp,%rdx,8),%xmm6
movupd 0x0(%rbp,%rdx,8),%xmm7
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm7,%xmm5
movupd %xmm4,-0x10(%rbp,%rdx,8)
movupd %xmm5,0x0(%rbp,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rax
jne    12c0 <xgerc+0x1c0>
cmp    %r15,%r11
je     11b0 <xgerc+0xb0>
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x40(%rsp),%rax
add    %rcx,%rax
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%r8,%r9,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  -0x8(%r8,%r14,8),%xmm3
movsd  %xmm3,-0x8(%r8,%r14,8)
add    $0x1,%r9
add    $0x1,%r14
cmp    %rax,%r14
jl     1320 <xgerc+0x220>
jmp    11b0 <xgerc+0xb0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @xgerc(i32 %0, i32 %1, double %2, i32 %3, double* nocapture %4, double* nocapture %5, i32 %6, i32 %7) {
%9 = fcmp une double %2, 0.000000e+00
%10 = icmp sgt i32 %1, 0
%11 = and i1 %9, %10
br i1 %11, label %12, label %116
12: ; preds = %8
%13 = sext i32 %3 to i64
%14 = sext i32 %6 to i64
%15 = sext i32 %7 to i64
%16 = sext i32 %0 to i64
%17 = zext i32 %1 to i64
%18 = icmp sgt i32 %0, 0
%19 = add nsw i64 %14, -1
%20 = add nsw i64 %14, %16
%21 = add nsw i64 %14, 1
%22 = xor i64 %14, -1
%23 = add nsw i64 %13, -1
%24 = getelementptr double, double* %5, i64 %23
%25 = getelementptr double, double* %5, i64 %13
%26 = add nsw i64 %14, %16
%27 = add nsw i64 %14, 1
br label %28
28: ; preds = %12, %112
%29 = phi i64 [ 0, %12 ], [ %113, %112 ]
%30 = phi i64 [ %14, %12 ], [ %114, %112 ]
%31 = mul i64 %29, %15
%32 = add i64 %26, %31
%33 = add i64 %27, %31
%34 = call i64 @llvm.smax.i64(i64 %32, i64 %33)
%35 = mul i64 %29, %15
%36 = add i64 %35, %14
%37 = sub i64 %34, %36
%38 = mul i64 %29, %15
%39 = add i64 %19, %38
%40 = getelementptr double, double* %5, i64 %39
%41 = getelementptr double, double* %5, i64 %30
%42 = add i64 %20, %38
%43 = add i64 %21, %38
%44 = call i64 @llvm.smax.i64(i64 %42, i64 %43)
%45 = mul i64 %29, %15
%46 = sub i64 %22, %45
%47 = add i64 %44, %46
%48 = getelementptr double, double* %41, i64 %47
%49 = getelementptr double, double* %25, i64 %47
%50 = getelementptr inbounds double, double* %4, i64 %29
%51 = load double, double* %50
%52 = fcmp une double %51, 0.000000e+00
br i1 %52, label %53, label %112
53: ; preds = %28
%54 = fmul double %51, %2
%55 = add nsw i64 %30, %16
br i1 %18, label %56, label %112
56: ; preds = %53
%57 = icmp ult i64 %37, 4
br i1 %57, label %96, label %58
58: ; preds = %56
%59 = icmp ult double* %40, %49
%60 = icmp ult double* %24, %48
%61 = and i1 %59, %60
br i1 %61, label %96, label %62
62: ; preds = %58
%63 = and i64 %37, -4
%64 = add i64 %30, %63
%65 = add i64 %63, %13
%66 = insertelement <2 x double> poison, double %54, i64 0
%67 = shufflevector <2 x double> %66, <2 x double> poison, <2 x i32> zeroinitializer
%68 = insertelement <2 x double> poison, double %54, i64 0
%69 = shufflevector <2 x double> %68, <2 x double> poison, <2 x i32> zeroinitializer
br label %70
70: ; preds = %70, %62
%71 = phi i64 [ 0, %62 ], [ %92, %70 ]
%72 = add i64 %30, %71
%73 = add i64 %71, %13
%74 = add nsw i64 %73, -1
%75 = getelementptr inbounds double, double* %5, i64 %74
%76 = bitcast double* %75 to <2 x double>*
%77 = load <2 x double>, <2 x double>* %76
%78 = getelementptr inbounds double, double* %75, i64 2
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = add nsw i64 %72, -1
%82 = getelementptr inbounds double, double* %5, i64 %81
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = getelementptr inbounds double, double* %82, i64 2
%86 = bitcast double* %85 to <2 x double>*
%87 = load <2 x double>, <2 x double>* %86
%88 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %77, <2 x double> %67, <2 x double> %84)
%89 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %80, <2 x double> %69, <2 x double> %87)
%90 = bitcast double* %82 to <2 x double>*
store <2 x double> %88, <2 x double>* %90
%91 = bitcast double* %85 to <2 x double>*
store <2 x double> %89, <2 x double>* %91
%92 = add nuw i64 %71, 4
%93 = icmp eq i64 %92, %63
br i1 %93, label %94, label %70
94: ; preds = %70
%95 = icmp eq i64 %37, %63
br i1 %95, label %112, label %96
96: ; preds = %58, %56, %94
%97 = phi i64 [ %30, %58 ], [ %30, %56 ], [ %64, %94 ]
%98 = phi i64 [ %13, %58 ], [ %13, %56 ], [ %65, %94 ]
br label %99
99: ; preds = %96, %99
%100 = phi i64 [ %110, %99 ], [ %97, %96 ]
%101 = phi i64 [ %109, %99 ], [ %98, %96 ]
%102 = add nsw i64 %101, -1
%103 = getelementptr inbounds double, double* %5, i64 %102
%104 = load double, double* %103
%105 = add nsw i64 %100, -1
%106 = getelementptr inbounds double, double* %5, i64 %105
%107 = load double, double* %106
%108 = tail call double @llvm.fmuladd.f64(double %104, double %54, double %107)
store double %108, double* %106
%109 = add nsw i64 %101, 1
%110 = add nsw i64 %100, 1
%111 = icmp slt i64 %110, %55
br i1 %111, label %99, label %112
112: ; preds = %99, %94, %53, %28
%113 = add nuw nsw i64 %29, 1
%114 = add i64 %30, %15
%115 = icmp eq i64 %113, %17
br i1 %115, label %116, label %28
116: ; preds = %112, %8
ret void
}