22 research outputs found
Numerical approximation of statistical solutions of scalar conservation laws
We propose efficient numerical algorithms for approximating statistical
solutions of scalar conservation laws. The proposed algorithms combine finite
volume spatio-temporal approximations with Monte Carlo and multi-level Monte
Carlo discretizations of the probability space. Both sets of methods are proved
to converge to the entropy statistical solution. We also prove that there is a
considerable gain in efficiency resulting from the multi-level Monte Carlo
method over the standard Monte Carlo method. Numerical experiments illustrating
the ability of both methods to accurately compute multi-point statistical
quantities of interest are also presented
On C*-algebras Related to the Roe algebra
Vi betrakter en gruppevirkning av G pĂĄ en mengde X. Opgaven omhandler C*-algebraer relatert til den uniforme Roe-algebraen. Vi bruker konstruksjoner fra kryssprodukter av C*-algebraer til ĂĄ konstruere en analog til den uniforme Roe-algebraen. Videre konstruerer vi en konkret representert C*-algebra som vil tilsvare den uniforme Roe-algebraen i tilfellet hvor G virker pĂĄ seg selv. Vi ser ogsĂĄ pĂĄ en vridd variant av kryssproduktet, og lager en generalisering av Roe-algebraen i denne retningen.
Vi betrakter spesielt hvordan hvordan egenskapene til den uniforme Roe-algebraen overføres til konstruksjonen med det vridde kryssproduktet, og ser på hvordan eksakthet av G manifesterer seg selv i den vridde konstruksjonen.
I tilfellet med den konkret representerte algebraen ser vi hvordan eksistens av Følnernet for virkningen av G på X overføres til denne algebraen.
Til slutt betrakter vi algebraen av nestenperiodiske funksjoner for ĂĄ gi en mulig fin underalgebra av Roe-algebraen. I det abelske tilfellet gir vi en karakterisering av en viktig klasse av underalgebraer av denne algebraen
Pseudo-Hamiltonian neural networks for learning partial differential equations
Pseudo-Hamiltonian neural networks (PHNN) were recently introduced for
learning dynamical systems that can be modelled by ordinary differential
equations. In this paper, we extend the method to partial differential
equations. The resulting model is comprised of up to three neural networks,
modelling terms representing conservation, dissipation and external forces, and
discrete convolution operators that can either be learned or be given as input.
We demonstrate numerically the superior performance of PHNN compared to a
baseline model that models the full dynamics by a single neural network.
Moreover, since the PHNN model consists of three parts with different physical
interpretations, these can be studied separately to gain insight into the
system, and the learned model is applicable also if external forces are removed
or changed.Comment: 33 pages, 14 figures; v2: minor changes to text, updated numerical
experiment
A Multi-level procedure for enhancing accuracy of machine learning algorithms
We propose a multi-level method to increase the accuracy of machine learning
algorithms for approximating observables in scientific computing, particularly
those that arise in systems modeled by differential equations. The algorithm
relies on judiciously combining a large number of computationally cheap
training data on coarse resolutions with a few expensive training samples on
fine grid resolutions. Theoretical arguments for lowering the generalization
error, based on reducing the variance of the underlying maps, are provided and
numerical evidence, indicating significant gains over underlying single-level
machine learning algorithms, are presented. Moreover, we also apply the
multi-level algorithm in the context of forward uncertainty quantification and
observe a considerable speed-up over competing algorithms
Iterative Surrogate Model Optimization (ISMO): An active learning algorithm for PDE constrained optimization with deep neural networks
We present a novel active learning algorithm, termed as iterative surrogate
model optimization (ISMO), for robust and efficient numerical approximation of
PDE constrained optimization problems. This algorithm is based on deep neural
networks and its key feature is the iterative selection of training data
through a feedback loop between deep neural networks and any underlying
standard optimization algorithm. Under suitable hypotheses, we show that the
resulting optimizers converge exponentially fast (and with exponentially
decaying variance), with respect to increasing number of training samples.
Numerical examples for optimal control, parameter identification and shape
optimization problems for PDEs are provided to validate the proposed theory and
to illustrate that ISMO significantly outperforms a standard deep neural
network based surrogate optimization algorithm
Constraint Preserving Mixers for the Quantum Approximate Optimization Algorithm
The quantum approximate optimization algorithm/quantum alternating operator ansatz (QAOA) is a heuristic to find approximate solutions of combinatorial optimization problems. Most of the literature is limited to quadratic problems without constraints. However, many practically relevant optimization problems do have (hard) constraints that need to be fulfilled. In this article, we present a framework for constructing mixing operators that restrict the evolution to a subspace of the full Hilbert space given by these constraints. We generalize the “XY”-mixer designed to preserve the subspace of “one-hot” states to the general case of subspaces given by a number of computational basis states. We expose the underlying mathematical structure which reveals more of how mixers work and how one can minimize their cost in terms of the number of CX gates, particularly when Trotterization is taken into account. Our analysis also leads to valid Trotterizations for an “XY”-mixer with fewer CX gates than is known to date. In view of practical implementations, we also describe algorithms for efficient decomposition into basis gates. Several examples of more general cases are presented and analyzed.publishedVersio
Mortality and pulmonary complications in patients undergoing surgery with perioperative SARS-CoV-2 infection: an international cohort study
Background: The impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) on postoperative recovery needs to be understood to inform clinical decision making during and after the COVID-19 pandemic. This study reports 30-day mortality and pulmonary complication rates in patients with perioperative SARS-CoV-2 infection. Methods: This international, multicentre, cohort study at 235 hospitals in 24 countries included all patients undergoing surgery who had SARS-CoV-2 infection confirmed within 7 days before or 30 days after surgery. The primary outcome measure was 30-day postoperative mortality and was assessed in all enrolled patients. The main secondary outcome measure was pulmonary complications, defined as pneumonia, acute respiratory distress syndrome, or unexpected postoperative ventilation. Findings: This analysis includes 1128 patients who had surgery between Jan 1 and March 31, 2020, of whom 835 (74·0%) had emergency surgery and 280 (24·8%) had elective surgery. SARS-CoV-2 infection was confirmed preoperatively in 294 (26·1%) patients. 30-day mortality was 23·8% (268 of 1128). Pulmonary complications occurred in 577 (51·2%) of 1128 patients; 30-day mortality in these patients was 38·0% (219 of 577), accounting for 81·7% (219 of 268) of all deaths. In adjusted analyses, 30-day mortality was associated with male sex (odds ratio 1·75 [95% CI 1·28–2·40], p\textless0·0001), age 70 years or older versus younger than 70 years (2·30 [1·65–3·22], p\textless0·0001), American Society of Anesthesiologists grades 3–5 versus grades 1–2 (2·35 [1·57–3·53], p\textless0·0001), malignant versus benign or obstetric diagnosis (1·55 [1·01–2·39], p=0·046), emergency versus elective surgery (1·67 [1·06–2·63], p=0·026), and major versus minor surgery (1·52 [1·01–2·31], p=0·047). Interpretation: Postoperative pulmonary complications occur in half of patients with perioperative SARS-CoV-2 infection and are associated with high mortality. Thresholds for surgery during the COVID-19 pandemic should be higher than during normal practice, particularly in men aged 70 years and older. Consideration should be given for postponing non-urgent procedures and promoting non-operative treatment to delay or avoid the need for surgery. Funding: National Institute for Health Research (NIHR), Association of Coloproctology of Great Britain and Ireland, Bowel and Cancer Research, Bowel Disease Research Foundation, Association of Upper Gastrointestinal Surgeons, British Association of Surgical Oncology, British Gynaecological Cancer Society, European Society of Coloproctology, NIHR Academy, Sarcoma UK, Vascular Society for Great Britain and Ireland, and Yorkshire Cancer Research
Computation of statistical solutions of hyperbolic systems of conservation laws
Statistical solutions are time-parameterized probability measures on spaces of integrable functions, that have been proposed recently as a framework for global solutions and uncertainty quantification for multi-dimensional hyperbolic system of conservation laws. By combining high-resolution finite volume methods with a Monte Carlo sampling procedure, we present a numerical algorithm to approximate statistical solutions. Under verifiable assumptions on the finite volume method, we prove that the approximations, generated by the proposed algorithm, converge in an appropriate topology to a statistical solution. Numerical experiments illustrating the convergence theory and revealing interesting properties of statistical solutions, are also presented.
We furthermore show that the multi-level Monte Carlo algorithm converges in the weak topology, and provide testable conditions for when the multi-level Monte Carlo algorithm outperforms the Monte Carlo algorithm.
Finally we present the Alsvinn simulator, a fast multi general purpose graphical processing unit (GPGPU) finite volume solver for hyperbolic conservation laws in multiple space dimensions. Alsvinn has native support for uncertainty quantifications, and exhibits excellent scaling on top tier compute clusters