C++ Динамические матрицы Умножая возвращаемый выпуск

Question

Я пытаюсь умножить на n ** x ** n динамические матрицы и вернуть результат. Это код для него:

long long int** Multiply(long long int** m1, long long int **m2) 
{ 
    static long long int** output; 
    output= new long long int* [k]; 

    for (int i=0; i<k; i++) 
     output[k]= new long long int [k]; 

    long long int cellRes= 0; 
    for (int i=0; i<k; i++) 
    { 
     for (int f=0; f<k; f++) 
     { 
      for (int j=0; j<k; j++) 
      { 
       cellRes+= m1[i][j]*m2[j][f]; 
      } 
      output[i][f]= cellRes; 
      cellRes=0; 
     } 
    } 

    return output; 
} 

Это, кажется, делает это хорошо при возвращении новой матрицы, но по какой-то причине программа сбой после выполнения Multiply() ... Даже если я создаю новый матрицы, выделите его и затем назначьте ему Multiply(), он сработает.

Я не могу понять, что я делаю неправильно. Есть идеи?

Answer 1

for (int i=0; i<k; i++) 
    output[k]= new long long int [k]; 

должно быть:

for (int i=0; i<k; i++) 
    output[i]= new long long int [k]; 
     //^^^^ output[k] is out of bound and you allocate space for each row