RK8Propagator

class ssapy.propagator.RK8Propagator(accel, h)[source][source]

Bases: RKPropagator

Runge-Kutta 8th order numerical integrator.

Parameters:
  • accel (ssapy.Accel) – Accel object containing the acceleration model by which to propagate.

  • h (float) – Step size in seconds. ~70s yields accuracy of ~1e-6 meters at GEO over a couple of days. ~20s yields accuracy of ~1e-5 meters at LEO over a few hours.

Attributes Summary

a

b8

c

Attributes Documentation

a = array([[ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 5.55555556e-02,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 2.08333333e-02,  6.25000000e-02,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 3.12500000e-02,  0.00000000e+00,  9.37500000e-02,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 3.12500000e-01,  0.00000000e+00, -1.17187500e+00,          1.17187500e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 3.75000000e-02,  0.00000000e+00,  0.00000000e+00,          1.87500000e-01,  1.50000000e-01,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 4.79101371e-02,  0.00000000e+00,  0.00000000e+00,          1.12248713e-01, -2.55056738e-02,  1.28468239e-02,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 1.69179898e-02,  0.00000000e+00,  0.00000000e+00,          3.87848278e-01,  3.59773699e-02,  1.96970214e-01,         -1.72713852e-01,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 6.90957534e-02,  0.00000000e+00,  0.00000000e+00,         -6.34247977e-01, -1.61197575e-01,  1.38650309e-01,          9.40928614e-01,  2.11636326e-01,  0.00000000e+00,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 1.83556997e-01,  0.00000000e+00,  0.00000000e+00,         -2.46876808e+00, -2.91286888e-01, -2.64730202e-02,          2.84783876e+00,  2.81387331e-01,  1.23744900e-01,          0.00000000e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [-1.21542482e+00,  0.00000000e+00,  0.00000000e+00,          1.66726087e+01,  9.15741828e-01, -6.05660580e+00,         -1.60035736e+01,  1.48493031e+01, -1.33715757e+01,          5.13418265e+00,  0.00000000e+00,  0.00000000e+00,          0.00000000e+00],        [ 2.58860916e-01,  0.00000000e+00,  0.00000000e+00,         -4.77448579e+00, -4.35093014e-01, -3.04948333e+00,          5.57792004e+00,  6.15583159e+00, -5.06210459e+00,          2.19392617e+00,  1.34627999e-01,  0.00000000e+00,          0.00000000e+00],        [ 8.22427600e-01,  0.00000000e+00,  0.00000000e+00,         -1.16586733e+01, -7.57622117e-01,  7.13973588e-01,          1.20757750e+01, -2.12765911e+00,  1.99016621e+00,         -2.34286472e-01,  1.75898578e-01,  0.00000000e+00,          0.00000000e+00]])[source]
b8 = array([ 0.04174749,  0.        ,  0.        ,  0.        ,  0.        ,        -0.05545233,  0.23931281,  0.70351067, -0.75975961,  0.66056303,         0.15818748, -0.23810954,  0.25      ])[source]
c = array([0.        , 0.05555556, 0.08333333, 0.125     , 0.3125    ,        0.375     , 0.1475    , 0.465     , 0.56486545, 0.65      ,        0.92465628, 1.        , 1.        ])[source]