| Commit message (Collapse) | Author | Age | Files | Lines |
... | |
|
|
|
| |
Investigate whether the vectorization will improve the speed.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Merge PR#3: Add new simulation component "extragalactic/pointsources", including SF, SB, RQ, FRI and FRII.
* extragalactic/pointsource: add point source simulation module
* extracgalactic/pointsource: add point source simulation module
* Add configurating spec to extragalgactic point sources.
* Modified some variables
* base.py: modified
* flux.py: modified
* fr1.py: modified
* fr2.py: modified
* pointsources.py: modified
* psparams.py: modified
* radioquiet.py: modified
* starforming.py: modified
* starbursting.py: modified
* Rewritten the comments.
* base.py: modified
* flux.py: modified
* fr1.py: modified
* fr2.py: modified
* Modified
* psparams.py: modified
* radioquiet.py: modified
* starforming.py: modified
* starbursting.py: modified
* Modified
* Modified
* Modified
* Modified
* Modified
* MOdified
* Modified
* Modified
* Modified
* Modified
* Modified
* Modified
* Modified
* Modified
* Modified
* Changed pointsource to pointsources
* Fixed some config keywords
* Fixed some config keywords
* Fixed some config keywords
* Fixed some config keywords
* Fixed some config keywords
* Fixed some config keywords
* Fixed some config keywords
* base.py:rewrited
* fr1.py:rewritten
* pointsources.py: rewritten
* radioquiet.py: rewritten
* Rewritten
* Rewritten
* base.py: modified
* fr1.py: modified
* fr2.py: modified
* radioquiet.py: modified
* starbursting.py: modified
* starforming.py: modified
* Fix conflicts
* fg21sim: fixed conflicts
* base.py: modified frequencies loading in _get_configs()
* Rewritten as forground.py
* fg21sim: fixed conflicts
* base.py: deteled loading for frequencies configurations.
* fr1.py: modified
* fr2.py: modified
* pointsources.py: modified
* radioquiet.py: modified
* starbursting.py: modified
* starforming.py: modified
* Add new methods to calculate Tb.
* Add new methods to calculate Tb.
* Add new methods to calculate Tb.
* Add new methods to calculate Tb.
* Add new methods to calculate Tb.
* Add new methods to calculate Tb.
* Deleted useless comments.
* Add pscomps to deal with multi-type PS problem.
* Add a new key.
* Fixed permission to 755.
* Rejusted PS subsections.
* Add methods to calcualte luminosity function and redshift distribution.
* Rejusted generation of samples redshift and luminosity.
* Fixed mistakes on FRII structure, added hotspots and offsets.
* Reajusted generation of samples redshift and luminosity.
* Readujsted generation of samples radii, redshifts and luminosity.
* Readujsted generation of samples radii, redshifts and luminosity.
* Fixed conflicts.
* Fixed conficts.
* Combined configurations of pointsources.
* Removed the older extragalactic configuration file.
* Fixed some mistakes.
* Fixed mistakes of drawing PS.
* Fixed mistakes of drawing PS.
* Fixed code style by pep8 checking.
* Fixed code style by pep8 checking.
* Fixed code style by pep8 checking.
* Fixed some coding style.
* Reconfigured default redshift interval.
* Reconfigured default redshift interval.
* Reconfigured default redshift interval.
* Fixed mistakes in method calc_single_Tb and changed resolution of grid.
* Fixed mistakes in method calc_single_Tb and changed resolution of grid.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Deleted astropy.units style code to accelerate.
* Fixed some mistakes.
* Fixed some mistakes.
* Changed dA from au.Mpc to float64.
* Fixed some mistakes.
* Fixed some mistakes.
* Fixed some mistakes.
* Reajusted grid resolution to generate discs.
* Reajusted grid resolution to generate discs.
* Reajusted loading strategy of parameter resolution.
* Reajusted code style of configuration loading.
* Reajusted code style of configuration loading.
* Reajusted code style of configuration loading.
* Reajusted code style of configuration loading.
* Reajusted code style of configuration loading.
|
|
|
|
|
|
|
| |
Import from ".foregrounds" involves a lot of imports of other
components, which cause a lot of numba JIT compilations.
Therefore, the simple "import fg21sim" causes a long lag for running the
scripts, e.g., bin/healpix2hpx.
|
|
|
|
|
| |
Also add a warning log message for "galactic/snr.py" and
"extragalactic/clusters.py".
|
|
|
|
| |
Also update the docstring and one comment.
|
| |
|
|
|
|
|
|
|
| |
* Fix the wrong output shape: wrongly swap the row and column
* Fix the area mapping method, which may give zero values when the index
of the input is integer, which cause problem to calculate the
overlapping areas.
|
|
|
|
| |
NOTE that the original optional paramter "rotation" is now required.
|
|
|
|
| |
Also remove the dependence on "scipy".
|
|
|
|
|
|
|
|
|
| |
The "rotate_center()" function rotates a gray-scale image by a given
angle about its center, using the area mapping method or simple nearest
neighbor.
This function is `numba` JIT-optimized, and can replace the
`scipy.ndimage.rotate()` function.
|
| |
|
|
|
|
|
|
|
|
|
|
| |
NOTE: However, there are two significant problems:
* Cannot figure out why ``nb.optional(nb.types.UniTuple(nb.int64, 2))``
does NOT work. Therefore, make ``shape`` as mandatory parameter
instead of optional.
* Cannot figure out multi-dispatch that allows both int and float types
for ``r``, ``c``, ``r_radius`` and ``c_radius``. Thus only support
the int type for the moment.
|
| |
|
|
|
|
|
| |
Remove the original "make_coordinate_grid()" function, and rename the
"make_coordinate_grid_fast()" to override the former one.
|
|
|
|
|
| |
Change the parameters of "make_coordinate_grid_fast()", make it share
the same parameter list as "make_coordinate_grid()".
|
|
|
|
|
| |
There are some limitations with the numba JIT, which are commented in
the code.
|
|
|
|
| |
Also fix a wrong notation in "random.py"
|
|
|
|
|
|
|
|
|
|
| |
To optimize "map_grid_to_healpix()" in `grid.py` which uses
"healpy.ang2pix()" using `numba.jit`, implement the latter's
functionality explicitly with `numba.jit` support.
The implementation simply mimic the corresponding functions
"ang2pix_ring_z_phi()" and "pix2ang_ring_z_phi()" in HEALPix's
"src/C/subs/chealpix.c". Thanks!
|
|
|
|
|
| |
Also place "_calc_hpx_row_idx()" before "_calc_hpx_indicies()", which is
required by `numba`.
|
|
|
|
|
|
|
|
|
| |
Avoid using "astropy.units" in "calc_Tb()", which significantly speeds
up this method.
NOTE:
whether the distance given by the simulation needs additional conversion
to derive the *luminosity distance* ??
|
|
|
|
|
|
|
|
| |
Optimize "make_coordinate_grid()" using "numba.jit", which improves the
speed a bit, since the other two functions (i.e., "make_grid_ellipse()",
and "map_grid_to_healpix()") are much slow and time-consuming.
TODO: optimize "make_grid_ellipse()" and "map_grid_to_healpix()"
|
| |
|
|
|
|
|
|
|
|
|
| |
* Optimize the "Fnu_to_Tb()" and "Sb_to_Tb()" functions by explicitly
calculating the conversions, avoiding the slow `astropy.unit`
conversions (which is rather slow). The new fast functions are named
as "Fnu_to_Tb_fast()" and "Sb_to_Tb_fast()".
* Optimize the new functions using `numba.jit` further.
* Add `numba` to the dependencies.
|
|
|
|
|
|
| |
* Fix several bugs, e.g., typo, wrong quantity operation, etc.
* Add new method "_calc_specindex()", also allow further improvement
* Improve docstrings, comments, and log messages.
|
| |
|
|
|
|
|
|
|
|
|
|
|
| |
Support simulate the radio emission from clusters of galaxies.
Currently only the radio halos is considered, while the radio relics
need more investigations.
Also, this simulation requires a cluster catalog derived from the Hubble
Volume Project simulation.
NOTE: not tested at the moment
|
|
|
|
|
| |
Currently implement "spherical_uniform()" to generate random points
uniformly distributed on the spherical surface.
|
|
|
|
|
|
|
| |
When two quantities of same unit class divides, to get the correct
numerical value, whether use `.decompose().value`, or just convert to
plain Python object, e.g., `float(...)`, which will fail if the input
quantity is not dimensionless.
|
| |
|
| |
|
| |
|
| |
|
|
|
|
| |
All HEALPix maps used are in *RING* ordering.
|
|
|
|
|
| |
* freefree.py: Add the missing "/" to the config "halphamap"
* synchrotron.py: Fix "comp" to have the correct value
|
|
|
|
|
|
| |
* Fix the wrong import path
* Add the missing "output_dir" property
* Tweak the comments/messages a bit
|
|
|
|
|
| |
This module provides a simple/easy-to-use interface to the simulation of
various supported foreground components.
|
|
|
|
|
|
| |
* Add the "preprocess()" and "postprocess()" for a better interface;
* Add "_make_filepath()" to compose the output filename;
* Reorder some methods for consistency.
|
|
|
|
|
|
|
| |
* Add the "preprocess()" and "postprocess()" for a better interface;
* Add "_make_filepath()" to compose the output filename;
* Reorder some methods for consistency;
* Some refactorings and cleanups.
|
|
|
|
|
|
|
| |
* Add "preprocess()" and "postprocess()" methods
* Add the "name" property
* Tweak the places of preparation procedures been invoked
* Reorder some methods
|
| |
|
| |
|
| |
|
|
|
|
| |
Barely test and seems OK.
|
|
|
|
|
|
|
| |
The "filter()" returns an iterator, that is consumed after the use,
therefore, the "specfiles" is used by the print under "DEBUG_FG21SIM".
Then "specfiles" is empty and thus prevents the reading of the
specification files.
|
| |
|
|
|
|
| |
Also update the general config for the new "galactic/snr" component.
|
|
|
|
|
|
|
|
| |
The necessary but missing functionalities to simulate the Galactic SNRs
emission maps are implemented, and this new emission component is ready
for testing.
Also fix a typo in "utils/grid.py".
|
|
|
|
|
|
|
|
| |
The "map_grid_to_healpix()" maps the generated coordinate grid to a
HEALPix map. Note that only effective HEALPix pixels are returned
instead of a full HEALPix map.
TODO: Try to avoid the explicit for loop to optimize the speed.
|
|
|
|
|
|
|
| |
* make_coordinate_grid(): Make a Cartesian coordinate grid according to
the specified center coordinate, size and grid resolution.
* make_grid_ellipse(): Also generate an (rotated) ellipse and place into
the coordinate grid.
|
|
|
|
| |
NOTE: very preliminary, and missing many necessary functions!
|