非平方矩阵乘法帮助C++

首先，很抱歉代码有点乱，因为我在处理不同的事情以使其正常工作。到目前为止，我的代码可以很好地乘以平方矩阵；然而，它在计算非平方矩阵方面存在困难。调试后，我的最佳猜测是如何重新调整向量的大小，并且存在导致程序崩溃的越界错误。如果有任何帮助，我们将不胜感激，我的代码应该能够在矩阵乘法规则中乘以任何向量大小。

我还想注意的是，这是一个硬件任务，所以我只限于如何构建我的代码，基本上只使用向量，不能编写自己的类等。…

#include <iostream>
#include <vector>
using namespace std;
void multiply_matrices(vector <vector<int> > matrix1,vector <vector<int> > matrix2, int cols, int rows2); 
void setMatrix(vector <vector<int> > &matrix, int rows, int cols);
int main()
{
    int rows, cols, rows2, cols2;    
    vector< vector<int> > matrix, matrix2;        
    cout<<"Please enter the number of Rows and Columns for your first Matrix."<<endl;
    cout<<"Rows: ";
    cin>>rows;
    cout<<"Columns: ";
    cin>>cols;
    matrix.resize(cols, vector<int>(rows,0));  //Saw this online so not sure how it works but it works, if i take out one i cant do row<column and vice versa
    matrix.resize(rows, vector<int>(cols,0));
    cout<<"Size has been declared, please enter data for your matrix"<<endl;
    setMatrix(matrix,rows,cols);
    cout<<"Second Matrix Automatically Set by Matrix Multiplication Rule"<<endl; //Just automatically sets second matrix as per Matrix Multiplication Rule
    rows2=cols;
    cols2=rows;
    cout<<"Second Matrix Size is: " << rows2 << " by " << cols2 << endl;
    matrix2.resize(cols2, vector<int>(rows2,0));
    matrix2.resize(rows2, vector<int>(cols2,0));
    setMatrix(matrix2,rows2,cols2);        
    cout<<"Multiplied Matrix is:"<<endl;
    multiply_matrices(matrix,matrix2,cols,rows2);
    system("PAUSE");
    return 0;
}
void setMatrix(vector <vector<int> > &matrix, int rows,int cols){
     int num;
     for(int i = 0; i < rows; i ++)
     {
        for (int j = 0; j < cols; j++)
        {
            cout << "Enter Value for Row " << (i+1) << " Column " << (j+1) << ": ";
            cin>>num;
            matrix[i][j]=num;            
        }        
        cout <<  endl;
    }
 /*for(int i = 0; i < rows; i ++)
    {
        for (int j = 0; j < cols; j++)
        {
            cout << matrix[i][j] << " ";
        }       
        cout <<  endl;
    }          
  */   
     }
void multiply_matrices(vector <vector<int> > matrix1,vector <vector<int> > matrix2, int cols, int rows2){
    vector< vector<int> > tempMatrix;
    int newrows=rows2;
    int newcols=cols;
    int sum;
    tempMatrix.resize(newcols, vector<int>(newrows,0));   //Resizing new matrix to proper size, so if it was (2x3)(3x2), new matrix is (3x3)
    for (int i = 0; i < newrows; i++)                    //This Works Fine for Square Matrixes but not for others, i have no clue how to fix it?
    {
        for (int j = 0; j < newcols; j++){
            //sum=0;    
            for (int u = 0; u < newcols; u++)
            {
                //sum+=matrix1[i][u] * matrix2[u][j];
                //tempMatrix[i][j]=sum;
                tempMatrix[i][j] += matrix1[i][u] * matrix2[u][j];
            }
        }
    }
    for(int i = 0; i < newrows; i ++)
    {
        for (int j = 0; j < newcols; j++)
        {
            cout << tempMatrix[i][j] << " ";
        }        
        cout <<  endl;
    }          
}