#include "Tools.cuh"
#include "defines.cuh"
#include "function_defines.cuh"


__host__ inline void sampler_PureRandom	( Sampler *sampler );
__host__ inline void sampler_Regular		( Sampler *sampler );
__host__ inline void sampler_Jittered		( Sampler *sampler );
__host__ inline void sampler_Nrooks		( Sampler *sampler );
__host__ inline void sampler_MultiJittered( Sampler *sampler );
__host__ inline void sampler_Hammersley	( Sampler *sampler );

__host__ inline
void GenerateSample(Sampler *sampler){
	switch( sampler->type ){
	case SAMPLER_PURERANDOM :
		sampler_PureRandom(sampler);
		break;
	case SAMPLER_REGULAR :
		sampler_Regular(sampler);
		break;
	case SAMPLER_JITTERED :
		sampler_Jittered(sampler);
		break;
	case SAMPLER_NROOKS :
		sampler_Nrooks(sampler);
		break;
	case SAMPLER_MULTIJITTERED :
		sampler_MultiJittered(sampler);
		break;
	case SAMPLER_HAMMERSLEY :
		sampler_Hammersley(sampler);
		break;
	default:
		break;
	}
}

__device__ inline
Point2D getSampleUnitSquare(Sampler *sampler){
	return sampler->sample[ sampler->count++ % SAMPLE_POOL_SIZE ];
}

__host__ inline
void sampler_PureRandom		( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Regular		( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Jittered		( Sampler *sampler  ){

	for(int j = 0 ; j < SAMPLE_POOL_SIZE / sampler->numSamples ; j ++){
		for( int i = 0 ; i < sampler->numSamples ; i ++ ){

			int dim = sqrtf( float(sampler->numSamples ) );

			int r = i / dim;
			int c = i % dim;
	
			float t = 1 / dim;

			sampler->sample[j * sampler->numSamples + i].x = float( c * t ) + floatRand()/dim;
			sampler->sample[j * sampler->numSamples + i].y = float( r * t ) + floatRand()/dim;
		}
	}
}

__host__ inline
void sampler_Nrooks			( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_MultiJittered	( Sampler *sampler ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Hammersley		( Sampler *sampler ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__device__ inline
Point3D MapSquareToHemiSphere(Point2D point, float exp){
	float cosPhi = cosf( point.x * 2.0 * PI );
	float sinPhi = sinf( point.x * 2.0 * PI );
	float cosTheta = powf( (1.0 - point.y ), 1.0 / ( exp + 1 ));
	float sinTheta = sqrtf( 1.0 - cosTheta * cosTheta );
	float pu = sinTheta * cosTheta;
	float pv = sinTheta * sinTheta;
	float pw = cosTheta;
	return  Point3D( pu , pv, pw );
}

__device__ inline
Point3D getSampleUnitHemiSphere(Sampler *sampler, float exp){
	return MapSquareToHemiSphere( sampler->sample[ sampler->count++ % SAMPLE_POOL_SIZE ],exp );
}
	.file	"tmpxft_0037ae13_00000000-6_Sampler.cudafe1.cpp"
	.text
#APP

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB2024:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE2024:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB2047:
	.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
.LFE2047:
	.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 <hip/hip_runtime.h>
#include "Tools.cuh"
#include "defines.cuh"
#include "function_defines.cuh"


__host__ inline void sampler_PureRandom	( Sampler *sampler );
__host__ inline void sampler_Regular		( Sampler *sampler );
__host__ inline void sampler_Jittered		( Sampler *sampler );
__host__ inline void sampler_Nrooks		( Sampler *sampler );
__host__ inline void sampler_MultiJittered( Sampler *sampler );
__host__ inline void sampler_Hammersley	( Sampler *sampler );

__host__ inline
void GenerateSample(Sampler *sampler){
	switch( sampler->type ){
	case SAMPLER_PURERANDOM :
		sampler_PureRandom(sampler);
		break;
	case SAMPLER_REGULAR :
		sampler_Regular(sampler);
		break;
	case SAMPLER_JITTERED :
		sampler_Jittered(sampler);
		break;
	case SAMPLER_NROOKS :
		sampler_Nrooks(sampler);
		break;
	case SAMPLER_MULTIJITTERED :
		sampler_MultiJittered(sampler);
		break;
	case SAMPLER_HAMMERSLEY :
		sampler_Hammersley(sampler);
		break;
	default:
		break;
	}
}

__device__ inline
Point2D getSampleUnitSquare(Sampler *sampler){
	return sampler->sample[ sampler->count++ % SAMPLE_POOL_SIZE ];
}

__host__ inline
void sampler_PureRandom		( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Regular		( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Jittered		( Sampler *sampler  ){

	for(int j = 0 ; j < SAMPLE_POOL_SIZE / sampler->numSamples ; j ++){
		for( int i = 0 ; i < sampler->numSamples ; i ++ ){

			int dim = sqrtf( float(sampler->numSamples ) );

			int r = i / dim;
			int c = i % dim;
	
			float t = 1 / dim;

			sampler->sample[j * sampler->numSamples + i].x = float( c * t ) + floatRand()/dim;
			sampler->sample[j * sampler->numSamples + i].y = float( r * t ) + floatRand()/dim;
		}
	}
}

__host__ inline
void sampler_Nrooks			( Sampler *sampler  ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_MultiJittered	( Sampler *sampler ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__host__ inline 
void sampler_Hammersley		( Sampler *sampler ){
	/*not yet..*/
	for( int i = 0 ; i < sampler->numSamples ; i ++ ){
		sampler->sample[i].x = 0;
		sampler->sample[i].y = 0;
	}
}

__device__ inline
Point3D MapSquareToHemiSphere(Point2D point, float exp){
	float cosPhi = cosf( point.x * 2.0 * PI );
	float sinPhi = sinf( point.x * 2.0 * PI );
	float cosTheta = powf( (1.0 - point.y ), 1.0 / ( exp + 1 ));
	float sinTheta = sqrtf( 1.0 - cosTheta * cosTheta );
	float pu = sinTheta * cosTheta;
	float pv = sinTheta * sinTheta;
	float pw = cosTheta;
	return  Point3D( pu , pv, pw );
}

__device__ inline
Point3D getSampleUnitHemiSphere(Sampler *sampler, float exp){
	return MapSquareToHemiSphere( sampler->sample[ sampler->count++ % SAMPLE_POOL_SIZE ],exp );
}
	.text
	.file	"Sampler.hip"
	.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
