#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 NAN;
void clear_dmatrix_nan(double **m, int nrl, int nrh, int ncl, int nch) {
  int i, j;
  for (i = nrl; i <= nrh; i++) {
    for (j = ncl; j <= nch; j++)
      m[i][j] = NAN;
  }
}


x86 assembly (clang, -O0):
<clear_dmatrix_nan>:
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    -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     117a <clear_dmatrix_nan+0x7a>
mov    -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1167 <clear_dmatrix_nan+0x67>
mov    0x2eb8(%rip),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112d <clear_dmatrix_nan+0x2d>
jmp    116c <clear_dmatrix_nan+0x6c>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    111b <clear_dmatrix_nan+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<clear_dmatrix_nan>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %edx,%esi
jle    1119 <clear_dmatrix_nan+0x19>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movslq %ecx,%r11
movslq %esi,%r15
add    $0x1,%edx
mov    %r8d,%r9d
sub    %r11d,%r9d
lea    (%r11,%r9,1),%rax
add    $0x1,%rax
mov    %rax,-0x8(%rsp)
lea    0x1(%r9),%rax
mov    %rax,-0x10(%rsp)
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%rsi
mov    %rsi,-0x20(%rsp)
shr    $0x2,%rsi
add    $0x1,%rsi
lea    0x1(%r8),%r12d
mov    %rax,-0x18(%rsp)
add    %r11,%rax
mov    %rax,-0x30(%rsp)
mov    %rsi,%rax
mov    %rsi,-0x28(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x38(%rsp)
lea    0x6(%r11),%rax
mov    %rax,-0x40(%rsp)
mov    0x2e73(%rip),%r13
jmp    119d <clear_dmatrix_nan+0x9d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r15
cmp    %r15d,%edx
je     110e <clear_dmatrix_nan+0xe>
cmp    %ecx,%r8d
jl     1190 <clear_dmatrix_nan+0x90>
mov    (%rdi,%r15,8),%rbp
mov    %r11,%rax
cmp    $0x3,%r9d
jb     1260 <clear_dmatrix_nan+0x160>
lea    0x0(,%r11,8),%rax
add    %rbp,%rax
lea    0x8(%r13),%rsi
cmp    %rsi,%rax
jae    11e3 <clear_dmatrix_nan+0xe3>
mov    -0x8(%rsp),%rax
lea    0x0(,%rax,8),%rsi
add    %rbp,%rsi
mov    %r11,%rax
cmp    %r13,%rsi
ja     1260 <clear_dmatrix_nan+0x160>
cmpq   $0x0,-0x20(%rsp)
je     12de <clear_dmatrix_nan+0x1de>
movsd  0x0(%r13),%xmm0
movlhps %xmm0,%xmm0
mov    -0x40(%rsp),%rax
lea    0x0(,%rax,8),%rsi
add    %rbp,%rsi
mov    -0x38(%rsp),%rbx
xor    %eax,%eax
nop
movups %xmm0,-0x30(%rsi,%rax,8)
movups %xmm0,-0x20(%rsi,%rax,8)
movups %xmm0,-0x10(%rsi,%rax,8)
movups %xmm0,(%rsi,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rbx
jne    1210 <clear_dmatrix_nan+0x110>
testb  $0x1,-0x28(%rsp)
je     124a <clear_dmatrix_nan+0x14a>
add    %r11,%rax
movsd  0x0(%r13),%xmm0
movlhps %xmm0,%xmm0
movups %xmm0,0x0(%rbp,%rax,8)
movups %xmm0,0x10(%rbp,%rax,8)
mov    -0x30(%rsp),%rax
mov    -0x18(%rsp),%rsi
cmp    %rsi,-0x10(%rsp)
je     1190 <clear_dmatrix_nan+0x90>
nop
mov    %r12d,%r10d
sub    %eax,%r10d
mov    %r8d,%r14d
sub    %eax,%r14d
and    $0x3,%r10d
je     1296 <clear_dmatrix_nan+0x196>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%r13),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x1,%rax
add    $0xffffffff,%r10d
jne    1280 <clear_dmatrix_nan+0x180>
cmp    $0x3,%r14d
jb     1190 <clear_dmatrix_nan+0x90>
movsd  0x0(%r13),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
movsd  0x0(%r13),%xmm0
movsd  %xmm0,0x8(%rbp,%rax,8)
movsd  0x0(%r13),%xmm0
movsd  %xmm0,0x10(%rbp,%rax,8)
movsd  0x0(%r13),%xmm0
movsd  %xmm0,0x18(%rbp,%rax,8)
add    $0x4,%rax
cmp    %eax,%r12d
jne    12a0 <clear_dmatrix_nan+0x1a0>
jmp    1190 <clear_dmatrix_nan+0x90>
xor    %eax,%eax
testb  $0x1,-0x28(%rsp)
jne    1234 <clear_dmatrix_nan+0x134>
jmp    124a <clear_dmatrix_nan+0x14a>

LLVM IR (-O0):
@NAN = global double 0.000000e+00
define void @clear_dmatrix_nan(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 i32
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 i32, i32* %7
store i32 %13, i32* %11
br label %14
14: ; preds = %38, %5
%15 = load i32, i32* %11
%16 = load i32, i32* %8
%17 = icmp sle i32 %15, %16
br i1 %17, label %18, label %41
18: ; preds = %14
%19 = load i32, i32* %9
store i32 %19, i32* %12
br label %20
20: ; preds = %34, %18
%21 = load i32, i32* %12
%22 = load i32, i32* %10
%23 = icmp sle i32 %21, %22
br i1 %23, label %24, label %37
24: ; preds = %20
%25 = load double, double* @NAN
%26 = load double**, double*** %6
%27 = load i32, i32* %11
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double*, double** %26, i64 %28
%30 = load double*, double** %29
%31 = load i32, i32* %12
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
store double %25, double* %33
br label %34
34: ; preds = %24
%35 = load i32, i32* %12
%36 = add nsw i32 %35, 1
store i32 %36, i32* %12
br label %20
37: ; preds = %20
br label %38
38: ; preds = %37
%39 = load i32, i32* %11
%40 = add nsw i32 %39, 1
store i32 %40, i32* %11
br label %14
41: ; preds = %14
ret void
}

LLVM IR (-O2):
@NAN = global double 0.000000e+00
define void @clear_dmatrix_nan(double** nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %1, %2
br i1 %6, label %124, label %7
7: ; preds = %5
%8 = icmp slt i32 %4, %3
%9 = sext i32 %3 to i64
%10 = add i32 %4, 1
%11 = sext i32 %1 to i64
%12 = add i32 %2, 1
%13 = sub i32 %4, %3
%14 = zext i32 %13 to i64
%15 = add nsw i64 %9, %14
%16 = add nsw i64 %15, 1
%17 = zext i32 %13 to i64
%18 = add nuw nsw i64 %17, 1
%19 = and i64 %18, 8589934588
%20 = add nsw i64 %19, -4
%21 = lshr exact i64 %20, 2
%22 = add nuw nsw i64 %21, 1
%23 = add i32 %4, 1
%24 = icmp ult i32 %13, 3
%25 = and i64 %18, 8589934588
%26 = add nsw i64 %25, %9
%27 = and i64 %22, 1
%28 = icmp eq i64 %20, 0
%29 = and i64 %22, 9223372036854775806
%30 = icmp eq i64 %27, 0
%31 = icmp eq i64 %18, %25
br label %32
32: ; preds = %7, %120
%33 = phi i64 [ %11, %7 ], [ %121, %120 ]
br i1 %8, label %120, label %34
34: ; preds = %32
%35 = getelementptr inbounds double*, double** %0, i64 %33
%36 = load double*, double** %35
br i1 %24, label %86, label %37
37: ; preds = %34
%38 = getelementptr double, double* %36, i64 %9
%39 = getelementptr double, double* %36, i64 %16
%40 = icmp ult double* %38, getelementptr inbounds (double, double* @NAN, i64 1)
%41 = icmp ugt double* %39, @NAN
%42 = and i1 %40, %41
br i1 %42, label %86, label %43
43: ; preds = %37
br i1 %28, label %72, label %44
44: ; preds = %43
%45 = load double, double* @NAN
%46 = insertelement <2 x double> poison, double %45, i64 0
%47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
%48 = insertelement <2 x double> poison, double %45, i64 0
%49 = shufflevector <2 x double> %48, <2 x double> poison, <2 x i32> zeroinitializer
%50 = load double, double* @NAN
%51 = insertelement <2 x double> poison, double %50, i64 0
%52 = shufflevector <2 x double> %51, <2 x double> poison, <2 x i32> zeroinitializer
%53 = insertelement <2 x double> poison, double %50, i64 0
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> zeroinitializer
br label %55
55: ; preds = %55, %44
%56 = phi i64 [ 0, %44 ], [ %69, %55 ]
%57 = phi i64 [ 0, %44 ], [ %70, %55 ]
%58 = add i64 %56, %9
%59 = getelementptr inbounds double, double* %36, i64 %58
%60 = bitcast double* %59 to <2 x double>*
store <2 x double> %47, <2 x double>* %60
%61 = getelementptr inbounds double, double* %59, i64 2
%62 = bitcast double* %61 to <2 x double>*
store <2 x double> %49, <2 x double>* %62
%63 = or i64 %56, 4
%64 = add i64 %63, %9
%65 = getelementptr inbounds double, double* %36, i64 %64
%66 = bitcast double* %65 to <2 x double>*
store <2 x double> %52, <2 x double>* %66
%67 = getelementptr inbounds double, double* %65, i64 2
%68 = bitcast double* %67 to <2 x double>*
store <2 x double> %54, <2 x double>* %68
%69 = add nuw i64 %56, 8
%70 = add i64 %57, 2
%71 = icmp eq i64 %70, %29
br i1 %71, label %72, label %55
72: ; preds = %55, %43
%73 = phi i64 [ 0, %43 ], [ %69, %55 ]
br i1 %30, label %85, label %74
74: ; preds = %72
%75 = add i64 %73, %9
%76 = load double, double* @NAN
%77 = insertelement <2 x double> poison, double %76, i64 0
%78 = shufflevector <2 x double> %77, <2 x double> poison, <2 x i32> zeroinitializer
%79 = insertelement <2 x double> poison, double %76, i64 0
%80 = shufflevector <2 x double> %79, <2 x double> poison, <2 x i32> zeroinitializer
%81 = getelementptr inbounds double, double* %36, i64 %75
%82 = bitcast double* %81 to <2 x double>*
store <2 x double> %78, <2 x double>* %82
%83 = getelementptr inbounds double, double* %81, i64 2
%84 = bitcast double* %83 to <2 x double>*
store <2 x double> %80, <2 x double>* %84
br label %85
85: ; preds = %72, %74
br i1 %31, label %120, label %86
86: ; preds = %37, %34, %85
%87 = phi i64 [ %9, %37 ], [ %9, %34 ], [ %26, %85 ]
%88 = trunc i64 %87 to i32
%89 = sub i32 %23, %88
%90 = sub i32 %4, %88
%91 = and i32 %89, 3
%92 = icmp eq i32 %91, 0
br i1 %92, label %101, label %93
93: ; preds = %86, %93
%94 = phi i64 [ %98, %93 ], [ %87, %86 ]
%95 = phi i32 [ %99, %93 ], [ 0, %86 ]
%96 = load double, double* @NAN
%97 = getelementptr inbounds double, double* %36, i64 %94
store double %96, double* %97
%98 = add nsw i64 %94, 1
%99 = add i32 %95, 1
%100 = icmp eq i32 %99, %91
br i1 %100, label %101, label %93
101: ; preds = %93, %86
%102 = phi i64 [ %87, %86 ], [ %98, %93 ]
%103 = icmp ult i32 %90, 3
br i1 %103, label %120, label %104
104: ; preds = %101, %104
%105 = phi i64 [ %117, %104 ], [ %102, %101 ]
%106 = load double, double* @NAN
%107 = getelementptr inbounds double, double* %36, i64 %105
store double %106, double* %107
%108 = add nsw i64 %105, 1
%109 = load double, double* @NAN
%110 = getelementptr inbounds double, double* %36, i64 %108
store double %109, double* %110
%111 = add nsw i64 %105, 2
%112 = load double, double* @NAN
%113 = getelementptr inbounds double, double* %36, i64 %111
store double %112, double* %113
%114 = add nsw i64 %105, 3
%115 = load double, double* @NAN
%116 = getelementptr inbounds double, double* %36, i64 %114
store double %115, double* %116
%117 = add nsw i64 %105, 4
%118 = trunc i64 %117 to i32
%119 = icmp eq i32 %10, %118
br i1 %119, label %120, label %104
120: ; preds = %101, %104, %85, %32
%121 = add nsw i64 %33, 1
%122 = trunc i64 %121 to i32
%123 = icmp eq i32 %12, %122
br i1 %123, label %124, label %32
124: ; preds = %120, %5
ret void
}