dss_handle_umfpack.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <memory>
 
#include <umfpack.h>
 
#include "dss_algorithm.hpp"
 
namespace kynema::dss {
template <>
class Handle<Algorithm::UMFPACK> {
    struct umfpackDssHandleType {
        void* Symbolic = nullptr;
        void* Numeric = nullptr;
        std::array<double, UMFPACK_CONTROL> Control{};
 
        umfpackDssHandleType() {
            umfpack_di_defaults(Control.data());
            Control[UMFPACK_ORDERING] = UMFPACK_ORDERING_BEST;
        }
 
        umfpackDssHandleType(umfpackDssHandleType&) = delete;
        void operator=(umfpackDssHandleType&) = delete;
        umfpackDssHandleType(umfpackDssHandleType&&) = delete;
        void operator=(umfpackDssHandleType&&) = delete;
 
        ~umfpackDssHandleType() {
            umfpack_di_free_symbolic(&Symbolic);
            umfpack_di_free_numeric(&Numeric);
        }
    };
    std::shared_ptr<umfpackDssHandleType> umfpack_dss_handle;
 
public:
    Handle() : umfpack_dss_handle(std::make_shared<umfpackDssHandleType>()) {}
 
    void*& get_symbolic() { return umfpack_dss_handle->Symbolic; }
 
    void*& get_numeric() { return umfpack_dss_handle->Numeric; }
 
    double* get_control() { return umfpack_dss_handle->Control.data(); }
};
 
}  // namespace kynema::dss