#include <stdio.h>
//#include <cuda.h>
#include <stdlib.h>

#define N 64
#define R 3

// Kernel definition
//__global__ void MatMul(int A[N][N], int B[N][N], int C[N][N])
//{
//	//finish matrix multiplication, each thread calculate one element of C
//	int mi = threadIdx.x + blockIdx.x*blockDim.x;
//	int mj = threadIdx.y + blockIdx.y*blockDim.y;
//	C[mi][mj] = 0;
//	for(int i = 0; i < N; i++){
//			C[mi][mj] += A[mi][i] * B[i][mj];
//	}
//}

int main() {

	int a_h[N], b_h[N];
	int (*a_d)[N], (*b_d)[N];

	int size = N*N*sizeof(int);

	////allocate the memory on the device
	//cudaMalloc((void*)&a_d, size);
	//cudaMalloc((void*)&b_d, size);

	//assign values to matrixes
	for(int i=0; i<N; i++)
	{
		a_h[i] = i;
		b_h[i] = 0;
	}


	//sequential implementation
	for(int i = 0; i < N; i++){
		for(int j = -R; j < R+1; j++){
			if(!(i+j < 0 || i+j >= N)){
				b_h[i] += a_h[i+j];
			}
		}
	}

	for(int i = 0; i < N; i++){
		printf("%d\n", b_h[i]);
	}

	//copy matrixes to the device
	//cudaMemcpy(a_d, a_h, size, cudaMemcpyHostToDevice);
	//cudaMemcpy(b_d, b_h, size, cudaMemcpyHostToDevice);
	//cudaMemcpy(c_d, c_h, size, cudaMemcpyHostToDevice);

	//// launch kernel on the device by combining two-dimensional blocks with two-dimensional threads
	////define numBlocks and threadsPerBlock
	//int nThreads = 4;
	//dim3 numBlocks(N/nThreads,N/nThreads);
	//dim3 numThreads(nThreads,nThreads);
	//MatMul<<<numBlocks, numThreads>>>(a_d, b_d, c_d);

	////copy results from the device to the host
	//cudaMemcpy(c_h, c_d, size, cudaMemcpyDeviceToHost);



	////print the results
	//for(int i=0;i<N;i++){
   	//	for(int j=0;j<N;j++){
        //		printf("%d ", c_h[i][j]);
    	//	}	
    	//	printf("\n");
	//}

	////free the memory
	//cudaFree(a_d); 
	//cudaFree(b_d); 
	//cudaFree(c_d);

	return 0;
}

	.file	"tmpxft_0033e04b_00000000-6_e3_sequential.cudafe1.cpp"
	.text
#APP

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB2030:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE2030:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.section	.rodata.str1.1,"aMS",@progbits,1
.LC0:
	.string	"%d\n"
	.section	.text.startup,"ax",@progbits
	.globl	main
	.type	main, @function
main:
.LFB2027:
	.cfi_startproc
	endbr64
	pushq	%r12
	.cfi_def_cfa_offset 16
	.cfi_offset 12, -16
	xorl	%ecx, %ecx
	pushq	%rbp
	.cfi_def_cfa_offset 24
	.cfi_offset 6, -24
	pushq	%rbx
	.cfi_def_cfa_offset 32
	.cfi_offset 3, -32
	subq	$528, %rsp
	.cfi_def_cfa_offset 560
	movq	%fs:40, %rax
	movq	%rax, 520(%rsp)
	xorl	%eax, %eax
.L3:
	leaq	8(%rsp), %rdx
	movl	%ecx, (%rdx,%rcx,4)
	incq	%rcx
	cmpq	$64, %rcx
	jne	.L3
	leaq	264(%rsp), %rbx
	xorl	%eax, %eax
	movq	%rbx, %rdi
	rep stosl
	xorl	%eax, %eax
.L4:
	movl	%eax, %esi
	movq	$-3, %rcx
.L6:
	leal	(%rsi,%rcx), %edi
	cmpl	$63, %edi
	ja	.L5
	movl	(%rdx,%rcx,4), %edi
	addl	%edi, (%rbx,%rax,4)
.L5:
	incq	%rcx
	cmpq	$4, %rcx
	jne	.L6
	incq	%rax
	addq	$4, %rdx
	cmpq	$64, %rax
	jne	.L4
	xorl	%ebp, %ebp
	leaq	.LC0(%rip), %r12
.L7:
	movl	(%rbx,%rbp,4), %edx
	movq	%r12, %rsi
	movl	$2, %edi
	xorl	%eax, %eax
	incq	%rbp
	call	__printf_chk@PLT
	cmpq	$64, %rbp
	jne	.L7
	movq	520(%rsp), %rax
	subq	%fs:40, %rax
	je	.L8
	call	__stack_chk_fail@PLT
.L8:
	addq	$528, %rsp
	.cfi_def_cfa_offset 32
	xorl	%eax, %eax
	popq	%rbx
	.cfi_def_cfa_offset 24
	popq	%rbp
	.cfi_def_cfa_offset 16
	popq	%r12
	.cfi_def_cfa_offset 8
	ret
	.cfi_endproc
.LFE2027:
	.size	main, .-main
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB2053:
	.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
.LFE2053:
	.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 <stdio.h>
//#include <cuda.h>
#include <stdlib.h>

#define N 64
#define R 3

// Kernel definition
//__global__ void MatMul(int A[N][N], int B[N][N], int C[N][N])
//{
//	//finish matrix multiplication, each thread calculate one element of C
//	int mi = threadIdx.x + blockIdx.x*blockDim.x;
//	int mj = threadIdx.y + blockIdx.y*blockDim.y;
//	C[mi][mj] = 0;
//	for(int i = 0; i < N; i++){
//			C[mi][mj] += A[mi][i] * B[i][mj];
//	}
//}

int main() {

	int a_h[N], b_h[N];
	int (*a_d)[N], (*b_d)[N];

	int size = N*N*sizeof(int);

	////allocate the memory on the device
	//cudaMalloc((void*)&a_d, size);
	//cudaMalloc((void*)&b_d, size);

	//assign values to matrixes
	for(int i=0; i<N; i++)
	{
		a_h[i] = i;
		b_h[i] = 0;
	}


	//sequential implementation
	for(int i = 0; i < N; i++){
		for(int j = -R; j < R+1; j++){
			if(!(i+j < 0 || i+j >= N)){
				b_h[i] += a_h[i+j];
			}
		}
	}

	for(int i = 0; i < N; i++){
		printf("%d\n", b_h[i]);
	}

	//copy matrixes to the device
	//cudaMemcpy(a_d, a_h, size, cudaMemcpyHostToDevice);
	//cudaMemcpy(b_d, b_h, size, cudaMemcpyHostToDevice);
	//cudaMemcpy(c_d, c_h, size, cudaMemcpyHostToDevice);

	//// launch kernel on the device by combining two-dimensional blocks with two-dimensional threads
	////define numBlocks and threadsPerBlock
	//int nThreads = 4;
	//dim3 numBlocks(N/nThreads,N/nThreads);
	//dim3 numThreads(nThreads,nThreads);
	//MatMul<<<numBlocks, numThreads>>>(a_d, b_d, c_d);

	////copy results from the device to the host
	//cudaMemcpy(c_h, c_d, size, cudaMemcpyDeviceToHost);



	////print the results
	//for(int i=0;i<N;i++){
   	//	for(int j=0;j<N;j++){
        //		printf("%d ", c_h[i][j]);
    	//	}	
    	//	printf("\n");
	//}

	////free the memory
	//cudaFree(a_d); 
	//cudaFree(b_d); 
	//cudaFree(c_d);

	return 0;
}

	.text
	.file	"e3_sequential.hip"
	.globl	main                            # -- Begin function main
	.type	main,@function
main:                                   # @main
	.cfi_startproc
# %bb.0:
	pushq	%rbx
	.cfi_def_cfa_offset 16
	subq	$512, %rsp                      # imm = 0x200
	.cfi_def_cfa_offset 528
	.cfi_offset %rbx, -16
	movq	%rsp, %rdi
	xorl	%ebx, %ebx
	movl	$256, %edx                      # imm = 0x100
	xorl	%esi, %esi
	callq	memset@PLT
.LBB0_1:                                # =>This Inner Loop Header: Depth=1
	movl	%ebx, 256(%rsp,%rbx,4)
	incq	%rbx
	cmpq	$64, %rbx
	jne	.LBB0_1
# %bb.2:                                # %.preheader24.preheader
	xorl	%eax, %eax
.LBB0_3:                                # %.preheader24
                                        # =>This Loop Header: Depth=1
                                        #     Child Loop BB0_4 Depth 2
	movl	$-3, %ecx
.LBB0_4:                                #   Parent Loop BB0_3 Depth=1
                                        # =>  This Inner Loop Header: Depth=2
	leal	(%rax,%rcx), %edx
	cmpl	$63, %edx
	ja	.LBB0_6
# %bb.5:                                #   in Loop: Header=BB0_4 Depth=2
	movl	%edx, %edx
	movl	256(%rsp,%rdx,4), %edx
	addl	%edx, (%rsp,%rax,4)
.LBB0_6:                                #   in Loop: Header=BB0_4 Depth=2
	incl	%ecx
	cmpl	$4, %ecx
	jne	.LBB0_4
# %bb.7:                                #   in Loop: Header=BB0_3 Depth=1
	incq	%rax
	cmpq	$64, %rax
	jne	.LBB0_3
# %bb.8:                                # %.preheader.preheader
	xorl	%ebx, %ebx
.LBB0_9:                                # %.preheader
                                        # =>This Inner Loop Header: Depth=1
	movl	(%rsp,%rbx,4), %esi
	movl	$.L.str, %edi
	xorl	%eax, %eax
	callq	printf
	incq	%rbx
	cmpq	$64, %rbx
	jne	.LBB0_9
# %bb.10:
	xorl	%eax, %eax
	addq	$512, %rsp                      # imm = 0x200
	.cfi_def_cfa_offset 16
	popq	%rbx
	.cfi_def_cfa_offset 8
	retq
.Lfunc_end0:
	.size	main, .Lfunc_end0-main
	.cfi_endproc
                                        # -- End function
	.type	.L.str,@object                  # @.str
	.section	.rodata.str1.1,"aMS",@progbits,1
.L.str:
	.asciz	"%d\n"
	.size	.L.str, 4

	.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
