C程序接收信号SIGSEGV,矩阵乘法中的分段错误?

时间:2017-05-22 17:58:02

标签: c matrix segmentation-fault malloc sigsegv

是的,我已经阅读了这样的问题,但是在我的问题中没有一个被证明是有用的。我的程序在另一个函数中运行完美的函数存在问题。 / p>

  

编程接收信号SIGSEGV,分段故障。 malloc.c中的_int_malloc(av = 0x7ffff7dd8e40,bytes = 32):4703 4703 malloc.c:没有这样的文件或目录。

这是我收到的错误。

这是我的计划:

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

void endl()
{
  printf("\n\n");
}

void loadVector(FILE *infile, double *vector, int col)
{
  int i;
  for(i=0;i<col;i++)
  {
    fscanf(infile,"%lf",&vector[i]);
  }
  printf("vector read...\n");
}

void printVec(double *vect, int run)
{
  int k;
  for(k=0;k<run;k++)
  {
   printf("%.2f ",vect[k]);
  }
}

double* duplic(double* matr, double *vect, double *sol,  int col, int  row)
{
  int k, l;
  for(k=0;k<row;k++)
  {
    for(l=0;l<col;l++)
    {
      sol[k] += matr[col*k+l] * vect[l];
    }
  }
  return sol;
}


int main(int argc, char **argv)
{
  endl();
  printf("=== ===");
  endl();

  if(argc!=3)
  {
    perror("not enough arguments");
    endl();
    return 1;
  }

  char ins[300];
  char *cut;
  int row=0, col=0, rnr = 1, index = 0;
  double temp;
  double *mat = NULL;
  double *vec = NULL;
  double *res = NULL;

  FILE *in;
  in = fopen(argv[1], "r");

  while(fgets(ins,sizeof(ins),in) != NULL)
  {
      col = 0;
      cut = strtok(ins," ");
      while(cut != NULL)
      {
        col++;
        sscanf(cut,"%lf",&temp);
        if(mat==NULL)
         {
            mat = malloc(sizeof(temp));
            *mat = temp;
         }
         else
         {
            rnr ++;
            mat = realloc(mat,sizeof(mat)*rnr);
            index = rnr - 1;
            *(mat+index)=temp;
         }
        cut = strtok(NULL," ");
    }
    row ++;
 }

fclose(in);

printf("Matrix read...");
endl();

printf("Printing matrix: ");
endl();

int i, j;
for(i=0;i<row;i++)
{
   for(j=0;j<col;j++)
   {
      printf("%.2f ",mat[col*i+j]);
   }
   printf("\n");
}
endl();

// printf("rows: %i\tcols: %i\n",row,col);

vec = malloc(sizeof(col));

FILE *inv;
inv = fopen(argv[2],"r");
loadVector(inv,vec,col);

endl();
printVec(vec,col);
endl();

res = malloc(sizeof(row)); 
res = duplic(mat,vec,res,col,row);

printf("Solution-vector calculated...");
endl();
printf("\nSolution vector:");
endl();

printVec(res,row);

endl();

free(mat);
free(vec);
free(res);
fclose(inv);

printf("Exitting\n");
return 0;
}

该程序是一个矩阵 - 向量乘法,它读取一个未知大小的矩阵和一个匹配的向量(,如果你能提供一个更好的方法来读取一个更好的代码的随机长度行,我'我很高兴)。 问题来自乘法代码 - 虽然它应该正常工作,因为它是我写的另一个程序,其中大小是修复。 所以它很好地读取矩阵和向量,然后给出段错误。

有什么想法吗?

1 个答案:

答案 0 :(得分:2)

您编写sizeof(col),但col是一个整数,这意味着您的vec数组将始终分配sizeof(int)个字节的内存(通常为4个字节但会取决于你的架构)。如果你想要一个col数为double的双数组,你想使用malloc(sizeof(double) * col)来分配正确的字节数。

我怀疑您的res数组发生了同样的问题。它可能在某些情况下工作,在这种情况下,你分配的空间之后的空间仍可用于写入,在这种情况下计算机不会抱怨,因为malloc的内存都在堆中。但是,如果碰巧在malloc&gt; ed块(太小)之后尝试写入受保护或保留的内存部分,操作系统将抛出SegFault。