isosurface.h

/* *************************************************************************
                          isosurface.h  -  description
                             -------------------
    begin                : Wed Oct 16 2002        
    copyright            : (C) 2002 by Micha Riser
    email                : mriser@gmx.net

    $Id: isosurface_8h-source.html,v 1.1 2005/02/04 21:02:27 micha Exp $
    
 ************************************************************************* */

/* *************************************************************************
 *                                                                         *
 *   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.                                   *
 *                                                                         *
 ************************************************************************* */

#ifndef ISOSURFACE_H
#define ISOSURFACE_H

#include "object3d.h"
#include "intersection.h"
#include "function.h"

class Isosurface: public Object3D  {

public: // Type definitions
    typedef Function<const Vector3,DBL> IsosurfaceFunction;

public: // Constructor and destructor

    Isosurface(IsosurfaceFunction* func, Object3D* cont, DBL maxd, bool doeval);
    
    ~Isosurface();
    
private: // Private constructor
    Isosurface(const Isosurface&);
    Isosurface& operator=(const Isosurface&);

public: // Public methods

    // Inherited from Object3D
    Object3D* duplicate() const;
    bool isInsideAt(const Vector3& locallocation) const;
    Intersection* intersectLocalLimitedTime(const Ray& r, DBL maxtime) const;

private: // Constants
    static const DBL NORMAL_E;
    static const int MAX_STEP_INCREASE = 1024;
#ifdef WRONG_INNER_CLASSES_ACCESS
public:
#endif
    static const DBL K_INCREASE_FACT;
    
private: // Static members
    static HitStatistics stat;
    static int isosurfacecount;

private: // Private attributes
#ifdef WRONG_INNER_CLASSES_ACCESS
public:
#endif

    Object3D* container;

    IsosurfaceFunction* f;

    mutable DBL K;
    mutable DBL missstepsize;
    bool eval;

    // max derivative statistics
    mutable MaxStatistics<DBL> maxdstat;

private: // Private inner classes

    class StackElement {
    public: // Public attributes
        DBL time;
        DBL value;
    public: // Constructors and destructor
        StackElement() {}
        StackElement(DBL t, DBL v): time(t), value(v) {}
    };

    class IsosurfaceIntersection: public Intersection {

    public: // Public attributes
        bool nointersection;

    public: // Constructor and destructor
        IsosurfaceIntersection(const Isosurface* object, const Ray& r, DBL tmax);
        ~IsosurfaceIntersection();

    public: // Public methods
        void next();
        DBL currentTime();
        Vector3 normal();
        
    private: // Private attributes
        const Isosurface* iso;
        StackElement prelast, last;
        DBL root;
        DBL klocal;
        DBL maxtime;
        DBL t_container_max;
        bool containerhit;
        Intersection* i_container;

    private: // Private methods
        void findExactRoot();
        void findRoot();
        DBL calcPolyRoot();

    };

};

#endif

---