gsl_vector.hpp

Go to the documentation of this file.

// ====================================================================
// This file is part of FlexibleSUSY.
//
// FlexibleSUSY is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published
// by the Free Software Foundation, either version 3 of the License,
// or (at your option) any later version.
//
// FlexibleSUSY is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FlexibleSUSY.  If not, see
// <http://www.gnu.org/licenses/>.
// ====================================================================
 
#ifndef GSL_VECTOR_H
#define GSL_VECTOR_H
 
#include "gsl.hpp"
#include <gsl/gsl_vector.h>
#include <initializer_list>
#include <iosfwd>
#include <cstddef>
 
namespace flexiblesusy {
 
class GSL_vector {
public:
   GSL_vector() = default;
   GSL_vector(std::size_t);
   GSL_vector(const gsl_vector*);        
   GSL_vector(const GSL_vector&);
   GSL_vector(GSL_vector&&) noexcept;
   GSL_vector(std::initializer_list<double>);
   ~GSL_vector() noexcept;
 
   GSL_vector& operator=(const GSL_vector&);
   GSL_vector& operator=(GSL_vector&&) noexcept;
   double& operator[](std::size_t);      
   double operator[](std::size_t) const; 
   double& operator()(std::size_t);      
   double operator()(std::size_t) const; 
 
   bool empty() const noexcept;             
   const gsl_vector* raw() const noexcept;  
   gsl_vector* raw() noexcept;              
   gsl::owner<gsl_vector>* release() noexcept; 
   void set_all(double) noexcept;           
   std::size_t size() const noexcept;       
 
private:
   gsl::owner<gsl_vector>* vec{nullptr};    
 
   void assign(const gsl_vector*);          
   void move_assign(GSL_vector&&) noexcept; 
   void range_check(std::size_t) const;
};
 
double* begin(GSL_vector&); 
double* end(GSL_vector&);   
 
const double* cbegin(const GSL_vector&); 
const double* cend(const GSL_vector&);   
 
bool is_finite(const GSL_vector&);
 
std::ostream& operator<<(std::ostream&, const GSL_vector&);
 
} // namespace flexiblesusy
 
#endif