동적 인 2 차원 배열을 정적 배열

Question

어떤 몸이라도 2d 동적 배열을 정적으로 변환하여 lapacke에서 사용할 수 있습니다. dgels 함수는 C에서 정적 행렬 만 가져옵니다. malloc을 사용할 때 정답을주지 않습니다. 이 it.thankyou와 함께 작동하도록하는 방법이 malloc을 사용할 수 있습니다

#include <stdio.h> 
#include <lapacke.h> 
#include <conio.h> 

int main (int argc, const char * argv[]) 
{ 
    /*double a[5][3] = {1,1,1,2,3,4,3,5,2,4,2,5,5,4,3};*/ 

    double b[5][2] = {-10,-3,12,14,14,12,16,16,18,16}; 
    lapack_int info,m,n,lda,ldb,nrhs; 
    int i,j; 
    double **a; 

    a=(double**)malloc(5* sizeof(double*)); 
    for (i=0;i<5;i++) 
    { 
     a[i]=(double*)malloc(3* sizeof(double)); 
    } 
    a[0][0]=1; 
    a[0][1]=1; 
    a[0][2]=1; 
    a[1][0]=2; 
    a[1][1]=3; 
    a[1][2]=4; 
    a[2][0]=3; 
    a[2][1]=5; 
    a[2][2]=2; 
    a[3][0]=4; 
    a[3][1]=2; 
    a[3][2]=5; 
    a[4][0]=5; 
    a[4][1]=4; 
    a[4][2]=3; 

    m = 5; 
    n = 3; 
    nrhs = 2; 
    lda = 3; 
    ldb = 2; 

    info = LAPACKE_dgels(LAPACK_ROW_MAJOR,'N',m,n,nrhs,*a,lda,*b,ldb); 

    for(i=0;i<n;i++) 
    { 
     for(j=0;j<nrhs;j++) 
     { 
      printf("%lf ",b[i][j]); 
     } 
     printf("\n"); 
    } 
    getch(); 
    return(info); 
} 

Answer 1

a가 아니다 2d 배열은 1d 배열을 분리하기위한 포인터의 배열입니다. *a을 LAPACKE_dgels에 전달하면 첫 번째 행에 대한 포인터 만 제공됩니다. 독립적으로 할당 된 이후에 다른 모든 행이 할당 된 위치를 알 수있는 방법이 없습니다. 전체 배열이 하나의 인접한 메모리 블록에 있어야합니다. a은 double**이 아닌 double*이어야하며 전달할 때 참조를 취소하지 마십시오. 행 또는 열의 주요 형식 (함수에 전달한 것)을 사용하여 2d 색인을 직접 1d 색인으로 전개해야합니다.

EDIT

다음 코드 m*ndouble의 방 평평한 1D 배열을 할당한다. 그런 다음 수식 row * n + col을 사용하여 2 차원 인덱스를 1 차원 행 - 주요 인덱스로 변환하여 배열을 채 웁니다. 열 주요 인덱스가 필요하면 col * m + row을 사용합니다.

#include <stdio.h> 
#include <lapacke.h> 
#include <conio.h> 

int main (int argc, const char * argv[]) 
{ 
    double b[5][2] = {-10,-3,12,14,14,12,16,16,18,16}; 
    lapack_int info,m,n,lda,ldb,nrhs; 
    int i,j; 
    double *a; 

    m = 5; 
    n = 3; 
    nrhs = 2; 
    lda = 3; 
    ldb = 2; 

    a = malloc(m * n * sizeof(double)); 
    a[0 * n + 0] = 1; 
    a[0 * n + 1] = 1; 
    a[0 * n + 2] = 1; 
    a[1 * n + 0] = 2; 
    a[1 * n + 1] = 3; 
    a[1 * n + 2] = 4; 
    a[2 * n + 0] = 3; 
    a[2 * n + 1] = 5; 
    a[2 * n + 2] = 2; 
    a[3 * n + 0] = 4; 
    a[3 * n + 1] = 2; 
    a[3 * n + 2] = 5; 
    a[4 * n + 0] = 5; 
    a[4 * n + 1] = 4; 
    a[4 * n + 2] = 3; 

    info = LAPACKE_dgels(LAPACK_ROW_MAJOR,'N',m,n,nrhs,a,lda,*b,ldb); 

    for(i=0;i<n;i++) 
    { 
     for(j=0;j<nrhs;j++) 
     { 
      printf("%lf ",b[i][j]); 
     } 
     printf("\n"); 
    } 
    getch(); 
    return(info); 
} 

Answer 2

나는 lapacke.dgels을 알고 있지만 변경하려고하지 않습니다

double **a; 
a=(double**)malloc(5* sizeof(double*)); 
for (i=0;i<5;i++) 
{ 
    a[i]=(double*)malloc(3* sizeof(double)); 
} 

에 :

double (*a)[3]; 
a = malloc(5 * 3 * sizeof(double));