Каждый столбец матрицы A сложить со столбцом k

#include <iostream>
#include <iomanip>
 
template<typename T>
class Matrix {
public:
    Matrix() : height(0), width(0), data(nullptr) {}
    Matrix(size_t height, size_t width) : height(height), width(width), data(new T[height * width]{}) {}
    Matrix(const Matrix<T> &other) : height(other.height), width(other.width), data(new T[height * width]) {
        for (size_t i = 0, size = height * width; i < size; ++i) {
            data[i] = other.data[i];
        }
    }
    Matrix(Matrix<T> &&other) noexcept : Matrix() {
        swap(*this, other);
    }
    Matrix<T> &operator=(Matrix<T> other) noexcept {
        swap(*this, other);
        return *this;
    }
    ~Matrix() {
        delete [] data;
    }
    friend void swap(Matrix<T> &a, Matrix<T> &b) {
        using std::swap;
        swap(a.height, b.height);
        swap(a.width, b.width);
        swap(a.data, b.data);
    }
 
    size_t getHeight() const {
        return height;
    }
 
    size_t getWidth() const {
        return width;
    }
 
    T &operator()(size_t row, size_t column) {
        return data[row * width + column];
    }
 
    const T &operator()(size_t row, size_t column) const {
        return data[row * width + column];
    }
private:
    size_t height;
    size_t width;
    T *data;
};
 
template<typename T>
std::istream &operator>>(std::istream &in, Matrix<T> &matrix) {
    for (size_t i = 0, height = matrix.getHeight(); i < height; ++i) {
        for (size_t j = 0, width = matrix.getWidth(); j < width; ++j) {
            in >> matrix(i, j);
        }
    }
    return in;
}
 
template<typename T>
std::ostream &operator<<(std::ostream &out, const Matrix<T> &matrix) {
    for (size_t i = 0, height = matrix.getHeight(); i < height; ++i) {
        for (size_t j = 0, width = matrix.getWidth(); j < width; ++j) {
            out << std::setw(7) << std::setprecision(2) << std::fixed << matrix(i, j);
        }
        out << "\n";
    }
    return out;
}
 
int main() {
    size_t n, m;
    std::cin >> n >> m;
 
    Matrix<double> a(n, m);
    std::cin >> a;
 
    std::cout << a << std::endl;
 
    size_t k;
    std::cin >> k;
    Matrix<double> b = a;
 
    for (size_t i = 0, height = a.getHeight(); i < height; ++i) {
        for (size_t j = 0, width = a.getWidth(); j < width; ++j) {
            b(i, j) += a(i, k);
        }
    }
    std::cout << b << std::endl;
 
    return 0;
}