jonpsy

from poliastro.iod.vallado import lambert as lambert_v
from poliastro.iod.izzo import lambert as lambert_i 
from astropy import units as u, constants as c


if __name__ == "__main__":
    k = (c.GM_earth.to(u.km ** 3/ u.s ** 2))
    r0 = ([10000., 0, 0] * u.km)
    rf = ([8000., -5000, 0] * u.km)
    tof = 2*u.hour

jonpsy

from poliastro.iod.vallado import lambert as vallado
from astropy import units as u, constants as c


if __name__ == "__main__":
    k = (c.GM_earth.to(u.km ** 3/ u.s ** 2))
    r0 = [10000., 0, 0] * u.km
    rf = [8000., -5000, 0] * u.km
    tof = 2 * u.hour
    numiter = 10**3

jonpsy

import numpy as np
from numpy import cross, pi
from numpy.linalg import norm

from poliastro.core.hyper import hyp2f1b
from poliastro.core.stumpff import c2, c3

from ._jit import jit

jonpsy

import numpy as np
from numpy import pi

from poliastro.core.hyper import hyp2f1b
from poliastro.core.stumpff import c2, c3
from poliastro.core.util import cross, norm

from ._jit import jit

jonpsy

jonpsy

jonpsy

jonpsy

jonpsy

  #include <bits/stdc++.h>
  #include <armadillo>

  template<typename MatType,
          typename... ArbitraryFunctionType,
          typename... CallbackTypes>
  typename MatType::elem_type SPEA2::Optimize(
      std::tuple<ArbitraryFunctionType...>& objectives,
      MatType& iterate,
      CallbackTypes&&... callbacks)

jonpsy

#include <bits/stdc++.h>
#include <armadillo>
/**
 * @brief Code heavily inspired from Pagmo. All credits to Dario Izzo.
 */

double archiveSize = 10;

// Possibly make it a inner class???
template<typename MatType>