CurvedSpaceIntegrator.hpp

#include "spacetime.h"

namespace Traum
{
        class   SchwarzschildGeodesic;
}

class   SPACETIME_API CurvedSpaceIntegrator
{
        Traum::SchwarzschildGeodesic*swg;

public:

        struct pt
        {
                double x,y,z;
        };


        CurvedSpaceIntegrator();
        ~CurvedSpaceIntegrator();

        void setMass(long double M) const;

        void setCenter(const pt&Center) const;


        int integrate(double time, const pt&O, const pt&d, double dist=-1, int maxsteps=5000);

        virtual int straight_ray(const pt&O, const pt&d) = 0;

        virtual int segment_intersect(const pt&where, const pt&dir, double segment_length) = 0;

        virtual int at_singularity() = 0;
};