Extrude.h

/*
Hyperspace Explorer - visualizing higher-dimensional geometry
Copyright (C) 2010  Lene Preuss <lene.preuss@gmail.com>

This program 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 2 of the License, or
(at your option) any later version.

This program 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 this program; if not, write to the Free Software Foundation, Inc.,
51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

*/

#ifndef EXTRUSION_H
#define EXTRUSION_H

#include "Rotope.h"


template <unsigned D>
    class extrude_base: public VertexData<D> {

    public:

        extrude_base(): VertexData<D>() { }
        extrude_base(const VertexData<D> &v): VertexData<D>(v) { }


        void extrude(unsigned d);

    };


template <unsigned D, unsigned Dmin, unsigned Dmax>
    class Extrude: public Extrude<D, Dmin, Dmax-1> {
    public:

        Extrude(): Extrude<D, Dmin, Dmax-1>() {
            extrude_base<D>::extrude(Dmax);
        }


        Extrude(const VertexData<D> &v): Extrude<D, Dmin, Dmax-1>(v) {
            extrude_base<D>::extrude(Dmax);
        }
    };


template <unsigned D, unsigned Dmin>
    class Extrude<D, Dmin, Dmin>: public extrude_base<D> {
    public:

        Extrude(): extrude_base<D>() {
            extrude_base<D>::extrude(Dmin);
        }


        Extrude(const VertexData<D> &v): extrude_base<D>(v) {
            extrude_base<D>::extrude(Dmin);
        }
    };


template <unsigned D> void extrude_base<D>::extrude(unsigned d) {
    if (d >= D) {
        throw std::logic_error(
           "extrude_base::extrude() called on a higher dimension"
           " than the vector space allows");
    }

    VecMath::Vector<D> x;
    x[d] = 1.;
    unsigned xsize = VertexData<D>::raw_vertices().size();
    for (unsigned i = 0; i < xsize; ++i) {
        VertexData<D>::raw_vertices()[i] -= x;
        VertexData<D>::raw_vertices().push_back(VertexData<D>::raw_vertices()[i]+x*2.);
    }

    if (d == 0) VertexData<D>::realm() = Realm(0);

    VertexData<D>::realm() = VertexData<D>::realm().extruded(xsize);

    VertexData<D>::dimension()++;   //  object is now one dimension higher
}

#endif