#include #include #include void initMatrix(size_t m, size_t n, double *A, ptrdiff_t incRow, ptrdiff_t incCol) { for (size_t i=0; i A\n"); dger(m, n, 1.5, y, incY, x, incX, A, incRowA, incColA); printMatrix(m, n, A, incRow, incCol); // 2) dger printf("1.5*x * y^T + A^T -> A^T\n"); dger(n, m, 1.5, x, incX, y, incY, A, incColA, incRowA); printMatrix(m, n, A, incRow, incCol); // Deallocate matrix A_ free(A_); }