#include <bits/stdc++.h>
#include <cuda.h>
#include <stdlib.h>

#define IFOR(v, s, e) for(int v = s; v < e; ++v)
#define UFOR(v, s, e) for(unsigned v = s; v < e; v++)

using namespace std;

class MatrixUtility
{
    public:
        void print1Dmat(double *arr, int m) {

            IFOR(i, 0, m)    
                cout << arr[i] << " ";
            cout << '\n';
        }

        void print2Dmat(double **arr, int m, int n) {
            IFOR(i, 0, m)
            {
                IFOR(j, 0, n)
                    printf("%0.9f ", arr[i][j]);
                cout << '\n';
            }
            cout << '\n';
        }

        void init_1D_mat(double *(&arr), int n) {
            arr = (double *)malloc(n * sizeof(double));
        }

        void init_2D_mat(double **(&arr), int row, int col) {
            arr = (double **)malloc(row * sizeof(double *));
            IFOR(i, 0, row)
                arr[i] = (double *)malloc(col * sizeof(double));
        }

        double **mat_add(double **(&A), double **(&B), int row, int col)
        {
            double **res = NULL;
            init_2D_mat(res, row, col);

            IFOR(i, 0, row)
                IFOR(j, 0, col)
                res[i][j] = A[i][j] + B[i][j];
            return res;
        }

        double **mat_multiply(double **(&a), double **(&b), int r1, int c1, int r2, int c2) {
            double **c = NULL;
            init_2D_mat(c, r1, c2);

            IFOR(i, 0, r1)
                IFOR(j, 0, c2)
                    IFOR(k, 0, c1)
                        c[i][j] = c[i][j] + a[i][k] * b[k][j];

            return c;
        }


        double *vector_add(double *(&a), double *(&b), int row) {
            double *add = NULL;
            init_1D_mat(add, row);

            IFOR(i, 0, row)
                add[i] = a[i] + b[i];

            return add;
        }

        double **add_2D_mat_1D_mat(double **a, double *b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    a[i][j] += b[j];
            return a;
        }

        double **diff_2D_mat_1D_mat(double **a, double *b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    a[i][j] -= b[i];
            return a;
        }


        double **scalar_add_2D_mat(double **mat, int scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] += scalar;
            return mat;
        }

        double **scalar_divide_2D_mat(double **mat, double scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] /= scalar;
            return mat;
        }

        double **mat_transpose(double **a, int r, int c) {
            double **trans;
            init_2D_mat(trans, c, r);
            IFOR(i, 0, r)
                IFOR(j, 0, c)
                    trans[i][j] = a[j][i];
            return trans;
        }

        double **scalar_multiply_2D_mat(double **mat, int scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] *= scalar;
            return mat;
        }


        double *scalar_divide_1D_mat(double *mat, int scalar, int r) {
            UFOR(i, 0, r)
                mat[i] /= scalar;
            return mat;
        }

        double *scalar_multiply_1D_mat(double *mat, int scalar, int r) {
            UFOR(i, 0, r)
            {
                mat[i] *= scalar;
            }
            return mat;
        }



        double *sum_across_2nd_dim(double **a, int r, int c) {
            double *sum;
            init_1D_mat(sum, r);
            UFOR(i, 0, r)
            {
                sum[i] = 0;
                UFOR(j, 0, c)
                    sum[i] += a[i][j];
            }
            return sum;
        }

        double **element_wise_multiply(double **a, double **b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                a[i][j] *= b[i][j];
            return a;
        }

};

	.file	"tmpxft_002abc55_00000000-6_matrix_utility.cudafe1.cpp"
	.text
#APP
	.globl _ZSt21ios_base_library_initv

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB10844:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE10844:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB10867:
	.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
.LFE10867:
	.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 <bits/stdc++.h>
#include <hip/hip_runtime.h>
#include <stdlib.h>

#define IFOR(v, s, e) for(int v = s; v < e; ++v)
#define UFOR(v, s, e) for(unsigned v = s; v < e; v++)

using namespace std;

class MatrixUtility
{
    public:
        void print1Dmat(double *arr, int m) {

            IFOR(i, 0, m)    
                cout << arr[i] << " ";
            cout << '\n';
        }

        void print2Dmat(double **arr, int m, int n) {
            IFOR(i, 0, m)
            {
                IFOR(j, 0, n)
                    printf("%0.9f ", arr[i][j]);
                cout << '\n';
            }
            cout << '\n';
        }

        void init_1D_mat(double *(&arr), int n) {
            arr = (double *)malloc(n * sizeof(double));
        }

        void init_2D_mat(double **(&arr), int row, int col) {
            arr = (double **)malloc(row * sizeof(double *));
            IFOR(i, 0, row)
                arr[i] = (double *)malloc(col * sizeof(double));
        }

        double **mat_add(double **(&A), double **(&B), int row, int col)
        {
            double **res = NULL;
            init_2D_mat(res, row, col);

            IFOR(i, 0, row)
                IFOR(j, 0, col)
                res[i][j] = A[i][j] + B[i][j];
            return res;
        }

        double **mat_multiply(double **(&a), double **(&b), int r1, int c1, int r2, int c2) {
            double **c = NULL;
            init_2D_mat(c, r1, c2);

            IFOR(i, 0, r1)
                IFOR(j, 0, c2)
                    IFOR(k, 0, c1)
                        c[i][j] = c[i][j] + a[i][k] * b[k][j];

            return c;
        }


        double *vector_add(double *(&a), double *(&b), int row) {
            double *add = NULL;
            init_1D_mat(add, row);

            IFOR(i, 0, row)
                add[i] = a[i] + b[i];

            return add;
        }

        double **add_2D_mat_1D_mat(double **a, double *b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    a[i][j] += b[j];
            return a;
        }

        double **diff_2D_mat_1D_mat(double **a, double *b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    a[i][j] -= b[i];
            return a;
        }


        double **scalar_add_2D_mat(double **mat, int scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] += scalar;
            return mat;
        }

        double **scalar_divide_2D_mat(double **mat, double scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] /= scalar;
            return mat;
        }

        double **mat_transpose(double **a, int r, int c) {
            double **trans;
            init_2D_mat(trans, c, r);
            IFOR(i, 0, r)
                IFOR(j, 0, c)
                    trans[i][j] = a[j][i];
            return trans;
        }

        double **scalar_multiply_2D_mat(double **mat, int scalar, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                    mat[i][j] *= scalar;
            return mat;
        }


        double *scalar_divide_1D_mat(double *mat, int scalar, int r) {
            UFOR(i, 0, r)
                mat[i] /= scalar;
            return mat;
        }

        double *scalar_multiply_1D_mat(double *mat, int scalar, int r) {
            UFOR(i, 0, r)
            {
                mat[i] *= scalar;
            }
            return mat;
        }



        double *sum_across_2nd_dim(double **a, int r, int c) {
            double *sum;
            init_1D_mat(sum, r);
            UFOR(i, 0, r)
            {
                sum[i] = 0;
                UFOR(j, 0, c)
                    sum[i] += a[i][j];
            }
            return sum;
        }

        double **element_wise_multiply(double **a, double **b, int r, int c) {
            UFOR(i, 0, r)
                UFOR(j, 0, c)
                a[i][j] *= b[i][j];
            return a;
        }

};

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