#include <math.h>
#include <iostream>
#include <array>
#include <cmath>
#include <cstdint>
#include "cuda_runtime.h"
#include <stdlib.h>
#include <cuda_runtime_api.h>
#include <cuda.h>

using namespace std;

template<int E, int M, int T, int P, int B = (1 << (E - 1)) - 1>
static inline __device__ uint64_t compress(float* xyz)
{
    static_assert(E + M + P + T == 64, "Invalid number of bits");
    
    const float rpis = 1.0/static_cast<float>(M_PI);
    const double told = 1e-13;

    float ts = atan2(static_cast<float>(xyz[1]),static_cast<float>(xyz[0]));
    uint64_t nt = round(static_cast<float>((1 << (T-1)) - 1)*(abs(ts)*rpis));
    if(ts >= 0) nt = nt | (1 << (T-1));

    double rd = sqrt((double)xyz[0]*xyz[0] + (double)xyz[1]*xyz[1] + (double)xyz[2]*xyz[2]);
    float rs = (float)rd;
    float pd = (rd < told) ? 0 : static_cast<float>(acos(static_cast<double>(xyz[2])/rd));
    uint64_t np = round(static_cast<float>((1 << P)-1)*pd*rpis);

    uint32_t* ptr = reinterpret_cast<uint32_t*>(&rs);
    uint64_t fre = (*ptr & 0x7f800000) >> 23;
    uint64_t frm = (*ptr & 0x007fffff) >> 23-M ;
    
    return nt | (np << (T)) | (frm << (T + P)) | ((fre - 127 + B) << (T + P + M));
    
}

template<int E, int M, int T, int P, int B = (1 << (E - 1)) - 1>
static inline __device__ void decompress(const uint64_t& C, float& x, float& y, float& z)
{
    static_assert(E + M + P + T == 64, "Invalid number of bits");
    
    const float pi = static_cast<float>(M_PI);
    
    uint32_t tm = (1 << T) - 1;
    uint32_t pm = (1 << P) - 1;
    uint32_t mm = (1 << M) - 1;
    uint32_t em = (1 << E) - 1;
    
    uint32_t td = C & tm;
    uint32_t pd = (C >> T) & pm;
    uint32_t frm = (C >> (T+P)) & mm;
    uint32_t fre = (C >> (T+P+M)) & em;
    uint32_t ri = ((fre - B + 127) << 23) | (frm << (23-M));
    
    float r = *reinterpret_cast<float*>(&ri);    
    float t = pi*(-1.0 + 2.0*(td >> (T-1)))*(td & (tm >> 1))/static_cast<float>((1 << (T-1))-1);    
    float p = pi*pd/static_cast<float>((1 << P)-1);
    
    x = r*cosf(t)*sinf(p);
    y = r*sinf(t)*sinf(p);
    z = r*cosf(p);
    
}

	.file	"tmpxft_002cbde5_00000000-6_compress.cudafe1.cpp"
	.text
#APP
	.globl _ZSt21ios_base_library_initv

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB3679:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE3679:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB3702:
	.cfi_startproc
	endbr64
	pushq	%rax
	.cfi_def_cfa_offset 16
	leaq	_ZL15__fatDeviceText(%rip), %rdi
	call	__cudaRegisterFatBinary@PLT
	movq	%rax, %rdi
	movq	%rax, _ZL20__cudaFatCubinHandle(%rip)
	call	__cudaRegisterFatBinaryEnd@PLT
	leaq	_ZL26__cudaUnregisterBinaryUtilv(%rip), %rdi
	popq	%rdx
	.cfi_def_cfa_offset 8
	jmp	atexit@PLT
	.cfi_endproc
.LFE3702:
	.size	_ZL24__sti____cudaRegisterAllv, .-_ZL24__sti____cudaRegisterAllv
	.section	.init_array,"aw"
	.align 8
	.quad	_ZL24__sti____cudaRegisterAllv
	.section	.nvFatBinSegment,"aw"
	.align 8
	.type	_ZL15__fatDeviceText, @object
	.size	_ZL15__fatDeviceText, 24
_ZL15__fatDeviceText:
	.long	1180844977
	.long	1
	.quad	fatbinData
	.quad	0
	.local	_ZL20__cudaFatCubinHandle
	.comm	_ZL20__cudaFatCubinHandle,8,8
	.ident	"GCC: (Ubuntu 13.3.0-6ubuntu2~24.04) 13.3.0"
	.section	.note.GNU-stack,"",@progbits
	.section	.note.gnu.property,"a"
	.align 8
	.long	1f - 0f
	.long	4f - 1f
	.long	5
0:
	.string	"GNU"
1:
	.align 8
	.long	0xc0000002
	.long	3f - 2f
2:
	.long	0x3
3:
	.align 8
4:

#include <math.h>
#include <iostream>
#include <array>
#include <cmath>
#include <cstdint>
#include "hip/hip_runtime.h"
#include <stdlib.h>
#include <hip/hip_runtime_api.h>


using namespace std;

template<int E, int M, int T, int P, int B = (1 << (E - 1)) - 1>
static inline __device__ uint64_t compress(float* xyz)
{
    static_assert(E + M + P + T == 64, "Invalid number of bits");
    
    const float rpis = 1.0/static_cast<float>(M_PI);
    const double told = 1e-13;

    float ts = atan2(static_cast<float>(xyz[1]),static_cast<float>(xyz[0]));
    uint64_t nt = round(static_cast<float>((1 << (T-1)) - 1)*(abs(ts)*rpis));
    if(ts >= 0) nt = nt | (1 << (T-1));

    double rd = sqrt((double)xyz[0]*xyz[0] + (double)xyz[1]*xyz[1] + (double)xyz[2]*xyz[2]);
    float rs = (float)rd;
    float pd = (rd < told) ? 0 : static_cast<float>(acos(static_cast<double>(xyz[2])/rd));
    uint64_t np = round(static_cast<float>((1 << P)-1)*pd*rpis);

    uint32_t* ptr = reinterpret_cast<uint32_t*>(&rs);
    uint64_t fre = (*ptr & 0x7f800000) >> 23;
    uint64_t frm = (*ptr & 0x007fffff) >> 23-M ;
    
    return nt | (np << (T)) | (frm << (T + P)) | ((fre - 127 + B) << (T + P + M));
    
}

template<int E, int M, int T, int P, int B = (1 << (E - 1)) - 1>
static inline __device__ void decompress(const uint64_t& C, float& x, float& y, float& z)
{
    static_assert(E + M + P + T == 64, "Invalid number of bits");
    
    const float pi = static_cast<float>(M_PI);
    
    uint32_t tm = (1 << T) - 1;
    uint32_t pm = (1 << P) - 1;
    uint32_t mm = (1 << M) - 1;
    uint32_t em = (1 << E) - 1;
    
    uint32_t td = C & tm;
    uint32_t pd = (C >> T) & pm;
    uint32_t frm = (C >> (T+P)) & mm;
    uint32_t fre = (C >> (T+P+M)) & em;
    uint32_t ri = ((fre - B + 127) << 23) | (frm << (23-M));
    
    float r = *reinterpret_cast<float*>(&ri);    
    float t = pi*(-1.0 + 2.0*(td >> (T-1)))*(td & (tm >> 1))/static_cast<float>((1 << (T-1))-1);    
    float p = pi*pd/static_cast<float>((1 << P)-1);
    
    x = r*cosf(t)*sinf(p);
    y = r*sinf(t)*sinf(p);
    z = r*cosf(p);
    
}

	.text
	.file	"compress.hip"
                                        # Start of file scope inline assembly
	.globl	_ZSt21ios_base_library_initv

                                        # End of file scope inline assembly
	.type	__hip_cuid_,@object             # @__hip_cuid_
	.bss
	.globl	__hip_cuid_
__hip_cuid_:
	.byte	0                               # 0x0
	.size	__hip_cuid_, 1

	.ident	"AMD clang version 18.0.0git (https://github.com/RadeonOpenCompute/llvm-project roc-6.3.2 25012 e5bf7e55c91490b07c49d8960fa7983d864936c4)"
	.section	".note.GNU-stack","",@progbits
	.addrsig
	.addrsig_sym __hip_cuid_

	.text
	.p2alignl 7, 3214868480
	.fill 96, 4, 3214868480
	.type	__hip_cuid_,@object             ; @__hip_cuid_
	.section	.bss,"aw",@nobits
	.globl	__hip_cuid_
__hip_cuid_:
	.byte	0                               ; 0x0
	.size	__hip_cuid_, 1

	.ident	"AMD clang version 18.0.0git (https://github.com/RadeonOpenCompute/llvm-project roc-6.3.2 25012 e5bf7e55c91490b07c49d8960fa7983d864936c4)"
	.section	".note.GNU-stack","",@progbits
	.addrsig
	.addrsig_sym __hip_cuid_
	.amdgcn_target "amdgcn-amd-amdhsa--gfx1100"
	.amdgpu_metadata
---
amdhsa.kernels:  []
amdhsa.target:   amdgcn-amd-amdhsa--gfx1100
amdhsa.version:
  - 1
  - 2
...

	.end_amdgpu_metadata
