numeric.h

// -*- mode: c++; tab-width: 4; indent-tabs-mode: t; eval: (progn (c-set-style "stroustrup") (c-set-offset 'innamespace 0)); -*-
// vi:set ts=4 sts=4 sw=4 noet :
// Copyright 2011, 2012, The TPIE development team
// 
// This file is part of TPIE.
// 
// TPIE is free software: you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License as published by the
// Free Software Foundation, either version 3 of the License, or (at your
// option) any later version.
// 
// TPIE 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 Lesser General Public
// License for more details.
// 
// You should have received a copy of the GNU Lesser General Public License
// along with TPIE.  If not, see <http://www.gnu.org/licenses/>

#ifndef __TPIE_PIPELINING_NUMERIC_H__
#define __TPIE_PIPELINING_NUMERIC_H__

#include <iostream>
#include <tpie/pipelining/node.h>
#include <tpie/pipelining/pipe_base.h>
#include <tpie/pipelining/factory_helpers.h>

namespace tpie {

namespace pipelining {

namespace bits {

template <typename dest_t>
class linear_t : public node {
public:
    typedef typename push_type<dest_t>::type item_type;

    inline linear_t(dest_t dest, item_type factor, item_type term) : factor(factor), term(term), dest(std::move(dest)) {
        add_push_destination(this->dest);
        set_name("Linear transform", PRIORITY_INSIGNIFICANT);
    }
    inline void push(const item_type & item) {
        dest.push(item*factor+term);
    }
private:
    item_type factor, term;
    dest_t dest;
};

template <typename dest_t>
class range_t : public node {
public:
    typedef typename push_type<dest_t>::type item_type;

    range_t(dest_t dest, item_type from, item_type to, item_type increment) : from(from), to(to), increment(increment), dest(std::move(dest)) {}

    void propagate() {
        stream_size_type items = (from - to) / increment;
        set_steps(items);
        forward("items", items);
    }
    
    virtual void go() override {
        for (item_type i=from; i < to; i += increment) {
            dest.push(i);
            step(1);
        }
    }

private:
    item_type from, to, increment;
    dest_t dest;
};

} // namespace bits

template <typename T>
inline pipe_middle<factory<bits::linear_t, T, T> >
linear(T factor, T term) {
    return factory<bits::linear_t, T, T>(factor, term);
}

template <typename T>
inline pipe_begin<factory<bits::range_t, T, T, T> >
range(T from, T to, T increment = 1) {
    return factory<bits::range_t, T, T, T>(from, to, increment);
}


} // namespace pipelining

} // namespace tpie

#endif