Cublas Matrix LU Разложение

Question

Cublas Matrix LU Разложение

У меня проблемы с вызовом dgetrf в cuda. Из того, что я нашел, я могу назвать только пакетную версию (http://docs.nvidia.com/cuda/cublas/#cublas-lt-t-gt-getrfbatched). Когда я вызываю его, я получаю значение ошибки, равное 7, и мне не удалось найти соответствующее перечисление для этого кода ошибки. Ниже мой код, любая помощь будет высоко ценится;

void cuda_matrix_inverse (int m, int n, double* a){

cublasHandle_t handle;
cublasStatus_t status;
double **devPtrA = 0;
double **devPtrA_dev = NULL;
int *d_pivot_array;
int *d_info_array;
int rowsA = m;
int colsA = n;
int matrixSizeA;
cudaError_t error;

fprintf(stderr,"starting cuda inverse\n");

error = cudaMalloc((void **)&d_pivot_array, sizeof(int));
if (error != cudaSuccess) fprintf(stderr,"\nError: %s\n",cudaGetErrorString(error));
error = cudaMalloc((void **)&d_info_array, sizeof(int));
if (error != cudaSuccess) fprintf(stderr,"\nError: %s\n",cudaGetErrorString(error));

fprintf(stderr,"malloced pivot and info\n");

status = cublasCreate(&handle);
if (status != CUBLAS_STATUS_SUCCESS) fprintf(stderr,"error %i\n",status);

matrixSizeA = rowsA * colsA;

devPtrA =(double **)malloc(1 * sizeof(*devPtrA));

fprintf(stderr,"malloced devPtrA\n");

error = cudaMalloc((void **)&devPtrA[0], matrixSizeA * sizeof(devPtrA[0][0]));
if (error != cudaSuccess) fprintf(stderr,"\nError: %s\n",cudaGetErrorString(error));

error = cudaMalloc((void **)&devPtrA_dev, 1 * sizeof(*devPtrA));
if (error != cudaSuccess) fprintf(stderr,"\nError: %s\n",cudaGetErrorString(error));

fprintf(stderr,"malloced device variables\n");

error = cudaMemcpy(devPtrA_dev, devPtrA, 1 * sizeof(*devPtrA), cudaMemcpyHostToDevice);
if (error != cudaSuccess) fprintf(stderr,"\nError: %s\n",cudaGetErrorString(error));

fprintf(stderr,"copied from devPtrA to d_devPtrA\n");

status = cublasSetMatrix(rowsA, colsA, sizeof(a[0]), a, rowsA, devPtrA[0], rowsA);
if (status != CUBLAS_STATUS_SUCCESS) fprintf(stderr,"error %i\n",status);status = cublasDgetrfBatched(handle, m, devPtrA_dev,m,d_pivot_array,d_info_array,1); //cannot get this to work
if (status != CUBLAS_STATUS_SUCCESS) fprintf(stderr,"error in dgetrf %i\n",status);fprintf(stderr,"done with cuda inverse\n");
}

4

c++cuda matrix-inverse

Решение

Другие решения

Других решений пока нет …

Источник

Accepted Answer

Код ошибки 7 означает кублас CUBLAS_STATUS_INVALID_VALUE, Инверсия матриц в кублах возможна только для квадратных матриц, поэтому я предполагаю, что m == n в твоем случае. Это, как говорится, функции cublas<t>getrfBatched требуется, чтобы сводный массив имел длину n для каждой матрицы, поэтому вы должны выделить d_pivot_array как:

error = cudaMalloc((void **)&d_pivot_array, n * sizeof(int));

Чтобы быть более общим, оно выделяется как:

error = cudaMalloc((void **)&d_pivot_array, n * batchSize * sizeof(int));

Вот код инверсии квадратной матрицы, который я написал при тестировании функций CUBLAS Функция ввода и вывода float Тип квадратных матриц размещается на устройстве.

3