布尔乘法答案

【问题标题】：Boolean Multiplication布尔乘法
【发布时间】：2014-02-21 01:15:40
【问题描述】：

我想使用 cuda 6 进行布尔乘法，但我无法以正确的方式进行。 B 是一个布尔对称矩阵，我必须做 B^n 布尔乘法。

我的 C++ 代码是：

for (m=0; m<dimension; m++) {
    for (n=0; n<dimension; n++) {
        for (k=0; k<dimension; k++) {
            temp=mat[m*dimension+k] && mat[k*dimension+n];
            B[m*dimension+n]= B[m*dimension+n] || temp;
        }
    }
}

我找到了一段代码，但我无法根据我的情况调整它。

__global__ void gpuMM(float *A, float *B, float *C, int N)
{
    int row = blockIdx.y*blockDim.y + threadIdx.y;
    int col = blockIdx.x*blockDim.x + threadIdx.x;

    float sum = 0.f;
    for (int n = 0; n < N; ++n)
         sum += A[row*N+n]*B[n*N+col];

    C[row*N+col] = sum;
}

【问题讨论】：

虽然它不会是一个高效的matrixMul，但它是直截了当的。将代码的最内层循环而不是内核中的循环。还将sum、row、col和N替换为temp、m、n和dimension。
你正在做的是正则乘法B = A * transpose(A)。 ` && ` 等价于* 而|| 等价于+。只需使用 cublas 将这两个矩阵相乘即可。如果 A 是布尔值，请先转换为浮点数。

标签： c++ c cuda parallel-processing multiplication

【解决方案1】：

在代码的第一部分使用内核的表示法，您会得到这个。

for (row=0; row<N; row++) {
    for (col=0; col<N; col++) {
        for (n=0; n<N; n++) {
            temp=mat[row*N+n] && mat[n*N+col];
            B[row*N+col] = B[row*N+col] || temp; 
        }       
    }
}

所以你的内核应该是这样的：

__global__ void gpu_booleanMM(char *mat, char *B, int N)
{
    int row = blockIdx.y*blockDim.y + threadIdx.y;
    int col = blockIdx.x*blockDim.x + threadIdx.x;

    for (int n=0; n<N; n++) {
        temp=mat[row*N+n] && mat[n*N+col];
        B[row*N+col] = B[row*N+col] || temp; 
    }   
}

我怀疑这是否非常有效，但这样的事情仍然应该给出正确的结果。

【讨论】：