search

bgmeulem / Luminet

This repo provides material for recreating the famous paper by Jean-Pierre Luminet (1979), simulating the first image of a Schwarzschild black hole.
MIT License
63 stars 5 forks source link

how to change inclination? #7

Open Specom opened 2 weeks ago

Specom commented 2 weeks ago

I really like this project. I run the code in anaconda, however, I somehow can not change the inclination. Do I need to run the code in jupiter notebook instead? If I rember right, the inclination from the famouse image was 10 degree. Could you give an example for the command line for this inclination?

Thanks a lot!

bgmeulem commented 6 days ago

Hi, thanks for asking. Normally, the inclination is set upon initializing either the black hole class or the isoradial class. If you want an inclination of 10 degrees, you can do something like 

bh = BlackHole(inclination=10, mass=1)
bh.plot_isoradials([10, 20, 30], [10, 20, 30])
# or to sample photons 
bh.sample_points(n_points=20000)

Let me know if that works for you