#include <cuda.h>

namespace cyclic_reduction{
/**
* Calculates:
*	alpha_i = -a_i/b_(i-2^(l-i))
*	where:
*		x = a_i
*		y = b_(i-2^(l-i))
*
*	beta_i = -c_i/b_(i+2^(l-i))
*	where:
*		x = c_i
*		y = b_(i+2^(l-i))	
**/
struct AlphaBeta{
	__host__ __device__
	double operator()(double x, double y){
		if(y == 0.00){
			return 0.00;
		} else{
			return (-x)/y;
		}
	}
};





struct FirstStepFunctor{
	
	__host__ __device__
	double operator()(double a, double b){
		if(b == 0.00){
			return 0.00;
		} else{
			return -a/b;
		}	
	}
};

struct SecondStepFunctor{
	__host__ __device__
	double operator()(double c, double b){
		if(b == 0.00){
			return 0.00;
		}else{
			return -c/b;
		}
	}
};	

}





/**
* Calulates:
*	a_i = alpha_u * a_(i-2^(l-1))_prime
*	where:
*		x,y = as in AlphaBeta
*		z = a_(i-2^(l-1))_prime
*
*	c_i = beta_i * c_(i+2^(l-1))_prime
*	where:
*		x,y = as in AlphaBeta
*		z = c_(i+2^(l-1))_prime
**/
struct AC{
	__host__ __device__
	double operator()(double x, double y, double z){
		cyclic_reduction::AlphaBeta alphaBeta;
		return alphaBeta(x,y)*z;
	}
};



/**
* Calculates:
*	b_i = b_i_prime + alpha_i * c_(i-2^(l-1))_prime + beta_i * a_(i+2^(l-1))_prime
* 	where:
*		alpha_i , beta_i are results from AlphaBeta
*		x = b_i_prime
* 		y = c_(i-2^(l-1))_prime
*		z = (a_(i+2^(l-1))_prime
*
* 	d_i = d_i_prime + alpha_i * d_(i-2^(l-1))_prime + beta_i * d_(i+2^(l-1))_prime
* 	where:
*		alpha_i , beta_i are results from AlphaBeta
*		x = d_i_prime
* 		y = d_(i-2^(l-1))_prime
*		z = (d_(i+2^(l-1))_prime
*
*		
**/
struct BD{
	__host__ __device__
	double operator()(double x, double alpha, double y, double beta, double z){
		double a = alpha*y;
		double b = beta*z;
		return x + a + b;
	}
};


	.file	"tmpxft_0033b20b_00000000-6_cu_cr_functors.cudafe1.cpp"
	.text
#APP

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB2015:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE2015:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB2038:
	.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
.LFE2038:
	.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>

namespace cyclic_reduction{
/**
* Calculates:
*	alpha_i = -a_i/b_(i-2^(l-i))
*	where:
*		x = a_i
*		y = b_(i-2^(l-i))
*
*	beta_i = -c_i/b_(i+2^(l-i))
*	where:
*		x = c_i
*		y = b_(i+2^(l-i))	
**/
struct AlphaBeta{
	__host__ __device__
	double operator()(double x, double y){
		if(y == 0.00){
			return 0.00;
		} else{
			return (-x)/y;
		}
	}
};





struct FirstStepFunctor{
	
	__host__ __device__
	double operator()(double a, double b){
		if(b == 0.00){
			return 0.00;
		} else{
			return -a/b;
		}	
	}
};

struct SecondStepFunctor{
	__host__ __device__
	double operator()(double c, double b){
		if(b == 0.00){
			return 0.00;
		}else{
			return -c/b;
		}
	}
};	

}





/**
* Calulates:
*	a_i = alpha_u * a_(i-2^(l-1))_prime
*	where:
*		x,y = as in AlphaBeta
*		z = a_(i-2^(l-1))_prime
*
*	c_i = beta_i * c_(i+2^(l-1))_prime
*	where:
*		x,y = as in AlphaBeta
*		z = c_(i+2^(l-1))_prime
**/
struct AC{
	__host__ __device__
	double operator()(double x, double y, double z){
		cyclic_reduction::AlphaBeta alphaBeta;
		return alphaBeta(x,y)*z;
	}
};



/**
* Calculates:
*	b_i = b_i_prime + alpha_i * c_(i-2^(l-1))_prime + beta_i * a_(i+2^(l-1))_prime
* 	where:
*		alpha_i , beta_i are results from AlphaBeta
*		x = b_i_prime
* 		y = c_(i-2^(l-1))_prime
*		z = (a_(i+2^(l-1))_prime
*
* 	d_i = d_i_prime + alpha_i * d_(i-2^(l-1))_prime + beta_i * d_(i+2^(l-1))_prime
* 	where:
*		alpha_i , beta_i are results from AlphaBeta
*		x = d_i_prime
* 		y = d_(i-2^(l-1))_prime
*		z = (d_(i+2^(l-1))_prime
*
*		
**/
struct BD{
	__host__ __device__
	double operator()(double x, double alpha, double y, double beta, double z){
		double a = alpha*y;
		double b = beta*z;
		return x + a + b;
	}
};


	.text
	.file	"cu_cr_functors.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
