分段错误，但数组对象未超出范围 (C++）

关于分段错误的常见原因的帖子，但我认为我在这里创建的内置数组对象（result）在我为其赋值时不会越界。
我认为这可能会对将来拥有数组不越界的人有所帮助，而且我也没有看到很多关于制作 2D 内置数组对象的东西 - 我见过的例子几乎完全是向量或 std：array 对象。

以下是可运行的相关代码：

矩阵.h

#ifndef MATRIX_H
#define MATRIX_H
#include <initializer_list>
using std::initializer_list;
typedef unsigned int uint;
class Matrix {
  public:
    Matrix(uint rows, uint cols);  
    ~Matrix();  
    Matrix add(double s) const;  
    const uint numRows() const;  
    const uint numCols() const;  
    double & at(uint row, uint col);  
    const double & at(uint row, uint col) const;  
  private:
    uint rows, cols;
    double ** matrix;
    void makeArray() {
      matrix = new double * [rows];
      for(uint i = 0; i < rows; ++i) {
        matrix[i] = new double [cols];
      }
    }
};
#endif  

矩阵.cpp

#include "matrix.h"
Matrix::Matrix(uint rows, uint cols) {
  //Make matrix of desired size
  this->rows = rows;
  this->cols = cols;
  makeArray();
  //Initialize all elements to 0
  for(uint i = 0; i < rows; ++i) {
    for(uint j = 0; j < cols; ++j) {
      this->matrix[i][j] = 0.0;
    }
  }
}
Matrix::~Matrix() {
  for(uint i = 0; i < numRows(); ++i) {
    delete[] matrix[i];
  }
  delete[] matrix;
}
const uint Matrix::numRows() const {
  return this->rows;
}
const uint Matrix::numCols() const {
  return this->cols;
}
double & Matrix::at(uint row, uint col) {
  return matrix[row][col];
}
const double & Matrix::at(uint row, uint col) const {
  return matrix[row][col];
}
Matrix Matrix::add(double s) const {
  uint r = this->numRows();
  uint c = this->numCols();
  Matrix * result;
  result = new Matrix(r, c);
  for(uint i = 0; i < r; ++i) {
    for(uint j = 0; j < c; ++j) {
      result->at(i,j) = (this->at(i,j)) + s;
    }
  }
  return * result;
}  

主.cpp

#include <iostream>
#include <cstdlib>
#include "matrix.h"
using namespace std;
typedef unsigned int uint;
int main() {
  Matrix * matrix;
  matrix = new Matrix(3, 2); //Works fine
  double scaler = 5;
  matrix->at(2,1) = 5.0; //Works fine
  Matrix r = matrix->add(scaler); //DOESN'T WORK
  return EXIT_SUCCESS;
}  

知道为什么add函数会导致分段错误吗？我用来填充结果矩阵对象的 for 循环并没有越界，而且我对C++不够熟悉，不知道还有什么可能导致它。
提前谢谢。