Question

Write you code in a way that you can compute any size matrix-vector product.

However, I specifically want you to confirm that the following product works correctly.

e output of the matrix-vector program should look like 2.00 1.00 0.00 0.00 0.00 1.00 3.00 4.00 1.00 2.00 1.00 0.00 0.00 1.00 0.00 1.00 2.00 1.00 0.00 1.00 4.00 4.00 3.00 0.00 0.00 1.00 2.00 1.00 1.00 0.00 e.00 0.00 1.00 2.00 1.00

Required to implement these functions: assignMat, assignVec, printMatVec, allocVec, and allocMat.

As well as these function prototypes:

double* mvp1(double* mat, double* vec, int n);

double* mvp2(double** mat, double* vec, int n);

Cannot take user input. Only has to work for a square matrix, so no need for rows/column variables. Just need variable n

Note one takes the matrix as an int* and the other takes it as an int**

Answer 1

code

solution

#include <stdio.h> //standard library #include <stdlib.h> double* mvp-1 (double* matrix, double* vector, int n); double* mvp_2 (double** matrix,double vector,int n); void free matrix (double **matrix, int row value, int col value) vold assign matrix (double **matrix, int row value, Int col value) void assign_vector (double *vector, int len) v id free vect r (double *vect r, int len); void print matrix_vector (double **matrix, double *vector, double res, int row_value, int col_value) double *all c ve ct r (int len) ; double *alloc matrix (int row value, int col_value) //Main int main (int argc, char *argv) [ int row value-5; int col value 5; double *"matrix! = alloc matrix (row value, col value) *matrix2 double -alloc vector(row value*col value) assign_matrix (matrixl, row_value, col_value); for (int i1-0; il<row value; i1++) for (int j1-0 j1<col_value; j1++) * (matrix2 + il*row value 1 matrix1 [i1] [j1;

double *vectoralloc vector (col value) assign_vector (vector, col _value) double *res1mvp 1 (matrix2, vector, col value) double *res2 mvp_2 (matrixl, vector, col_value); print matrix_vector (matrixl, vector, resl, row _value, col value); printf("n\n") print matrix_vector (matrixl, vector, res2, row _value, col_value); free_ matrix (matrixl, row value, col_value) free vector (vector, col value) free vector (matrix2, row value*col value); ree vector (resl, col value); free vector (res2, col value); return 0; double alloc matrix (int row value, int col value) double **matrix - (double**)malloc (row value sizeof (double*)) for (int 1-0 ; 1くrow value ; 1++) matrix [i] = (double)malloc (col value*sizeof (double)) return matrix; ˇ ld assign matrix (double **matrix, İnt row value, İnt col value) double data;

for (int i-0; i< row_value; i++) for (int j-o; j<col _value; j++) scanf ("Slf", &data) matrix[i]data; v İd free matrix (double **matrix, int r w value, int c 1 value) for (int i-0; i< row_value; i++) free (matrix[i]); free (matrix); double *all c vect r (int ïen1) double *vector (double+)malloc (len1+sizeof (double)) return vector; void assign_vector (double *vectorl, int len1) double haia; scanf ("Slf", &data)

vectorl[ij data; void free_vector (double vectorl, int len1) free (vector1); void print matrix vector (double matrix, double double vector, res, int row_value, int col_value) int vi = 0; for (int i-0 i<row_value; i++) for (int j 0 j<col_value; j++) printf("%.21f", matrix [i] [j]); printf(" " if (vi<col value) printf("%" 21f vi++; vector [vi]); ", if (i--(int) row value/2) printf("-%.21f\n", res[i]); else printf(" %.21f\n", res[i]);

double* mvp_2 (double** matrix, double *vector, int n) double *res = (double*) malloc (n*sizeof (double) ) ; for (int 1-0 ; ¡<n ; i++) res[i] -0; for (int i-o; i<n; i++) for (int j=0; j++) J<n; res [i] += matrix [i] [j] *vector [j]; return res double* mvp_1 (double* matrix, double* vector,int n) double (d uble*)malloc (n*sizeof (double)) ; *res = for (int i-0; i<n; i++) res[i] for (int 1-0 ; ¡<n ; i++) 0; r ( int j=0 ; jen; j++) f res [i] += * (matrix + i"n + j ) *vector [j ] ; return res;

//output

2 1 0 0 o 1 2 1 0 o 0 1 2 1 o 0 0 1 2 1 0 0 0 1 2 2.00 1.00 0.00 0.00 0.001.00 3.00 1.00 2.00 1.00 0.00 0.001.00 4.00 0.00 1.00 2.00 1.00 0.00 1.00 4.00 0.00 0.00 1.00 2.00 1.00 1.00 4.00 0.00 0.00 0.00 1.00 2.00 1.00 3.00 2.00 1.00 0.00 0.00 0.001.00 3.00 1.00 2.00 1.00 0.00 0.001.00 4.00 0.00 1.00 2.00 1.00 0.00 1.00 4.00 0.00 0.00 1.00 2.00 1.00 1.00 4.00 0.00 0.00 0.00 1.00 2.00 1.00 3.00

//copyable code

#include <stdio.h>

//standard library

#include <stdlib.h>

double* mvp_1(double* matrix, double* vector,int n);

double* mvp_2(double** matrix,double*vector,int n);

void free_matrix(double **matrix, int row_value, int col_value);

void assign_matrix(double **matrix, int row_value, int col_value);

void assign_vector(double *vector, int len);

void free_vector(double *vector, int len);

void print_matrix_vector(double **matrix, double *vector, double *res, int row_value, int col_value);

double *alloc_vector(int len);

double **alloc_matrix(int row_value, int col_value);

//Main

int main(int argc, char **argv){

    int row_value=5;

    int col_value=5;

    double **matrix1 = alloc_matrix(row_value, col_value);

    double *matrix2 = alloc_vector(row_value*col_value);

    assign_matrix(matrix1, row_value, col_value);

    for(int i1=0; i1<row_value; i1++)

    {

        for(int j1=0; j1<col_value; j1++)

        {

            *(matrix2 + i1*row_value + j1) = matrix1[i1][j1];

      }

    }

    double *vector = alloc_vector(col_value);

    assign_vector(vector, col_value);

    double *res1 = mvp_1(matrix2, vector, col_value);

    double *res2 = mvp_2(matrix1, vector, col_value);

    print_matrix_vector(matrix1, vector, res1, row_value, col_value);

    printf("\n\n");

    print_matrix_vector(matrix1, vector, res2, row_value, col_value);

    free_matrix(matrix1, row_value, col_value);

    free_vector(vector, col_value);

    free_vector(matrix2, row_value*col_value);

    free_vector(res1, col_value);

    free_vector(res2, col_value);

    return 0;

}

double **alloc_matrix(int row_value, int col_value)

{

    double **matrix = (double**)malloc(row_value*sizeof(double*));

    for(int i=0; i<row_value; i++) matrix[i] = (double*)malloc(col_value*sizeof(double));

    return matrix;

}

void assign_matrix(double **matrix, int row_value, int col_value)

{

    double data;

    for(int i=0; i< row_value; i++)

    {

        for(int j=0; j<col_value; j++)

        {

            scanf("%lf", &data);

            matrix[i][j] = data;

        }

    }

}

void free_matrix(double **matrix, int row_value, int col_value)

{

    for(int i=0; i< row_value; i++)

    {

        free(matrix[i]);

    }

    free(matrix);

}

double *alloc_vector(int len1)

{

    double *vector = (double*)malloc(len1*sizeof(double));

    return vector;

}

void assign_vector(double *vector1, int len1)

{

    double data;

    for(int i=0; i< len1; i++)

    {

        scanf("%lf", &data);

        vector1[i] = data;

    }

}

void free_vector(double *vector1, int len1)

{

    free(vector1);

}

void print_matrix_vector(double **matrix, double *vector, double *res, int row_value, int col_value)

{

    int vi = 0;

    for(int i=0; i<row_value; i++)

    {

        for(int j=0; j<col_value; j++)

        {

            printf("%.2lf ", matrix[i][j]);

        }

        printf(" ");

        if(vi<col_value)

        {

            printf("%.2lf ", vector[vi]);

            vi++;

        }

        if(i==(int)row_value/2)

            printf(" = %.2lf\n", res[i]);

        else printf(" %.2lf\n", res[i]);

    }

}

double* mvp_2(double** matrix, double *vector,int n)

{

    double *res = (double*)malloc(n*sizeof(double));

    for(int i=0; i<n; i++) res[i] = 0;

    for(int i=0; i<n; i++)

    {

        for(int j=0; j<n; j++)

        {

            res[i] += matrix[i][j]*vector[j];

        }

    }

    return res;

}

double* mvp_1(double* matrix, double* vector,int n)

{

    double *res = (double*)malloc(n*sizeof(double));

    for(int i=0; i<n; i++) res[i] = 0;

    for(int i=0; i<n; i++)

    {

        for(int j=0; j<n; j++)

        {

            res[i] += *(matrix + i*n + j)*vector[j];

        }

    }

    return res;

}