c - 使用双指针创建函数进行矩阵运算

我正在尝试创建一个库，其中包括一些函数，如创建矩阵、执行加法、sub、转置和反转矩阵，我需要使用双指针
一开始，我写这段代码是为了分配矩阵，但它似乎不起作用，我不知道问题出在哪里

#include <stdio.h>
#include <math.h>
#include <stdlib.h>

static double P[4][4]={ { 1,   0,   0,   0},
                        { 0,   1,   0,   0},
                        { 0,   0,   1,   0},
                        { 0,   0,   0,   1}
                      };
double **P_M;
void show_matrix(int n,int m,double **matrix)
{
    int i,j;
    printf("\n The matrix is:\n");
    for (i=0;i<n;i++)
    {
        for (j=0;j<m;j++);
        printf(" \t",&matrix[i][j]);
        printf("\n");
    }
}

double matrix( int n, int m, double **matrix)
{
    int row;
    /*  allocate N 'rows'. */
    matrix = malloc( sizeof( double* ) * n );
    /*  for each row, allocate M actual doubles. */
    for( row = 0; row < n; row++ )
    matrix[ row ] = malloc( sizeof( double ) * m );

}

void main()
{
    int i, j;
    matrix(4,4,P_M);
    for(i=1; i<5; i++)
            for(j=1; j<5; j++)
                P_M[i][j] = P[i-1][j-1];
    //show_matrix(4,4,P_M);

}

最佳答案

很多问题。
超出界限-因为索引是从零开始的。
printf(" \t",&matrix[i][j]);>printf("%lf \t",matrix[i][j]);
double matrix( int n, int m, double **matrix)>double **matrix( int n, int m, double ***matrix)以及函数内部的适当更改（如果需要）。否则就无效。称之为return *martix;
可能还有更多我没注意到的。***指针很愚蠢，不需要将地址传递给指针。

double **matrix(int n, int m)
{
    int row;
    double **array;
    /*  allocate N 'rows'. */
    if (!(array = malloc(sizeof(double*) * n)))
    {
        return NULL;
    }
    /*  for each row, allocate M actual doubles. */
    for (row = 0; row < n; row++)
        if (!(array[row] = malloc(sizeof(double) * m)))
        {
            //do something if malloc failed - for example free already allocated space.
            return NULL;
        }
    return array;
}

总的来说
matrix(4,4,&P_M);

关于c - 使用双指针创建函数进行矩阵运算，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/50315160/