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;
};