паралелльное определение детерминанта

float det(float* *m,int n)
{
if(n==1)
return m[0][0];
float s=0;
for(int i=0; i<n; i++)
{
float k=m[0][i]*det(minor(i),n-1);
if(i%2==0)
s+=k;
else s-=k;
}
return s;
}
 
float* *minor(float* *m,int n, int i)//Возвращает матрицу-минор к элементу (0,i)
{
float* *min=new float*[n-1];
for(int j=0; j<n-1; j++)
{
min[j]=new float[n-1];
for(int k=0; k<n-1; k++)
{
if(k<i)
min[j][k]=m[j+1][k];
else if(k>i)
min[j][k]=m[j+1][k+1];
}
}
return min;
}

float* *l=new float*[n];
float* *u=new float*[n];
for(int i=0; i<n; i++)
{
    u[i]=new float[n];
    l[i]=new float[n];
    for(int j=0; j<n; j++)
    {
        u[i][j]=0;
        l[i][j]=0;
    }
}
for(int i=0; i<n; i++)
{
    u[i][i]=1;
    for(int j=i; j<n; j++)
    {
        l[j][i]=a[j][i];
        for(int k=0; k<i; k++)
            l[j][i]=l[j][i]-l[j][k]*u[k][i];
    }
    for(int j=i+1; j<n; j++)
    {
        u[i][j]=a[i][j];
        for(int k=0; k<i; k++)
            u[i][j]=u[i][j]-l[i][k]*u[k][j];
        u[i][j]=u[i][j]/l[i][i];
    }

float det(float* *m1, int n, float tolerance)//tolerance -- точность
{
if(n==1)
return m1[0][0];
float* *m=new float*[n];
for(int i=0; i<n; i++)
{
m[i]=new float[n];
for(int j=0; j<n; j++)
m[i][j]=m1[i][j];
}//Копируем матрицу, чтобы не изменять то, что мы передавали функцию
if(abs(m[0][0])>tolerance)//Не 0 на главной диаг.
{
for(int i=1; i<n; i++)
{
float k=m[i][0]/m[0][0];
for(int j=0; j<n; j++)
m[i][j]-=k*m[0][j];
}
float* *m2=new float*[n-1];
for(int i=0; i<n-1; i++)
{
m2[i]=new float[n-1];
for(int j=0; j<n-1; j++)
m2[i][j]=m[i+1][j+1];
}
float s=m[0][0]*det(m2,n-1,tolerance);
delete m;
delete m2;
return s;
}
else//0 на главной
{
int flag=false;//Есть ли не 0 в первом столбце
for(int i=1; i<n; i++)
if(abs(m[i][0])>tolerance)//Нашли такую строку
{
if(i%2==0)
flag=-1;//Она четная
else
flag=1;
float tmp;//буффер
for(int j=0; j<n; j++)
{
tmp=m[0][j];
m[0][j]=m[i][j];
m[i][j]=tmp;
}//Поставили ее наверх -- теперь m[0][0]!=0
break;
}
if(flag==0)
return 0;
return flag*det(m,n,tolerance);
}
}