uhrwecker / GRDonuts

Toolkit for creating plots and computing stable equipotential surfaces of fluids orbiting a compact stellar object (see Zarotti: Relativistic Hydrodynamics))
MIT License
0 stars 0 forks source link

Implement the search for closed configurations #18

Closed uhrwecker closed 4 years ago

uhrwecker commented 4 years ago

Implement a function to the tuility class of the potential in potential.py which allows the user to specify a range of one or many parameters, which then prompts a list of possible combinations when the potential is closed. Possibly print it on console/log it/save it to .csv file.

uhrwecker commented 4 years ago

Code idea/snippet:

def try_stable_point(q, l): try: k = QMetric(q=q, l=l, cont_without_eq=True) k.compute_w() print('!!! Stable at q={} and l={} !!!'.format(q, l)) except: pass

print('Not stable at q={}.'.format(q))

for q in np.linspace(1, 2, num=100): for l in np.linspace(1, 10, num=100): try_stable_point(q, l)

uhrwecker commented 4 years ago

Done. Merging.