Source author record

Thomas Schulthess

Thomas Schulthess appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el Distributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Productive Performance Engineering for Weather and Climate Modeling with Python

Earth system models are developed with a tight coupling to target hardware, often containing specialized code predicated on processor characteristics. This coupling stems from using imperative languages that hard-code computation schedules and layout. We present a detailed account of optimizing the Finite Volume Cubed-Sphere Dynamical Core (FV3), improving productivity and performance. By using a declarative Python-embedded stencil domain-specific language and data-centric optimization, we abstract hardware-specific details and define a semi-automated workflow for analyzing and optimizing weather and climate applications. The workflow utilizes both local and full-program optimization, as well as user-guided fine-tuning. To prune the infeasible global optimization space, we automatically utilize repeating code motifs via a novel transfer tuning approach. On the Piz Daint supercomputer, we scale to 2,400 GPUs, achieving speedups of up to 3.92x over the tuned production implementation at a fraction of the original code.

preprint2014arXiv

Two-particle correlations in a dynamic cluster approximation with continuous momentum dependence: Superconductivity in the 2D Hubbard model

The DCA$^+$ algortihm was recently introduced to extend the dynamic cluster approximation (DCA) with a continuous lattice self-energy in order to achieve better convergence with cluster size. Here we extend the DCA$^+$ algorithm to the calculation of two-particle correlation functions by introducing irreducible vertex functions with continuous momentum dependence consistent with the DCA$^+$ self-energy. This enables a significantly more controlled and reliable study of phase transitions than with the DCA. We test the new method by calculating the superconducting transition temperature $T_{c}$ in the attractive Hubbard model and show that it reproduces previous high-precision determinantal quantum Monte Carlo results. We then calculate $T_c$ in the doped repulsive Hubbard model, for which previous DCA calculations could only access the weak-coupling ($U=4t$) regime for large clusters. We show that the new algorithm provides access to much larger clusters and delivers asymptotically converged results for $T_c$ for both the weak ($U=4t$) and intermediate ($U=7t$) coupling regimes, and thereby enables the accurate determination of the exact infinite cluster size result.

preprint2013arXiv

The Continuous-Pole-Expansion method to obtain spectra of electronic lattice models

We present a new algorithm to analytically continue the self-energy of quantum many-body systems from Matsubara frequencies to the real axis. The method allows straightforward, unambiguous computation of electronic spectra for lattice models of strongly correlated systems from self-energy data that has been collected with state-of-the are continuous time solvers within dynamical mean field simulations. Using well-known analytical properties of the self-energy, the analytic continuation is cast into a constrained minimization problem that can be formulated as a quadratic programmable optimization with linear constraints. The algorithm is validated against exactly solvable finite size problems, showing that all features of the spectral function near the Femi level are very well reproduced and coarse features are reproduced for all energies. The method is applied to two well known lattice problems, the two-dimensional Hubbard model at half filling where the momentum dependence of the gap formation is studied, as well as a multi-band model of NiO, for which the spectral function can be directly compared to experiment. Agreement with results published results is very good.

preprint2011arXiv

Sub-matrix updates for the Continuous-Time Auxiliary Field algorithm

We present a sub-matrix update algorithm for the continuous-time auxiliary field method that allows the simulation of large lattice and impurity problems. The algorithm takes optimal advantage of modern CPU architectures by consistently using matrix instead of vector operations, resulting in a speedup of a factor of $\approx 8$ and thereby allowing access to larger systems and lower temperature. We illustrate the power of our algorithm at the example of a cluster dynamical mean field simulation of the Néel transition in the three-dimensional Hubbard model, where we show momentum dependent self-energies for clusters with up to 100 sites.