#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 *c_natgridd(int, double *, double *, double *, int, int, double *,
                   double *, int *);
int /*<<< orphan*/ free(double *);
void natgridd_(int *n, double *x, double *y, double *z, int *nxg, int *nyg,
               double *xg, double *yg, double *zg, int *ier) {
  double *zar;
  int nn, mm;
  zar = c_natgridd(*n, x, y, z, *nxg, *nyg, xg, yg, ier);
  if (*ier)
    return;
  for (mm = 0; mm < *nxg; mm++) {
    for (nn = 0; nn < *nyg; nn++) {
      *(zg + nn * (*nxg) + mm) = zar[mm * (*nyg) + nn];
    }
  }
  free(zar);
  return;
}


x86 assembly (clang, -O0):
<natgridd_>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%edi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    (%rax),%r8d
mov    -0x30(%rbp),%rax
mov    (%rax),%r9d
mov    0x10(%rbp),%r11
mov    0x18(%rbp),%r10
mov    0x28(%rbp),%rax
mov    %r11,(%rsp)
mov    %r10,0x8(%rsp)
mov    %rax,0x10(%rsp)
call   1040 <c_natgridd@plt>
mov    %rax,-0x38(%rbp)
mov    0x28(%rbp),%rax
cmpl   $0x0,(%rax)
je     11a5 <natgridd_+0x85>
jmp    1235 <natgridd_+0x115>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jge    122c <natgridd_+0x10c>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x30(%rbp),%rcx
cmp    (%rcx),%eax
jge    1219 <natgridd_+0xf9>
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x30(%rbp),%rdx
imul   (%rdx),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x28(%rbp),%rdx
imul   (%rdx),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11c2 <natgridd_+0xa2>
jmp    121e <natgridd_+0xfe>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11ac <natgridd_+0x8c>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<natgridd_>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9,%r14
mov    %r8,%r15
mov    0x80(%rsp),%r12
movaps 0x70(%rsp),%xmm0
mov    0x88(%rsp),%rbx
mov    (%rdi),%edi
mov    (%r8),%r8d
mov    (%r9),%r9d
mov    %rbx,0x10(%rsp)
movups %xmm0,(%rsp)
call   1040 <c_natgridd@plt>
cmpl   $0x0,(%rbx)
je     1173 <natgridd_+0x53>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rax,0x18(%rsp)
mov    (%r15),%r10d
test   %r10d,%r10d
jle    12ab <natgridd_+0x18b>
mov    (%r14),%eax
movslq %eax,%rcx
mov    %rax,0x30(%rsp)
add    $0xffffffffffffffff,%rax
mov    %rax,0x20(%rsp)
mov    %ecx,%esi
and    $0x3,%esi
mov    %ecx,%r11d
and    $0xfffffffc,%r11d
mov    0x18(%rsp),%r9
mov    %r9,%rdx
add    $0x18,%rdx
shl    $0x3,%rcx
mov    %rcx,0x28(%rsp)
mov    %r10,%r14
shl    $0x4,%r14
mov    %r10,%r15
shl    $0x5,%r15
lea    0x0(,%r10,8),%rdi
lea    (%rdi,%rdi,2),%rcx
xor    %r13d,%r13d
jmp    11fc <natgridd_+0xdc>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r13
mov    0x28(%rsp),%rax
add    %rax,%rdx
add    $0x8,%r12
add    %rax,%r9
cmp    %r10,%r13
je     12ab <natgridd_+0x18b>
cmpl   $0x0,0x30(%rsp)
jle    11e0 <natgridd_+0xc0>
cmpq   $0x3,0x20(%rsp)
jae    1220 <natgridd_+0x100>
mov    %r12,%rax
xor    %r8d,%r8d
jmp    126f <natgridd_+0x14f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r12,%rax
xor    %r8d,%r8d
cs nopw 0x0(%rax,%rax,1)
movsd  -0x18(%rdx,%r8,8),%xmm0
movsd  %xmm0,(%r12)
movsd  -0x10(%rdx,%r8,8),%xmm0
movsd  %xmm0,(%r12,%r10,8)
movsd  -0x8(%rdx,%r8,8),%xmm0
movsd  %xmm0,(%r12,%r14,1)
movsd  (%rdx,%r8,8),%xmm0
movsd  %xmm0,(%r12,%rcx,1)
add    $0x4,%r8
add    %r15,%r12
cmp    %r8,%r11
jne    1230 <natgridd_+0x110>
test   %rsi,%rsi
mov    %rax,%r12
je     11e0 <natgridd_+0xc0>
mov    %r10,%rax
imul   %r8,%rax
lea    (%r12,%rax,8),%rax
lea    (%r9,%r8,8),%rbp
xor    %ebx,%ebx
nopl   0x0(%rax)
movsd  0x0(%rbp,%rbx,8),%xmm0
movsd  %xmm0,(%rax)
add    $0x1,%rbx
add    %rdi,%rax
cmp    %rbx,%rsi
jne    1290 <natgridd_+0x170>
jmp    11e0 <natgridd_+0xc0>
mov    0x18(%rsp),%rdi
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

LLVM IR (-O2):
define void @natgridd_(i32* nocapture %0, double* %1, double* %2, double* %3, i32* nocapture %4, i32* nocapture %5, double* %6, double* %7, double* nocapture writeonly %8, i32* %9) {
%11 = load i32, i32* %0
%12 = load i32, i32* %4
%13 = load i32, i32* %5
%14 = tail call double* @c_natgridd(i32 %11, double* %1, double* %2, double* %3, i32 %12, i32 %13, double* %6, double* %7, i32* %9)
%15 = load i32, i32* %9
%16 = icmp eq i32 %15, 0
br i1 %16, label %17, label %84
17: ; preds = %10
%18 = load i32, i32* %4
%19 = icmp sgt i32 %18, 0
br i1 %19, label %20, label %82
20: ; preds = %17
%21 = load i32, i32* %5
%22 = icmp sgt i32 %21, 0
%23 = zext i32 %18 to i64
%24 = sext i32 %21 to i64
%25 = zext i32 %18 to i64
%26 = zext i32 %21 to i64
%27 = add nsw i64 %26, -1
%28 = and i64 %26, 3
%29 = icmp ult i64 %27, 3
%30 = and i64 %26, 4294967292
%31 = icmp eq i64 %28, 0
br label %32
32: ; preds = %20, %79
%33 = phi i64 [ 0, %20 ], [ %80, %79 ]
br i1 %22, label %34, label %79
34: ; preds = %32
%35 = mul nsw i64 %33, %24
%36 = getelementptr inbounds double, double* %8, i64 %33
br i1 %29, label %66, label %37
37: ; preds = %34, %37
%38 = phi i64 [ %63, %37 ], [ 0, %34 ]
%39 = phi i64 [ %64, %37 ], [ 0, %34 ]
%40 = add nsw i64 %35, %38
%41 = getelementptr inbounds double, double* %14, i64 %40
%42 = load double, double* %41
%43 = mul nsw i64 %38, %23
%44 = getelementptr inbounds double, double* %36, i64 %43
store double %42, double* %44
%45 = or i64 %38, 1
%46 = add nsw i64 %35, %45
%47 = getelementptr inbounds double, double* %14, i64 %46
%48 = load double, double* %47
%49 = mul nsw i64 %45, %23
%50 = getelementptr inbounds double, double* %36, i64 %49
store double %48, double* %50
%51 = or i64 %38, 2
%52 = add nsw i64 %35, %51
%53 = getelementptr inbounds double, double* %14, i64 %52
%54 = load double, double* %53
%55 = mul nsw i64 %51, %23
%56 = getelementptr inbounds double, double* %36, i64 %55
store double %54, double* %56
%57 = or i64 %38, 3
%58 = add nsw i64 %35, %57
%59 = getelementptr inbounds double, double* %14, i64 %58
%60 = load double, double* %59
%61 = mul nsw i64 %57, %23
%62 = getelementptr inbounds double, double* %36, i64 %61
store double %60, double* %62
%63 = add nuw nsw i64 %38, 4
%64 = add i64 %39, 4
%65 = icmp eq i64 %64, %30
br i1 %65, label %66, label %37
66: ; preds = %37, %34
%67 = phi i64 [ 0, %34 ], [ %63, %37 ]
br i1 %31, label %79, label %68
68: ; preds = %66, %68
%69 = phi i64 [ %76, %68 ], [ %67, %66 ]
%70 = phi i64 [ %77, %68 ], [ 0, %66 ]
%71 = add nsw i64 %35, %69
%72 = getelementptr inbounds double, double* %14, i64 %71
%73 = load double, double* %72
%74 = mul nsw i64 %69, %23
%75 = getelementptr inbounds double, double* %36, i64 %74
store double %73, double* %75
%76 = add nuw nsw i64 %69, 1
%77 = add i64 %70, 1
%78 = icmp eq i64 %77, %28
br i1 %78, label %79, label %68
79: ; preds = %66, %68, %32
%80 = add nuw nsw i64 %33, 1
%81 = icmp eq i64 %80, %25
br i1 %81, label %82, label %32
82: ; preds = %79, %17
%83 = tail call i32 @free(double* %14)
br label %84
84: ; preds = %10, %82
ret void
}