double maxerror(float *A, float *B, int n) {
  int i, j;
  double error = 0.0;
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      double diff = (A[i * n + j] - B[i * n + j]) / A[i * n + j];
      if (diff < 0)
        diff = -diff;
      if (diff > error)
        error = diff;
    }
  }
  return error;
}


x86 assembly (clang, -O0):
<maxerror>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11f7 <maxerror+0xf7>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e4 <maxerror+0xe4>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
divss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    11b7 <maxerror+0xb7>
movsd  -0x30(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11d1 <maxerror+0xd1>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11d6 <maxerror+0xd6>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1131 <maxerror+0x31>
jmp    11e9 <maxerror+0xe9>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111e <maxerror+0x1e>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<maxerror>:
test   %edx,%edx
jle    1174 <maxerror+0x74>
mov    %edx,%eax
lea    0x0(,%rax,4),%r8
xorpd  %xmm0,%xmm0
xor    %edx,%edx
xorps  %xmm1,%xmm1
movaps 0xee2(%rip),%xmm2
jmp    112f <maxerror+0x2f>
add    $0x1,%rdx
add    %r8,%rdi
add    %r8,%rsi
cmp    %rax,%rdx
je     1173 <maxerror+0x73>
xor    %ecx,%ecx
movapd %xmm0,%xmm3
jmp    1151 <maxerror+0x51>
nopw   0x0(%rax,%rax,1)
maxsd  %xmm3,%xmm0
add    $0x1,%rcx
movapd %xmm0,%xmm3
cmp    %rcx,%rax
je     1120 <maxerror+0x20>
movss  (%rdi,%rcx,4),%xmm0
movaps %xmm0,%xmm4
subss  (%rsi,%rcx,4),%xmm4
divss  %xmm0,%xmm4
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
ucomiss %xmm4,%xmm1
jbe    1140 <maxerror+0x40>
xorps  %xmm2,%xmm0
jmp    1140 <maxerror+0x40>
ret
xorps  %xmm0,%xmm0
ret

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

LLVM IR (-O2):
define double @maxerror(float* nocapture %0, float* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %33
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = zext i32 %2 to i64
br label %8
8: ; preds = %30, %5
%9 = phi i64 [ 0, %5 ], [ %31, %30 ]
%10 = phi double [ 0.000000e+00, %5 ], [ %27, %30 ]
%11 = mul nsw i64 %9, %6
br label %12
12: ; preds = %8, %12
%13 = phi i64 [ 0, %8 ], [ %28, %12 ]
%14 = phi double [ %10, %8 ], [ %27, %12 ]
%15 = add nuw nsw i64 %13, %11
%16 = getelementptr inbounds float, float* %0, i64 %15
%17 = load float, float* %16
%18 = getelementptr inbounds float, float* %1, i64 %15
%19 = load float, float* %18
%20 = fsub float %17, %19
%21 = fdiv float %20, %17
%22 = fpext float %21 to double
%23 = fcmp olt float %21, 0.000000e+00
%24 = fneg double %22
%25 = select i1 %23, double %24, double %22
%26 = fcmp ogt double %25, %14
%27 = select i1 %26, double %25, double %14
%28 = add nuw nsw i64 %13, 1
%29 = icmp eq i64 %28, %7
br i1 %29, label %30, label %12
30: ; preds = %12
%31 = add nuw nsw i64 %9, 1
%32 = icmp eq i64 %31, %7
br i1 %32, label %33, label %8
33: ; preds = %30, %3
%34 = phi double [ 0.000000e+00, %3 ], [ %27, %30 ]
ret double %34
}