astropy VS NumPy

Astropy [0] lives at the heart of most work. It has a Python interface, often backed by Fortran and C++ extension modules. If you use Astropy, you're indirectly using libraries like ERFA [6] and cfitsio [7] which are in C/Fortran.

I personally end up doing a lot of work that uses the HEALPix sky tesselation, so I use healpy [2] as well.

Openorb is perhaps a good example of a pure-Fortran package that I use quite. frequently for orbit propagation [3].

In C, there's Rebound [4] (for N-body simulations) and ASSIST [5] (which extends Rebound to use JPL's pre-calculated positions of major perturbers, and expands the force model to account for general relativity).

There are many more, these are just ones that come to mind from frequent usage in the last few months.

[0] https://www.astropy.org/

Users interested in a broader range of astronomical tools beyond coordinate transformations may be interested in https://www.astropy.org/ and its affiliated packages.

I'm working on a pure Rust approximation of astropy. Up til now, I was able to recreate the intent by looking at an external API, but I'm moving on to functionality that I don't understand enough to implement without basically copying the code. Astropy uses the BSD-3 license, and it wraps the ERFA library which uses a custom license. My project currently uses the MIT license. My PR is here - my question is have I attributed everything correctly, or is there anything I need to change for everything to be above-board?

For folks who might be interested in astronomical calculations but who don't want to roll their own library, astropy (https://www.astropy.org/) is widely used by professional astronomers.

I agree with the other commenter. Check out their github. If you’re looking to build your skills long term (and have some experience with python) it’s worth checking out astropy and their fits file handling routines.

In NumPy with @, dot() or matmul():

In NumPy with * or multiply(). ` or multiply()` can multiply 0D or more D arrays by element-wise multiplication.

Data science projects often use numpy. However, numpy objects are not JSON-serializable and therefore require conversion to standard python objects in order to be saved:

numpy: A library for scientific computing in Python

Python has become a preferred language for data analysis due to its simplicity and robust library ecosystem. Among these, NumPy stands out with its efficient handling of numerical data. Let’s say you’re working with numbers for large data sets—something Python’s native data structures may find challenging. That’s where NumPy arrays come into play, making numerical computations seamless and speedy.

I was looking at all the big projects like numpy, pytorch, flask, etc.

I wish the NumPy devs would more thoroughly consider adding full fluent API support, e.g. x.sqrt().ceil(). [Issue #24081]