Скалярное произведение

#include <iostream>
#include <vector>
 
using namespace std;
using matrix = vector<vector<double>>;
 
matrix add(matrix A, matrix B) {
    matrix C;
    C.resize(A.size());
    for (int i = 0; i< A.size(); i++) {
        C[i].resize(B.size());
        for (int j = 0; j < B.size(); j++) {
            C[i][j] = A[i][j] + B[i][j];
        }
    }
    return C;
}
 
matrix multiple(matrix A, matrix B)
{
    matrix C;
    C.reserve(100);
    C.resize(B.size());
    for (int i = 0; i < A.size(); i++) {
        C[i].resize(B.size());
        for (int j = 0; j < B.size(); j++) {
            for (int k = 0; k < B.size(); k++)
                C[i][j] += A[i][k] * B[k][j];
        }
    }
    return C;
}
 
void main() {
    matrix A = { {1,2,3}, {1,2,1}, {3,2,0} };
    matrix B = { {4,1,2},{0,4,3},{1,1,1} };
    matrix r = { {-0.7f, 1.3, 0.2} };
    matrix q = { { -1.6f, 0.8, 1.1} };
    matrix p = { {0.1, 1.7, -1.5} };
 
 
    matrix r_q = add(r, q);
 
    for (int i = 0; i < r_q.size(); i++) {
        for (int j = 0; j < r_q.size(); j++) {
            cout << r_q[i][j] << "\t";
        }
        cout << "\n";
    }
 
    matrix a_a = multiple(A, A);
    matrix a_a_p = multiple(a_a,p);
 
    getchar();
}