Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Sep 28 22:27
    mikeheddes closed #109
  • Sep 28 22:27
    mikeheddes commented #109
  • Sep 28 21:18
    mikeheddes opened #109
  • Sep 19 10:50
    Sceki synchronize #108
  • Sep 19 00:00
    Sceki edited #108
  • Sep 18 23:59
    Sceki opened #108
  • Sep 18 23:41
  • Sep 10 14:09
    BookOwl opened #107
  • Sep 02 13:14

    darioizzo on master

    added spherical harmonics files new names duplicate and 25 more (compare)

  • Sep 02 13:14
    darioizzo closed #101
  • Sep 02 13:14

    darioizzo on master

    fix(juice): working multi-objec… fix(juice): unpack time only wh… Merge pull request #106 from mi… (compare)

  • Sep 02 13:14
    darioizzo closed #106
  • Sep 02 09:33
    mikeheddes review_requested #106
  • Aug 29 13:34
    mikeheddes synchronize #106
  • Aug 29 12:51
    mikeheddes opened #106
  • Aug 29 12:29
  • Aug 25 23:08
  • Aug 24 23:39
    cpsancha commented #95
  • Aug 24 23:37
    sork-dev commented #95
  • Aug 13 12:56
Jorge Martinez
@jorgepiloto
I expected that procedure when trying to work with those C subroutines. Anyway, thank you for your fast answer @darioizzo :smile:
Dario Izzo
@darioizzo
@matrixbot try to use floats when constructing l ... let me know
matrixbot
@matrixbot
astrojuanlu Dario Izzo (Gitter): I copy-pasted l = lambert_problem([1,0,0],[0,1,0],5 * pi / 2. ) from the docstring :)
Dario Izzo
@darioizzo
Can you use floats instead? @matrixbot
matrixbot
@matrixbot

astrojuanlu ```

l = lambert_problem([1., 0., 0.], [0., 1., 0.], 5 * pi / 2.)
l
Lambert's problem:
mu = 1
r1 = [1, 0, 0]
r2 = [0, 1, 0]
Time of flight: 7.8539816339745

chord = 1.4142135623731
semiperimeter = 1.7071067811865
lambda = 0.41421356237309
non dimensional time of flight = 4.9798223260205

Maximum number of revolutions: 1
Solutions:
0 revs, Iters: 3, x: -0.58644325543794, a: 1.3009812610796
v1= [0.90230487720848529, 0.64590659483860502, 0] v2= [-0.64590659483860502, -0.90230487720848529, 0]
1 revs, left. Iters: 4, x: -0.11188837995933, a: 0.86437450327173
v1= [0.45213333668550759, 0.79916800651739384, 0] v2= [-0.79916800651739384, -0.45213333668550759, 0]
1 revs, right. Iters: 4, a: 0.38268343236509, a: 1
v1= [1.0257123289739564e-16, 0.99999999999999978, 0] v2= [-0.99999999999999978, -1.0257123289739564e-16, 0]

l.get_iters()
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
TypeError: No to_python (by-value) converter found for C++ type: std::vector<int, std::allocator<int> >

```

astrojuanlu and another attempt:

astrojuanlu ```

l = lambert_problem([1.0001, 0.0001, 0.0001], [0.0001, 1.0001, 0.00001], 5 * pi / 2.)
l
Lambert's problem:
mu = 1
r1 = [1.0001, 0.0001, 0.0001]
r2 = [0.0001, 1.0001, 1.0000000000000001e-05]
Time of flight: 7.8539816339745

chord = 1.4142135652369
semiperimeter = 1.7072067901427
lambda = 0.41427213666497
non dimensional time of flight = 4.9793847520179

Maximum number of revolutions: 1
Solutions:
0 revs, Iters: 3, x: -0.58641808361264, a: 1.3009989341246
v1= [0.90219671389791434, 0.64591093990950621, 9.6661755009850417e-05] v2= [-0.64591093991027215, -0.90219672155874819, -7.3596035427494358e-05]
1 revs, left. Iters: 4, x: -0.11177712481321, a: 0.86440335923796
v1= [0.4519570329407957, 0.79916676280745691, 5.3173611312322879e-05] v2= [-0.79916676280807564, -0.45195703913200069, -8.4422487062228236e-05]
1 revs, right. Iters: 4, a: 0.38258941011596, a: 0.99997428825661
v1= [2.5732593140190884e-05, 0.99988713704973886, 9.9904477876137778e-06] v2= [-0.99988713705023335, -2.5737541224102957e-05, -9.997797423818353e-05]

l.get_iters()
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
TypeError: No to_python (by-value) converter found for C++ type: std::vector<int, std::allocator<int> >

``

matrixbot
@matrixbot
astrojuanlu (ooooops sorry for the broken formatting, I will upload a gist instead)
Jorge Martinez
@jorgepiloto
Hi all, I was wondering if anyone of you has ever heard about Gauss' algorithm for Kepler's problem. I found this algorithm in Battin 1999. However, and even after a hard debugging phase I was not able to make it totally work. Still doing a long research and looking for more information about it on different astrodynamicists software groups :rocket:
astrojuanlu 😉
Jorge Martinez
@jorgepiloto
Oh, sorry I meant the "Extended Gauss Method" algorithm :sweat_smile:
matrixbot
@matrixbot
astrojuanlu oh, no idea! which section of Battin? you might also want to ask in #openastrodynamics:matrix.org
Giacomo Acciarini
@Sceki
Hello, I have tried to generate 20 individuals for the juice_mo problem, but it does not seem to work.. in particular, when doing:
import pykep as pk
from pygmo import *
prob=pk.trajopt.gym.juice_mo
pop=population(prob, 20)
Jupyter notebook crashes
with a malloc error that I cannot really decipher:
python(23217,0x11b68f5c0) malloc: *** error for object 0x7fa0113dc408: pointer being freed was not allocated
python(23217,0x11b68f5c0) malloc: *** set a breakpoint in malloc_error_break to debug
Do you know what could cause the issue?
Ops, I didn't see PR #106 , maybe that fixes the issue
Giacomo Acciarini
@Sceki
indeed, synching the local repo with the one upstream fixed the issue..
sorry for the noise :)
Giacomo Acciarini
@Sceki
I tried to make doctest on the sphinx folder, but I get the following error:
Running Sphinx v1.6.7
making output directory...

Configuration error:
There is a programable error in your configuration file:

Traceback (most recent call last):
  File "/Users/giacomoacciarini/anaconda3/lib/python3.7/site-packages/sphinx/config.py", line 157, in __init__
    execfile_(filename, config)
  File "/Users/giacomoacciarini/anaconda3/lib/python3.7/site-packages/sphinx/util/pycompat.py", line 150, in execfile_
    exec_(code, _globals)
  File "conf.py", line 92, in <module>
    from jupyter_sphinx_theme import *
ModuleNotFoundError: No module named 'jupyter_sphinx_theme'

make: *** [doctest] Error 1
the tutorial is ready (i.e., I made the .ipynb file and I added the relevant line in the index.rst).. I could also simply add it via PR, if this error is some configuration one
Dario Izzo
@darioizzo
@Sceki you need to install the jupyter_sphinx_theme
pip install jupyter_sphinx_theme
maybe you also need other deps .. if so just install them via pip.
Giacomo Acciarini
@Sceki
okay!
what about the .html files? how can I produce them locally?
Dario Izzo
@darioizzo
I saw your PR and the notebook, ... in general its good, when possible to avoid just one huge code segment
And break it apart in a few bits commenting inbetween (if possible)
Giacomo Acciarini
@Sceki
I thought it would have been a bit ugly but there are two nested for loops
Dario Izzo
@darioizzo
sure ... maybe inavoidable ...

what about the .html files? how can I produce them locally?

in what sense?

make html
Giacomo Acciarini
@Sceki
okay!
I used to do sh compile.sh in pagmo
Dario Izzo
@darioizzo
oh right ... you can also make the script and add it to the PR
so that its uniform acreoss the projects
do not use doctest as I do not think we use it in pykep
Giacomo Acciarini
@Sceki
interesting.. I am gonna check what's inside the compile.sh file of pagmo
Dario Izzo
@darioizzo
probably just make html :)
Giacomo Acciarini
@Sceki
btw yesterday I tried to install jupyter_sphinx_theme via conda but it was not there
I will do it via pip
Dario Izzo
@darioizzo
no its in pip
for some reason
Giacomo Acciarini
@Sceki
I was afraid of mixing up the envs ahah
Dario Izzo
@darioizzo
true, in this case its allowed