は行列から特定の列を抽出し、Cの新しい行列として格納します

Question

行列は動的配列、つまりdouble * inputMatとして格納されています。私は、特定の列を抽出できる行と列の数を知っています。今問題は、私は抽出され、別の動的配列に格納される列のセットを持っているということです。これを行う方法？

double *extractMatrix(double *inputMat,int rows, int *columnIndex, int columnTotal) 
{  
    double *outputMat=malloc(sizeof(double)*rows*columnTotal); 
    for(int i=0; i<columnTotal; i++) 
     memcpy(outputMat, &inputMat[rows*columnIndex[i]],rows*sizeof(double)); 
    return outputMat; 
} 

columnIndexには、行列から抽出する列のインデックスが含まれています。 ColumnTotalはcolumnIndex配列のサイズです。しかし、これはinputMatの特定の1つの列だけをoutputMatにコピーし、おそらく上書きされます。私はcolumnIndexのすべての列の完全な配列を必要とします。私はlapackとBLASライブラリを扱っています。これを行うための組み込み方法がある場合は、共有してください。

Answer 1

それは2次元配列であるかのようにあなたの基本的な目標は、インデックスにダブル（を倍増するために、実際にポインタ）のご配列です。これを行うには、特定の列columnIndexを抽出し、その列を構成する値を保持するメモリブロック（rowsの値）を動的に割り当てて、新しく割り当てられたブロックへのポインタを返します。

あなたのアプローチは正しい方向にあります。インデックスはちょうどオフです。 forループ自体でインデックスを処理する方がはるかに簡単です。基本的なアプローチは次のとおりです。

int n = 0; 
for (int i = colindex; i < rows * cols; i += cols) 
    output[n++] = input[i]; 

（つまり、列内n < INT_MAX値をとる - 必要に応じて調整）

は、あなたがやろうとしているように見えるものを、あなたが行うことができない小さな例をまとめます次のようなもの：

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

int *getcol (const int *a, const int rows, const int col, const int ncols) 
{ 
    int sz = rows * ncols, 
     n = 0, /* output index */ 
     *out = malloc (rows * sizeof *a); 

    if (out) 
     for (int i = col; i < sz; i += ncols) /* index with loop vars */ 
      out[n++] = a[i];     /* assign column values */ 

    return out; 
} 

int main (void) { 

    int a[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9 }, 
     rows = 3,     /* simulate a 3x3 with 9 values */ 
     cols = sizeof a/(rows * sizeof *a), 
     colidx = 1,     /* zero based index of wanted column */ 
     *output = NULL; 

    if ((output = getcol (a, rows, colidx, cols))) { 
     for (int i = 0; i < rows; i++) 
      printf (" %2d\n", output[i]); 
     free (output); /* don't forget to free memory */ 
    } 

    return 0; 
} 

（注：関数のパラメータは、あなたが列挙されている順序になっている - しかし、短い名前を使用して好ましくはあなたがように順番を入れ替えますが最後ですが、それはあなた次第です。あなたが望む列はポインタとして渡される必要はなく、単純割り当てが動作するmemcpyの必要はありません。また、あなたがしたい列のインデックスはゼロベースのインデックス）

使用例/出力として渡される）9つの値からシミュレートされた3×3アレイの2列目にある

$ ./bin/array_idx_1d_as_2d 
    2 
    5 
    8 

（

ご不明な点がありましたらお気軽にお問い合わせください。