interval.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_INTERVAL_HH
#define DUNE_FUFEM_INTERVAL_HH
 
#include <array>
#include <iostream>
#include <algorithm>
 
template <class field_type>
class Interval
{
public:
    field_type& operator[](int i)
    {
        return data_[i];
    }
    
    const field_type& operator[](int i) const
    {
        return data_[i];
    }
    
    field_type projectIn(const field_type& x) const
    {
        return std::max(std::min(x,data_[1]), data_[0]);
    }
    
    field_type projectFromBelow(const field_type& x) const
    {
        return std::max(x,data_[0]);
    };
 
    field_type projectFromAbove(const field_type& x) const
    {
        return std::min(x,data_[1]);
    };
 
    bool containsZero(const field_type& safety) const
    {
        return (data_[0] <= safety) and (-safety <= data_[1]);
    };
 
private:
    
    std::array<field_type,2> data_;
};
 
template <class field_type>
inline std::ostream& operator<< (std::ostream& s, const Interval<field_type>& interval)
{
    s << "[" << interval[0] << ", " << interval[1] << "]";
    return s;
}
  
#endif