ASTROPHOT: fitting everything everywhere all at once in astronomical images

We present ASTROPHOT, a fast, powerful, and user-friendly PYTHON based astronomical image photometry solver. ASTROPHOT incorporates automatic differentiation and graphics processing unit (GPU), or parallel central processing unit (CPU), acceleration, powered by the machine learning library PYTORCH. Everything: ASTROPHOT can fit models for sky, stars, galaxies, point spread functions (PSFs), and more in a principled χ² forward optimization, recovering Bayesian posterior information and covariance of all parameters. Everywhere: ASTROPHOT can optimize forward models on CPU or GPU; across images that are large, multiband, multi-epoch, rotated, dithered, and more. All at once: The models are optimized together, thus handling overlapping objects and including the covariance between parameters (including PSF and galaxy parameters). A number of optimization algorithms are available including Levenberg-Marquardt, Gradient descent, and No-U-Turn Markov chain Monte Carlo sampling. With an object-oriented user interface, ASTROPHOT makes it easy to quickly extract detailed information from complex astronomical data for individual images or large survey programs. This paper outlines novel features of the ASTROPHOT code and compares it to other popular astronomical image modelling software. ASTROPHOT is open-source, fully PYTHON based, and freely accessible at https://github.com/Autostronomy/AstroPhot .