Sponge.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.

*/

#if(!defined SPONGE_H)
#define SPONGE_H

#define USE_ALT_SPONGE 0

#include "ObjectImplementations.h"



class Sponge: public Hypercube {
public:

    Sponge (unsigned _level = 1, unsigned _distance = 2, double _rad = 0.8,
            VecMath::Vector<4> _center = VecMath::Vector<4>(0., 0.,0., 0.));
    virtual ~Sponge ();
    virtual std::string getFunctionName() const { return "4D Menger Sponge"; }

    virtual void Transform (const VecMath::Rotation<4> &R,
                            const VecMath::Vector<4> &T,
                            const VecMath::Vector<4> &scale);
    virtual void Project (double ScrW, double CamW, bool DepthCue4D);
    virtual void Draw (UI::View *view);

    virtual void SetParameters(const ParameterMap &parms) {
        std::cerr << "Sponge::SetParameters(" << parms.toString() << ")\n";
#       if 1
            for (ParameterMap::const_iterator i = parms.begin();
                 i != parms.end(); ++i) {
                if (i->second->getName() == "Level")
                    Level = i->second->toUnsigned();
                if (i->second->getName() == "Distance")
                    distance = i->second->toInt();
                if (i->second->getName() == "Size")
                    rad = i->second->toDouble();
            }
#       else
            setParameter(parms, this->Phase, "Phase");
#       endif
        }

    virtual void ReInit (double, double, double,
                         double, double, double,
                         double, double, double) {
        List.clear();
        Object::ReInit(0,0,0,0,0,0,0,0,0);
    }

    virtual std::string description () {
        std::ostringstream out;
        out << "Sponge (level = " << Level << ")" << std::ends;
        return out.str ();
    }

protected:
    virtual void Initialize();
    virtual unsigned long MemRequired (unsigned);
    unsigned Level;                 


    std::vector<Hypercube *> List;

    int distance;                   
    double rad;                     
    VecMath::Vector<4> center;      
};

namespace {
    Displayable *createSponge() { return new Sponge(); }
    const bool registeredS = TheFunctionFactory::Instance().registerFunction(createSponge, "Object");
}


class AltSponge: public Hypercube {
public:

  AltSponge (unsigned _level = 1, unsigned _distance = 2, double _rad = 0.8,
             VecMath::Vector<4> _center = VecMath::Vector<4>(0., 0.,0., 0.));
  virtual ~AltSponge () { }
  virtual std::string getFunctionName() const { return "4D Menger Sponge (version)"; }

  virtual void SetParameters(const ParameterMap &parms) {
    std::cerr << "AltSponge::SetParameters(" << parms.toString() << ")\n";
#       if 1
    for (ParameterMap::const_iterator i = parms.begin();
         i != parms.end(); ++i) {
      if (i->second->getName() == "Level")
        Level = i->second->toUnsigned();
      if (i->second->getName() == "Distance")
        distance = i->second->toInt();
      if (i->second->getName() == "Size")
        rad = i->second->toDouble();
     }
#       else
     setParameter(parms, this->Phase, "Phase");
#       endif
  }

  virtual std::string description () {
    std::ostringstream out;
    out << "Alternative Sponge (level = " << Level << ")" << std::ends;
    return out.str ();
  }

protected:
  virtual void Initialize();
  virtual unsigned long MemRequired (unsigned);

  void reduceVertices();
  void removeDuplicateSurfaces();

  void renumberSurfaces(unsigned original_vertex, unsigned duplicate_vertex);
  
  void X_push_back(const VecMath::Vector<4> &x);

  unsigned Level;                 
  int distance;                   
        double rad;                     
        VecMath::Vector<4> center;      
};

#if USE_ALT_SPONGE
namespace {
  Displayable *createAltSponge() { return new AltSponge(); }
  const bool registeredAS = TheFunctionFactory::Instance().registerFunction(createAltSponge, "Object");
}
#endif

#endif