We present an algorithm for efficient evaluation of Boys functions $F_0,\dots,F_{k_\mathrm{max}}$ tailored to modern computing architectures, in particular graphical processing units (GPUs), where maximum throughput is high and data movement is costly. The method combines rational minimax approximations with upward and downward recurrence relations. The non-negative real axis is partitioned into three regions, $[0,\infty\rangle = A\cup B\cup C$, where regions $A$ and $B$ are treated using rational minimax approximations and region $C$ by an asymptotic approximation. This formulation avoids lookup tables and irregular memory access, making it well suited hardware with high maximum throughput and low latency. The rational minimax coefficients are generated using the rational Remez algorithm. For a target maximum absolute error of $\varepsilon_\mathrm{tol} = 5\cdot10^{-14}$, the corresponding approximation regions and coefficients for Boys functions $F_0,\dots,F_{32}$ are provided in the appendix.
A fast approximation to the Boys functions (related to the lower incomplete gamma function of half-integer parameter) by a single closed-form analytical expression for all argument values have been developed and tested. Besides the exponential function needed anyway for downward recursion, it uses a small number of addition, multiplication, division, and square root operations, and thus is straightforward to vectorize.
We present a new fast algorithm for computing the Boys function using a nonlinear approximation of the integrand via exponentials. The resulting algorithms evaluate the Boys function with real and complex valued arguments and are competitive with previously developed algorithms for the same purpose.
Based on the crystal structures of the previously proposed low energy phosphorene allotropes η-P and θ-P (Nano. Lett. 2015, 15, 3557), we propose five new structural stable phosphorene boys (XX-XY or XY-XY) and girls (XX-XX) through gene (XY from η-P and XX from θ-P) recombination methods. All of these five new phosphorene allotropes are obviously different from their parents, showing very different and fascinating two-dimensional patterns between each other. The dynamical stabilities of these phosphorene allotropes are confirmed positive and some of them are confirmed energetically more favorable than their parents (η-P and θ-P). Especially, the XX-XX type girl G1-P is confirmed energetically more favorable than all the previously proposed phosphorene allotropes, including black phosphorene (α-P, ACS Nano, 2014, 8, 4033) and blue phosphorene (β-P, Phys. Rev. Lett. 2014, 112, 176802), which is expected to be synthesized in future experiment through vapor deposition. Our results show that such a new promising phosphorene allotrope G1-P is an excellent candidate for potential applications in nano-electronics according to its middle band gap about 1.58 eV from DFT-HSE06 calculation.
This article describes Boy's surface in a nice way that does not make many demands on three-dimensional visualization. The article includes a kit that you can print out onto card stock and assemble with scissors and tape.
We prove Hodge-theoretic formulas for the $\mathbb{Q}$-factoriality defect of a normal projective variety, and for the local analytic $\mathbb{Q}$-factoriality defect of an analytic germ of a normal variety. These formulas lead to consequences ranging from a local analytic version of Samuel's conjecture to a characterization of projective rational homology threefolds with rational singularities, or the invariance of Hodge-Du Bois numbers under flops of projective threefolds.
We construct some version of the trace morphism between the Du Bois complexes, with applications towards the behavior of the local cohomological dimension and some Hodge theoretic aspects of singularities under finite morphisms.
This paper resolves a question of Huneke and Watanabe by proving a sharp upper bound for the multiplicity of Du Bois singularities: at a point of a $d$-dimensional variety with Du Bois singularities and embedding dimension $e$, the multiplicity is at most $\binom{e}{d}$. Additionally, the result recovers the previously known upper bound for the multiplicity of rational singularities.
We extend the notions of higher Du Bois and higher rational singularities to pairs in the sense of the minimal model program. We extend numerous results to these higher pairs, including Bertini type theorems, stability under finite maps and that m-rational pairs are m-Du Bois. We prove these using a generalized Kovács-Schwede-type injectivity theorem for pairs, the main technical result of this paper.
A partial answer is given to a question raised by Kovács and Taji in arxiv:2307.07192, namely that the relative Du Bois complex of a family parametrized by a non-singular curve commutes with base change to a general point on the base. It is also shown that this property usually fails for special points.
We study the relationship between higher Du Bois singularities and $K$-regularity, a notion that measures the $\mathbb{A}^1$-invariance of the algebraic $K$-groups. Building on this relationship, we establish a strengthened form of Vorst's conjecture for local complete intersections in characteristic zero. Our work also provides tools to construct new examples that illustrate various phenomena in the study of $K$-regularity. The main inputs for our results are vanishing theorems for the Du Bois complexes.
We study the behavior of Du~Bois singularities under base change and fiber products. For embeddable varieties in characteristic zero, we show that Du~Bois singularities descend from any field extension. We also prove that the product of a variety with rational singularities and a variety with Du~Bois singularities again has Du~Bois singularities, and we discuss the case of products of curves.
We prove an injectivity theorem for the cohomology of the Du Bois complexes of varieties with isolated singularities. We use this to deduce vanishing statements for the cohomologies of higher Du Bois complexes of such varieties. Besides some extensions and conjectures in the non-isolated case, we also provide analogues for intersection complexes.
We compute the Du Bois complexes of abstract cones over singular varieties, and use this to describe the local cohomological dimension and the non-positive K-groups of such cones.
We analyze the computational complexity of several popular video games released for the Nintendo Game Boy video game console. We analyze the complexity of generalized versions of four popular Game Boy games: Donkey Kong, Wario Land, Harvest Moon GB, and Mole Mania. We provide original proofs showing that these games are \textbf{NP}-hard. Our proofs rely on Karp reductions from four of Karp's original 21 \textbf{NP}-complete problems: \textsc{Sat}, \textsc{3-Cnf-Sat}, \textsc{Hamiltonian Cycle}, and \textsc{Knapsack}. We also discuss proofs easily derived from known results demonstrating the \textbf{NP}-hardness of Lock `n' Chase and The Lion King.
Diffusion generative models unlock new possibilities for inverse problems as they allow for the incorporation of strong empirical priors in scientific inference. Recently, diffusion models are repurposed for solving inverse problems using Gaussian approximations to conditional densities of the reverse process via Tweedie's formula to parameterise the mean, complemented with various heuristics. To address various challenges arising from these approximations, we leverage higher order information using Tweedie's formula and obtain a statistically principled approximation. We further provide a theoretical guarantee specifically for posterior sampling which can lead to a better theoretical understanding of diffusion-based conditional sampling. Finally, we illustrate the empirical effectiveness of our approach for general linear inverse problems on toy synthetic examples as well as image restoration. We show that our method (i) removes any time-dependent step-size hyperparameters required by earlier methods, (ii) brings stability and better sample quality across multiple noise levels, (iii) is the only method that works in a stable way with variance exploding (VE) forward processes as op
Let $X$ be a normal complex variety and $π:\tilde X \to X$ a resolution of singularities. We show that the inclusion morphism $π_*Ω_{\tilde X}^p\hookrightarrow Ω_X^{[p]}$ is an isomorphism for $p < \mathrm{codim}_X(X_{\mathrm{sing}})$ when $X$ has du Bois singularities, giving an improvement on Flenner's criterion for arbitrary singularities. We also study the $k$-du Bois definition from the perspective of holomorphic extension and compare how different restrictions on $\mathscr H^0(\underline Ω_X^p)$ affect the singularities of $X$, where $\underlineΩ_X^p$ is the $p^{th}$-graded piece of the du Bois complex.
We establish a characterization of the Du Bois complex of a reduced pair $(X,Z)$ when $X\smallsetminus Z$ has rational singularities. As an application, when $X$ has normal Du Bois singularities and $Z$ is the locus of non-rational singularities of $X$, holomorphic $p$-forms on the smooth locus of $X$ extend regularly to forms on a resolution of singularities for $p\le\mathrm{codim}_X Z-1$, and to forms with log poles over $Z$ for $p\ge\mathrm{codim}_X Z$. If $X$ is not necessarily Du Bois, then $p$-forms extend regularly for $p\le\mathrm{codim}_X Z-2$. This is a generalization of the theorems of Flenner, Greb-Kebekus-Kovács-Peternell, and Kebekus-Schnell on extending holomorphic (log) forms. A by-product of our methods is a new proof of the theorem of Kollár-Kovács that log canonical singularities are Du Bois. We also show that the Proj of the log canonical ring of a log canonical pair is Du Bois if this ring is finitely generated. The proofs are based on Saito's theory of mixed Hodge modules.
We introduce new notions of $k$-Du Bois and $k$-rational singularities, extending the previous definitions in the case of local complete intersections (lci), to include natural examples outside of this setting. We study the stability of these notions under general hyperplane sections and show that varieties with $k$-rational singularities are $k$-Du Bois, extending previous results in [MP22b] and [FL22b] in the lci and the isolated singularities cases. In the process, we identify the aspects of the theory that depend only on the vanishing of higher cohomologies of Du Bois complexes (or related objects), and not on the behaviour of the Kähler differentials.
Bayesian optimization (BO) has proven to be an effective paradigm for the global optimization of expensive-to-sample systems. One of the main advantages of BO is its use of Gaussian processes (GPs) to characterize model uncertainty which can be leveraged to guide the learning and search process. However, BO typically treats systems as black-boxes and this limits the ability to exploit structural knowledge (e.g., physics and sparse interconnections). Composite functions of the form $f(x, y(x))$, wherein GP modeling is shifted from the performance function $f$ to an intermediate function $y$, offer an avenue for exploiting structural knowledge. However, the use of composite functions in a BO framework is complicated by the need to generate a probability density for $f$ from the Gaussian density of $y$ calculated by the GP (e.g., when $f$ is nonlinear it is not possible to obtain a closed-form expression). Previous work has handled this issue using sampling techniques; these are easy to implement and flexible but are computationally intensive. In this work, we introduce a new paradigm which allows for the efficient use of composite functions in BO; this uses adaptive linearizations of