#include <iostream>
#include <LatAnalyze/Io.hpp>
#include <LatAnalyze/Mat.hpp>
#include <LatAnalyze/Math.hpp>

using namespace std;
using namespace Latan;

int main(void)
{
    DMat A(2, 3), B(3, 2), C = DMat::Random(6, 6);
    const string fileName = "exMat.h5";
    
    A << 1, 2, 3,
         4, 5, 6;
    
    B << 1.0, 2.5,
         4.5, 8.9,
         1.2, 3.5;
    
    // read
    cout << "A=\n"       << A                         << '\n' << endl;
    cout << "B=\n"       << B                         << '\n' << endl;
    cout << "A*B=\n"     << A*B                       << '\n' << endl;
    cout << "cos(A)=\n"  << A.unaryExpr(StdMath::cos) << '\n' << endl;
    cout << "C=\n"       << C                         << '\n' << endl;
    cout << "inv(C)=\n"  << C.inverse()               << '\n' << endl;
    cout << "pinv(C)=\n" << C.pInverse()              << '\n' << endl;
    
    // write
    cout << "-- saving and loading A*B using '" + fileName + "'..." << endl;
    Io::save(A*B, fileName, File::Mode::write);

    DMat D = Io::load<DMat>(fileName);
    cout << D << endl;

    return EXIT_SUCCESS;
}