double reduceElement(double **A, int i, int j, int n, int l) {
  int k;
  double sum;
  sum = A[i][j];
  for (k = 1; k < l; k++)
    sum -= A[i][k] * A[k][j];
  return sum;
}


x86 assembly (clang, -O0):
<reduceElement>:
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)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11a3 <reduceElement+0xa3>
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1136 <reduceElement+0x36>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<reduceElement>:
movslq %esi,%rax
mov    (%rdi,%rax,8),%rcx
movslq %edx,%rax
movsd  (%rcx,%rax,8),%xmm0
cmp    $0x2,%r8d
jl     1190 <reduceElement+0x90>
mov    %r8d,%r9d
add    $0xffffffffffffffff,%r9
cmp    $0x2,%r8d
jne    1129 <reduceElement+0x29>
mov    $0x1,%esi
jmp    1178 <reduceElement+0x78>
mov    %r9,%r8
and    $0xfffffffffffffffe,%r8
neg    %r8
mov    $0x1,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,8),%xmm1
movsd  0x8(%rcx,%rsi,8),%xmm2
mov    (%rdi,%rsi,8),%rdx
mulsd  (%rdx,%rax,8),%xmm1
mov    0x8(%rdi,%rsi,8),%rdx
subsd  %xmm1,%xmm0
mulsd  (%rdx,%rax,8),%xmm2
subsd  %xmm2,%xmm0
lea    (%r8,%rsi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rsi
cmp    $0x1,%rdx
jne    1140 <reduceElement+0x40>
test   $0x1,%r9b
je     1190 <reduceElement+0x90>
movsd  (%rcx,%rsi,8),%xmm1
mov    (%rdi,%rsi,8),%rcx
mulsd  (%rcx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
ret

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

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