#pragma once

#include <iostream>
#include <sstream>
#include <stdexcept>


template <typename T, int TD_T, int BLOCK_ROWS>
__global__ void transposeNoBankConflicts(T *odata, T *idata, int width, int height)
{
	__shared__ T tile[TD_T][TD_T+1];

	int xIndex = blockIdx.x * TD_T + threadIdx.x;
	int yIndex = blockIdx.y * TD_T + threadIdx.y;  

	int index_in = xIndex + yIndex * width;

	xIndex = blockIdx.y * TD_T + threadIdx.x;
	yIndex = blockIdx.x * TD_T + threadIdx.y;

	int index_out = xIndex + (yIndex) * height;

	if ( xIndex < height && yIndex < width)
	{
		for (int i=0; i<TD_T; i+=BLOCK_ROWS) 
		{
			tile[threadIdx.y + i][threadIdx.x] = idata[ index_in + i * width];
		}	
 
		__syncthreads();

		for ( int i = 0; i < TD_T; i += BLOCK_ROWS) 
		{
			odata[ index_out + i * height] = tile[ threadIdx.x ][ threadIdx.y + i ];
		}
	}
}

template <typename T>
void transpose_2( size_t height, size_t width, T* idata, T* odata, int tile_dimension,cudaStream_t stream )
{
	int gridx=width / tile_dimension;
	if	( width % tile_dimension != 0)
	{
		++ gridx;
		std::stringstream ss;
		ss << "Transpose 2: Width " << width << " is not divisible by tile dimension: " << tile_dimension << ". Aborting\n";
		throw std::runtime_error( ss.str() );
	}
  
	int gridy = height/tile_dimension;
  	if ( height % tile_dimension != 0)
	{
		++ gridy;
		std::stringstream ss;
		ss << "Transpose 2: Height " << height<< " is not divisible by tile dimension: " << tile_dimension << ". Aborting\n";
		throw std::runtime_error( ss.str() );
	}
  
	dim3 grid( gridx, gridy), 
		threads( tile_dimension, tile_dimension );
	switch (tile_dimension)
	{
		case 2:
			transposeNoBankConflicts<T,2,2><<<grid, threads,0,stream>>>(odata, idata, width, height);
			break;
		case 4:
			transposeNoBankConflicts<T,4,4><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
		case 8:
			transposeNoBankConflicts<T,8,8><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
		case 16:
			transposeNoBankConflicts<T,16,16><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
/*		case 24:
			transposeNoBankConflicts<T,24,24><<<grid, threads>>>(odata, idata, width, height);
			break;
		case 32:
			transposeNoBankConflicts<T,32,32><<<grid, threads>>>(odata, idata, width, height);
			break;
		case 64:
			transposeNoBankConflicts<T,64,64><<<grid, threads>>>(odata, idata, width, height);
		case 128:
			transposeNoBankConflicts<T,128,128><<<grid, threads>>>(odata, idata, width, height);
*/		default:
			std::cerr << "Tile Dimension: " << tile_dimension << " not supported. Aborting\n";
			exit( -1 );
	}
}


	.file	"tmpxft_00228aec_00000000-6_transpose_kernel.cudafe1.cpp"
	.text
#APP
	.globl _ZSt21ios_base_library_initv

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

#pragma once
#include <hip/hip_runtime.h>


#include <iostream>
#include <sstream>
#include <stdexcept>


template <typename T, int TD_T, int BLOCK_ROWS>
__global__ void transposeNoBankConflicts(T *odata, T *idata, int width, int height)
{
	__shared__ T tile[TD_T][TD_T+1];

	int xIndex = blockIdx.x * TD_T + threadIdx.x;
	int yIndex = blockIdx.y * TD_T + threadIdx.y;  

	int index_in = xIndex + yIndex * width;

	xIndex = blockIdx.y * TD_T + threadIdx.x;
	yIndex = blockIdx.x * TD_T + threadIdx.y;

	int index_out = xIndex + (yIndex) * height;

	if ( xIndex < height && yIndex < width)
	{
		for (int i=0; i<TD_T; i+=BLOCK_ROWS) 
		{
			tile[threadIdx.y + i][threadIdx.x] = idata[ index_in + i * width];
		}	
 
		__syncthreads();

		for ( int i = 0; i < TD_T; i += BLOCK_ROWS) 
		{
			odata[ index_out + i * height] = tile[ threadIdx.x ][ threadIdx.y + i ];
		}
	}
}

template <typename T>
void transpose_2( size_t height, size_t width, T* idata, T* odata, int tile_dimension,hipStream_t stream )
{
	int gridx=width / tile_dimension;
	if	( width % tile_dimension != 0)
	{
		++ gridx;
		std::stringstream ss;
		ss << "Transpose 2: Width " << width << " is not divisible by tile dimension: " << tile_dimension << ". Aborting\n";
		throw std::runtime_error( ss.str() );
	}
  
	int gridy = height/tile_dimension;
  	if ( height % tile_dimension != 0)
	{
		++ gridy;
		std::stringstream ss;
		ss << "Transpose 2: Height " << height<< " is not divisible by tile dimension: " << tile_dimension << ". Aborting\n";
		throw std::runtime_error( ss.str() );
	}
  
	dim3 grid( gridx, gridy), 
		threads( tile_dimension, tile_dimension );
	switch (tile_dimension)
	{
		case 2:
			transposeNoBankConflicts<T,2,2><<<grid, threads,0,stream>>>(odata, idata, width, height);
			break;
		case 4:
			transposeNoBankConflicts<T,4,4><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
		case 8:
			transposeNoBankConflicts<T,8,8><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
		case 16:
			transposeNoBankConflicts<T,16,16><<<grid, threads,0,stream>>>(odata, idata, width, height);
		break;
/*		case 24:
			transposeNoBankConflicts<T,24,24><<<grid, threads>>>(odata, idata, width, height);
			break;
		case 32:
			transposeNoBankConflicts<T,32,32><<<grid, threads>>>(odata, idata, width, height);
			break;
		case 64:
			transposeNoBankConflicts<T,64,64><<<grid, threads>>>(odata, idata, width, height);
		case 128:
			transposeNoBankConflicts<T,128,128><<<grid, threads>>>(odata, idata, width, height);
*/		default:
			std::cerr << "Tile Dimension: " << tile_dimension << " not supported. Aborting\n";
			exit( -1 );
	}
}


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