diff --git a/.gitignore b/.gitignore index 18901a88..bb0d57e0 100644 --- a/.gitignore +++ b/.gitignore @@ -65,3 +65,5 @@ target/ .idea/ .idea/* .vscode/settings.json + +sedkit/data/models/atmospheric/spexprismlibrary/index.p diff --git a/CITATION.cff b/CITATION.cff new file mode 100644 index 00000000..1012d77e --- /dev/null +++ b/CITATION.cff @@ -0,0 +1,17 @@ +authors: +- family-names: Filippazzo + given-names: Joseph + orcid: https://orcid.org/0000-0002-0201-8306 +- family-names: Alejandro Merchan + given-names: Sherelyn + orcid: https://orcid.org/0000-0003-0548-0093 +- family-names: Cruz + given-names: Kelle + orcid: https://orcid.org/0000-0002-1821-0650 +cff-version: 1.2.0 +doi: 10.5281/zenodo.14014658 +repository-code: https://github.com/BDNYC/sedkit/ +title: 'sedkit' +type: software +message: "If you use this software, please cite the Zenodo record" +license: BSD-3-Clause \ No newline at end of file diff --git a/LICENSE b/LICENSE index 42b35567..c1e4dd5a 100644 --- a/LICENSE +++ b/LICENSE @@ -1,21 +1,11 @@ -The MIT License (MIT) +Copyright (c) 2024, Joe Filippazzo, Sherelyn Alejandro Merchan, Kelle Cruz -Copyright (c) 2019 Joe Filippazzo +Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met: -Permission is hereby granted, free of charge, to any person obtaining a copy -of this software and associated documentation files (the "Software"), to deal -in the Software without restriction, including without limitation the rights -to use, copy, modify, merge, publish, distribute, sublicense, and/or sell -copies of the Software, and to permit persons to whom the Software is -furnished to do so, subject to the following conditions: +Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer. -The above copyright notice and this permission notice shall be included in all -copies or substantial portions of the Software. +Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution. -THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR -IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, -FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE -AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER -LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, -OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE -SOFTWARE. +Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission. + +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. \ No newline at end of file diff --git a/README.md b/README.md index 7ba62bec..28f9a01c 100644 --- a/README.md +++ b/README.md @@ -3,27 +3,27 @@ [![Build Status](https://github.com/hover2pi/sedkit/actions/workflows/ci.yml/badge.svg)](https://github.com/hover2pi/sedkit/actions/workflows/ci.yml) [![Coverage Status](https://coveralls.io/repos/github/hover2pi/sedkit/badge.svg?branch=master&service=github)](https://coveralls.io/github/hover2pi/sedkit?branch=master) [![Documentation Status](https://readthedocs.org/projects/sedkit/badge/?version=latest)](https://sedkit.readthedocs.io/en/latest/?badge=latest) +[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.14014658.svg)](https://doi.org/10.5281/zenodo.14014658) + ## Introduction -`sedkit` is a collection of pure Python modules for simple SED construction and analysis. Users can create individual SEDs or SED catalogs from spectra and/or photometry and calculate fundamental parameters (fbol, Mbol, Lbol, Teff, mass, log(g)) using the methods presented in [Filippazzo et al. (2015)](http://adsabs.harvard.edu/abs/2015ApJ...810..158F). +`sedkit` is a collection of pure Python modules for simple SED construction and analysis. Users can create individual SEDs or SED catalogs from spectra and/or photometry and calculate fundamental parameters (fbol, Mbol, Lbol, Teff, mass, log(g)) using the methods presented in [Filippazzo et al. (2015)](http://adsabs.harvard.edu/abs/2015ApJ...810..158F) and described in the documentation: https://sedkit.readthedocs.io/. + ## Installation -Install via PyPI with +Install in an environment with Python 3.11 or 3.12. ``` pip install sedkit ``` +### Developer Instructions -or via `conda` with +If you plan to modify and/or contribute to the code, clone the repository and then install an editable version of the package and the requirements for running the tests. ``` -git clone https://github.com/hover2pi/sedkit.git -cd sedkit -conda env create -f env/environment-3.11.yml --force -conda activate sedkit -python setup.py install +pip install -e ".[test]" ``` ## Demo @@ -58,7 +58,7 @@ Spectrum arrays or ASCII/FITS files can also be added to the SED data. from pkg_resources import resource_filename spec_file = resource_filename('sedkit', 'data/Trappist-1_NIR.fits') import astropy.units as u -trap1.add_spectrum_file(spec_file, wave_units=u.um, flux_units=u.erg/u.s/q.cm**2/u.AA) +trap1.add_spectrum_file(spec_file, wave_units=u.um, flux_units=u.erg/u.s/u.cm**2/u.AA) ``` Other data which may affect the calculated and inferred fundamantal parameters can be set at any time. @@ -201,20 +201,19 @@ trap1.evo_model = 'DUSTY00' trap1.infer_mass() ``` - - -A variety of atmospheric model grids can be fit to the data with or without MCMC analysis, + + -And any arbitrary atlas of models can be applied as well. +Any arbitrary atlas of models (stored in `data\models\atmospheric`) can be used to fit the data with or without MCMC analysis. ```python from sedkit import SpexPrismLibrary -trap1.fit_modelgrid(SpexPrismLibrary()) +trap1.fit_modelgrid(SpexPrismLibrary(), mcmc=True) ``` Inspect the SED at any time with the interactive plotting method. @@ -223,28 +222,22 @@ Inspect the SED at any time with the interactive plotting method. trap1.plot(integral=True, best_fit=True) ``` - + References for all data can be accessed via the `refs` attribute. Entire catalogs of `SED` objects can also be created and their properties can be arbitrarily compared and analyzed with the `sedkit.catalog.Catalog()` object. - + -Please read the full documentation for details on this functionality and much more. +Please read the [full documentation](https://sedkit.readthedocs.io/en/latest/) for details on this functionality and much more. -## Documentation +## Documentation and Notebooks -Full documentation for the latest build can be found on [ReadTheDocs](https://sedkit.readthedocs.io/en/latest/). +Full documentation is available: https://sedkit.readthedocs.io/. The package also contains detailed Jupyter notebooks highlighting the core functionality of its primary classes, including -- [sedkit.spectrum.Spectrum](https://github.com/hover2pi/sedkit/blob/master/sedkit/notebooks/working_with_spectra.ipynb) -- [sedkit.sed.SED](https://github.com/hover2pi/sedkit/blob/master/sedkit/notebooks/create_sed.ipynb) -- [sedkit.catalog.Catalog](https://github.com/hover2pi/sedkit/blob/master/sedkit/notebooks/create_catalog.ipynb) - -If you use or reference this software, please cite [Filippazzo et al. (submitted to PASP)]() - -## Licensed - -This project is Copyright (c) Joe Filippazzo and licensed under the terms of the BSD 3-Clause license. See the licenses folder for more information. +- [sedkit.spectrum.Spectrum](sedkit/notebooks/working_with_spectra.ipynb) +- [sedkit.sed.SED](sedkit/notebooks/create_sed.ipynb) +- [sedkit.catalog.Catalog](sedkit/notebooks/create_catalog.ipynb) diff --git a/pyproject.toml b/pyproject.toml index 20d61799..a9a602e7 100644 --- a/pyproject.toml +++ b/pyproject.toml @@ -29,7 +29,7 @@ dependencies = [ "dill>=0.3.4", "dustmaps>=1.0.9", "emcee>=3.1.1", - "numpy>=1.25.1", + "numpy>=2.1.0", "pandas>=1.3.5", "scipy>=1.8.0", "svo-filters>=0.4.4", diff --git a/sedkit/data/figures/Lbol_v_SpT.png b/sedkit/data/figures/Lbol_v_SpT.png index e69de29b..92fd5786 100644 Binary files a/sedkit/data/figures/Lbol_v_SpT.png and b/sedkit/data/figures/Lbol_v_SpT.png differ diff --git a/sedkit/data/figures/Lbol_v_mass.png b/sedkit/data/figures/Lbol_v_mass.png index e69de29b..f39bc31c 100644 Binary files a/sedkit/data/figures/Lbol_v_mass.png and b/sedkit/data/figures/Lbol_v_mass.png differ diff --git a/sedkit/data/figures/sed_plot.png b/sedkit/data/figures/sed_plot.png index e69de29b..aa07352d 100644 Binary files a/sedkit/data/figures/sed_plot.png and b/sedkit/data/figures/sed_plot.png differ diff --git a/sedkit/tests/test_isochrone.py b/sedkit/tests/test_isochrone.py index aec4c105..e9185c20 100644 --- a/sedkit/tests/test_isochrone.py +++ b/sedkit/tests/test_isochrone.py @@ -26,7 +26,7 @@ (-4, 4 * q.Gyr, 'Lbol', 'logg', 5.395, 5.36, 5.43), # No xparam and yparam uncertaintiesd # Young age ((-4, 0.1), (0.4 * q.Gyr, 0.1 * q.Gyr), 'Lbol', 'mass', 0.0515, 0.045, 0.055), # mass with uncertainties - ((-4, 0.1), (0.4 * q.Gyr, 0.1 * q.Gyr), 'Lbol', 'logg', 5.08, 5.01, 5.15), # logg with uncertainties + ((-4, 0.1), (0.4 * q.Gyr, 0.1 * q.Gyr), 'Lbol', 'logg', 5.085, 5.01, 5.15), # logg with uncertainties ((-4, 0.1), (0.4 * q.Gyr, 0.1 * q.Gyr), 'Lbol', 'radius', 0.105, 0.10, 0.11), # radius with uncertainties # Old age ((-4, 0.1), (9 * q.Gyr, 0.1 * q.Gyr), 'Lbol', 'mass', 0.074, 0.070, 0.080), # mass with uncertainties diff --git a/source/index.rst b/source/index.rst index 7ef5f53b..18d7cca5 100644 --- a/source/index.rst +++ b/source/index.rst @@ -16,15 +16,6 @@ Install via PyPI with pip install sedkit -or via ``conda`` with - -.. code:: - - git clone https://github.com/hover2pi/sedkit.git - cd sedkit - conda env create -f environment.yml --force - conda activate sedkit - python setup.py install Contents ======== diff --git a/test.txt b/test.txt deleted file mode 100644 index 64ca14b4..00000000 --- a/test.txt +++ /dev/null @@ -1,233 +0,0 @@ -# New_Spectrum -# Wavelength [um], Flux Density [erg / (Angstrom s cm2)] -# Foo -6.000000238418579102e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.017467379570007324e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.034934520721435547e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.052401661872863770e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.069868803024291992e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.087336540222167969e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.104803681373596191e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.122270822525024414e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.139737963676452637e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.157205104827880859e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.174672245979309082e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.192139983177185059e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.209607124328613281e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.227074265480041504e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.244541406631469727e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.262008547782897949e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.279475688934326172e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.296943426132202148e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.314410567283630371e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.331877708435058594e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.349344849586486816e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.366811990737915039e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.384279727935791016e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.401746869087219238e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.419214010238647461e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.436681151390075684e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.454148292541503906e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.471615433692932129e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.489083170890808105e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.506550312042236328e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.524017453193664551e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.541484594345092773e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.558951735496520996e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.576419472694396973e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.593886613845825195e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.611353754997253418e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.628820896148681641e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.646288037300109863e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.663755178451538086e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.681222915649414062e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.698690056800842285e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.716157197952270508e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.733624339103698730e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.751091480255126953e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.768559217453002930e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.786026358604431152e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.803493499755859375e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.820960640907287598e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.838427782058715820e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.855894923210144043e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.873362660408020020e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.890829801559448242e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.908296942710876465e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.925764083862304688e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.943231225013732910e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.960698962211608887e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.978166103363037109e-01 1.000000000000000000e+00 9.999999776482582092e-03 -6.995633244514465332e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.013100385665893555e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.030567526817321777e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.048034667968750000e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.065502405166625977e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.082969546318054199e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.100436687469482422e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.117903828620910645e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.135370969772338867e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.152838706970214844e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.170305848121643066e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.187772989273071289e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.205240130424499512e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.222707271575927734e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.240174412727355957e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.257642149925231934e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.275109291076660156e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.292576432228088379e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.310043573379516602e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.327510714530944824e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.344978451728820801e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.362445592880249023e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.379912734031677246e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.397379875183105469e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.414847016334533691e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.432314157485961914e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.449781894683837891e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.467249035835266113e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.484716176986694336e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.502183318138122559e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.519650459289550781e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.537118196487426758e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.554585337638854980e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.572052478790283203e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.589519619941711426e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.606986761093139648e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.624453902244567871e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.641921639442443848e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.659388780593872070e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.676855921745300293e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.694323062896728516e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.711790204048156738e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.729257345199584961e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.746725082397460938e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.764192223548889160e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.781659364700317383e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.799126505851745605e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.816593647003173828e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.834061384201049805e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.851528525352478027e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.868995666503906250e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.886462807655334473e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.903929948806762695e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.921397089958190918e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.938864827156066895e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.956331968307495117e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.973799109458923340e-01 1.000000000000000000e+00 9.999999776482582092e-03 -7.991266250610351562e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.008733391761779785e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.026201128959655762e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.043668270111083984e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.061135411262512207e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.078602552413940430e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.096069693565368652e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.113536834716796875e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.131004571914672852e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.148471713066101074e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.165938854217529297e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.183405995368957520e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.200873136520385742e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.218340873718261719e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.235808014869689941e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.253275156021118164e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.270742297172546387e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.288209438323974609e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.305676579475402832e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.323144316673278809e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.340611457824707031e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.358078598976135254e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.375545740127563477e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.393012881278991699e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.410480618476867676e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.427947759628295898e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.445414900779724121e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.462882041931152344e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.480349183082580566e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.497816324234008789e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.515284061431884766e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.532751202583312988e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.550218343734741211e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.567685484886169434e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.585152626037597656e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.602620363235473633e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.620087504386901855e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.637554645538330078e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.655021786689758301e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.672488927841186523e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.689956068992614746e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.707423806190490723e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.724890947341918945e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.742358088493347168e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.759825229644775391e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.777292370796203613e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.794760107994079590e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.812227249145507812e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.829694390296936035e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.847161531448364258e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.864628672599792480e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.882095813751220703e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.899563550949096680e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.917030692100524902e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.934497833251953125e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.951964974403381348e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.969432115554809570e-01 1.000000000000000000e+00 9.999999776482582092e-03 -8.986899852752685547e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.004366993904113770e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.021834135055541992e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.039301276206970215e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.056768417358398438e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.074235558509826660e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.091703295707702637e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.109170436859130859e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.126637578010559082e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.144104719161987305e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.161571860313415527e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.179039597511291504e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.196506738662719727e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.213973879814147949e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.231441020965576172e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.248908162117004395e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.266375303268432617e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.283843040466308594e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.301310181617736816e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.318777322769165039e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.336244463920593262e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.353711605072021484e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.371178746223449707e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.388646483421325684e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.406113624572753906e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.423580765724182129e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.441047906875610352e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.458515048027038574e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.475982785224914551e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.493449926376342773e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.510917067527770996e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.528384208679199219e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.545851349830627441e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.563318490982055664e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.580786228179931641e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.598253369331359863e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.615720510482788086e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.633187651634216309e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.650654792785644531e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.668122529983520508e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.685589671134948730e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.703056812286376953e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.720523953437805176e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.737991094589233398e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.755458235740661621e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.772925972938537598e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.790393114089965820e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.807860255241394043e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.825327396392822266e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.842794537544250488e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.860262274742126465e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.877729415893554688e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.895196557044982910e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.912663698196411133e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.930130839347839355e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.947597980499267578e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.965065717697143555e-01 1.000000000000000000e+00 9.999999776482582092e-03 -9.982532858848571777e-01 1.000000000000000000e+00 9.999999776482582092e-03 -1.000000000000000000e+00 1.000000000000000000e+00 9.999999776482582092e-03