Search results for "Mathematical software"
showing 10 items of 60 documents
Array programming with NumPy.
2020
Array programming provides a powerful, compact and expressive syntax for accessing, manipulating and operating on data in vectors, matrices and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It has an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, materials science, engineering, finance and economics. For example, in astronomy, NumPy was an important part of the software stack used in the discovery of gravitational waves1 and in the first imaging of a black hole2. Here we review how a few fundamental array concepts lead to a simple and powerful programmi…
RationalizeRoots: Software Package for the Rationalization of Square Roots
2019
The computation of Feynman integrals often involves square roots. One way to obtain a solution in terms of multiple polylogarithms is to rationalize these square roots by a suitable variable change. We present a program that can be used to find such transformations. After an introduction to the theoretical background, we explain in detail how to use the program in practice.
Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements
2014
We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time.
An HLLC Riemann solver for resistive relativistic magnetohydrodynamics
2017
We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through one- and two-dimensional test problems.
Factorization of denominators in integration-by-parts reductions
2020
We present a Mathematica package which finds a basis of master integrals for the Feynman integral reduction. In this basis the dependence on the dimensional regularization in the denominators factorizes in kinematic independent polynomials.
The dyon charge in noncommutative gauge theories
2007
We present an explicit classical dyon solution for the noncommutative version of the Yang-Mills-Higgs model (in the Prasad-Sommerfield limit) with a tehta term. We show that the relation between classical electric and magnetic charges also holds in noncommutative space. Extending the Noether approach to the case of a noncommutative gauge theory, we analyze the effect of CP violation at the quantum level, induced both by the theta term and by noncommutativity and we prove that the Witten effect formula for the dyon charge remains the same as in ordinary space.
Subleading Regge limit from a soft anomalous dimension
2018
Wilson lines capture important features of scattering amplitudes, for example soft effects relevant for infrared divergences, and the Regge limit. Beyond the leading power approximation, corrections to the eikonal picture have to be taken into account. In this paper, we study such corrections in a model of massive scattering amplitudes in N = 4 super Yang-Mills, in the planar limit, where the mass is generated through a Higgs mechanism. Using known three-loop analytic expressions for the scattering amplitude, we find that the first power suppressed term has a very simple form, equal to a single power law. We propose that its exponent is governed by the anomalous dimension of a Wilson loop w…
A new method for computing one-loop integrals
1994
We present a new program package for calculating one-loop Feynman integrals, based on a new method avoiding Feynman parametrization and the contraction due to Passarino and Veltman. The package is calculating one-, two- and three-point functions both algebraically and numerically to all tensor cases. This program is written as a package for Maple. An additional Mathematica version is planned later.
oneloop 2.0 — A program package calculating one-loop integrals
1997
We present an improved version of our program package oneloop which -- written as a package for MAPLE -- solves one-loop Feynman integrals. The package is calculating one-, two- and three-point functions both algebraically and numerically to any tensor rank. In addition to the original version oneloop 2.0 also calculates infrared divergent integrals. Higher powers of propagator terms and the $O(\eps)$ parts relevant for two-loop calculations are now supported.
A Projected Algebraic Multigrid Method for Linear Complementarity Problems
2011
We present an algebraic version of an iterative multigrid method for obstacle problems, called projected algebraic multigrid (PAMG) here. We show that classical AMG algorithms can easily be extended to deal with this kind of problem. This paves the way for efficient multigrid solution of obstacle problems with partial differential equations arising, for example, in financial engineering.