EigenMatrix.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "EigenMatrix.h"
 
#include <fstream>
 
namespace MathLib
{

void EigenMatrix::write(const std::string& filename) const
{
    std::ofstream of(filename);
    if (of)
    {
        write(of);
    }
}

void EigenMatrix::write(std::ostream& os) const
{
    for (int k = 0; k < mat_.outerSize(); ++k)
    {
        for (RawMatrixType::InnerIterator it(mat_, k); it; ++it)
        {
            os << it.row() << " " << it.col() << " " << it.value() << "\n";
        }
    }
    os << std::endl;
}
 
void EigenMatrix::addToDiagonal(double const value)
{
    // Add value to all diagonal entries
    for (IndexType i = 0; i < mat_.rows(); ++i)
    {
        mat_.coeffRef(i, i) += value;
    }
}
}  // namespace MathLib