Найти определитель матрицы на графическом процессоер

[Nik_Nez]

Необходимо переписать код для нахождения определителя матрицы на CPU,
в код для GPU на CUDA.
Код для нахождения матрицы есть, а вот как его переписать на CUDA не могу разобраться. Помогите пожалуйста кто чем может.

PHP код:

#include "stdafx.h"
#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#include <time.h>
typedef double *A;
int main()
{     int  i,a=0,j,n=4,k, matrix_size= 42;
    //A *mas;
//double *mas= new double;
    double mas[500][500];//={{0.0,5.0,6.0,2.0},{9.0,0.0,7.0,3.0,},{8.0,1.0,0.0,5.0},{7.0,9.0,9.0,0.0}};
    double b,det=1;
    clock_t start, finish;
double duration;

printf("razmer massiva=");
scanf_s("%d",&n);


for(i=0;i<n;i++){
for(j=0;j<n;j++)
{mas[i][j]=rand()%9 + 1
;
}}

start = clock();



for (j = 0; j < n-1; j++){
for(i=j+1;i<n;i++){

if(mas[i][j]==0)
continue;
b=mas[i][j]/mas[j][j];
for(int r=0;r<n;r++){
mas[i][r]=mas[i][r]-(mas[j][r])*b;}}
for(i=0;i<n;i++){
if(mas[i][i]==0){
int q=i;
for(i=0;i<n;i++)
{
if
    (mas[i][q]!=0)
{k=i;
a++;
break;}
}
for(j=0;j<n;j++)
{b=mas[q][j];
mas[q][j]=mas[k][j];
mas[k][j]=b;}
}} 
}
for (i = 0; i < n; i++){
det*=mas[i][i];
}
a=a%2;
if(a!=0)
{det=det*(-1);}
finish = clock();
for(i=0;i<n;i++){
for (j=0;j<n;j++){
printf_s("%5.3f",mas[i][j]);
if (j==n-1)
printf("\n");}}
duration=(finish-start);
printf("\nopredelitel=%5.0f\n",det);
printf("t=%2.2f\n",duration);

_getch();
}