SimpleMatrixVectorProvider.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 <map>
#include <memory>
 
#include "MatrixProvider.h"
#include "VectorProvider.h"
 
namespace NumLib
{
class SimpleMatrixVectorProvider final : public MatrixProvider,
                                         public VectorProvider
{
public:
    SimpleMatrixVectorProvider() = default;
 
    // no copies
    SimpleMatrixVectorProvider(SimpleMatrixVectorProvider const&) = delete;
    SimpleMatrixVectorProvider& operator=(SimpleMatrixVectorProvider const&) =
        delete;
 
    GlobalVector& getVector(std::size_t& id) override;
 
    GlobalVector& getVector(GlobalVector const& x) override;
    GlobalVector& getVector(GlobalVector const& x, std::size_t& id) override;
 
    GlobalVector& getVector(MathLib::MatrixSpecifications const& ms) override;
    GlobalVector& getVector(MathLib::MatrixSpecifications const& ms,
                            std::size_t& id) override;
 
    void releaseVector(GlobalVector const& x) override;
 
    GlobalMatrix& getMatrix(std::size_t& id) override;
 
    GlobalMatrix& getMatrix(MathLib::MatrixSpecifications const& ms,
                            std::size_t& id) override;
 
    void releaseMatrix(GlobalMatrix const& A) override;
 
    void clear();
 
    ~SimpleMatrixVectorProvider() override;
 
private:
    template <typename... Args>
    std::pair<GlobalMatrix*, bool> getMatrix_(std::size_t& id, Args&&... args);
 
    template <typename... Args>
    std::pair<GlobalVector*, bool> getVector_(std::size_t& id, Args&&... args);
 
    // returns a pair with the pointer to the matrix/vector and
    // a boolean indicating if a new object has been built (then true else
    // false)
    template <typename MatVec, typename... Args>
    std::pair<MatVec*, bool> get_(std::size_t& id,
                                  std::map<MatVec*, std::size_t>& used_map,
                                  Args&&... args);
 
    std::size_t _next_id = 1;
 
    std::map<GlobalMatrix*, std::size_t> _used_matrices;
 
    std::map<GlobalVector*, std::size_t> _used_vectors;
};
 
}  // namespace NumLib