2

我想将矩阵和向量相乘，我已经编写了将矩阵和向量存储在 malloc 数组中的函数。对于这个函数，我需要先使用 malloc 创建另一个数组来存储我的答案。然后进行计算（http://www.facstaff.bucknell.edu/mastascu/elessonsHTML/Circuit/MatVecMultiply.htm）

#include <stdlib.h>
/* Multiply a matrix by a vector. */
double *matrix_vector_multiply(int rows, int cols,
               double **mat, double *vec){

 // creating an vecotr to hold the answer first

 double *ans= malloc(rows* sizeof(double));

 // do the multiplication:
 mulitply **mat and *vec (mat = the matrix and vec is the vector)

 for (rows=0; rows< ; rows++)
    for (cols=0; cols< ; cols++)
        ans[rows] = ans[rows] + vec[rows][cols] * mat[rows];    

 //not sure if it is right

 // then store the answer back to the ans array            


}

主要功能：

double *matrix_vector_multiply(int rows, int cols,
               double **mat, double *vec);

int main(){
  double *answer = matrix_vector_multiply(rows, cols, matrix, vector);

  printf("Answer vector: \n");
  print_vector(rows, answer);
  return 0;
 }

不知道如何用指针进行这种乘法然后将其存储回来..任何帮助将不胜感激！谢谢！

编辑：乘法功能：

#include <stdlib.h>
/* Multiply a matrix by a vector. */
double *matrix_vector_multiply(int rows, int cols,
               double **mat, double *vec){

double *ans = malloc(rows * sizeof (double));   
int i;  
for (i=0; i<rows; rows++)
    for (i=0; i<cols; cols++)
        ans[rows] = ans[rows] + vec[rows][cols] * mat[rows];    

return ans;

}

但我在第 12 行遇到错误，下标值既不是数组也不是指针

4

1 回答 1

2

您的功能有几个错误：

  1. 在 for 循环中迭代矩阵的变量不应该是您传递给函数的参数。尝试这样的事情： for(y=0;y<rows;y++)

  2. 你必须交换vec并mat在 for 循环中

  3. 您必须将答案向量初始化为0（另一个 for 循环）

  4. return ans;您必须在乘法结束时返回答案 ( )

希望有帮助，扬

于 2012-10-23T06:36:31.127 回答