We consider the solution of the large-scale nonsymmetric algebraic Riccati equation XCX − XD − AX + B = 0, with M ≡ [D,−C;−B,A] ∈ R(n1+n2)×(n1+n2) being a nonsingular M-matrix. In addition, A and D are sparselike, with the products A−1u, A−⊤u, D−1v, and D−⊤v computable in O(n) complexity (with n = max{n1, n2}), for some vectors u and v, and B, C are low ranked. The structure-preserving doubling algorithms (SDA) by Guo, Lin, and Xu [Numer. Math., 103 (2006), pp. 392–412] is adapted, with the appropriate applications of the Sherman–Morrison– Woodbury formula and the sparse-plus-low-rank representations of various iterates. The resulting large-scale doubling algorithm has an O(n) computational complexity and memory requirement per iteration and converges essentially quadratically. A detailed error analysis, on the effects of truncation of iterates with an explicit forward error bound for the approximate solution from the SDA, and some numerical results will be presented.

Let $D={\mathbb H} \setminus \cup_{k=1}^N C_k$ be a standard slit domain where $\mathbb H$ is the upper half plane and $C_k$, $1\leq k\leq N$, are mutually disjoint horizontal line segments in $\HH$. Given a Jordan arc $\gamma\subset D$ starting at $\partial {\mathbb H},$ let $g_t$ be the unique conformal map from $D\setminus \gamma[0,t]$ onto a standard slit domain $D_t$ satisfying the hydrodynamic normalization. We prove that $g_t$ satisfies an ODE with the kernel on its righthand side being the complex Poisson kernel of the Brownian motion with darning (BMD) for $D_t$, generalizing the chordal Loewner equation for the simply connected domain $D={\mathbb H}.$ Such a generalization has been obtained by Y. Komatu in the case of circularly slit annuli and by R. O. Bauer and R. M. Friedrich in the present chordal case, but only in the sense of the left derivative in $t$. We establish the differentiability of $g_t$ in $t$ to make the equation a genuine ODE. To this end, we first derive the continuity of $g_t(z)$ in $t$ with a certain uniformity in $z$ from a probabilistic expression of $\Im g_t(z)$ in terms of the BMD for $D$, which is then combined with a Lipschitz continuity of the complex Poisson kernel under the perturbation of standard slit domains to get the desired differentiability.

Skew-orthogonal polynomials (SOPs) arise in the study of the <i>n</i>-point distribution function for orthogonal and symplectic random matrix ensembles. Motivated by the average of characteristic polynomials of the Bures random matrix ensemble studied in Forrester and Kieburg (Commun Math Phys 342(1):151187, 2016), we propose the concept of <i>partial-skew-orthogonal polynomials</i> (PSOPs) as a modification of the SOPs, and then the PSOPs with a variety of special skew-symmetric kernels and weight functions are addressed. By considering appropriate deformations of the weight functions, we derive nine integrable lattices in different dimensions. As a consequence, the tau-functions for these systems are shown to be expressed in terms of Pfaffians and the wave vectors PSOPs. In fact, the tau-functions also admit the multiple integral representations. Among these integrable lattices, some of them are

We study the existence of special Lagrangian submanifolds of log Calabi–Yau manifolds equipped with the complete Ricci-flat Kähler metric constructed by Tian and Yau. We prove that if X is a Tian–Yau manifold and if the compact Calabi–Yau manifold at infinity admits a single special Lagrangian, then X admits infinitely many disjoint special Lagrangians. In complex dimension 2, we prove that if Y is a del Pezzo surface or a rational elliptic surface and D∈|-KY| is a smooth divisor with D^2=d , then X=Y\D admits a special Lagrangian torus fibration, as conjectured by Strominger–Yau–Zaslow and Auroux. In fact, we show that X admits twin special Lagrangian fibrations, confirming a prediction of Leung and Yau. In the special case that Y is a rational elliptic surface or Y=P^2, we identify the singular fibers for generic data, thereby confirming two conjectures of Auroux. Finally, we prove that after a hyper-Kähler rotation, X can be compactified to the complement of a Kodaira type I_d fiber appearing as a singular fiber in a rational elliptic surface π: Y →P^1.

Blocking is commonly used in randomized experiments to increase efficiency of estimation. A generalization of blocking removes allocations with imbalance in covariate distributions between treated and control units, and then randomizes within the remaining set of allocations with balance. This idea of rerandomization was formalized by Morgan and Rubin (Annals of Statistics, 2012, 40, 1263–1282), who suggested using Mahalanobis distance between treated and control covariate means as the criterion for removing unbalanced allocations. Kallus (Journal of the Royal Statistical Society, Series B: Statistical Methodology, 2018, 80, 85–112) proposed reducing the set of balanced allocations to the minimum. Here we discuss the implication of such an ‘optimal’ rerandomization design for inferences to the units in the sample and to the population from which the units in the sample were randomly drawn. We argue that, in general, it is a bad idea to seek the optimal design for an inference because that inference typically only reflects uncertainty from the random sampling of units, which is usually hypothetical, and not the randomization of units to treatment versus control.