void crossxy(double **x, double *y, int n, int k, double *xpy) {
  int i, j;
  double *xrow, yi;
  for (j = 0; j < k; j++) {
    xpy[j] = 0.0;
  }
  for (i = 0; i < n; i++) {
    xrow = x[i];
    yi = y[i];
    for (j = 0; j < k; j++) {
      xpy[j] += xrow[j] * y[i];
    }
  }
  return;
}


x86 assembly (clang, -O0):
<crossxy>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1147 <crossxy+0x47>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111d <crossxy+0x1d>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e4 <crossxy+0xe4>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11d1 <crossxy+0xd1>
mov    -0x30(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1183 <crossxy+0x83>
jmp    11d6 <crossxy+0xd6>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    114e <crossxy+0x4e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<crossxy>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8,%rbx
mov    %ecx,%r14d
mov    %edx,%r12d
mov    %rsi,%rbp
mov    %rdi,%r15
mov    %ecx,%r13d
test   %ecx,%ecx
jle    1146 <crossxy+0x36>
lea    0x0(,%r13,8),%rdx
mov    %rbx,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %r12d,%r12d
jle    12eb <crossxy+0x1db>
mov    %r12d,%r11d
lea    (%rbx,%r13,8),%rax
mov    %rax,0x8(%rsp)
mov    %r13d,%eax
and    $0xfffffffc,%eax
mov    %rax,(%rsp)
mov    %r13,%rax
neg    %rax
mov    %rax,0x10(%rsp)
lea    0x0(,%r13,8),%rdi
and    $0xffffffffffffffe0,%rdi
xor    %edx,%edx
jmp    118d <crossxy+0x7d>
add    $0x1,%rdx
cmp    %r11,%rdx
je     12eb <crossxy+0x1db>
test   %r14d,%r14d
jle    1180 <crossxy+0x70>
mov    (%r15,%rdx,8),%rcx
cmp    $0x4,%r14d
jae    11b0 <crossxy+0xa0>
xor    %eax,%eax
jmp    1270 <crossxy+0x160>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r15,%r9
lea    0x8(,%rdx,8),%rax
add    %rbp,%rax
lea    (%rcx,%r13,8),%rsi
cmp    %rbx,%rsi
seta   %r8b
mov    0x8(%rsp),%rsi
cmp    %rsi,%rcx
setb   %r15b
cmp    %rbx,%rax
seta   %r10b
lea    0x0(,%rdx,8),%rax
add    %rbp,%rax
cmp    %rsi,%rax
setb   %r12b
test   %r15b,%r8b
jne    125c <crossxy+0x14c>
mov    $0x0,%eax
and    %r12b,%r10b
mov    %r9,%r15
jne    1270 <crossxy+0x160>
movsd  0x0(%rbp,%rdx,8),%xmm0
unpcklpd %xmm0,%xmm0
xor    %eax,%eax
nopl   0x0(%rax)
movupd (%rcx,%rax,1),%xmm1
movupd 0x10(%rcx,%rax,1),%xmm2
movupd (%rbx,%rax,1),%xmm3
movupd 0x10(%rbx,%rax,1),%xmm4
mulpd  %xmm0,%xmm1
addpd  %xmm3,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm1,(%rbx,%rax,1)
movupd %xmm2,0x10(%rbx,%rax,1)
add    $0x20,%rax
cmp    %rax,%rdi
jne    1210 <crossxy+0x100>
mov    (%rsp),%rsi
mov    %rsi,%rax
cmp    %r13,%rsi
je     1180 <crossxy+0x70>
jmp    1270 <crossxy+0x160>
xor    %eax,%eax
mov    %r9,%r15
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rax,%rsi
test   $0x1,%r13b
je     1295 <crossxy+0x185>
movsd  (%rcx,%rax,8),%xmm0
mulsd  0x0(%rbp,%rdx,8),%xmm0
addsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
mov    %rax,%rsi
or     $0x1,%rsi
not    %rax
cmp    0x10(%rsp),%rax
je     1180 <crossxy+0x70>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx,%rsi,8),%xmm0
mulsd  0x0(%rbp,%rdx,8),%xmm0
addsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,(%rbx,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm0
mulsd  0x0(%rbp,%rdx,8),%xmm0
addsd  0x8(%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%r13
jne    12b0 <crossxy+0x1a0>
jmp    1180 <crossxy+0x70>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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