In this paper, we address the problem of constructing $C^2$ cubic spline functions on a given arbitrary triangulation $\mathcal{T}$. To this end, we endow every triangle of $\mathcal{T}$ with a Wang-Shi macro-structure. The $C^2$ cubic space on such a refined triangulation has a stable dimension and optimal approximation power. Moreover, any spline function in such space can be locally built on each of the macro-triangles independently via Hermite interpolation. We provide a simplex spline basis for the space of $C^2$ cubics defined on a single macro-triangle which behaves like a Bernstein/B-spline basis over the triangle. The basis functions inherit recurrence relations and differentiation formulas from the simplex spline construction, they form a nonnegative partition of unity, they admit simple conditions for $C^2$ joins across the edges of neighboring triangles, and they enjoy a Marsden-like identity. Also, there is a single control net to facilitate control and early visualization of a spline function over the macro-triangle. Thanks to these properties, the complex geometry of the Wang-Shi macro-structure is transparent to the user. Stable global bases for the full space of $C^2$ cubics on the Wang-Shi refined triangulation $\mathcal{T}$ are deduced from the local simplex spline basis by extending the concept of minimal determining sets.
This paper defines analysis-suitable T-splines for arbitrary degree (including even and mixed degrees) and arbitrary dimension. We generalize the concept of anchor elements known from the two-dimensional setting, extend two existing concepts of analysis-suitability and justify their sufficiency for linear independence of the T-spline basis. Finally, we propose a local refinement scheme for multivariate T-splines that allows anisotropic refinement and preserves weak geometric analysis-suitability.
For a fixed type of Petri nets $\tau$, \textsc{$\tau$-Synthesis} is the task of finding for a given transition system $A$ a Petri net $N$ of type $\tau$ ($\tau$-net, for short) whose reachability graph is isomorphic to $A$ if there is one. The decision version of this search problem is called \textsc{$\tau$-Solvability}. If an input $A$ allows a positive decision, then it is called $\tau$-solvable and a sought net $N$ $\tau$-solves $A$. As a well known fact, $A$ is $\tau$-solvable if and only if it has the so-called $\tau$-\emph{event state separation property} ($\tau$-ESSP, for short) and the $\tau$-\emph{state separation property} ($\tau$-SSP, for short). The question whether $A$ has the $\tau$-ESSP or the $\tau$-SSP defines also decision problems. In this paper, for all $b\in \mathbb{N}$, we completely characterize the computational complexity of \textsc{$\tau$-Solvability}, \textsc{$\tau$-ESSP} and \textsc{$\tau$-SSP} for the types of pure $b$-bounded Place/Transition-nets, the $b$-bounded Place/Transition-nets and their corresponding $\mathbb{Z}_{b+1}$-extensions.
How efficiently can we find an unknown graph using distance queries between its vertices? We assume that the unknown graph is connected, unweighted, and has bounded degree. The goal is to find every edge in the graph. This problem admits a reconstruction algorithm based on multi-phase Voronoi-cell decomposition and using $\tilde O(n^{3/2})$ distance queries. In our work, we analyze a simple reconstruction algorithm. We show that, on random $\Delta$-regular graphs, our algorithm uses $\tilde O(n)$ distance queries. As by-products, we can reconstruct those graphs using $O(\log^2 n)$ queries to an all-distances oracle or $\tilde O(n)$ queries to a betweenness oracle, and we bound the metric dimension of those graphs by $\log^2 n$. Our reconstruction algorithm has a very simple structure, and is highly parallelizable. On general graphs of bounded degree, our reconstruction algorithm has subquadratic query complexity.
We show that reconstructing a curve in $\mathbb{R}^d$ for $d\geq 2$ from a $0.66$-sample is always possible using an algorithm similar to the classical NN-Crust algorithm. Previously, this was only known to be possible for $0.47$-samples in $\mathbb{R}^2$ and $\frac{1}{3}$-samples in $\mathbb{R}^d$ for $d\geq 3$. In addition, we show that there is not always a unique way to reconstruct a curve from a $0.72$-sample; this was previously only known for $1$-samples. We also extend this non-uniqueness result to hypersurfaces in all higher dimensions.
We establish improved uniform error bounds on time-splitting methods for the long-time dynamics of the Dirac equation with small electromagnetic potentials characterized by a dimensionless parameter $\varepsilon\in (0, 1]$ representing the amplitude of the potentials. We begin with a semi-discritization of the Dirac equation in time by a time-splitting method, and then followed by a full-discretization in space by the Fourier pseudospectral method. Employing the unitary flow property of the second-order time-splitting method for the Dirac equation, we prove uniform error bounds at $C(t)\tau^2$ and $C(t)(h^m+\tau^2)$ for the semi-discretization and full-discretization, respectively, for any time $t\in[0,T_\varepsilon]$ with $T_\varepsilon = T/\varepsilon$ for $T > 0$, which are uniformly for $\varepsilon \in (0, 1]$, where $\tau$ is the time step, $h$ is the mesh size, $m\geq 2$ depends on the regularity of the solution, and $C(t) = C_0 + C_1\varepsilon t\le C_0+C_1T$ grows at most linearly with respect to $t$ with $C_0\ge0$ and $C_1>0$ two constants independent of $t$, $h$, $\tau$ and $\varepsilon$. Then by adopting the regularity compensation oscillation (RCO) technique which controls the high frequency modes by the regularity of the solution and low frequency modes by phase cancellation and energy method, we establish improved uniform error bounds at $O(\varepsilon\tau^2)$ and $O(h^m +\varepsilon\tau^2)$ for the semi-discretization and full-discretization, respectively, up to the long-time $T_\varepsilon$. Numerical results are reported to confirm our error bounds and to demonstrate that they are sharp. Comparisons on the accuracy of different time discretizations for the Dirac equation are also provided.
This work is devoted to design and study efficient and accurate numerical schemes to approximate a chemo-attraction model with consumption effects, which is a nonlinear parabolic system for two variables; the cell density and the concentration of the chemical signal that the cell feel attracted to. We present several finite element schemes to approximate the system, detailing the main properties of each of them, such as conservation of cells, energy-stability and approximated positivity. Moreover, we carry out several numerical simulations to study the efficiency of each of the schemes and to compare them with others classical schemes.
The mean width of a convex body is the average distance between parallel supporting hyperplanes when the normal direction is chosen uniformly over the sphere. The Simplex Mean Width Conjecture (SMWC) is a longstanding open problem that says the regular simplex has maximum mean width of all simplices contained in the unit ball and is unique up to isometry. We give a self contained proof of the SMWC in $d$ dimensions. The main idea is that when discussing mean width, $d+1$ vertices $v_i\in\mathbb{S}^{d-1}$ naturally divide $\mathbb{S}^{d-1}$ into $d+1$ Voronoi cells and conversely any partition of $\mathbb{S}^{d-1}$ points to selecting the centroids of regions as vertices. We will show that these two conditions are enough to ensure that a simplex with maximum mean width is a regular simplex.
The present paper focuses on the construction of a set of submatrices of a circulant matrix such that it is a smaller set to verify that the circulant matrix is an MDS (maximum distance separable) one, comparing to the complete set of square submatrices needed in general case. The general MDS verification method requires to test for singular submatrices: if at least one square submatrix is singular the matrix is not MDS. However, the complexity of the general method dramatically increases for matrices of a greater dimension. We develop an algorithm that constructs a smaller subset of submatrices thanks to a simple structure of circulant matrices. The algorithm proposed in the paper reduces the size of the testing set by approximately two matrix orders.
Linear kinetic transport equations play a critical role in optical tomography, radiative transfer and neutron transport. The fundamental difficulty hampering their efficient and accurate numerical resolution lies in the high dimensionality of the physical and velocity/angular variables and the fact that the problem is multiscale in nature. Leveraging the existence of a hidden low-rank structure hinted by the diffusive limit, in this work, we design and test the angular-space reduced order model for the linear radiative transfer equation, the first such effort based on the celebrated reduced basis method (RBM). Our method is built upon a high-fidelity solver employing the discrete ordinates method in the angular space, an asymptotic preserving upwind discontinuous Galerkin method for the physical space, and an efficient synthetic accelerated source iteration for the resulting linear system. Addressing the challenge of the parameter values (or angular directions) being coupled through an integration operator, the first novel ingredient of our method is an iterative procedure where the macroscopic density is constructed from the RBM snapshots, treated explicitly and allowing a transport sweep, and then updated afterwards. A greedy algorithm can then proceed to adaptively select the representative samples in the angular space and form a surrogate solution space. The second novelty is a least-squares density reconstruction strategy, at each of the relevant physical locations, enabling the robust and accurate integration over an arbitrarily unstructured set of angular samples toward the macroscopic density. Numerical experiments indicate that our method is effective for computational cost reduction in a variety of regimes.
De-homogenization is becoming an effective method to significantly expedite the design of high-resolution multiscale structures, but existing methods have thus far been confined to simple static compliance minimization. There are two critical challenges to be addressed in accommodating general cases: enabling the design of unit-cell orientation and using free-form microstructures. In this paper, we propose a data-driven de-homogenization method that allows effective design of the unit-cell orientation angles and conformal mapping of spatially varying, complex microstructures. We devise a parameterized microstructure composed of rods in different directions to provide more diversity in stiffness while retaining geometrical simplicity. The microstructural geometry-property relationship is then surrogated by a neural network to avoid costly homogenization. A Cartesian representation of the unit-cell orientation is incorporated into homogenization-based optimization to design the angles. Corresponding high-resolution multiscale structures are obtained from the homogenization-based designs through a conformal mapping constructed with sawtooth function fields. This allows us to assemble complex microstructures with an oriented and compatible tiling pattern, while preserving the local homogenized properties. To demonstrate our method with a specific application, we optimize the frequency response of structures under harmonic excitations within a given frequency range. It is the first time that a sawtooth function is applied in a de-homogenization framework for complex design scenarios beyond static compliance minimization. The examples illustrate that multiscale structures can be generated with high efficiency and much better dynamic performance compared with the macroscale-only optimization. Beyond frequency response design, our proposed framework can be applied to other general problems.