testGravitationalPotential.h

Go to the documentation of this file.

#ifndef __TEST_GRAVITATIONAL_POTENTIAL_H
#define __TEST_GRAVITATIONAL_POTENTIAL_H

#include "EccentricOrbit.h"
#include "TidalPotential.h"
#include "TidalPotentialExpansion.h"

#include "../shared/Common.h"

#include <iostream>
#include <iomanip>
#include <sstream>

namespace Evolve {
    class test_GravitationalPotential : public Test::Suite {
    private:
        void print_orbit(const EccentricOrbit &orbit) const;

        void print_tidal_potential(
            double primary_mass,

            double secondary_mass,

            double semimajor,

            double eccentricity,

            double inclination,

            double arg_of_periapsis,

            const Eigen::Matrix<long double, 3, 1> &position,

            unsigned expansion_order
        ) const;

        double abs_precision(
            const Eigen::Matrix<long double, 3, 1> &position,

            const EccentricOrbit &orbit,

            double expansion_precision
        ) const;

        void test_single_point(
            TidalPotentialExpansion &approx_potential,

            const TidalPotential &exact_potential,

            const Eigen::Matrix<long double, 3, 1> &position
        );

        void test_system(
            double primary_mass,

            double secondary_mass,

            double semimajor,

            double eccentricity,

            double inclination,

            double arg_of_periapsis
        );

    protected:
        void setup() {};

        void tear_down() {};

    public:
        test_GravitationalPotential();

        void test_expansion();
    }; //End test_GravitationalPotential class.
} //End Evolve namespace.

#endif