void clear_dmatrix(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] = 0.0;
  }
}


x86 assembly (clang, -O0):
<clear_dmatrix>:
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     1172 <clear_dmatrix+0x72>
mov    -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     115f <clear_dmatrix+0x5f>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112d <clear_dmatrix+0x2d>
jmp    1164 <clear_dmatrix+0x64>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    111b <clear_dmatrix+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<clear_dmatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %edx,%eax
sub    %esi,%eax
jl     11a7 <clear_dmatrix+0x97>
mov    %ecx,%r15d
mov    %rdi,%r14
mov    %eax,0x8(%rsp)
movslq %ecx,%rcx
mov    %r8d,%eax
sub    %ecx,%eax
lea    0x8(,%rax,8),%rbp
movslq %esi,%r12
add    $0x1,%edx
mov    %edx,0xc(%rsp)
mov    %edx,%r13d
sub    %esi,%r13d
and    $0x3,%r13d
mov    %r8d,0x4(%rsp)
je     1199 <clear_dmatrix+0x89>
mov    %rcx,%rbx
lea    0x0(,%rcx,8),%rax
mov    %rax,0x10(%rsp)
jmp    117a <clear_dmatrix+0x6a>
nop
add    $0x1,%r12
add    $0xffffffff,%r13d
je     119c <clear_dmatrix+0x8c>
cmp    %r15d,%r8d
jl     1170 <clear_dmatrix+0x60>
mov    (%r14,%r12,8),%rdi
add    0x10(%rsp),%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1030 <memset@plt>
mov    0x4(%rsp),%r8d
jmp    1170 <clear_dmatrix+0x60>
mov    %rcx,%rbx
cmpl   $0x3,0x8(%rsp)
mov    0xc(%rsp),%eax
jae    11b6 <clear_dmatrix+0xa6>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %eax,%r13d
shl    $0x3,%rbx
jmp    11c9 <clear_dmatrix+0xb9>
nop
add    $0x4,%r12
cmp    %r12d,%r13d
je     11a7 <clear_dmatrix+0x97>
cmp    %r15d,%r8d
jl     11c0 <clear_dmatrix+0xb0>
mov    (%r14,%r12,8),%rdi
add    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1030 <memset@plt>
mov    0x8(%r14,%r12,8),%rdi
add    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1030 <memset@plt>
mov    0x10(%r14,%r12,8),%rdi
add    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1030 <memset@plt>
mov    0x18(%r14,%r12,8),%rdi
add    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1030 <memset@plt>
mov    0x4(%rsp),%r8d
jmp    11c0 <clear_dmatrix+0xb0>

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

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