range_view.hpp

Go to the documentation of this file.

#ifndef RANGE_VIEW_HPP
#define RANGE_VIEW_HPP
#include "space_settings.hpp"

// Create range
inline View<double*, HostType> create_range_view(std::string name, double x_min, double x_max, int n){
    View<double*, HostType> range(NoInit(name), n);
    double dx = (x_max - x_min)/(n-1);
    for(int i = 0; i<n; i++){
        range(i) = x_min + i*dx;
    }
    // To ensure to machine precision
    range(n-1) = x_max;

    return range;
}

#endif