Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

one loop is superseded with fortran subroutine #1

Open
wants to merge 97 commits into
base: master
Choose a base branch
from

Conversation

belliavesha
Copy link
Owner

one loop is superseded with fortran subroutine

belliavesha and others added 30 commits April 13, 2018 18:51
one loop is superseded with fortran subroutine
uploaded
Mainly, the code of the angles transformation, and time lags and other connected things was implemented


in the Redistribution matrix an annoying doneness string added
and other small corrections
Some bug fixed, maybe
some stuff with output files changed
Thomson scattering model implemented.
bug in the computation of redistribution matrix is fixed. again.
small changes
fix the typos
reductant term in the Schwarzschild routine is deleted. But It does not seem to make the weather
polarization angle transformation description
Geometry of the problem
Making OblattePol more article like
High energies cut off in Thomson approximation
added the cubic spline along the mu axis instead of just linear interpolation
Several bugs, noted by Tuomo. 
- extra "wt" term deleted from Schwarzschild
- missing 2 added to the R derivative
- way of phase interpolation is changed
Fit PA with emcee and plot tools
rho is the radius of the spot
also a parameter 'antipodal' for two antipodal spot or only one spot if False
oblateness factor added in the dS term
Tuomo Salmi and others added 30 commits July 31, 2020 01:10
bug fix fix
merge with sigma.jl,
bug fixes
If we want to stay within 1% difference from "maximum values" scenario (250 energy points / 20 mu-s and optical depth 60), we should stick to 150 points, 9 Mu-s and optical depth 50.
ScatterNum is reduced to 12, originally 18. Calculation time is still enormous though
T_e comes from the cluster task index
grid is the one suggested by Tuomo, small ("test") one 
new HDU for each combination of T_bb and tau_T, new fits file for each T_e. this can be rearranged within the code or later on
all great but still 8 hours of calculations for 1 T_e-T_bb-tau_T combination with 12 orders, 150 NEnergies, 9 NMu and tau 50
probably CRM can further be accelerated
…at is used in ixpe_analysis at the moment (with a fixed Thomson atmosphere), and new options for analytical light bending / time delay formulas
Latest pulse profile modelling codes and plot tools updated
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants