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


x86 assembly (clang, -O0):
<reduceElement2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%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)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1142 <reduceElement2+0x42>
mov    -0xc(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    1148 <reduceElement2+0x48>
mov    -0x10(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     11c5 <reduceElement2+0xc5>
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%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    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1158 <reduceElement2+0x58>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<reduceElement2>:
movslq %esi,%rsi
mov    (%rdi,%rsi,8),%rcx
movslq %edx,%rax
movsd  (%rcx,%rax,8),%xmm0
cmp    %eax,%esi
cmovge %eax,%esi
cmp    $0x2,%esi
jl     1190 <reduceElement2+0x90>
mov    %esi,%r8d
add    $0xffffffffffffffff,%r8
cmp    $0x2,%esi
jne    112c <reduceElement2+0x2c>
mov    $0x1,%esi
jmp    1178 <reduceElement2+0x78>
mov    %r8,%r9
and    $0xfffffffffffffffe,%r9
neg    %r9
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    (%r9,%rsi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rsi
cmp    $0x1,%rdx
jne    1140 <reduceElement2+0x40>
test   $0x1,%r8b
je     1190 <reduceElement2+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 @reduceElement2(double** %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca double**
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
store double** %0, double*** %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%12 = load double**, double*** %5
%13 = load i32, i32* %6
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double*, double** %12, i64 %14
%16 = load double*, double** %15
%17 = load i32, i32* %7
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %16, i64 %18
%20 = load double, double* %19
store double %20, double* %11
%21 = load i32, i32* %6
%22 = load i32, i32* %7
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %26
24: ; preds = %4
%25 = load i32, i32* %6
br label %28
26: ; preds = %4
%27 = load i32, i32* %7
br label %28
28: ; preds = %26, %24
%29 = phi i32 [ %25, %24 ], [ %27, %26 ]
%30 = sub nsw i32 %29, 1
store i32 %30, i32* %10
store i32 1, i32* %9
br label %31
31: ; preds = %57, %28
%32 = load i32, i32* %9
%33 = load i32, i32* %10
%34 = icmp sle i32 %32, %33
br i1 %34, label %35, label %60
35: ; preds = %31
%36 = load double**, double*** %5
%37 = load i32, i32* %6
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double*, double** %36, i64 %38
%40 = load double*, double** %39
%41 = load i32, i32* %9
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
%44 = load double, double* %43
%45 = load double**, double*** %5
%46 = load i32, i32* %9
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double*, double** %45, i64 %47
%49 = load double*, double** %48
%50 = load i32, i32* %7
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = load double, double* %11
%55 = fneg double %44
%56 = call double @llvm.fmuladd.f64(double %55, double %53, double %54)
store double %56, double* %11
br label %57
57: ; preds = %35
%58 = load i32, i32* %9
%59 = add nsw i32 %58, 1
store i32 %59, i32* %9
br label %31
60: ; preds = %31
%61 = load double, double* %11
ret double %61
}

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