hermitian.hpp

Go to the documentation of this file.

//=======================================================================
// Copyright (c) 2014-2023 Baptiste Wicht
// Distributed under the terms of the MIT License.
// (See accompanying file LICENSE or copy at
//  http://opensource.org/licenses/MIT)
//=======================================================================

#pragma once

#include <exception>

#include "etl/adapters/adapter.hpp"             // The adapter base class
#include "etl/adapters/hermitian_reference.hpp" // The reference proxy

namespace etl {

struct hermitian_exception : std::exception {
    const char* what() const noexcept override {
        return "Invalid assignment to a hermitian matrix";
    }
};

template <adaptable Matrix>
struct hermitian_matrix final : adapter<Matrix>, iterable<const hermitian_matrix<Matrix>> {
    using matrix_t  = Matrix;                   
    using expr_t    = matrix_t;                 
    using this_type = hermitian_matrix<Matrix>; 

    static constexpr size_t n_dimensions = etl_traits<matrix_t>::dimensions();  
    static constexpr order storage_order = etl_traits<matrix_t>::storage_order; 
    static constexpr size_t alignment    = matrix_t::alignment;                 

    using value_type        = value_t<matrix_t>; 
    using memory_type       = value_type*;       
    using const_memory_type = const value_type*; 

    using iterator       = typename matrix_t::const_iterator; 
    using const_iterator = typename matrix_t::const_iterator; 

    using base_type = adapter<Matrix>; 

    template <typename V = default_vec>
    using vec_type = typename V::template vec_type<value_type>;

    using base_type::value;

public:
    hermitian_matrix() noexcept : base_type() {
        //Nothing else to init
    }

    explicit hermitian_matrix(value_type value) noexcept : base_type(value) {
        //Nothing else to init
    }

    explicit hermitian_matrix(size_t dim) noexcept : base_type(dim) {
        //Nothing else to init
    }

    hermitian_matrix(size_t dim, value_type value) noexcept : base_type(dim, value) {
        //Nothing else to init
    }

    hermitian_matrix(const hermitian_matrix& rhs) = default;

    hermitian_matrix& operator=(const hermitian_matrix& rhs) = default;

    hermitian_matrix(hermitian_matrix&& rhs) noexcept = default;

    hermitian_matrix& operator=(hermitian_matrix&& rhs) noexcept = default;

    template <convertible_expr<value_type> E>
    hermitian_matrix& operator=(E&& e) noexcept(false) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(e)) {
            throw hermitian_exception();
        }

        validate_assign(*this, e);

        // Avoid aliasing issues
        if constexpr (!decay_traits<E>::is_linear) {
            if (e.alias(*this)) {
                // Create a temporary to hold the result
                this_type tmp(*this);

                // Assign the expression to the temporary
                tmp = e;

                // Assign the temporary to this matrix
                *this = tmp;
            } else {
                e.assign_to(*this);
            }
        } else {
            // Direct assignment of the expression into this matrix
            e.assign_to(*this);
        }

        return *this;
    }

    hermitian_matrix& operator+=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_add_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator+=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_add_to(*this);
        return *this;
    }

    hermitian_matrix& operator-=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_sub_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator-=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_sub_to(*this);
        return *this;
    }

    hermitian_matrix& operator*=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_mul_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator*=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_mul_to(*this);
        return *this;
    }

    hermitian_matrix& operator>>=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_mul_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator>>=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_mul_to(*this);
        return *this;
    }

    hermitian_matrix& operator/=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_div_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator/=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_div_to(*this);
        return *this;
    }

    hermitian_matrix& operator%=(const value_type& rhs) noexcept {
        etl::scalar<value_type>(rhs).assign_mod_to(*this);
        return *this;
    }

    template <etl_expr R>
    hermitian_matrix& operator%=(const R & rhs) {
        // Make sure the other matrix is hermitian
        if (!is_hermitian(rhs)) {
            throw hermitian_exception();
        }

        validate_expression(*this, rhs);
        rhs.assign_mod_to(*this);
        return *this;
    }

    herm_detail::hermitian_reference<matrix_t> operator()(size_t i, size_t j) noexcept {
        return {value, i, j};
    }

    using base_type::operator();
};

template <typename Matrix>
struct etl_traits<hermitian_matrix<Matrix>> : wrapper_traits<hermitian_matrix<Matrix>> {};

} //end of namespace etl