#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double const fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, double, double, int,
                           int);
int /*<<< orphan*/ stderr;
void lm_check_inverse(double *const A, double *const IA, const int n,
                      const double Tol, int *const failure) {
  for (int j = 0; j < n; j++) {
    for (int i = 0; i < n; i++) {
      double s = 0;
      for (int k = 0; k < n; k++)
        s += IA[j * n + k] * A[k * n + i];
      double expected = (i == j) ? 1. : 0.;
      if (fabs(s - expected) > Tol) {
        fprintf(stderr, "lm_check_inverse found %g instead of %g at %i,%i \n",
                s, expected, j, i);
        *failure = 22;
        return;
      }
    }
  }
}


x86 assembly (clang, -O0):
<lm_check_inverse>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1274 <lm_check_inverse+0x154>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1261 <lm_check_inverse+0x141>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c9 <lm_check_inverse+0xa9>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1171 <lm_check_inverse+0x51>
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
movsd  0xe29(%rip),%xmm0
xorps  %xmm1,%xmm1
movsd  %xmm1,-0x58(%rbp)
cmp    %ecx,%eax
movsd  %xmm0,-0x50(%rbp)
je     11f6 <lm_check_inverse+0xd6>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
call   1030 <fabs@plt>
ucomisd -0x20(%rbp),%xmm0
jbe    124e <lm_check_inverse+0x12e>
mov    0x2dbf(%rip),%rax
mov    (%rax),%edi
movsd  -0x38(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mov    -0x2c(%rbp),%edx
mov    -0x30(%rbp),%ecx
lea    0xdce(%rip),%rsi
call   1040 <fprintf@plt>
mov    -0x28(%rbp),%rax
movl   $0x16,(%rax)
jmp    1274 <lm_check_inverse+0x154>
jmp    1253 <lm_check_inverse+0x133>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1156 <lm_check_inverse+0x36>
jmp    1266 <lm_check_inverse+0x146>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1143 <lm_check_inverse+0x23>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lm_check_inverse>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rcx,0x20(%rsp)
mov    %rsi,(%rsp)
test   %edx,%edx
jle    12b9 <lm_check_inverse+0x1a9>
movapd %xmm0,%xmm3
mov    %edx,%r8d
mov    %rdi,%r10
mov    %edx,%r13d
mov    %r13d,%r14d
and    $0xfffffffe,%r14d
mov    (%rsp),%rax
lea    0x8(%rax),%r11
mov    %r13,%r15
shl    $0x4,%r15
xor    %ebx,%ebx
movsd  0xea4(%rip),%xmm4
movapd 0xeac(%rip),%xmm5
movsd  %xmm0,0x18(%rsp)
mov    %rdi,0x10(%rsp)
mov    %r11,0x8(%rsp)
jmp    11df <lm_check_inverse+0xcf>
cs nopw 0x0(%rax,%rax,1)
mov    0x2e59(%rip),%rax
mov    (%rax),%edi
lea    0xe90(%rip),%rsi
mov    %ebx,%edx
mov    %r8d,%ebp
call   1030 <fprintf@plt>
movapd 0xe6e(%rip),%xmm5
movsd  0xe56(%rip),%xmm4
mov    0x8(%rsp),%r11
mov    0x10(%rsp),%r10
mov    %ebp,%r8d
movsd  0x18(%rsp),%xmm3
mov    0x20(%rsp),%rax
movl   $0x16,(%rax)
add    $0x1,%rbx
cmp    %r13,%rbx
sete   %al
or     %al,%r12b
test   $0x1,%r12b
jne    12b9 <lm_check_inverse+0x1a9>
mov    %rbx,%r9
imul   %r13,%r9
lea    (%r11,%r9,8),%rdx
mov    $0x1,%r12b
mov    %r10,%rbp
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmp    $0x1,%r8d
jne    1210 <lm_check_inverse+0x100>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
jmp    124d <lm_check_inverse+0x13d>
xchg   %ax,%ax
xorpd  %xmm1,%xmm1
mov    %rbp,%rsi
xor    %edi,%edi
nopl   0x0(%rax)
movsd  -0x8(%rdx,%rdi,8),%xmm2
mulsd  (%rsi),%xmm2
movsd  (%rdx,%rdi,8),%xmm0
addsd  %xmm1,%xmm2
mulsd  (%rsi,%r13,8),%xmm0
addsd  %xmm2,%xmm0
add    $0x2,%rdi
add    %r15,%rsi
movapd %xmm0,%xmm1
cmp    %rdi,%r14
jne    1220 <lm_check_inverse+0x110>
test   $0x1,%r13b
je     1271 <lm_check_inverse+0x161>
lea    (%rdi,%r9,1),%rsi
mov    (%rsp),%rax
movsd  (%rax,%rsi,8),%xmm1
imul   %r13,%rdi
add    %rcx,%rdi
mulsd  (%r10,%rdi,8),%xmm1
addsd  %xmm1,%xmm0
movapd %xmm4,%xmm1
cmp    %rbx,%rcx
je     127e <lm_check_inverse+0x16e>
xorpd  %xmm1,%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
andpd  %xmm5,%xmm2
ucomisd %xmm3,%xmm2
ja     1180 <lm_check_inverse+0x70>
add    $0x1,%rcx
cmp    %r13,%rcx
setb   %r12b
add    $0x8,%rbp
cmp    %r13,%rcx
jne    1200 <lm_check_inverse+0xf0>
add    $0x1,%rbx
cmp    %r13,%rbx
jne    11df <lm_check_inverse+0xcf>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [51 x i8] c"lm_check_inverse found %g instead of %g at %i,%i \0A\00"
define void @lm_check_inverse(double* %0, double* %1, i32 %2, double %3, i32* %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca double
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca i32
%15 = alloca double
store double* %0, double** %6
store double* %1, double** %7
store i32 %2, i32* %8
store double %3, double* %9
store i32* %4, i32** %10
store i32 0, i32* %11
br label %16
16: ; preds = %79, %5
%17 = load i32, i32* %11
%18 = load i32, i32* %8
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %82
20: ; preds = %16
store i32 0, i32* %12
br label %21
21: ; preds = %75, %20
%22 = load i32, i32* %12
%23 = load i32, i32* %8
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %78
25: ; preds = %21
store double 0.000000e+00, double* %13
store i32 0, i32* %14
br label %26
26: ; preds = %51, %25
%27 = load i32, i32* %14
%28 = load i32, i32* %8
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %54
30: ; preds = %26
%31 = load double*, double** %7
%32 = load i32, i32* %11
%33 = load i32, i32* %8
%34 = mul nsw i32 %32, %33
%35 = load i32, i32* %14
%36 = add nsw i32 %34, %35
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %31, i64 %37
%39 = load double, double* %38
%40 = load double*, double** %6
%41 = load i32, i32* %14
%42 = load i32, i32* %8
%43 = mul nsw i32 %41, %42
%44 = load i32, i32* %12
%45 = add nsw i32 %43, %44
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %40, i64 %46
%48 = load double, double* %47
%49 = load double, double* %13
%50 = call double @llvm.fmuladd.f64(double %39, double %48, double %49)
store double %50, double* %13
br label %51
51: ; preds = %30
%52 = load i32, i32* %14
%53 = add nsw i32 %52, 1
store i32 %53, i32* %14
br label %26
54: ; preds = %26
%55 = load i32, i32* %12
%56 = load i32, i32* %11
%57 = icmp eq i32 %55, %56
%58 = zext i1 %57 to i64
%59 = select i1 %57, double 1.000000e+00, double 0.000000e+00
store double %59, double* %15
%60 = load double, double* %13
%61 = load double, double* %15
%62 = fsub double %60, %61
%63 = call double @fabs(double %62)
%64 = load double, double* %9
%65 = fcmp ogt double %63, %64
br i1 %65, label %66, label %74
66: ; preds = %54
%67 = load i32, i32* @stderr
%68 = load double, double* %13
%69 = load double, double* %15
%70 = load i32, i32* %11
%71 = load i32, i32* %12
%72 = call i32 @fprintf(i32 %67, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str, i64 0, i64 0), double %68, double %69, i32 %70, i32 %71)
%73 = load i32*, i32** %10
store i32 22, i32* %73
br label %82
74: ; preds = %54
br label %75
75: ; preds = %74
%76 = load i32, i32* %12
%77 = add nsw i32 %76, 1
store i32 %77, i32* %12
br label %21
78: ; preds = %21
br label %79
79: ; preds = %78
%80 = load i32, i32* %11
%81 = add nsw i32 %80, 1
store i32 %81, i32* %11
br label %16
82: ; preds = %66, %16
ret void
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [51 x i8] c"lm_check_inverse found %g instead of %g at %i,%i \0A\00"
define void @lm_check_inverse(double* nocapture %0, double* nocapture %1, i32 %2, double %3, i32* nocapture writeonly %4) {
%6 = icmp sgt i32 %2, 0
br i1 %6, label %7, label %81
7: ; preds = %5
%8 = zext i32 %2 to i64
%9 = zext i32 %2 to i64
%10 = and i64 %9, 1
%11 = icmp eq i32 %2, 1
%12 = and i64 %9, 4294967294
%13 = icmp eq i64 %10, 0
br label %14
14: ; preds = %79, %7
%15 = phi i64 [ 0, %7 ], [ %80, %79 ]
%16 = mul nsw i64 %15, %8
br label %17
17: ; preds = %72, %14
%18 = phi i64 [ 0, %14 ], [ %73, %72 ]
%19 = phi i1 [ true, %14 ], [ %74, %72 ]
br i1 %11, label %20, label %40
20: ; preds = %40, %17
%21 = phi double [ undef, %17 ], [ %60, %40 ]
%22 = phi i64 [ 0, %17 ], [ %61, %40 ]
%23 = phi double [ 0.000000e+00, %17 ], [ %60, %40 ]
br i1 %13, label %33, label %24
24: ; preds = %20
%25 = add nsw i64 %22, %16
%26 = getelementptr inbounds double, double* %1, i64 %25
%27 = load double, double* %26
%28 = mul nsw i64 %22, %8
%29 = add nuw nsw i64 %28, %18
%30 = getelementptr inbounds double, double* %0, i64 %29
%31 = load double, double* %30
%32 = tail call double @llvm.fmuladd.f64(double %27, double %31, double %23)
br label %33
33: ; preds = %20, %24
%34 = phi double [ %21, %20 ], [ %32, %24 ]
%35 = icmp eq i64 %18, %15
%36 = select i1 %35, double 1.000000e+00, double 0.000000e+00
%37 = fsub double %34, %36
%38 = tail call double @llvm.fabs.f64(double %37)
%39 = fcmp ogt double %38, %3
br i1 %39, label %64, label %72
40: ; preds = %17, %40
%41 = phi i64 [ %61, %40 ], [ 0, %17 ]
%42 = phi double [ %60, %40 ], [ 0.000000e+00, %17 ]
%43 = phi i64 [ %62, %40 ], [ 0, %17 ]
%44 = add nsw i64 %41, %16
%45 = getelementptr inbounds double, double* %1, i64 %44
%46 = load double, double* %45
%47 = mul nsw i64 %41, %8
%48 = add nuw nsw i64 %47, %18
%49 = getelementptr inbounds double, double* %0, i64 %48
%50 = load double, double* %49
%51 = tail call double @llvm.fmuladd.f64(double %46, double %50, double %42)
%52 = or i64 %41, 1
%53 = add nsw i64 %52, %16
%54 = getelementptr inbounds double, double* %1, i64 %53
%55 = load double, double* %54
%56 = mul nsw i64 %52, %8
%57 = add nuw nsw i64 %56, %18
%58 = getelementptr inbounds double, double* %0, i64 %57
%59 = load double, double* %58
%60 = tail call double @llvm.fmuladd.f64(double %55, double %59, double %51)
%61 = add nuw nsw i64 %41, 2
%62 = add i64 %43, 2
%63 = icmp eq i64 %62, %12
br i1 %63, label %20, label %40
64: ; preds = %33
%65 = trunc i64 %18 to i32
%66 = load i32, i32* @stderr
%67 = trunc i64 %15 to i32
%68 = tail call i32 @fprintf(i32 %66, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str, i64 0, i64 0), double %34, double %36, i32 %67, i32 %65)
store i32 22, i32* %4
%69 = add nuw nsw i64 %15, 1
%70 = icmp eq i64 %69, %9
%71 = select i1 %19, i1 true, i1 %70
br i1 %71, label %81, label %79
72: ; preds = %33
%73 = add nuw nsw i64 %18, 1
%74 = icmp ult i64 %73, %8
%75 = icmp eq i64 %73, %9
br i1 %75, label %76, label %17
76: ; preds = %72
%77 = add nuw nsw i64 %15, 1
%78 = icmp eq i64 %77, %9
br i1 %78, label %81, label %79
79: ; preds = %76, %64
%80 = phi i64 [ %77, %76 ], [ %69, %64 ]
br label %14
81: ; preds = %64, %76, %5
ret void
}