#include<iostream>
#include<fstream>
#include <cmath>
using namespace std;

class Matrix {
private:
void method();

private:
double** matr; //matrix
int rown; //rown
int column; //column

std::string m_name;
public:
Matrix(const char *name)
{
    m_name = name;

    matr = NULL;
    rown = column = 0;

    cout << __FILE__ << ":" << __LINE__ << " CONSTRUCTOR " << m_name.data() << " " << matr << endl;
}

void allocate(int in_rown, int in_column)
{
    rown = in_rown;
    column = in_column;

    matr = new double*[rown];
    for (int i = 0; i < rown; i++) {
        matr[i] = new double[column];
    }

    cout << __FILE__ << ":" << __LINE__ << " ALLOCATOR " << m_name.data() << " " << matr << endl;
}

Matrix operator*(Matrix & mult);
void read_matrix(const char*filename);
void write_matrix(const char*filename) const;

~Matrix()
{
    cout << __FILE__ << ":" << __LINE__ << " DESTRUCTOR " << m_name.data() << " " << matr << endl;

    if (matr != NULL) {
        for (int i = 0; i < rown; i++) {
            delete[] matr[i];
        }

        delete[] matr;
    }
}
};

Matrix Matrix::operator*(Matrix & mult) {
Matrix result("result");


result.allocate(rown, mult.column);

for (int i = 0; i < result.rown; i++)
{
    for (int k = 0; k < result.column; k++)
    {
        result.matr[i][k] = 0;
        for (int j = 0; j < result.column; j++)
        {
            result.matr[i][k] += matr[i][j] * mult.matr[j][k];
        }
    }
}

return result;
}

void Matrix::read_matrix(const char* filename) {
ifstream in(filename);

in >> rown;
in >> column;

allocate(rown, column);

//Считаем матрицу из файла
for (int i = 0; i < rown; i++)
    for (int j = 0; j < column; j++)
        in >> matr[i][j];

in.close();
return;
}

void Matrix::write_matrix(const char* filename) const {
ofstream on(filename);

on << rown << 't' << column << endl;

//запишем матрицу в файл
for (int i = 0;i < rown;i++) {
    for (int k = 0;k < column;k++) {
        on << matr[i][k] << 't';
    }
    on << endl;
}

on.close();
return;
}

int main(int argc, char*argv[])
{
if (argc == 4) {

    Matrix input1("input1");
    Matrix input2("input2");
    input1.read_matrix(argv[1]);
    input2.read_matrix(argv[2]);

    Matrix a("a");
    a = input1*input2;
    /*Matrix b;*/
    a.write_matrix(argv[3]);
}
system("pause");
return 0;
 }

a = input1*input2;