On probabilistic results for the discrepancy of a hybrid-Monte Carlo sequence

In many applications it has been observed that hybrid-Monte Carlo sequences perform better than Monte Carlo and quasi-Monte Carlo sequences, especially in difficult problems. For a mixed s$s$-dimensional sequence m$m$, whose elements are vectors obtained by concatenating d$d$-dimensional vectors from a low-discrepancy sequence q$q$ with (s−d)$(s-d)$-dimensional random vectors, probabilistic upper bounds for its star discrepancy have been provided. In a paper of G. Ökten, B. Tuffin and V. Burago [G. Ökten, B. Tuffin, V. Burago, J. Complexity 22 (2006), 435–458] it was shown that for arbitrary ε>0$\epsilon >0$ the difference of the star discrepancies of the first N$N$ points of m$m$ and q$q$ is bounded by ε$\epsilon$ with probability at least 1−2exp(−ε²N/2)$1-2exp(-{\epsilon }^{2}N/2)$ for N$N$ sufficiently large. The authors did not study how large N$N$ actually has to be and if and how this actually depends on the parameters s$s$ and ε$\epsilon$. In this note we derive a lower bound for N$N$, which significantly depends on s$s$ and ε$\epsilon$. Furthermore, we provide a probabilistic bound for the difference of the star discrepancies of the first N$N$ points of m$m$ and q$q$, which holds without any restrictions on N$N$. In this sense it improves on the bound of Ökten, Tuffin and Burago and is more helpful in practice, especially for small sample sizes N$N$. We compare this bound to other known bounds.     

Critical exponents for non-simultaneous blow-up in a localized parabolic system

The paper deals with the radially symmetric solutions of _ut=Δu+u^m(x,t)vⁿ(0,t)$u_t=\Delta u+u^m(x,t)v^n(0,t)$, _vt=Δv+u^p(0,t)v^q(x,t)$v_t=\Delta v+u^p(0,t)v^q(x,t)$, subject to null Dirichlet boundary conditions. For the blow-up classical solutions, we propose the critical exponents for non-simultaneous blow-up by determining the complete and optimal classification for all the non-negative exponents: (i) There exist initial data such that u$u$ (v$v$) blows up alone if and only if m>p+1$m>p+1$ (q>n+1$q>n+1$), which means that any blow-up is simultaneous if and only if m≤p+1$m\le p+1$, q≤n+1$q\le n+1$. (ii) Any blow-up is u$u$ (v$v$) blowing up with v$v$ (u$u$) remaining bounded if and only if m>p+1$m>p+1$, q≤n+1$q\le n+1$ (m≤p+1$m\le p+1$, q>n+1$q>n+1$). (iii) Both non-simultaneous and simultaneous blow-up may occur if and only if m>p+1$m>p+1$, q>n+1$q>n+1$. Moreover, we consider the blow-up rate and set estimates which were not obtained in the previously known work for the same model.     

Quadrature rule for Abel’s equations: Uniformly approximating fractional derivatives

An automatic quadrature method is presented for approximating fractional derivative D^qf(x)$D^{q}f\left(x\right)$ of a given function f(x)$f\left(x\right)$, which is defined by an indefinite integral involving f(x)$f\left(x\right)$. The present method interpolates f(x)$f\left(x\right)$ in terms of the Chebyshev polynomials in the range [0, 1] to approximate the fractional derivative D^qf(x)$D^{q}f\left(x\right)$ uniformly for 0≤x≤1$0\le x\le 1$, namely the error is bounded independently of x$x$. Some numerical examples demonstrate the performance of the present automatic method.     

An approximation method for strictly pseudocontractive mappings

Let K$K$ be a closed convex subset of a q$q$-uniformly smooth separable Banach space, T:K→K$T:K\to K$ a strictly pseudocontractive mapping, and f:K→K$f:K\to K$ an L$L$-Lispschitzian strongly pseudocontractive mapping. For any t∈(0,1)$t\in (0,1)$, let _xt$x_t$ be the unique fixed point of tf+(1-t)T$\mathit{tf}+(1-t)T$. We prove that if T$T$ has a fixed point, then {_xt}$\{x_t\}$ converges to a fixed point of T$T$ as t$t$ approaches to 0.     

Asymptotic average shadowing property on compact metric spaces

We prove that if for a continuous map f$f$ on a compact metric space X$X$, the chain recurrent set, R(f)$R\left(f\right)$ has more than one chain component, then f$f$ does not satisfy the asymptotic average shadowing property. We also show that if a continuous map f$f$ on a compact metric space X$X$ has the asymptotic average shadowing property and if A$A$ is an attractor for f$f$, then A$A$ is the single attractor for f$f$ and we have A=R(f)$A=R\left(f\right)$. We also study diffeomorphisms with asymptotic average shadowing property and prove that if M$M$ is a compact manifold which is not finite with dimM=2$dimM=2$, then the C¹$C^1$ interior of the set of all C¹$C^1$ diffeomorphisms with the asymptotic average shadowing property is characterized by the set of Ω$\Omega$-stable diffeomorphisms.     

Feedback vertex set on graphs of low clique-width

The Feedback Vertex Set problem asks whether a graph contains q$q$ vertices meeting all its cycles. This is not a local property, in the sense that we cannot check if q$q$ vertices meet all cycles by looking only at their neighbors. Dynamic programming algorithms for problems based on non-local properties are usually more complicated. In this paper, given a graph G$G$ of clique-width cw$cw$ and a cw$cw$-expression of G$G$, we solve the Minimum Feedback Vertex Set problem in time O(n²2^O(cwlogcw))$O\left(n^2{2}^{O(cwlogcw)}\right)$. Our algorithm applies dynamic programming on a so-called k$k$-module decomposition of a graph, as defined by Rao (2008) [29], which is easily derivable from ak$k$-expression of the graph. The related notion of module-width of a graph is tightly linked to both clique-width and NLC-width, and in this paper we give an alternative equivalent characterization of module-width.     

Modules satisfying the weak Nakayama property

Let R$R$ be a commutative ring with identity. We will say that an R$R$-module M$M$ satisfies the weak Nakayama property, if IM=M$IM=M$, where I$I$ is an ideal of R$R$, implies that for any x∈M$x\in M$ there exists a∈I$a\in I$ such that (a−1)x=0$(a-1)x=0$. In this paper, we will study modules satisfying the weak Nakayama property. It is proved that if R$R$ is a local ring, then R$R$ is a Max ring if and only if J(R)$J\left(R\right)$, the Jacobson radical of R$R$, is T$T$-nilpotent if and only if every R$R$-module satisfies the weak Nakayama property.     

Convergence of invariant measures for singular stochastic diffusion equations

It is proved that the solutions to the singular stochastic p$p$-Laplace equation, p∈(1,2)$p\in (1,2)$ and the solutions to the stochastic fast diffusion equation with nonlinearity parameter r∈(0,1)$r\in (0,1)$ on a bounded open domain Λ⊂R^d$\Lambda \subset {\mathbb{R}}^d$ with Dirichlet boundary conditions are continuous in mean, uniformly in time, with respect to the parameters p$p$ and r$r$ respectively (in the Hilbert spaces L²(Λ)$L^2\left(\Lambda \right)$, H⁻¹(Λ)$H^{-1}\left(\Lambda \right)$ respectively). The highly singular limit case p=1$p=1$ is treated with the help of stochastic evolution variational inequalities, where P$\mathbb{P}$-a.s. convergence, uniformly in time, is established.     

Blow-up and propagation of disturbances for fast diffusion equations

This paper is concerned with the Cauchy problem for the fast diffusion equation _ut−Δu^m=αu^_p1$u_t-\Delta u^m=\alpha u^{p_1}$ in R^N${\mathbb{R}}^N$ (N≥1$N\ge 1$), where m∈(0,1)$m\in (0,1)$, _p1>1$p_1>1$ and α>0$\alpha >0$. The initial condition _u0$u_0$ is assumed to be continuous, nonnegative and bounded. Using a technique of subsolutions, we set up sufficient conditions on the initial value _u0$u_0$ so that u(t,x)$u(t,x)$ blows up in finite time, and we show how to get estimates on the profile of u(t,x)$u(t,x)$ for small enough values of t>0$t>0$.     

Iterative approximation of solutions of nonlinear equations of Hammerstein type

Suppose X$X$ is a real q$q$-uniformly smooth Banach space and F,K:X→X$F,K:X\to X$ are Lipschitz ?$?$-strongly accretive maps with D(K)=F(X)=X$D\left(K\right)=F\left(X\right)=X$. Let u^∗$u^{\ast }$ denote the unique solution of the Hammerstein equation u+KFu=0$u+KFu=0$. An iteration process recently introduced by Chidume and Zegeye is shown to converge strongly to u^∗$u^{\ast }$. No invertibility assumption is imposed on K$K$ and the operators K$K$ and F$F$ need not be defined on compact subsets of X$X$. Furthermore, our new technique of proof is of independent interest. Finally, some interesting open questions are included.     

Approximation of solutions of Hammerstein equations with bounded strongly accretive nonlinear operators

Suppose X$X$ is a real q$q$-uniformly smooth Banach space and F,K:X→X$F,K:X\to X$ are bounded strongly accretive maps with D(K)=F(X)=X$D\left(K\right)=F\left(X\right)=X$. Let u^∗$u^{\ast }$ denote the unique solution of the Hammerstein equation u+KFu=0$u+KFu=0$. A new explicit coupled iteration process is shown to converge strongly to u^∗$u^{\ast }$. No invertibility assumption is imposed on K$K$ and the operators K$K$ and F$F$ need not be defined on compact subsets of X$X$. Furthermore, our new technique of proof is of independent interest. Finally, some interesting open questions are included.     

Stable categories of higher preprojective algebras

We introduce (n+1)$(n+1)$-preprojective algebras of algebras of global dimension n$n$. We show that if an algebra is n$n$-representation-finite then its (n+1)$(n+1)$-preprojective algebra is self-injective. In this situation, we show that the stable module category of the (n+1)$(n+1)$-preprojective algebra is (n+1)$(n+1)$-Calabi–Yau, and, more precisely, it is the (n+1)$(n+1)$-Amiot cluster category of the stable n$n$-Auslander algebra of the original algebra. In particular this stable category contains an (n+1)$(n+1)$-cluster tilting object. We show that even if the (n+1)$(n+1)$-preprojective algebra is not self-injective, under certain assumptions (which are always satisfied for n∈{1,2}$n\in \{1,2\}$) the results above still hold for the stable category of Cohen–Macaulay modules.     

On the length of an external branch in the Beta-coalescent

In this paper, we consider Beta(2−α,α)$(2-\alpha ,\alpha )$ (with 1<α<2$1<\alpha <2$) and related Λ$\Lambda$-coalescents. If T⁽ⁿ⁾$T^{\left(n\right)}$ denotes the length of a randomly chosen external branch of the n$n$-coalescent, we prove the convergence of n^α−1T⁽ⁿ⁾$n^{\alpha -1}{T}^{\left(n\right)}$ when n$n$ tends to ∞$\infty$, and give the limit. To this aim, we give asymptotics for the number σ⁽ⁿ⁾${\sigma }^{\left(n\right)}$ of collisions which occur in the n$n$-coalescent until the end of the chosen external branch, and for the block counting process associated with the n$n$-coalescent.     

Unilateral small deviations of processes related to the fractional Brownian motion

Let x(s)$x\left(s\right)$, s∈R^d$s\in R^d$ be a Gaussian self-similar random process of index H$H$. We consider the problem of log-asymptotics for the probability _pT$p_T$ that x(s)$x\left(s\right)$, x(0)=0$x\left(0\right)=0$ does not exceed a fixed level in a star-shaped expanding domain T⋅Δ$T\cdot \Delta$ as T→∞$T\to \infty$. We solve the problem of the existence of the limit, θ?lim(−log_pT)/(logT)^D$\theta ?lim(-log{p}_T)/{(logT)}^D$, T→∞$T\to \infty$, for the fractional Brownian sheet x(s)$x\left(s\right)$, s∈[0,T]²$s\in {[0,T]}^2$ when D=2$D=2$, and we estimate θ$\theta$ for the integrated fractional Brownian motion when D=1$D=1$.     

On testing Hamiltonicity of graphs

Let us fix a function f(n)=o(nlnn)$f\left(n\right)=o(nlnn)$ and real numbers 0≤α<β≤1$0\le \alpha <\beta \le 1$. We present a polynomial time algorithm which, given a directed graph G$G$ with n$n$ vertices, decides either that one can add at most βn$\beta n$ new edges to G$G$ so that G$G$ acquires a Hamiltonian circuit or that one cannot add αn$\alpha n$ or fewer new edges to G$G$ so that G$G$ acquires at least e^−f(n)n!$e^{-f\left(n\right)}n!$ Hamiltonian circuits, or both.     

Parametric inference for discretely observed multidimensional diffusions with small diffusion coefficient

We consider a multidimensional diffusion X$X$ with drift coefficient b(α,_Xt)$b(\alpha ,X_t)$ and diffusion coefficient ?σ(β,_Xt)$?\sigma (\beta ,X_t)$. The diffusion sample path is discretely observed at times _tk=kΔ$t_k=k\Delta$ for k=1…n$k=1\dots n$ on a fixed interval [0,T]$[0,T]$. We study minimum contrast estimators derived from the Gaussian process approximating X$X$ for small ?$?$. We obtain consistent and asymptotically normal estimators of α$\alpha$ for fixed Δ$\Delta$ and ?→0$?\to 0$ and of (α,β)$(\alpha ,\beta )$ for Δ→0$\Delta \to 0$ and ?→0$?\to 0$ without any condition linking ?$?$ and Δ$\Delta$. We compare the estimators obtained with various methods and for various magnitudes of Δ$\Delta$ and ?$?$ based on simulation studies. Finally, we investigate the interest of using such methods in an epidemiological framework.     

Discrete Tomography and plane partitions

A plane partition   is a p×q$p\times q$ matrix A=(_aij)$A=(a_{ij})$, where 1?i?p$1?i?p$ and 1?j?q$1?j?q$, with non-negative integer entries, and whose rows and columns are weakly decreasing. From a geometric point of view plane partitions are equivalent to pyramids  , subsets of the integer lattice Z³${\mathbb{Z}}^3$ which play an important role in Discrete Tomography. As a consequence, some typical problems concerning the tomography of discrete lattice sets can be rephrased and considered via plane partitions. In this paper we focus on some of them. In particular, we get a necessary and sufficient condition for additivity, a canonical procedure for checking the existence of (weakly) bad configurations, and an algorithm which constructs minimal pyramids (with respect to the number of levels) with assigned projection of a bad configurations.