

__global__ void variance2(double *g_idata, double *g_odata, double* reduce_sum, double sum_divider){

	double mean = reduce_sum[0] / sum_divider;	
	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// each thread loads one element from global to shared mem
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*blockDim.x + threadIdx.x;
	sdata[tid] = pow(g_idata[i] - mean, 2.0);
	//sdata[tid] = g_idata[i] - mean;
	//sdata[tid] *= sdata[tid];
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] += sdata[tid+s];
		}
		__syncthreads();
	}

	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}

// HALVE THE NUMBER OF BLOCKS, AND REPLACE SINGLE LOAD
__global__ void variance3(double *g_idata, double *g_odata, double* reduce_sum, double sum_divider){

	/*
	jeli zmniejszamy liczb blokw o poow to tablica wynikowa te powinna by krtsza poow
	*/
	
	double mean = reduce_sum[0] / sum_divider;	
	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// perform first level of reduction
	// reading from global memory, writing to shared memory
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*(blockDim.x*2) + threadIdx.x;
	sdata[tid] = pow(g_idata[i] - mean, 2.0) + pow(g_idata[i+blockDim.x] - mean, 2.0);
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] += sdata[tid+s];
		}
		__syncthreads();
	}

	// poczekaj ...
	// konsekwencj tego, e zmniejszamy liczb blokw o poow jest to, e mamy krtsz o poow tablic wynikow
	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}


__global__ void norm_fro(double *g_idata, double *g_odata){

	/*
	function performs operation: sqrt(sum(diag(v'*v)))
	*/

	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// each thread loads one element from global to shared mem
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*blockDim.x + threadIdx.x;
	//sdata[tid] = pow(g_idata[i], 2.0);
	sdata[tid] = g_idata[i]*g_idata[i];
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] += sdata[tid+s];
		}
		__syncthreads();
	}

	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}

__global__ void norm_fro2(double *g_idata1, double *g_idata2, double *g_odata){

	/*
	function performs operation: sqrt(sum(diag(v'*v))) for v = g_idata1 - g_idata2
	*/

	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// each thread loads one element from global to shared mem
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*blockDim.x + threadIdx.x;
	sdata[tid] = pow(g_idata1[i] - g_idata2[i], 2.0);
	//sdata[tid] *= sdata[tid];
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] += sdata[tid+s];
		}
		__syncthreads();
	}

	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}

__global__ void max_value(double *g_idata, double *g_odata){
	
	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// each thread loads one element from global to shared mem
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*blockDim.x + threadIdx.x;

	sdata[tid] = g_idata[i];
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] = sdata[tid]*(sdata[tid]>=sdata[tid+s]) + sdata[tid+s]*(sdata[tid]<sdata[tid+s]);

			// a moe w sdata[tid+s] zapisywa min ?? - do przemylenia
		}
		__syncthreads();
	}

	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}

__global__ void min_value(double *g_idata, double *g_odata){
	
	// Sequential addressing is conflict free
	extern __shared__ double sdata[];

	// each thread loads one element from global to shared mem
	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*blockDim.x + threadIdx.x;

	sdata[tid] = g_idata[i];
	__syncthreads();

	for (unsigned int s=blockDim.x/2; s>0; s>>=1){
		if (tid < s){
			sdata[tid] = sdata[tid]*(sdata[tid]<=sdata[tid+s]) + sdata[tid+s]*(sdata[tid]>sdata[tid+s]);

			// a moe w sdata[tid+s] zapisywa min ?? - do przemylenia
		}
		__syncthreads();
	}

	if (tid==0) g_odata[blockIdx.x] = sdata[0];
}

	.file	"tmpxft_00263274_00000000-6_reductionMixed.cudafe1.cpp"
	.text
#APP

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB2010:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE2010:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.globl	_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d
	.type	_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d, @function
_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d:
.LFB2032:
	.cfi_startproc
	endbr64
	subq	$152, %rsp
	.cfi_def_cfa_offset 160
	movq	%rdi, 24(%rsp)
	leaq	48(%rsp), %rcx
	leaq	56(%rsp), %rdi
	movq	%rsi, 16(%rsp)
	leaq	68(%rsp), %rsi
	movq	%rdx, 8(%rsp)
	leaq	40(%rsp), %rdx
	movsd	%xmm0, (%rsp)
	movq	%fs:40, %rax
	movq	%rax, 136(%rsp)
	xorl	%eax, %eax
	leaq	24(%rsp), %rax
	movl	$1, 64(%rsp)
	movq	%rax, 104(%rsp)
	leaq	16(%rsp), %rax
	movq	%rax, 112(%rsp)
	leaq	8(%rsp), %rax
	movq	%rax, 120(%rsp)
	movq	%rsp, %rax
	movq	%rax, 128(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 56(%rsp)
	movq	%rax, 68(%rsp)
	movl	$1, 76(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L2
	pushq	48(%rsp)
	.cfi_def_cfa_offset 168
	leaq	_Z9variance2PdS_S_d(%rip), %rdi
	pushq	48(%rsp)
	.cfi_def_cfa_offset 176
	movq	84(%rsp), %rcx
	movl	92(%rsp), %r8d
	movq	72(%rsp), %rsi
	movl	80(%rsp), %edx
	leaq	120(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 168
	popq	%rdx
	.cfi_def_cfa_offset 160
.L2:
	movq	136(%rsp), %rax
	subq	%fs:40, %rax
	je	.L4
	call	__stack_chk_fail@PLT
.L4:
	addq	$152, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2032:
	.size	_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d, .-_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d
	.globl	_Z9variance2PdS_S_d
	.type	_Z9variance2PdS_S_d, @function
_Z9variance2PdS_S_d:
.LFB2033:
	.cfi_startproc
	endbr64
	jmp	_Z33__device_stub__Z9variance2PdS_S_dPdS_S_d
	.cfi_endproc
.LFE2033:
	.size	_Z9variance2PdS_S_d, .-_Z9variance2PdS_S_d
	.globl	_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d
	.type	_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d, @function
_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d:
.LFB2034:
	.cfi_startproc
	endbr64
	subq	$152, %rsp
	.cfi_def_cfa_offset 160
	movq	%rdi, 24(%rsp)
	leaq	48(%rsp), %rcx
	leaq	56(%rsp), %rdi
	movq	%rsi, 16(%rsp)
	leaq	68(%rsp), %rsi
	movq	%rdx, 8(%rsp)
	leaq	40(%rsp), %rdx
	movsd	%xmm0, (%rsp)
	movq	%fs:40, %rax
	movq	%rax, 136(%rsp)
	xorl	%eax, %eax
	leaq	24(%rsp), %rax
	movl	$1, 64(%rsp)
	movq	%rax, 104(%rsp)
	leaq	16(%rsp), %rax
	movq	%rax, 112(%rsp)
	leaq	8(%rsp), %rax
	movq	%rax, 120(%rsp)
	movq	%rsp, %rax
	movq	%rax, 128(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 56(%rsp)
	movq	%rax, 68(%rsp)
	movl	$1, 76(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L8
	pushq	48(%rsp)
	.cfi_def_cfa_offset 168
	leaq	_Z9variance3PdS_S_d(%rip), %rdi
	pushq	48(%rsp)
	.cfi_def_cfa_offset 176
	movq	84(%rsp), %rcx
	movl	92(%rsp), %r8d
	movq	72(%rsp), %rsi
	movl	80(%rsp), %edx
	leaq	120(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 168
	popq	%rdx
	.cfi_def_cfa_offset 160
.L8:
	movq	136(%rsp), %rax
	subq	%fs:40, %rax
	je	.L10
	call	__stack_chk_fail@PLT
.L10:
	addq	$152, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2034:
	.size	_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d, .-_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d
	.globl	_Z9variance3PdS_S_d
	.type	_Z9variance3PdS_S_d, @function
_Z9variance3PdS_S_d:
.LFB2035:
	.cfi_startproc
	endbr64
	jmp	_Z33__device_stub__Z9variance3PdS_S_dPdS_S_d
	.cfi_endproc
.LFE2035:
	.size	_Z9variance3PdS_S_d, .-_Z9variance3PdS_S_d
	.globl	_Z29__device_stub__Z8norm_froPdS_PdS_
	.type	_Z29__device_stub__Z8norm_froPdS_PdS_, @function
_Z29__device_stub__Z8norm_froPdS_PdS_:
.LFB2036:
	.cfi_startproc
	endbr64
	subq	$120, %rsp
	.cfi_def_cfa_offset 128
	movq	%rdi, 8(%rsp)
	leaq	32(%rsp), %rcx
	leaq	24(%rsp), %rdx
	movq	%rsi, (%rsp)
	leaq	40(%rsp), %rdi
	leaq	52(%rsp), %rsi
	movq	%fs:40, %rax
	movq	%rax, 104(%rsp)
	xorl	%eax, %eax
	leaq	8(%rsp), %rax
	movl	$1, 48(%rsp)
	movq	%rax, 88(%rsp)
	movq	%rsp, %rax
	movq	%rax, 96(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 40(%rsp)
	movq	%rax, 52(%rsp)
	movl	$1, 60(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L13
	pushq	32(%rsp)
	.cfi_def_cfa_offset 136
	leaq	_Z8norm_froPdS_(%rip), %rdi
	pushq	32(%rsp)
	.cfi_def_cfa_offset 144
	movq	68(%rsp), %rcx
	movl	76(%rsp), %r8d
	movq	56(%rsp), %rsi
	movl	64(%rsp), %edx
	leaq	104(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 136
	popq	%rdx
	.cfi_def_cfa_offset 128
.L13:
	movq	104(%rsp), %rax
	subq	%fs:40, %rax
	je	.L15
	call	__stack_chk_fail@PLT
.L15:
	addq	$120, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2036:
	.size	_Z29__device_stub__Z8norm_froPdS_PdS_, .-_Z29__device_stub__Z8norm_froPdS_PdS_
	.globl	_Z8norm_froPdS_
	.type	_Z8norm_froPdS_, @function
_Z8norm_froPdS_:
.LFB2037:
	.cfi_startproc
	endbr64
	jmp	_Z29__device_stub__Z8norm_froPdS_PdS_
	.cfi_endproc
.LFE2037:
	.size	_Z8norm_froPdS_, .-_Z8norm_froPdS_
	.globl	_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_
	.type	_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_, @function
_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_:
.LFB2038:
	.cfi_startproc
	endbr64
	subq	$136, %rsp
	.cfi_def_cfa_offset 144
	movq	%rdi, 24(%rsp)
	leaq	40(%rsp), %rcx
	leaq	48(%rsp), %rdi
	movq	%rsi, 16(%rsp)
	leaq	60(%rsp), %rsi
	movq	%rdx, 8(%rsp)
	leaq	32(%rsp), %rdx
	movq	%fs:40, %rax
	movq	%rax, 120(%rsp)
	xorl	%eax, %eax
	leaq	24(%rsp), %rax
	movl	$1, 56(%rsp)
	movq	%rax, 96(%rsp)
	leaq	16(%rsp), %rax
	movq	%rax, 104(%rsp)
	leaq	8(%rsp), %rax
	movq	%rax, 112(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 48(%rsp)
	movq	%rax, 60(%rsp)
	movl	$1, 68(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L18
	pushq	40(%rsp)
	.cfi_def_cfa_offset 152
	leaq	_Z9norm_fro2PdS_S_(%rip), %rdi
	pushq	40(%rsp)
	.cfi_def_cfa_offset 160
	movq	76(%rsp), %rcx
	movl	84(%rsp), %r8d
	movq	64(%rsp), %rsi
	movl	72(%rsp), %edx
	leaq	112(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 152
	popq	%rdx
	.cfi_def_cfa_offset 144
.L18:
	movq	120(%rsp), %rax
	subq	%fs:40, %rax
	je	.L20
	call	__stack_chk_fail@PLT
.L20:
	addq	$136, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2038:
	.size	_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_, .-_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_
	.globl	_Z9norm_fro2PdS_S_
	.type	_Z9norm_fro2PdS_S_, @function
_Z9norm_fro2PdS_S_:
.LFB2039:
	.cfi_startproc
	endbr64
	jmp	_Z32__device_stub__Z9norm_fro2PdS_S_PdS_S_
	.cfi_endproc
.LFE2039:
	.size	_Z9norm_fro2PdS_S_, .-_Z9norm_fro2PdS_S_
	.globl	_Z30__device_stub__Z9max_valuePdS_PdS_
	.type	_Z30__device_stub__Z9max_valuePdS_PdS_, @function
_Z30__device_stub__Z9max_valuePdS_PdS_:
.LFB2040:
	.cfi_startproc
	endbr64
	subq	$120, %rsp
	.cfi_def_cfa_offset 128
	movq	%rdi, 8(%rsp)
	leaq	32(%rsp), %rcx
	leaq	24(%rsp), %rdx
	movq	%rsi, (%rsp)
	leaq	40(%rsp), %rdi
	leaq	52(%rsp), %rsi
	movq	%fs:40, %rax
	movq	%rax, 104(%rsp)
	xorl	%eax, %eax
	leaq	8(%rsp), %rax
	movl	$1, 48(%rsp)
	movq	%rax, 88(%rsp)
	movq	%rsp, %rax
	movq	%rax, 96(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 40(%rsp)
	movq	%rax, 52(%rsp)
	movl	$1, 60(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L23
	pushq	32(%rsp)
	.cfi_def_cfa_offset 136
	leaq	_Z9max_valuePdS_(%rip), %rdi
	pushq	32(%rsp)
	.cfi_def_cfa_offset 144
	movq	68(%rsp), %rcx
	movl	76(%rsp), %r8d
	movq	56(%rsp), %rsi
	movl	64(%rsp), %edx
	leaq	104(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 136
	popq	%rdx
	.cfi_def_cfa_offset 128
.L23:
	movq	104(%rsp), %rax
	subq	%fs:40, %rax
	je	.L25
	call	__stack_chk_fail@PLT
.L25:
	addq	$120, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2040:
	.size	_Z30__device_stub__Z9max_valuePdS_PdS_, .-_Z30__device_stub__Z9max_valuePdS_PdS_
	.globl	_Z9max_valuePdS_
	.type	_Z9max_valuePdS_, @function
_Z9max_valuePdS_:
.LFB2041:
	.cfi_startproc
	endbr64
	jmp	_Z30__device_stub__Z9max_valuePdS_PdS_
	.cfi_endproc
.LFE2041:
	.size	_Z9max_valuePdS_, .-_Z9max_valuePdS_
	.globl	_Z30__device_stub__Z9min_valuePdS_PdS_
	.type	_Z30__device_stub__Z9min_valuePdS_PdS_, @function
_Z30__device_stub__Z9min_valuePdS_PdS_:
.LFB2042:
	.cfi_startproc
	endbr64
	subq	$120, %rsp
	.cfi_def_cfa_offset 128
	movq	%rdi, 8(%rsp)
	leaq	32(%rsp), %rcx
	leaq	24(%rsp), %rdx
	movq	%rsi, (%rsp)
	leaq	40(%rsp), %rdi
	leaq	52(%rsp), %rsi
	movq	%fs:40, %rax
	movq	%rax, 104(%rsp)
	xorl	%eax, %eax
	leaq	8(%rsp), %rax
	movl	$1, 48(%rsp)
	movq	%rax, 88(%rsp)
	movq	%rsp, %rax
	movq	%rax, 96(%rsp)
	movabsq	$4294967297, %rax
	movq	%rax, 40(%rsp)
	movq	%rax, 52(%rsp)
	movl	$1, 60(%rsp)
	call	__cudaPopCallConfiguration@PLT
	testl	%eax, %eax
	jne	.L28
	pushq	32(%rsp)
	.cfi_def_cfa_offset 136
	leaq	_Z9min_valuePdS_(%rip), %rdi
	pushq	32(%rsp)
	.cfi_def_cfa_offset 144
	movq	68(%rsp), %rcx
	movl	76(%rsp), %r8d
	movq	56(%rsp), %rsi
	movl	64(%rsp), %edx
	leaq	104(%rsp), %r9
	call	cudaLaunchKernel@PLT
	popq	%rax
	.cfi_def_cfa_offset 136
	popq	%rdx
	.cfi_def_cfa_offset 128
.L28:
	movq	104(%rsp), %rax
	subq	%fs:40, %rax
	je	.L30
	call	__stack_chk_fail@PLT
.L30:
	addq	$120, %rsp
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2042:
	.size	_Z30__device_stub__Z9min_valuePdS_PdS_, .-_Z30__device_stub__Z9min_valuePdS_PdS_
	.globl	_Z9min_valuePdS_
	.type	_Z9min_valuePdS_, @function
_Z9min_valuePdS_:
.LFB2043:
	.cfi_startproc
	endbr64
	jmp	_Z30__device_stub__Z9min_valuePdS_PdS_
	.cfi_endproc
.LFE2043:
	.size	_Z9min_valuePdS_, .-_Z9min_valuePdS_
	.section	.rodata.str1.1,"aMS",@progbits,1
.LC0:
	.string	"_Z9min_valuePdS_"
.LC1:
	.string	"_Z9max_valuePdS_"
.LC2:
	.string	"_Z9norm_fro2PdS_S_"
.LC3:
	.string	"_Z8norm_froPdS_"
.LC4:
	.string	"_Z9variance3PdS_S_d"
.LC5:
	.string	"_Z9variance2PdS_S_d"
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB2045:
	.cfi_startproc
	endbr64
	pushq	%rbx
	.cfi_def_cfa_offset 16
	.cfi_offset 3, -16
	leaq	_ZL15__fatDeviceText(%rip), %rdi
	call	__cudaRegisterFatBinary@PLT
	pushq	$0
	.cfi_def_cfa_offset 24
	xorl	%r9d, %r9d
	orl	$-1, %r8d
	pushq	$0
	.cfi_def_cfa_offset 32
	leaq	.LC0(%rip), %rdx
	movq	%rax, %rdi
	movq	%rax, %rbx
	pushq	$0
	.cfi_def_cfa_offset 40
	movq	%rdx, %rcx
	leaq	_Z9min_valuePdS_(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 48
	movq	%rax, _ZL20__cudaFatCubinHandle(%rip)
	call	__cudaRegisterFunction@PLT
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	movq	%rbx, %rdi
	xorl	%r9d, %r9d
	pushq	$0
	.cfi_def_cfa_offset 24
	leaq	.LC1(%rip), %rdx
	orl	$-1, %r8d
	leaq	_Z9max_valuePdS_(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 32
	movq	%rdx, %rcx
	pushq	$0
	.cfi_def_cfa_offset 40
	pushq	$0
	.cfi_def_cfa_offset 48
	call	__cudaRegisterFunction@PLT
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	movq	%rbx, %rdi
	xorl	%r9d, %r9d
	pushq	$0
	.cfi_def_cfa_offset 24
	leaq	.LC2(%rip), %rdx
	orl	$-1, %r8d
	leaq	_Z9norm_fro2PdS_S_(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 32
	movq	%rdx, %rcx
	pushq	$0
	.cfi_def_cfa_offset 40
	pushq	$0
	.cfi_def_cfa_offset 48
	call	__cudaRegisterFunction@PLT
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	movq	%rbx, %rdi
	xorl	%r9d, %r9d
	pushq	$0
	.cfi_def_cfa_offset 24
	leaq	.LC3(%rip), %rdx
	orl	$-1, %r8d
	leaq	_Z8norm_froPdS_(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 32
	movq	%rdx, %rcx
	pushq	$0
	.cfi_def_cfa_offset 40
	pushq	$0
	.cfi_def_cfa_offset 48
	call	__cudaRegisterFunction@PLT
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	movq	%rbx, %rdi
	xorl	%r9d, %r9d
	pushq	$0
	.cfi_def_cfa_offset 24
	leaq	.LC4(%rip), %rdx
	orl	$-1, %r8d
	leaq	_Z9variance3PdS_S_d(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 32
	movq	%rdx, %rcx
	pushq	$0
	.cfi_def_cfa_offset 40
	pushq	$0
	.cfi_def_cfa_offset 48
	call	__cudaRegisterFunction@PLT
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	movq	%rbx, %rdi
	xorl	%r9d, %r9d
	pushq	$0
	.cfi_def_cfa_offset 24
	leaq	.LC5(%rip), %rdx
	orl	$-1, %r8d
	leaq	_Z9variance2PdS_S_d(%rip), %rsi
	pushq	$0
	.cfi_def_cfa_offset 32
	movq	%rdx, %rcx
	pushq	$0
	.cfi_def_cfa_offset 40
	pushq	$0
	.cfi_def_cfa_offset 48
	call	__cudaRegisterFunction@PLT
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	addq	$32, %rsp
	.cfi_def_cfa_offset 16
	call	__cudaRegisterFatBinaryEnd@PLT
	leaq	_ZL26__cudaUnregisterBinaryUtilv(%rip), %rdi
	popq	%rbx
	.cfi_def_cfa_offset 8
	jmp	atexit@PLT
	.cfi_endproc
.LFE2045:
	.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>

//#ifndef VARIANCE_CUDA_H
//#define VARIANCE_CUDA_H
__global__ void variance2(double *g_idata, double *g_odata, double* reduce_sum, double sum_divider);
__global__ void norm_fro(double *g_idata, double *g_odata);
__global__ void norm_fro2(double *g_idata1, double *g_idata2, double *g_odata);

__global__ void min_value(double *g_idata, double *g_odata);
__global__ void max_value(double *g_idata, double *g_odata);

template <unsigned int blockSize> __global__ void variance6(double *g_idata, double* g_odata, unsigned int n, double* reduce_sum, double sum_divider);

template <unsigned int blockSize>
__global__ void variance6(double *g_idata, double* g_odata, unsigned int n, double* reduce_sum, double sum_divider){

	double mean = reduce_sum[0] / sum_divider;	
	extern __shared__ double sdata[];

	unsigned int tid = threadIdx.x;
	unsigned int i = blockIdx.x*(blockSize*2) + tid;
	unsigned int gridSize = blockSize*2*gridDim.x;
	sdata[tid] = 0;

	while (i < n) { sdata[tid] += pow(g_idata[i] - mean, 2.0) + pow(g_idata[i+blockSize] - mean, 2.0);  i += gridSize; }
	__syncthreads();

	if (blockSize >= 512) { if (tid < 256) { sdata[tid] += sdata[tid + 256]; } __syncthreads(); }
	if (blockSize >= 256) { if (tid < 128) { sdata[tid] += sdata[tid + 128]; } __syncthreads(); }
	if (blockSize >= 128) { if (tid <   64) { sdata[tid] += sdata[tid +   64]; } __syncthreads(); }

	if (tid < 32) {
		if (blockSize >=  64) sdata[tid] += sdata[tid + 32];
		if (blockSize >=  32) sdata[tid] += sdata[tid + 16];
		if (blockSize >=  16) sdata[tid] += sdata[tid +  8];
		if (blockSize >=    8) sdata[tid] += sdata[tid +  4];
		if (blockSize >=    4) sdata[tid] += sdata[tid +  2];
		if (blockSize >=    2) sdata[tid] += sdata[tid +  1];
	}

	if (tid == 0) g_odata[blockIdx.x] = sdata[0];
}

//#endif
	.text
	.file	"reductionMixed.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
