Source author record

Jochen Schütz

Jochen Schütz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

10 published item(s)

preprint2026arXiv

TrackHHL: The 1-Bit Quantum Filter for particle trajectory reconstruction

The transition to the High-Luminosity Large Hadron Collider (HL-LHC) presents a computational challenge where particle reconstruction complexity may outpace classical computing resources. While quantum computing offers potential speedups, standard algorithms like Harrow-Hassidim-Lloyd (HHL) require prohibitive circuit depths for near-term hardware. Here, we introduce the 1-Bit Quantum Filter, a domain-specific adaptation of HHL that reformulates tracking from matrix inversion to binary ground-state filtering. By replacing high-precision phase estimation with a single-ancilla spectral threshold and exploiting the Hamiltonian's sparsity, we achieve an asymptotic gate complexity of $O(\sqrt{N} \log N)$, given Hamiltonian dimension $N$. We validate this approach by simulating LHCb Vertex Locator events with a toy model, and benchmark performance using the noise models of Quantinuum H2 trapped-ion and IBM Heron superconducting processors. This work establishes a resource-efficient track reconstruction method capable of solving realistic event topologies on noise-free simulators and smaller tracking scenarios within the current constraints of the Noisy Intermediate Scale Quantum (NISQ) era.

preprint2021arXiv

Parallel-in-time high-order multiderivative IMEX solvers

In this work, we present a novel class of parallelizable high-order time integration schemes for the approximate solution of additive ODEs. The methods achieve high order through a combination of a suitable quadrature formula involving multiple derivatives of the ODE's right-hand side and a predictor-corrector ansatz. The latter approach is designed in such a way that parallelism in time is made possible. We present thorough analysis as well as numerical results that showcase scaling opportunities of methods from this class of solvers.

preprint2020arXiv

An asymptotic preserving semi-implicit multiderivative solver

In this work we construct a multiderivative implicit-explicit (IMEX) scheme for a class of stiff ordinary differential equations. Our solver is high-order accurate and has an asymptotic preserving (AP) property. The proposed method is based upon a two-derivative backward Taylor series base solver, which we show has an AP property. Higher order accuracies are found by iterating the result over a high-order multiderivative interpolant of the right hand side function, which we again prove has an AP property. Theoretical results showcasing the asymptotic consistency as well as the high-order accuracy of the solver are presented. In addition, an extension of the solver to an arbitrarily split right hand side function is also offered. Numerical results for a collection of standard test cases from the literature are presented that support the theoretical findings of the paper.

preprint2016arXiv

Implicit multistage two-derivative discontinuous Galerkin schemes for viscous conservation laws

In this paper we apply implicit two-derivative multistage time integrators to viscous conservation laws in one and two dimensions. The one dimensional solver discretizes space with the classical discontinuous Galerkin (DG) method, and the two dimensional solver uses a hybridized discontinuous Galerkin (HDG) spatial discretization for efficiency. We propose methods that permit us to construct implicit solvers using each of these spatial discretizations, wherein a chief difficulty is how to handle the higher derivatives in time. The end result is that the multiderivative time integrator allows us to obtain high-order accuracy in time while keeping the number of implicit stages at a minimum. We show numerical results validating and comparing methods.

preprint2014arXiv

A Comparison of Hybridized and Standard DG Methods for Target-Based hp-Adaptive Simulation of Compressible Flow

We present a comparison between hybridized and non-hybridized discontinuous Galerkin methods in the context of target-based hp-adaptation for compressible flow problems. The aim is to provide a critical assessment of the computational efficiency of hybridized DG methods. Hybridization of finite element discretizations has the main advantage, that the resulting set of algebraic equations has globally coupled degrees of freedom only on the skeleton of the computational mesh. Consequently, solving for these degrees of freedom involves the solution of a potentially much smaller system. This not only reduces storage requirements, but also allows for a faster solution with iterative solvers. Using a discrete-adjoint approach, sensitivities with respect to output functionals are computed to drive the adaptation. From the error distribution given by the adjoint-based error estimator, h- or p-refinement is chosen based on the smoothness of the solution which can be quantified by properly-chosen smoothness indicators. Numerical results are shown for subsonic, transonic, and supersonic flow around the NACA0012 airfoil. hp-adaptation proves to be superior to pure h-adaptation if discontinuous or singular flow features are involved. In all cases, a higher polynomial degree turns out to be beneficial. We show that for polynomial degree of approximation p=2 and higher, and for a broad range of test cases, HDG performs better than DG in terms of runtime and memory requirements.

preprint2014arXiv

A note on adjoint error estimation for one-dimensional stationary balance laws with shocks

We consider one-dimensional steady-state balance laws with discontinuous solutions. Giles and Pierce realized that a shock leads to a new term in the adjoint error representation for target functionals.This term disappears if and only if the adjoint solution satisfies an internal boundary condition. Curiously, most computer codes implementing adjoint error estimation ignore the new term in the functional, as well as the internal adjoint boundary condition. The purpose of this note is to justify this omission as follows: if one represents the exact forward and adjoint solutions as vanishing viscosity limits of the corresponding viscous problems, then the internal boundary condition is naturally satisfied in the limit.

preprint2014arXiv

A temporally adaptive hybridized discontinuous Galerkin method for time-dependent compressible flows

The potential of the hybridized discontinuous Galerkin (HDG) method has been recognized for the computation of stationary flows. Extending the method to time-dependent problems can, e.g., be done by backward difference formulae (BDF) or diagonally implicit Runge-Kutta (DIRK) methods. In this work, we investigate the use of embedded DIRK methods in an HDG solver, including the use of adaptive time-step control. Numerical results demonstrate the performance of the method for both linear and nonlinear (systems of) time-dependent convection-diffusion equations.

preprint2014arXiv

An asymptotic preserving method for linear systems of balance laws based on Galerkin's method

We apply the concept of Asymptotic Preserving (AP) schemes to the linearized p-system and discretize the resulting elliptic equation using standard continuous Finite Elements instead of Finite Differences. The fully discrete method is analyzed with respect to consistency, and we compare it numerically with more traditional methods such as Implicit Euler's method. Numerical results indicate that the AP method is indeed superior to more traditional methods.

preprint2014arXiv

Flux Splitting for stiff equations: A notion on stability

For low Mach number flows, there is a strong recent interest in the development and analysis of IMEX (implicit/explicit) schemes, which rely on a splitting of the convective flux into stiff and nonstiff parts. A key ingredient of the analysis is the so-called Asymptotic Preserving (AP) property, which guarantees uniform consistency and stability as the Mach number goes to zero. While many authors have focussed on asymptotic consistency, we study asymptotic stability in this paper: does an IMEX scheme allow for a CFL number which is independent of the Mach number? We derive a stability criterion for a general linear hyperbolic system. In the decisive eigenvalue analysis, the advective term, the upwind diffusion and a quadratic term stemming from the truncation in time all interact in a subtle way. As an application, we show that a new class of splittings based on characteristic decomposition, for which the commutator vanishes, avoids the deterioration of the time step which has sometimes been observed in the literature.

preprint2013arXiv

Adjoint-Based Error Estimation and Mesh Adaptation for Hybridized Discontinuous Galerkin Methods

We present a robust and efficient target-based mesh adaptation methodology, building on hybridized discontinuous Galerkin schemes for (nonlinear) convection-diffusion problems, including the compressible Euler and Navier-Stokes equations. Hybridization of finite element discretizations has the main advantage, that the resulting set of algebraic equations has globally coupled degrees of freedom only on the skeleton of the computational mesh. Consequently, solving for these degrees of freedom involves the solution of a potentially much smaller system. This not only reduces storage requirements, but also allows for a faster solution with iterative solvers. The mesh adaptation is driven by an error estimate obtained via a discrete adjoint approach. Furthermore, the computed target functional can be corrected with this error estimate to obtain an even more accurate value. The aim of this paper is twofold: Firstly, to show the superiority of adjoint-based mesh adaptation over uniform and residual-based mesh refinement, and secondly to investigate the efficiency of the global error estimate.