#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 */
int /*<<< orphan*/ memset(double *, int /*<<< orphan*/, int const);
void bmgs_zero(double *a, const int n[3], const int c[3], const int s[3]) {
  a += c[2] + (c[1] + c[0] * n[1]) * n[2];
  for (int i0 = 0; i0 < s[0]; i0++) {
    for (int i1 = 0; i1 < s[1]; i1++) {
      memset(a, 0, s[2] * sizeof(double));
      a += n[2];
    }
    a += n[2] * (n[1] - s[1]);
  }
}


x86 assembly (clang, -O0):
<bmgs_zero>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%ecx
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x18(%rbp),%rdx
mov    (%rdx),%edx
mov    -0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
add    %edx,%eax
mov    -0x10(%rbp),%rdx
imul   0x8(%rdx),%eax
add    %eax,%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    1207 <bmgs_zero+0xf7>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    11cf <bmgs_zero+0xbf>
mov    -0x8(%rbp),%rdi
mov    -0x20(%rbp),%rax
movslq 0x8(%rax),%rax
shl    $0x3,%rax
mov    %eax,%edx
xor    %esi,%esi
call   1030 <memset@plt>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    117f <bmgs_zero+0x6f>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%ecx
mov    -0x10(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x20(%rbp),%rdx
sub    0x4(%rdx),%eax
imul   %eax,%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1169 <bmgs_zero+0x59>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bmgs_zero>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %rcx,%r15
mov    (%rcx),%ecx
test   %ecx,%ecx
jle    119d <bmgs_zero+0x8d>
mov    %rsi,%r14
mov    0x4(%rsi),%esi
mov    (%rdx),%ebp
imul   %esi,%ebp
add    0x4(%rdx),%ebp
mov    0x8(%r14),%eax
imul   %eax,%ebp
movslq 0x8(%rdx),%rdx
movslq %ebp,%rbp
add    %rdx,%rbp
lea    (%rdi,%rbp,8),%rbx
mov    0x4(%r15),%edx
xor    %r12d,%r12d
jmp    1167 <bmgs_zero+0x57>
nopl   (%rax)
mov    %esi,%edi
sub    %edx,%edi
imul   %eax,%edi
movslq %edi,%rdi
lea    (%rbx,%rdi,8),%rbx
add    $0x1,%r12d
cmp    %ecx,%r12d
jge    119d <bmgs_zero+0x8d>
test   %edx,%edx
jle    1150 <bmgs_zero+0x40>
xor    %ebp,%ebp
nopl   (%rax)
mov    0x8(%r15),%edx
shl    $0x3,%edx
mov    %rbx,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
movslq 0x8(%r14),%rax
lea    (%rbx,%rax,8),%rbx
add    $0x1,%ebp
mov    0x4(%r15),%edx
cmp    %edx,%ebp
jl     1170 <bmgs_zero+0x60>
mov    0x4(%r14),%esi
mov    (%r15),%ecx
jmp    1150 <bmgs_zero+0x40>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @bmgs_zero(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
store double* %0, double** %5
store i32* %1, i32** %6
store i32* %2, i32** %7
store i32* %3, i32** %8
%11 = load i32*, i32** %7
%12 = getelementptr inbounds i32, i32* %11, i64 2
%13 = load i32, i32* %12
%14 = load i32*, i32** %7
%15 = getelementptr inbounds i32, i32* %14, i64 1
%16 = load i32, i32* %15
%17 = load i32*, i32** %7
%18 = getelementptr inbounds i32, i32* %17, i64 0
%19 = load i32, i32* %18
%20 = load i32*, i32** %6
%21 = getelementptr inbounds i32, i32* %20, i64 1
%22 = load i32, i32* %21
%23 = mul nsw i32 %19, %22
%24 = add nsw i32 %16, %23
%25 = load i32*, i32** %6
%26 = getelementptr inbounds i32, i32* %25, i64 2
%27 = load i32, i32* %26
%28 = mul nsw i32 %24, %27
%29 = add nsw i32 %13, %28
%30 = load double*, double** %5
%31 = sext i32 %29 to i64
%32 = getelementptr inbounds double, double* %30, i64 %31
store double* %32, double** %5
store i32 0, i32* %9
br label %33
33: ; preds = %79, %4
%34 = load i32, i32* %9
%35 = load i32*, i32** %8
%36 = getelementptr inbounds i32, i32* %35, i64 0
%37 = load i32, i32* %36
%38 = icmp slt i32 %34, %37
br i1 %38, label %39, label %82
39: ; preds = %33
store i32 0, i32* %10
br label %40
40: ; preds = %61, %39
%41 = load i32, i32* %10
%42 = load i32*, i32** %8
%43 = getelementptr inbounds i32, i32* %42, i64 1
%44 = load i32, i32* %43
%45 = icmp slt i32 %41, %44
br i1 %45, label %46, label %64
46: ; preds = %40
%47 = load double*, double** %5
%48 = load i32*, i32** %8
%49 = getelementptr inbounds i32, i32* %48, i64 2
%50 = load i32, i32* %49
%51 = sext i32 %50 to i64
%52 = mul i64 %51, 8
%53 = trunc i64 %52 to i32
%54 = call i32 @memset(double* %47, i32 0, i32 %53)
%55 = load i32*, i32** %6
%56 = getelementptr inbounds i32, i32* %55, i64 2
%57 = load i32, i32* %56
%58 = load double*, double** %5
%59 = sext i32 %57 to i64
%60 = getelementptr inbounds double, double* %58, i64 %59
store double* %60, double** %5
br label %61
61: ; preds = %46
%62 = load i32, i32* %10
%63 = add nsw i32 %62, 1
store i32 %63, i32* %10
br label %40
64: ; preds = %40
%65 = load i32*, i32** %6
%66 = getelementptr inbounds i32, i32* %65, i64 2
%67 = load i32, i32* %66
%68 = load i32*, i32** %6
%69 = getelementptr inbounds i32, i32* %68, i64 1
%70 = load i32, i32* %69
%71 = load i32*, i32** %8
%72 = getelementptr inbounds i32, i32* %71, i64 1
%73 = load i32, i32* %72
%74 = sub nsw i32 %70, %73
%75 = mul nsw i32 %67, %74
%76 = load double*, double** %5
%77 = sext i32 %75 to i64
%78 = getelementptr inbounds double, double* %76, i64 %77
store double* %78, double** %5
br label %79
79: ; preds = %64
%80 = load i32, i32* %9
%81 = add nsw i32 %80, 1
store i32 %81, i32* %9
br label %33
82: ; preds = %33
ret void
}

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