#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 */
float *c_natgrids(int, float *, float *, float *, int, int, float *, float *,
                  int *);
int /*<<< orphan*/ free(float *);
void natgrids_(int *n, float *x, float *y, float *z, int *nxg, int *nyg,
               float *xg, float *yg, float *zg, int *ier) {
  float *zar;
  int nn, mm;
  zar = c_natgrids(*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):
<natgrids_>:
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_natgrids@plt>
mov    %rax,-0x38(%rbp)
mov    0x28(%rbp),%rax
cmpl   $0x0,(%rax)
je     11a5 <natgrids_+0x85>
jmp    1235 <natgrids_+0x115>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jge    122c <natgrids_+0x10c>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x30(%rbp),%rcx
cmp    (%rcx),%eax
jge    1219 <natgrids_+0xf9>
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x30(%rbp),%rdx
imul   (%rdx),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x28(%rbp),%rdx
imul   (%rdx),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11c2 <natgrids_+0xa2>
jmp    121e <natgrids_+0xfe>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11ac <natgrids_+0x8c>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<natgrids_>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9,%r14
mov    %r8,%r15
mov    0x70(%rsp),%r12
movaps 0x60(%rsp),%xmm0
mov    0x78(%rsp),%rbx
mov    (%rdi),%edi
mov    (%r8),%r8d
mov    (%r9),%r9d
mov    %rbx,0x10(%rsp)
movups %xmm0,(%rsp)
call   1040 <c_natgrids@plt>
cmpl   $0x0,(%rbx)
je     116d <natgrids_+0x4d>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rax,0x18(%rsp)
mov    (%r15),%ecx
test   %ecx,%ecx
jle    128b <natgrids_+0x16b>
mov    (%r14),%r9d
movslq %r9d,%r10
mov    %r9,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,0x20(%rsp)
mov    %r10d,%esi
and    $0x3,%esi
mov    %r10d,%r11d
and    $0xfffffffc,%r11d
mov    0x18(%rsp),%r13
mov    %r13,%rdx
add    $0xc,%rdx
shl    $0x2,%r10
mov    %rcx,%r14
shl    $0x4,%r14
lea    0x0(,%rcx,4),%rdi
lea    (%rdi,%rdi,2),%r8
xor    %r15d,%r15d
jmp    11e7 <natgrids_+0xc7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r15
add    %r10,%rdx
add    $0x4,%r12
add    %r10,%r13
cmp    %rcx,%r15
je     128b <natgrids_+0x16b>
test   %r9d,%r9d
jle    11d0 <natgrids_+0xb0>
cmpq   $0x3,0x20(%rsp)
jae    1200 <natgrids_+0xe0>
mov    %r12,%rax
xor    %ebp,%ebp
jmp    124b <natgrids_+0x12b>
nopl   0x0(%rax,%rax,1)
mov    %r12,%rax
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
movss  -0xc(%rdx,%rbp,4),%xmm0
movss  %xmm0,(%r12)
movss  -0x8(%rdx,%rbp,4),%xmm0
movss  %xmm0,(%r12,%rcx,4)
movss  -0x4(%rdx,%rbp,4),%xmm0
movss  %xmm0,(%r12,%rcx,8)
movss  (%rdx,%rbp,4),%xmm0
movss  %xmm0,(%r12,%r8,1)
add    $0x4,%rbp
add    %r14,%r12
cmp    %rbp,%r11
jne    1210 <natgrids_+0xf0>
test   %rsi,%rsi
mov    %rax,%r12
je     11d0 <natgrids_+0xb0>
mov    %rcx,%rax
imul   %rbp,%rax
lea    (%r12,%rax,4),%rbx
lea    0x0(,%rbp,4),%rbp
add    %r13,%rbp
xor    %eax,%eax
nop
movss  0x0(%rbp,%rax,4),%xmm0
movss  %xmm0,(%rbx)
add    $0x1,%rax
add    %rdi,%rbx
cmp    %rax,%rsi
jne    1270 <natgrids_+0x150>
jmp    11d0 <natgrids_+0xb0>
mov    0x18(%rsp),%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @natgrids_(i32* %0, float* %1, float* %2, float* %3, i32* %4, i32* %5, float* %6, float* %7, float* %8, i32* %9) {
%11 = alloca i32*
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i32*
%16 = alloca i32*
%17 = alloca float*
%18 = alloca float*
%19 = alloca float*
%20 = alloca i32*
%21 = alloca float*
%22 = alloca i32
%23 = alloca i32
store i32* %0, i32** %11
store float* %1, float** %12
store float* %2, float** %13
store float* %3, float** %14
store i32* %4, i32** %15
store i32* %5, i32** %16
store float* %6, float** %17
store float* %7, float** %18
store float* %8, float** %19
store i32* %9, i32** %20
%24 = load i32*, i32** %11
%25 = load i32, i32* %24
%26 = load float*, float** %12
%27 = load float*, float** %13
%28 = load float*, float** %14
%29 = load i32*, i32** %15
%30 = load i32, i32* %29
%31 = load i32*, i32** %16
%32 = load i32, i32* %31
%33 = load float*, float** %17
%34 = load float*, float** %18
%35 = load i32*, i32** %20
%36 = call float* @c_natgrids(i32 %25, float* %26, float* %27, float* %28, i32 %30, i32 %32, float* %33, float* %34, i32* %35)
store float* %36, float** %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 float*, float** %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 float, float* %54, i64 %61
%63 = load float, float* %62
%64 = load float*, float** %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 float, float* %64, i64 %69
%71 = load i32, i32* %23
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds float, float* %70, i64 %72
store float %63, float* %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 float*, float** %21
%83 = call i32 @free(float* %82)
br label %84
84: ; preds = %81, %40
ret void
}

LLVM IR (-O2):
define void @natgrids_(i32* nocapture %0, float* %1, float* %2, float* %3, i32* nocapture %4, i32* nocapture %5, float* %6, float* %7, float* nocapture writeonly %8, i32* %9) {
%11 = load i32, i32* %0
%12 = load i32, i32* %4
%13 = load i32, i32* %5
%14 = tail call float* @c_natgrids(i32 %11, float* %1, float* %2, float* %3, i32 %12, i32 %13, float* %6, float* %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 float, float* %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 float, float* %14, i64 %40
%42 = load float, float* %41
%43 = mul nsw i64 %38, %23
%44 = getelementptr inbounds float, float* %36, i64 %43
store float %42, float* %44
%45 = or i64 %38, 1
%46 = add nsw i64 %35, %45
%47 = getelementptr inbounds float, float* %14, i64 %46
%48 = load float, float* %47
%49 = mul nsw i64 %45, %23
%50 = getelementptr inbounds float, float* %36, i64 %49
store float %48, float* %50
%51 = or i64 %38, 2
%52 = add nsw i64 %35, %51
%53 = getelementptr inbounds float, float* %14, i64 %52
%54 = load float, float* %53
%55 = mul nsw i64 %51, %23
%56 = getelementptr inbounds float, float* %36, i64 %55
store float %54, float* %56
%57 = or i64 %38, 3
%58 = add nsw i64 %35, %57
%59 = getelementptr inbounds float, float* %14, i64 %58
%60 = load float, float* %59
%61 = mul nsw i64 %57, %23
%62 = getelementptr inbounds float, float* %36, i64 %61
store float %60, float* %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 float, float* %14, i64 %71
%73 = load float, float* %72
%74 = mul nsw i64 %69, %23
%75 = getelementptr inbounds float, float* %36, i64 %74
store float %73, float* %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(float* %14)
br label %84
84: ; preds = %10, %82
ret void
}