[交流] 【求助】如何调用此函数分解具体矩阵？ 已有1人参与

int choldc(float **a, int n, float p[])
{
        int i, j, k;
        float sum = 0.0F;
        for (i=1; i<=n; i++)
        {
                for (j=i; j<=n; j++)
                {
                        for (sum=a[i][j], k=i-1; k>=1; k--)
                        {
                                sum -= a[i][k]*a[j][k];
                        }
                        if (i == j)
                        {
                                if (sum <= 0.0)
                                {
                                        return 0;
                                }
                                p[i]=sqrt(sum);
                        }
                        else
                        {
                                a[j][i]=sum/p[i];
                        }
                }
        }
        return 1;
}