BasicMathPlugin.h

Go to the documentation of this file.

 
 
#pragma once
 
#include "AlignmentSubsystemForMathPlugins.h"
#include "ConvexHull.h"
 
#include <gsl/gsl_matrix.h>
 
namespace INDI
{
namespace AlignmentSubsystem
{
class BasicMathPlugin : public AlignmentSubsystemForMathPlugins
{
    public:
        BasicMathPlugin();
 
        virtual ~BasicMathPlugin();
 
        virtual bool Initialise(InMemoryDatabase *pInMemoryDatabase) override;
 
        virtual bool TransformCelestialToTelescopeJD(double RightAscension, double Declination,
                double JulianDate,
                TelescopeDirectionVector &ApparentTelescopeDirectionVector) override;
 
        virtual bool TransformTelescopeToCelestialJD(const TelescopeDirectionVector &ApparentTelescopeDirectionVector,
                double &RightAscension, double &Declination, double JulianDate) override;
 
        virtual bool TransformCelestialToTelescope(const double RightAscension, const double Declination,
                double JulianOffset,
                TelescopeDirectionVector &ApparentTelescopeDirectionVector) override;
 
        virtual bool TransformTelescopeToCelestial(const TelescopeDirectionVector &ApparentTelescopeDirectionVector,
                double &RightAscension, double &Declination, double JulianOffset = 0) override;
 
    protected:
        virtual void
        CalculateTransformMatrices(const TelescopeDirectionVector &Alpha1, const TelescopeDirectionVector &Alpha2,
                                   const TelescopeDirectionVector &Alpha3, const TelescopeDirectionVector &Beta1,
                                   const TelescopeDirectionVector &Beta2, const TelescopeDirectionVector &Beta3,
                                   gsl_matrix *pAlphaToBeta, gsl_matrix *pBetaToAlpha) = 0;
 
        void Dump3(const char *Label, gsl_vector *pVector);
 
        void Dump3x3(const char *Label, gsl_matrix *pMatrix);
 
        double Matrix3x3Determinant(gsl_matrix *pMatrix);
 
        bool MatrixInvert3x3(gsl_matrix *pInput, gsl_matrix *pInversion);
 
        void MatrixMatrixMultiply(gsl_matrix *pA, gsl_matrix *pB, gsl_matrix *pC);
 
        void MatrixVectorMultiply(gsl_matrix *pA, gsl_vector *pB, gsl_vector *pC);
 
        bool RayTriangleIntersection(TelescopeDirectionVector &Ray, TelescopeDirectionVector &TriangleVertex1,
                                     TelescopeDirectionVector &TriangleVertex2, TelescopeDirectionVector &TriangleVertex3);
 
        // Transformation matrixes for 1, 2 and 3 sync points case
        gsl_matrix *pActualToApparentTransform;
        gsl_matrix *pApparentToActualTransform;
 
        // Convex hulls for 4+ sync points case
        ConvexHull ActualConvexHull;
        ConvexHull ApparentConvexHull;
        // Actual direction cosines for the 4+ case
        std::vector<TelescopeDirectionVector> ActualDirectionCosines;
};
 
} // namespace AlignmentSubsystem
} // namespace INDI