From 7612417495606c492582137f2e1de4aba2a142d2 Mon Sep 17 00:00:00 2001 From: Alejandro Cardenas Avendano Date: Wed, 24 Aug 2022 10:23:42 -0400 Subject: [PATCH] More comments added and csc notation and safe value --- aart_func/ellipi_f.py | 7 ++++--- params.py | 16 +++++++++------- 2 files changed, 13 insertions(+), 10 deletions(-) diff --git a/aart_func/ellipi_f.py b/aart_func/ellipi_f.py index 57dce22..48d9a52 100644 --- a/aart_func/ellipi_f.py +++ b/aart_func/ellipi_f.py @@ -1,9 +1,10 @@ from aart_func import * from params import * -def ellipi(n,phi,m): - csc=1/np.sin(phi)**2 - return 1/3*n*elliprj(csc-1,csc-m,csc,csc-n)+ellipf(phi,m) +#Legendre’s Integrals as Symmetric Integrals https://dlmf.nist.gov/19.25#i +def ellipi(n,phi,m): + csc2=1/(np.sin(phi)**2+1e-16) + return 1/3*n*elliprj(csc2-1,csc2-m,csc2,csc2-n)+ellipf(phi,m) def ellippi_lim(n,phi,m): #take arrays n, m, phi can either be a float or an array diff --git a/params.py b/params.py index 3defc3c..b368886 100644 --- a/params.py +++ b/params.py @@ -7,9 +7,6 @@ spin_case=0.94 #Observer's inclination i_case=17 -#Disk's inclination -i_disk=90 - # Distance to M87 in meters dM87=5.214795112e23 @@ -58,15 +55,16 @@ #Number of snapshots in that range snapshots=12 -# SU's parameters for the envelope -# Just used for the profiles computed within AART isco = rms(spin_case) - +# SU's parameters for the envelope +# Just used for the analytical profiles gammap=-3/2 mup=1-sqrt(1-spin_case**2) sigmap=1/2 -speed_p=1 +# Sarting point of the profile's cutoff cutoff_p=limits-1 +#Speed/slope of the cutoff +speed_p=1 #The power of the redshift factor gfactor=3 @@ -80,7 +78,11 @@ #For the parallel generation of images (movies) nthreads=4 +#Observer's inclination in radians thetao=i_case*np.pi/180 +#Disk's inclination +#Current version is just implemented for equatorial models +i_disk=90 thetad=i_disk*np.pi/180 Gc=6.67e-11 # G constant [m^3 kg^-1 s^-2]