We introduce three representation formulas for the fractional $p$-Laplace operator in the whole range of parameters $0<s<1$ and $1<p<\infty$. Note that for $p\ne 2$ this a nonlinear operator. The first representation is based on a splitting procedure that combines a renormalized nonlinearity with the linear heat semigroup. The second adapts the nonlinearity to the Caffarelli-Silvestre linear extension technique. The third one is the corresponding nonlinear version of the Balakrishnan formula. We also discuss the correct choice of the constant of the fractional $p$-Laplace operator in order to have continuous dependence as $p\to 2$ and $s \to 0^+, 1^-$. A number of consequences and proposals are derived. Thus, we propose a natural spectral-type operator in domains, different from the standard restriction of the fractional $p$-Laplace operator acting on the whole space. We also propose numerical schemes, a new definition of the fractional $p$-Laplacian on manifolds, as well as alternative characterizations of the $W^{s,p}(\mathbb{R}^n)$ seminorms.
We prove that a formula predicted on the basis of non-rigorous physics arguments [Zdeborova and Krzakala: Phys. Rev. E (2007)] provides a lower bound on the chromatic number of sparse random graphs. The proof is based on the interpolation method from mathematical physics. In the case of random regular graphs the lower bound can be expressed algebraically, while in the case of the binomial random we obtain a variational formula. As an application we calculate improved explicit lower bounds on the chromatic number of random graphs for small (average) degrees. Additionally, show how asymptotic formulas for large degrees that were previously obtained by lengthy and complicated combinatorial arguments can be re-derived easily from these new results.
We introduce a physically relevant stochastic representation of the rotating shallow water equations. The derivation relies mainly on a stochastic transport principle and on a decomposition of the fluid flow into a large-scale component and a noise term that models the unresolved flow components. As for the classical (deterministic) system, this scheme, referred to as modelling under location uncertainty (LU), conserves the global energy of any realization and provides the possibility to generate an ensemble of physically relevant random simulations with a good trade-off between the model error representation and the ensemble's spread. To maintain numerically the energy conservation feature, we combine an energy (in space) preserving discretization of the underlying deterministic model with approximations of the stochastic terms that are based on standard finite volume/difference operators. The LU derivation, built from the very same conservation principles as the usual geophysical models, together with the numerical scheme proposed can be directly used in existing dynamical cores of global numerical weather prediction models. The capabilities of the proposed framework is demonstrated for an inviscid test case on the f-plane and for a barotropically unstable jet on the sphere.
In this paper, we discuss DNA codes that are cyclic or quasi-cyclic over $\Z_{4}+\omega \Z_{4}$, where $\omega^{2}=2+2\omega$ along with methods to construct these with combinatorial constraints. We also generalize results obtained for the ring $\Z_{4}+\omega \Z_{4}$, where $\omega^{2}=2+2\omega$, and some other rings to the sixteen rings $R_{\theta}=\Z_{4}+\omega \Z_{4}$, where $\omega^{2}=\theta\in \Z_{4}+\omega \Z_{4}$, using the generalized Gau map and Gau distance in \cite{3}. We determine a relationship between the Gau distance and Hamming distance for linear codes over the sixteen rings $R_{\theta}$ which enables us to attain an upper boundary for the Gau distance of free codes that are self-dual over the rings $R_{\theta}$.
In this work, we study the propagation of solitons in lossy optical fibers. The main objective of this work is to study the loss of energy of the soliton wave during propagation and then to evaluate the impact of this loss on the transmission of the soliton signal. In this context, a numerical scheme was developed to solve a system of complex partial differential equations (CPDE) that describes the propagation of solitons in optical fibers with loss and nonlinear amplification mechanisms. The numerical procedure is based on the mathematical theory of Taylor series of complex functions. We adapted the Finite Difference Method (FDM) to approximate derivatives of complex functions. Then, we solve the algebraic system resulting from the discretization, implicitly, through the relaxation Gauss-Seidel method (RGSM). The numerical study of CPDE system with linear and cubic attenuation showed that soliton waves undergo attenuation, dispersion, and oscillation effects. On the other hand, we find that by considering the nonlinear term (cubic term) as an optical amplification, it is possible to partially compensate for the attenuation of the optical signal. Finally, we show that a gain of 9% triples the propagation distance of the fundamental soliton wave, when the dissipation rate is 1%.
Given independent identically-distributed samples from a one-dimensional distribution, IAs are formed by partitioning samples into pairs, triplets, or nth-order groupings and retaining the median of those groupings that are approximately equal. A new statistical method, Independent Approximates (IAs), is defined and proven to enable closed-form estimation of the parameters of heavy-tailed distributions. The pdf of the IAs is proven to be the normalized nth-power of the original density. From this property, heavy-tailed distributions are proven to have well-defined means for their IA pairs, finite second moments for their IA triplets, and a finite, well-defined (n-1)th-moment for the nth-grouping. Estimation of the location, scale, and shape (inverse of degree of freedom) of the generalized Pareto and Student's t distributions are possible via a system of three equations. Performance analysis of the IA estimation methodology is conducted for the Student's t distribution using between 1000 to 100,000 samples. Closed-form estimates of the location and scale are determined from the mean of the IA pairs and the variance of the IA triplets, respectively. For the Student's t distribution, the geometric mean of the original samples provides a third equation to determine the shape, though its nonlinear solution requires an iterative solver. With 10,000 samples the relative bias of the parameter estimates is less than 0.01 and the relative precision is less than +/-0.1. The theoretical precision is finite for a limited range of the shape but can be extended by using higher-order groupings for a given moment.
A new numerical approximation method for a class of Gaussian random fields on compact Riemannian manifolds is introduced. This class of random fields is characterized by the Laplace--Beltrami operator on the manifold. A Galerkin approximation is combined with a polynomial approximation using Chebyshev series. This so-called Galerkin--Chebyshev approximation scheme yields efficient and generic sampling algorithms for Gaussian random fields on manifolds. Strong and weak orders of convergence for the Galerkin approximation and strong convergence orders for the Galerkin--Chebyshev approximation are shown and confirmed through numerical experiments.
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.
We present symbolic and numerical methods for computing Poisson brackets on the spaces of measures with positive densities of the plane, the 2-torus, and the 2-sphere. We apply our methods to compute symplectic areas of finite regions for the case of the 2-sphere, including an explicit example for Gaussian measures with positive densities.
We discuss approaches to computing eigenfunctions of the Ornstein--Uhlenbeck (OU) operator in more than two dimensions. While the spectrum of the OU operator and theoretical properties of its eigenfunctions have been well characterized in previous research, the practical computation of general eigenfunctions has not been resolved. We review special cases for which the eigenfunctions can be expressed exactly in terms of commonly used orthogonal polynomials. Then we present a tractable approach for computing the eigenfunctions in general cases and comment on its dimension dependence.
We show that it is provable in PA that there is an arithmetically definable sequence $\{\phi_{n}:n \in \omega\}$ of $\Pi^{0}_{2}$-sentences, such that - PRA+$\{\phi_{n}:n \in \omega\}$ is $\Pi^{0}_{2}$-sound and $\Pi^{0}_{1}$-complete - the length of $\phi_{n}$ is bounded above by a polynomial function of $n$ with positive leading coefficient - PRA+$\phi_{n+1}$ always proves 1-consistency of PRA+$\phi_{n}$. One has that the growth in logical strength is in some sense "as fast as possible", manifested in the fact that the total general recursive functions whose totality is asserted by the true $\Pi^{0}_{2}$-sentences in the sequence are cofinal growth-rate-wise in the set of all total general recursive functions. We then develop an argument which makes use of a sequence of sentences constructed by an application of the diagonal lemma, which are generalisations in a broad sense of Hugh Woodin's "Tower of Hanoi" construction as outlined in his essay "Tower of Hanoi" in Chapter 18 of the anthology "Truth in Mathematics". The argument establishes the result that it is provable in PA that $P \neq NP$. We indicate how to pull the argument all the way down into EFA.