#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 M_PI;
int cos(double);
double log(scalar_t__);
scalar_t__ tan(double);
void latlon2tile(double lat, double lon, int zoom, unsigned int *x,
                 unsigned int *y) {
  double lat_rad = lat * M_PI / 180;
  unsigned long long n = 1LL << zoom;
  long long llx = n * ((lon + 180) / 360);
  long long lly = n * (1 - (log(tan(lat_rad) + 1 / cos(lat_rad)) / M_PI)) / 2;
  if (lat >= 85.0511) {
    lly = 0;
  }
  if (lat <= -85.0511) {
    lly = n - 1;
  }
  if (llx < 0) {
    llx = 0;
  }
  if (lly < 0) {
    lly = 0;
  }
  if (llx >= n) {
    llx = n - 1;
  }
  if (lly >= n) {
    lly = n - 1;
  }
  *x = llx;
  *y = lly;
}


x86 assembly (clang, -O0):
<latlon2tile>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    0x2e8f(%rip),%rax
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xeaf(%rip),%xmm4
divsd  %xmm4,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%cl
mov    $0x1,%eax
shl    %cl,%rax
mov    %rax,-0x38(%rbp)
movq   -0x38(%rbp),%xmm0
movaps 0xea3(%rip),%xmm2
punpckldq %xmm2,%xmm0
movapd 0xea7(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm3
unpckhpd %xmm0,%xmm0
addsd  %xmm3,%xmm0
movsd  -0x10(%rbp),%xmm3
addsd  %xmm4,%xmm3
movsd  0xe67(%rip),%xmm4
divsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x40(%rbp)
movq   -0x38(%rbp),%xmm0
punpckldq %xmm2,%xmm0
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1030 <tan@plt>
mov    %rax,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <cos@plt>
mov    -0x58(%rbp),%rdi
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
cltq
add    %rax,%rdi
call   1050 <log@plt>
movaps %xmm0,%xmm2
movsd  -0x50(%rbp),%xmm0
mov    0x2dc2(%rip),%rax
divsd  (%rax),%xmm2
movsd  0xdde(%rip),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdc6(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  0xda4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     126e <latlon2tile+0x13e>
movq   $0x0,-0x48(%rbp)
movsd  0xdb2(%rip),%xmm0
ucomisd -0x8(%rbp),%xmm0
jb     128d <latlon2tile+0x15d>
mov    -0x38(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x48(%rbp)
cmpq   $0x0,-0x40(%rbp)
jge    12a0 <latlon2tile+0x170>
movq   $0x0,-0x40(%rbp)
cmpq   $0x0,-0x48(%rbp)
jge    12b3 <latlon2tile+0x183>
movq   $0x0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x38(%rbp),%rax
jb     12cd <latlon2tile+0x19d>
mov    -0x38(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x38(%rbp),%rax
jb     12e7 <latlon2tile+0x1b7>
mov    -0x38(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x48(%rbp),%rax
mov    %eax,%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<latlon2tile>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rdx,0x18(%rsp)
mov    %rsi,%r15
mov    %edi,%ecx
movsd  %xmm0,0x10(%rsp)
mov    0x2e93(%rip),%rbp
movsd  0x0(%rbp),%xmm2
mulsd  %xmm0,%xmm2
movsd  0xe9a(%rip),%xmm4
divsd  %xmm4,%xmm2
movapd %xmm2,%xmm0
movsd  %xmm2,0x8(%rsp)
mov    $0x1,%r13d
shl    %cl,%r13
movq   %r13,%xmm2
punpckldq 0xea6(%rip),%xmm2
subpd  0xeae(%rip),%xmm2
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
addsd  %xmm4,%xmm1
divsd  0xe62(%rip),%xmm1
addsd  %xmm2,%xmm3
movapd %xmm3,0x20(%rsp)
mulsd  %xmm3,%xmm1
cvttsd2si %xmm1,%rbx
call   1030 <tan@plt>
mov    %rax,%r12
movsd  0x8(%rsp),%xmm0
call   1040 <cos@plt>
lea    0x1(%rax),%ecx
xor    %r14d,%r14d
cmp    $0x3,%ecx
cmovae %r14d,%eax
movslq %eax,%rdi
add    %r12,%rdi
call   1050 <log@plt>
divsd  0x0(%rbp),%xmm0
movsd  0xe1f(%rip),%xmm1
subsd  %xmm0,%xmm1
mulsd  0x20(%rsp),%xmm1
mulsd  0xe15(%rip),%xmm1
cvttsd2si %xmm1,%rax
movsd  0x10(%rsp),%xmm1
ucomisd 0xe0a(%rip),%xmm1
cmovae %r14,%rax
lea    -0x1(%r13),%rcx
movsd  0xe02(%rip),%xmm0
ucomisd %xmm1,%xmm0
cmovae %rcx,%rax
test   %rbx,%rbx
cmovs  %r14,%rbx
test   %rax,%rax
cmovs  %r14,%rax
cmp    %r13,%rbx
cmovae %ecx,%ebx
cmp    %r13,%rax
mov    %ebx,(%r15)
cmovae %ecx,%eax
mov    0x18(%rsp),%rcx
mov    %eax,(%rcx)
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define void @latlon2tile(double %0, double %1, i32 %2, i32* %3, i32* %4) {
%6 = alloca double
%7 = alloca double
%8 = alloca i32
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca double
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
store double %0, double* %6
store double %1, double* %7
store i32 %2, i32* %8
store i32* %3, i32** %9
store i32* %4, i32** %10
%15 = load double, double* %6
%16 = load double, double* @M_PI
%17 = fmul double %15, %16
%18 = fdiv double %17, 1.800000e+02
store double %18, double* %11
%19 = load i32, i32* %8
%20 = zext i32 %19 to i64
%21 = shl i64 1, %20
store i64 %21, i64* %12
%22 = load i64, i64* %12
%23 = uitofp i64 %22 to double
%24 = load double, double* %7
%25 = fadd double %24, 1.800000e+02
%26 = fdiv double %25, 3.600000e+02
%27 = fmul double %23, %26
%28 = fptosi double %27 to i64
store i64 %28, i64* %13
%29 = load i64, i64* %12
%30 = uitofp i64 %29 to double
%31 = load double, double* %11
%32 = call i64 @tan(double %31)
%33 = load double, double* %11
%34 = call i32 @cos(double %33)
%35 = sdiv i32 1, %34
%36 = sext i32 %35 to i64
%37 = add nsw i64 %32, %36
%38 = call double @log(i64 %37)
%39 = load double, double* @M_PI
%40 = fdiv double %38, %39
%41 = fsub double 1.000000e+00, %40
%42 = fmul double %30, %41
%43 = fdiv double %42, 2.000000e+00
%44 = fptosi double %43 to i64
store i64 %44, i64* %14
%45 = load double, double* %6
%46 = fcmp oge double %45, 8.505110e+01
br i1 %46, label %47, label %48
47: ; preds = %5
store i64 0, i64* %14
br label %48
48: ; preds = %47, %5
%49 = load double, double* %6
%50 = fcmp ole double %49, -8.505110e+01
br i1 %50, label %51, label %54
51: ; preds = %48
%52 = load i64, i64* %12
%53 = sub i64 %52, 1
store i64 %53, i64* %14
br label %54
54: ; preds = %51, %48
%55 = load i64, i64* %13
%56 = icmp slt i64 %55, 0
br i1 %56, label %57, label %58
57: ; preds = %54
store i64 0, i64* %13
br label %58
58: ; preds = %57, %54
%59 = load i64, i64* %14
%60 = icmp slt i64 %59, 0
br i1 %60, label %61, label %62
61: ; preds = %58
store i64 0, i64* %14
br label %62
62: ; preds = %61, %58
%63 = load i64, i64* %13
%64 = load i64, i64* %12
%65 = icmp uge i64 %63, %64
br i1 %65, label %66, label %69
66: ; preds = %62
%67 = load i64, i64* %12
%68 = sub i64 %67, 1
store i64 %68, i64* %13
br label %69
69: ; preds = %66, %62
%70 = load i64, i64* %14
%71 = load i64, i64* %12
%72 = icmp uge i64 %70, %71
br i1 %72, label %73, label %76
73: ; preds = %69
%74 = load i64, i64* %12
%75 = sub i64 %74, 1
store i64 %75, i64* %14
br label %76
76: ; preds = %73, %69
%77 = load i64, i64* %13
%78 = trunc i64 %77 to i32
%79 = load i32*, i32** %9
store i32 %78, i32* %79
%80 = load i64, i64* %14
%81 = trunc i64 %80 to i32
%82 = load i32*, i32** %10
store i32 %81, i32* %82
ret void
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define void @latlon2tile(double %0, double %1, i32 %2, i32* nocapture writeonly %3, i32* nocapture writeonly %4) {
%6 = load double, double* @M_PI
%7 = fmul double %6, %0
%8 = fdiv double %7, 1.800000e+02
%9 = zext i32 %2 to i64
%10 = shl nuw i64 1, %9
%11 = uitofp i64 %10 to double
%12 = fadd double %1, 1.800000e+02
%13 = fdiv double %12, 3.600000e+02
%14 = fmul double %13, %11
%15 = fptosi double %14 to i64
%16 = tail call i64 @tan(double %8)
%17 = tail call i32 @cos(double %8)
%18 = add i32 %17, 1
%19 = icmp ult i32 %18, 3
%20 = select i1 %19, i32 %17, i32 0
%21 = sext i32 %20 to i64
%22 = add nsw i64 %16, %21
%23 = tail call double @log(i64 %22)
%24 = load double, double* @M_PI
%25 = fdiv double %23, %24
%26 = fsub double 1.000000e+00, %25
%27 = fmul double %26, %11
%28 = fmul double %27, 5.000000e-01
%29 = fptosi double %28 to i64
%30 = fcmp ult double %0, 8.505110e+01
%31 = select i1 %30, i64 %29, i64 0
%32 = fcmp ugt double %0, -8.505110e+01
%33 = add i64 %10, -1
%34 = select i1 %32, i64 %31, i64 %33
%35 = icmp sgt i64 %15, 0
%36 = select i1 %35, i64 %15, i64 0
%37 = icmp sgt i64 %34, 0
%38 = select i1 %37, i64 %34, i64 0
%39 = icmp ult i64 %36, %10
%40 = select i1 %39, i64 %36, i64 %33
%41 = icmp ult i64 %38, %10
%42 = select i1 %41, i64 %38, i64 %33
%43 = trunc i64 %40 to i32
store i32 %43, i32* %3
%44 = trunc i64 %42 to i32
store i32 %44, i32* %4
ret void
}