c++ - 从 gsl_matrix 中删除列
问题描述
我想删除 a 的一列gsl_matrix
(第i
-th 列),然后将其内容复制到另一个gsl_matrix
. 这是我的代码:
#include <stdio.h>
#include <iostream>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
int remove_col (int K,
int N,
int i,//between range 1 to N
gsl_matrix *Sn, //Kx(N-1)
gsl_matrix *Z)
{
gsl_matrix_view Z_view;
gsl_matrix_view Sn_view;
if (i==0){
Z_view = gsl_matrix_submatrix (Z, 0, 1, K, N-1);
gsl_matrix_memcpy (Sn, &Z_view.matrix);
}else{
Z_view = gsl_matrix_submatrix (Z, 0, 0, K, i);
Sn_view = gsl_matrix_submatrix (Sn, 0, 0, K, i);
gsl_matrix_memcpy (&Sn_view.matrix, &Z_view.matrix);
Z_view = gsl_matrix_submatrix (Z, 0, i, K, N-i);
Sn_view = gsl_matrix_submatrix (Sn, 0, i-1, K, N-i);
gsl_matrix_memcpy (&Sn_view.matrix, &Z_view.matrix);
}
for ( int row = 0; row < K; ++row ) {
for ( int col = 0; col < N-1; ++col ) {
printf( "\t%3.1f", gsl_matrix_get( Sn , row, col));
}
printf( "\n" );
};
}
int main () {
double Zx[]={1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1,
0, 0, 1, 1, 1, 1,
0, 1, 0, 1, 1, 1,
1, 1, 0, 1, 1, 0,
1, 1, 1, 1, 1, 1,
1, 1, 0, 1, 0, 0,
1, 1, 1, 0, 1, 0,
1, 1, 0, 1, 0, 0,
0, 1, 1, 1, 0, 1};
int r=10;
int c=6;
gsl_matrix_view Z = gsl_matrix_view_array (Zx, r, c);
gsl_matrix *R= gsl_matrix_calloc(c,r);
for ( int row = 0; row < r; ++row ) {
for ( int col = 0; col < c; ++col ) {
gsl_matrix_set( R, col, row, gsl_matrix_get( &Z.matrix , row, col));
}
}
gsl_matrix *Sn = gsl_matrix_alloc(c,r-1);
remove_col (c,r,1, Sn, R);
for (int row=0; row<c; ++row)
for (int col=0; col<r-1; ++col)
printf(col==r-2?"%6.3f\n":"%6.3f ",gsl_matrix_get(Sn,row,col));
gsl_matrix_free(R);
gsl_matrix_free(Sn);
}
return 0;
}
我收到此错误:
gsl: submatrix_source.c:37: ERROR: first dimension overflows matrix
Default GSL error handler invoked.
Aborted (core dumped)
我的方法正确吗?有什么建议吗?
解决方案
此版本有效:
#include <stdio.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
int remove_col (int K,
int N,
int i,//between range 1 to N
gsl_matrix *Sn, //Kx(N-1)
gsl_matrix *Z) {
gsl_matrix_view Z_view;
gsl_matrix_view Sn_view;
if (i==0){
Z_view = gsl_matrix_submatrix (Z, 0, 1, K, N-1);
gsl_matrix_memcpy (Sn, &Z_view.matrix);
} else{
Z_view = gsl_matrix_submatrix (Z, 0, 0, K, i);
Sn_view = gsl_matrix_submatrix (Sn, 0, 0, K, i);
gsl_matrix_memcpy (&Sn_view.matrix, &Z_view.matrix);
Z_view = gsl_matrix_submatrix (Z, 0, i, K, N-i);
Sn_view = gsl_matrix_submatrix (Sn, 0, i-1, K, N-i);
gsl_matrix_memcpy (&Sn_view.matrix, &Z_view.matrix);
}
}
int main () {
double Zx[]={1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1,
0, 0, 1, 1, 1, 1,
0, 1, 0, 1, 1, 1,
1, 1, 0, 1, 1, 0,
1, 1, 1, 1, 1, 1,
1, 1, 0, 1, 0, 0,
1, 1, 1, 0, 1, 0,
1, 1, 0, 1, 0, 0,
0, 1, 1, 1, 0, 1};
int r=10;
int c=6;
gsl_matrix_view Z = gsl_matrix_view_array (Zx, r, c);
gsl_matrix *R= gsl_matrix_calloc(c,r);
for ( int row = 0; row < r; ++row ) {
for ( int col = 0; col < c; ++col ) {
gsl_matrix_set( R, col, row, gsl_matrix_get( &Z.matrix , row, col));
}
}
for (int row=0; row<r; ++row){
for (int col=0; col<c; ++col)
printf(col==r-2?"%6.3f\n":"%6.3f ",gsl_matrix_get(R,col,row));
printf("\n");
}
gsl_matrix *Sn = gsl_matrix_alloc(c,r-1);
remove_col (c,r,1, Sn, R);
for (int row=0; row<r-1; ++row){
for (int col=0; col<c; ++col)
printf(col==r-2?"%6.3f\n":"%6.3f ",gsl_matrix_get(Sn,col,row));
printf("\n");
}
gsl_matrix_free(R);
gsl_matrix_free(Sn);
return 0;
}
推荐阅读
- python - Python比较不同的时间值以获得以分钟为单位的时间增量
- python - 如何计算我输入的一系列数字的中位数?
- c# - Azure blobstorage sas 密钥在 Angular 应用程序中不安全
- sql - 对 RDBMS (SQL) 数据库进行分片
- xamarin.forms - 通过 USB 对笔记本电脑不可见的外部存储文件
- r - 在列表上应用应用函数
- database-design - Mysql 数据库设计销售报告
- asp.net - 如何生成包含登录用户 ID 的连接字符串
- c - 使用宏时出现 gcc 错误:预期为 ',' 或 ';'
- sql - ForEach ADO 枚举器中的 ssis 截断表