AssembledMatrixCache.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <memory>
 
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
#include "MathLib/LinAlg/MatrixVectorTraits.h"
 
namespace ProcessLib
{
struct AssembledMatrixCache final
{
    std::tuple<GlobalMatrix&, GlobalMatrix&, GlobalVector&> MKb() const
    {
        return std::tie(*M_, *K_, *b_);
    }

    bool hasMKb() const { return M_ && K_ && b_; }

    void storeMKb(GlobalMatrix const& M,
                  GlobalMatrix const& K,
                  GlobalVector const& b)
    {
        M_ = MathLib::MatrixVectorTraits<GlobalMatrix>::newInstance(M);
        K_ = MathLib::MatrixVectorTraits<GlobalMatrix>::newInstance(K);
        b_ = MathLib::MatrixVectorTraits<GlobalVector>::newInstance(b);
    }
 
private:
    std::unique_ptr<GlobalMatrix> M_;
    std::unique_ptr<GlobalMatrix> K_;
    std::unique_ptr<GlobalVector> b_;
};
}  // namespace ProcessLib