


/*!
 * Compute the Pearson correlation of a cluster in a pairwise data array.
 *
 * @param x
 * @param y
 * @param labels
 * @param sampleSize
 * @param cluster
 * @param minSamples
 */
__device__
float Pearson_computeCluster(
    const float *x,
    const float *y,
    const char *labels,
    int sampleSize,
    char cluster,
    int minSamples)
{
    // compute intermediate sums
    int n = 0;
    float sumx = 0;
    float sumy = 0;
    float sumx2 = 0;
    float sumy2 = 0;
    float sumxy = 0;

    for ( int i = 0; i < sampleSize; ++i )
    {
        if ( labels[i] == cluster )
        {
            float x_i = x[i];
            float y_i = y[i];

            sumx += x_i;
            sumy += y_i;
            sumx2 += x_i * x_i;
            sumy2 += y_i * y_i;
            sumxy += x_i * y_i;

            ++n;
        }
    }

    // compute correlation only if there are enough samples
    float result = NAN;

    if ( n >= minSamples )
    {
        result = (n*sumxy - sumx*sumy) / sqrt((n*sumx2 - sumx*sumx) * (n*sumy2 - sumy*sumy));
    }

    return result;
}



/*!
 * Compute the correlation of each cluster in a pairwise data array. The data array
 * should only contain the clean samples that were extracted from the expression
 * matrix, while the labels should contain all samples.
 *
 * @param x
 * @param y
 * @param sampleSize
 * @param clusterSize
 * @param labels
 * @param minSamples
 * @param correlations
 */
__device__
void Pearson_compute(
    const float *x,
    const float *y,
    int sampleSize,
    char clusterSize,
    const char *labels,
    int minSamples,
    float *correlations)
{
    for ( char k = 0; k < clusterSize; ++k )
    {
        correlations[k] = Pearson_computeCluster(
            x, y,
            labels,
            sampleSize,
            k,
            minSamples
        );
    }
}

	.file	"tmpxft_0039e712_00000000-6_pearson.cudafe1.cpp"
	.text
#APP

#NO_APP
	.type	_ZL26__cudaUnregisterBinaryUtilv, @function
_ZL26__cudaUnregisterBinaryUtilv:
.LFB2012:
	.cfi_startproc
	endbr64
	movq	_ZL20__cudaFatCubinHandle(%rip), %rdi
	jmp	__cudaUnregisterFatBinary@PLT
	.cfi_endproc
.LFE2012:
	.size	_ZL26__cudaUnregisterBinaryUtilv, .-_ZL26__cudaUnregisterBinaryUtilv
	.globl	_Z22Pearson_computeClusterPKfS0_PKcici
	.type	_Z22Pearson_computeClusterPKfS0_PKcici, @function
_Z22Pearson_computeClusterPKfS0_PKcici:
.LFB2008:
	.cfi_startproc
	endbr64
	pushq	%rax
	.cfi_def_cfa_offset 16
	popq	%rax
	.cfi_def_cfa_offset 8
	subq	$24, %rsp
	.cfi_def_cfa_offset 32
	movl	$1, 12(%rsp)
	movl	12(%rsp), %edi
	call	exit@PLT
	.cfi_endproc
.LFE2008:
	.size	_Z22Pearson_computeClusterPKfS0_PKcici, .-_Z22Pearson_computeClusterPKfS0_PKcici
	.globl	_Z15Pearson_computePKfS0_icPKciPf
	.type	_Z15Pearson_computePKfS0_icPKciPf, @function
_Z15Pearson_computePKfS0_icPKciPf:
.LFB2009:
	.cfi_startproc
	endbr64
	pushq	%rax
	.cfi_def_cfa_offset 16
	popq	%rax
	.cfi_def_cfa_offset 8
	subq	$24, %rsp
	.cfi_def_cfa_offset 32
	movl	$1, 12(%rsp)
	movl	12(%rsp), %edi
	call	exit@PLT
	.cfi_endproc
.LFE2009:
	.size	_Z15Pearson_computePKfS0_icPKciPf, .-_Z15Pearson_computePKfS0_icPKciPf
	.section	.text.startup,"ax",@progbits
	.type	_ZL24__sti____cudaRegisterAllv, @function
_ZL24__sti____cudaRegisterAllv:
.LFB2035:
	.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
.LFE2035:
	.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>



/*!
 * Compute the Pearson correlation of a cluster in a pairwise data array.
 *
 * @param x
 * @param y
 * @param labels
 * @param sampleSize
 * @param cluster
 * @param minSamples
 */
__device__
float Pearson_computeCluster(
    const float *x,
    const float *y,
    const char *labels,
    int sampleSize,
    char cluster,
    int minSamples)
{
    // compute intermediate sums
    int n = 0;
    float sumx = 0;
    float sumy = 0;
    float sumx2 = 0;
    float sumy2 = 0;
    float sumxy = 0;

    for ( int i = 0; i < sampleSize; ++i )
    {
        if ( labels[i] == cluster )
        {
            float x_i = x[i];
            float y_i = y[i];

            sumx += x_i;
            sumy += y_i;
            sumx2 += x_i * x_i;
            sumy2 += y_i * y_i;
            sumxy += x_i * y_i;

            ++n;
        }
    }

    // compute correlation only if there are enough samples
    float result = NAN;

    if ( n >= minSamples )
    {
        result = (n*sumxy - sumx*sumy) / sqrt((n*sumx2 - sumx*sumx) * (n*sumy2 - sumy*sumy));
    }

    return result;
}



/*!
 * Compute the correlation of each cluster in a pairwise data array. The data array
 * should only contain the clean samples that were extracted from the expression
 * matrix, while the labels should contain all samples.
 *
 * @param x
 * @param y
 * @param sampleSize
 * @param clusterSize
 * @param labels
 * @param minSamples
 * @param correlations
 */
__device__
void Pearson_compute(
    const float *x,
    const float *y,
    int sampleSize,
    char clusterSize,
    const char *labels,
    int minSamples,
    float *correlations)
{
    for ( char k = 0; k < clusterSize; ++k )
    {
        correlations[k] = Pearson_computeCluster(
            x, y,
            labels,
            sampleSize,
            k,
            minSamples
        );
    }
}

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