Related papers: Fast computation of temperature and polarization c…
We have developed a fast method for predicting the angular power spectrum, C_l, of the cosmic microwave background given cosmological parameters and a primordial power spectrum of perturbations. After pre--computing the radiation…
Matrix multiplication is a foundational operation in scientific computing and machine learning, yet its computational complexity makes it a significant bottleneck for large-scale applications. The shift to parallel architectures, primarily…
A very common problem in science is the numerical diagonalization of symmetric or hermitian 3x3 matrices. Since standard "black box" packages may be too inefficient if the number of matrices is large, we study several alternatives. We…
This paper generalises the hybrid power spectrum estimator developed in Efstathiou (2004a) to the estimation of polarization power spectra of the cosmic microwave background radiation. The hybrid power spectrum estimator is unbiased and we…
An estimation of the sky signal from streams of Time Ordered Data (TOD) acquired by Cosmic Microwave Background (\cmb) experiments is one of the most important steps in the context of \cmb data analysis referred to as the map-making…
Accurate and efficient thermal dynamics models of permanent magnet synchronous motors are vital to efficient thermal management strategies. Physics-informed methods combine model-based and data-driven methods, offering greater flexibility…
Observations of the polarization of the cosmic microwave backround (CMB) have the potential to place much tighter constraints on cosmological parameters than observations of the fluctuations in temperature alone. We discuss using CMB…
Signal processing in wireless communications, such as precoding, detection, and channel estimation, are basically about solving inverse matrix problems, which, however, are slow and inefficient in conventional digital computers, thus…
We propose an alternative formalism to simulate CMB temperature maps in $\Lambda$CDM universes with nontrivial spatial topologies. This formalism avoids the need to explicitly compute the eigenmodes of the Laplacian operator in the spatial…
The free-electron population during the reionized epoch rescatters CMB temperature quadrupole and generates a now well-known polarization signal at large angular scales. While this contribution has been detected in the…
The SLUSCHI (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) automated package, with interface to the first-principles code VASP (Vienna Ab initio Simulation Package), was developed by us for efficiently determining…
Recently, the polarimetricmethod for thermally-induced polarization changes driven power losses (TIPCL) mitigation in complex laser systems has been developed. However, the final optimization relied on the four-parameter numerical process.…
Magnetic particle imaging (MPI) data is commonly reconstructed using a system matrix acquired in a time-consuming calibration measurement. The calibration approach has the important advantage over model-based reconstruction that it takes…
The implicit 2D3V particle-in-cell (PIC) code developed to study the interaction of ultrashort pulse lasers with matter [G. M. Petrov and J. Davis, Computer Phys. Comm. 179, 868 (2008); Phys. Plasmas 18, 073102 (2011)] has been parallelized…
In the realm of Large Language Model (LLM) inference, the inherent structure of transformer models coupled with the multi-GPU tensor parallelism strategy leads to a sequential execution of computation and communication. This results in…
We present a method for fast optimal estimation of the temperature angular power spectrum from observations of the cosmic microwave background. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from the posterior…
The new generation of CMB polarization experiments will reach limits of sensitivity never achieved before to detect the primordial B-mode signal. However, all these efforts will be futile if we lack a tight control of systematics. Here, we…
We propose a parallel version of the cross interpolation algorithm and apply it to calculate high-dimensional integrals motivated by Ising model in quantum physics. In contrast to mainstream approaches, such as Monte Carlo and quasi Monte…
In this work I present a technique of construction and fast evaluation of a family of cubic polynomials for analytic smoothing and graphical rendering of particles trajectories for flows in a generic geometry. The principal result of the…
The Preconditioned Conjugate Gradient (PCG) method is widely used for solving linear systems of equations with sparse matrices. A recent version of PCG, Pipelined PCG, eliminates the dependencies in the computations of the PCG algorithm so…