LegacyBandedMatrix.hpp

Go to the documentation of this file.

 
//
// Copyright (c) 2021-2023, University of Colorado Denver. All rights reserved.
//
// This file is part of <T>LAPACK.
// <T>LAPACK is free software: you can redistribute it and/or modify it under
// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
 
#ifndef TLAPACK_LEGACY_BANDED_HH
#define TLAPACK_LEGACY_BANDED_HH
 
#include <cassert>
 
#include "tlapack/base/exceptionHandling.hpp"
 
namespace tlapack {
 
template <typename T, class idx_t = std::size_t>
struct LegacyBandedMatrix {
    idx_t m, n, kl, ku;  
    T* ptr;              
 
    constexpr const T& operator()(idx_t i, idx_t j) const noexcept
    {
        assert(i >= 0);
        assert(i < m);
        assert(j >= 0);
        assert(j < n);
        assert(j <= i + ku);
        assert(i <= j + kl);
        return ptr[(ku + i) + j * (ku + kl)];
    }
 
    constexpr T& operator()(idx_t i, idx_t j) noexcept
    {
        assert(i >= 0);
        assert(i < m);
        assert(j >= 0);
        assert(j < n);
        assert(j <= i + ku);
        assert(i <= j + kl);
        return ptr[(ku + i) + j * (ku + kl)];
    }
 
    constexpr LegacyBandedMatrix(idx_t m, idx_t n, idx_t kl, idx_t ku, T* ptr)
        : m(m), n(n), kl(kl), ku(ku), ptr(ptr)
    {
        tlapack_check(m >= 0);
        tlapack_check(n >= 0);
        tlapack_check((kl + 1 <= m || m == 0) && (ku + 1 <= n || n == 0));
    }
};
 
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_BANDED_HH