部分行列の抜き出し(Ｃ言語)
部分行列の抜き出し(Ｃ言語)
ある行列から、任意の部分行列を抜き出すような関数を作りたいのですが、わからないので手助けをお願いしたいです。
たとえば4行4列の行列Ａがあった場合、その2行2列、2行3列、3行2列、3行3列目の
4つの(2行2列の)要素を抜き出すような関数です。
枠組みは出来ているので載せておきます。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
typedef struct matrix{
double *data;
int size[2];
}MATRIX;
MATRIX create_matrix(int m,int l,int a); //行列の作成
void print_matrix(MATRIX A); //行列の表示
MATRIX submatrix(MATRIX matrix_A,int m,int n,int l,int k); //任意の部分行列の抜き出し
int main(int argc,char **argv)
{
//構造体の初期化
MATRIX matrix_A={NULL,NULL,NULL};
MATRIX matrix_B={NULL,NULL,NULL};
MATRIX matrix_C={NULL,NULL,NULL};
int m=0; //行列の行
int l=0; //行列の列
int a=1; //要素
double det=0; //行列計算
//行列Aの作成
m=4;//行
l=4;//列
matrix_A = create_matrix( m, l, a);
//関数
matrix_C = submatrix(matrix_A,2,2,3,3);
//行列結果出力
printf("\n\tmatrix input");
print_matrix(matrix_A);
printf("\n\tmatrix output");
print_matrix(matrix_C);
//行列の開放
free(matrix_A.data); matrix_A.data=NULL;
free(matrix_C.data); matrix_C.data=NULL;
printf("終了！\n\n");
return 0;
}
MATRIX create_matrix(int m,int l,int a)
{
MATRIX matrix_A={NULL,NULL,NULL};
int i=0;
//サイズの確保 行;size[0] 列;size[1]
matrix_A.size[0]=m;
matrix_A.size[1]=l;
//メモリの確保
matrix_A.data=(double *)malloc(sizeof(double)*matrix_A.size[0]*matrix_A.size[1]);
if(matrix_A.data == NULL){
printf("メモリ確保失敗!![matrix_A]\n");
exit(1);
}
//行列Aの作成
for(i=0; i<((matrix_A.size[0])*(matrix_A.size[1])) ;i++){
*(matrix_A.data+i)=a+i;
}
return(matrix_A);
}
/*行列の表示*/
void print_matrix(MATRIX matrix_A)
{
int i=0;
int j; //行
int k; //列
printf("\n行列の表示\n");
//行列Aのプリント
for(j=0;j<matrix_A.size[0];j++){
for(k=0;k<matrix_A.size[1];k++){
printf("%f\t",*(matrix_A.data+k*matrix_A.size[0]+j));
}
printf("\n");
}
return;
}
MATRIX submatrix(MATRIX matrix_A,int m,int n,int l,int k)
{
MATRIX matrix_C = {NULL,NULL,NULL};
int i,j;
//サイズの確保 行;size[0] 列;size[1]
matrix_C.size[0]=l-m+1;
matrix_C.size[1]=k-n+1;
//メモリの確保
matrix_C.data=(double *)malloc(sizeof(double)*matrix_C.size[0]*matrix_C.size[1]);
if(matrix_C.data == NULL){
printf("メモリ確保失敗!![matrix_C]\n");
exit(1);
}
//任意の部分行列の抜き出し
(ここがわかりません。)
}
よろしくお願いします。
補足
ありがとうございます 0,0,0,・・・,0,0 1,0,0,・・・,0,0 0,1,0,・・・,0,0 0,0,1,・・・,0,0 中略 0,0,0,・・・,1,0 0,0,0,・・・,0,1 のことです ０は横ベクトルです 「最初の１行目が全部0で、２行目からは(I-1)次の単位行列」 の場合 「２行目からは(I-1)次の単位行列」 としたかったのですが言葉について暗黙の約束を期待しているようで躊躇しました