Search

International audience; A database of ten type 1 diabetes patients wearing a continuous glucose monitoring device has enabled to record their blood glucose continuous variations every minute all day long during fourteen consecutive days. These recordings represent, for each patient, a time series consisting of 1 value of glycaemia per minute during 24 hours and 14 days, i.e., 20,160 data point. Thus, while using numerical methods, these time series have been anonymously analyzed. Nevertheless, because of the stochastic inputs induced by daily activities of any human being, it has not been possible to discriminate chaos from noise. So, we have decided to keep only the 14 nights of these ten patients. Then, the determination of the time delay and embedding dimension according to the delay coordinate embedding method has allowed us to estimate for each patient the correlation dimension and the maximal Lyapunov exponent. This has led us to show that type 1 diabetes could indeed be a chaotic phenomenon. Once this result has been confirmed by the determinism test, we have computed the Lyapunov time and found that the limit of predictability of this phenomenon is nearly equal to half the 90-minutes sleep-dream cycle. We hope that our results will prove to be useful to characterize and predict blood glucose variations.

We deal with the obstacle problem for the porous medium equation in the slow diffusion regime m > 1 . Our main interest is to treat fairly irregular obstacles assuming only boundedness and lower semicontinuity. In particular, the considered obstacles are not regular enough to work with the classical notion of variational solutions, and a different approach is needed. We prove the existence of a solution in the sense of the minimal supersolution lying above the obstacle. As a consequence, we can show that non-negative weak supersolutions to the porous medium equation can be approximated by a sequence of supersolutions which are bounded away from zero.

We study ideal lattices in R 2 coming from real quadratic fields, and give an explicit method for computing all well-rounded twists of any such ideal lattice. We apply this to ideal lattices coming from Markoff numbers to construct infinite families of non-equivalent planar lattices with good sphere-packing radius and good minimum product distance. We also provide a complete classification of all real quadratic fields such that the orthogonal lattice is the only well-rounded twist of the lattice corresponding to the ring of integers.

The Contextuality-by-Default theory is illustrated on contextuality analysis of the idealized double-slit experiment. The experiment is described by a system of contextually labeled binary random variables each of which answers the question: has the particle hit the detector, having passed through a given slit (left or right) in a given state (open or closed)? This system of random variables is a cyclic system of rank 4, formally the same as the system describing the EPR/Bell paradigm with signaling. Unlike the latter, however, the system describing the double-slit experiment is always noncontextual, i.e., the context-dependence in it is entirely explainable in terms of direct influences of contexts (closed-open arrangements of the slits) upon the marginal distributions of the random variables involved. The analysis presented is entirely within the framework of abstract classical probability theory (with contextually labeled random variables). The only physical constraint used in the analysis is that a particle cannot pass through a closed slit. The noncontextuality of the double-slit system does not generalize to systems describing experiments with more than two slits: in an abstract triple-slit system, almost any set of observable detection probabilities is compatible with both a contextual scenario and a noncontextual scenario of the particle passing though various combinations of open and closed slits (although the issue of physical realizability of these scenarios remains open)., Comment: Entropy 2018, 20, 278; doi:10.3390/e20040278

The eigenproblem for thin shells of revolution under uncertainty in material parameters is discussed. Here the focus is on the smallest eigenpairs. Shells of revolution have natural eigenclusters due to symmetries, moreover, the eigenpairs depend on a deterministic parameter, the dimensionless thickness. The stochastic subspace iteration algorithms presented here are capable of resolving the smallest eigenclusters. In the case of random material parameters, it is possible that the eigenmodes cross in the stochastic parameter space. This interesting phenomenon is demonstrated via numerical experiments. Finally, the effect of the chosen material model on the asymptotics in relation to the deterministic parameter is shown to be negligible.

For systems with hidden attractors and unstable equilibria, the property that hidden attractors are not connected with unstable equilibria is now accepted as one of their main characteristics. To the best of our knowledge this property has not been explored using realtime interactive three-dimensions graphics. Aided by advanced computer graphic analysis, in this paper, we explore this characteristic of a particular nonlinear system with very rich and unusual dynamics, the Rabinovich-Fabrikant system. It is shown that there exists a neighborhood of one of the unstable equilibria within which the initial conditions do not lead to the considered hidden chaotic attractor, but to one of the stable equilibria or are divergent. The trajectories starting from any neighborhood of the other unstable equlibria are attracted either by the stable equilibria, or are divergent, Comment: International Journal of Bifurcation and Chaos, accepted 2018

This paper considers approximate smoothing for discretely observed non-linear stochastic differential equations. The problem is tackled by developing methods for linearising stochastic differential equations with respect to an arbitrary Gaussian process. Two methods are developed based on (1) taking the limit of statistical linear regression of the discretised process and (2) minimising an upper bound to a cost functional. Their difference is manifested in the diffusion of the approximate processes. This in turn gives novel derivations of pre-existing Gaussian smoothers when Method 1 is used and a new class of Gaussian smoothers when Method 2 is used. Furthermore, based on the aforementioned development the iterative Gaussian smoothers in discrete-time are generalised to the continuous-time setting by iteratively re-linearising the stochastic differential equation with respect to the current Gaussian process approximation to the smoothed process. The method is verified in two challenging tracking problems, a reentry problem and a radar tracked coordinated turn model with state dependent diffusion. The results show that the method has competitive estimation accuracy with state-of-the-art smoothers.

Let $\Omega$ be an internal chord-arc Jordan domain and $\varphi:\mathbb S\rightarrow\partial\Omega$ be a homeomorphism. We show that $\varphi$ has finite dyadic energy if and only if $\varphi$ has a diffeomorphic extension $h: \mathbb D\rightarrow \Omega$ which has finite energy., Comment: 19 pages, 1 figure

We consider so-called regular invertible Gaussian Volterra processes and derive a formula for their prediction laws. Examples of such processes include the fractional Brownian motions and the mixed fractional Brownian motions. As an application, we consider conditional-mean hedging under transaction costs in Black-Scholes type pricing models where the Brownian motion is replaced with a more general regular invertible Gaussian Volterra process., Comment: arXiv admin note: text overlap with arXiv:1706.01534

In this article we study the large $N$ asymptotics of complex moments of the absolute value of the characteristic polynomial of a $N\times N$ complex Ginibre random matrix with the characteristic polynomial evaluated at a point in the unit disk. More precisely, we calculate the large $N$ asymptotics of $\mathbb{E}|\det(G_N-z)|^{\gamma}$, where $G_N$ is a $N\times N$ matrix whose entries are i.i.d and distributed as $N^{-1/2}Z$, $Z$ being a standard complex Gaussian, $\mathrm{Re}(\gamma)>-2$, and $|z, Comment: Result and proof extended to include complex-valued gamma; some reference added

This work reformulates the complete electrode model of electrical impedance tomography in order to enable more efficient numerical solution. The model traditionally assumes constant contact conductances on all electrodes, which leads to a discontinuous Robin boundary condition since the gaps between the electrodes can be described by vanishing conductance. As a consequence, the regularity of the electromagnetic potential is limited to less than two square-integrable weak derivatives, which negatively affects the convergence of, e.g., the finite element method. In this paper, a smoothened model for the boundary conductance is proposed, and the unique solvability and improved regularity of the ensuing boundary value problem are proven. Numerical experiments demonstrate that the proposed model is both computationally feasible and also compatible with real-world measurements. In particular, the new model allows faster convergence of the finite element method., Comment: 23 pages, 8 figures

The second author introduced with I. T\"orm\"a a two-player word-building game [Playing with Subshifts, Fund. Inform. 132 (2014), 131--152]. The game has a predetermined (possibly finite) choice sequence $\alpha_1$, $\alpha_2$, $\ldots$ of integers such that on round $n$ the player $A$ chooses a subset $S_n$ of size $\alpha_n$ of some fixed finite alphabet and the player $B$ picks a letter from the set $S_n$. The outcome is determined by whether the word obtained by concatenating the letters $B$ picked lies in a prescribed target set $X$ (a win for player $A$) or not (a win for player $B$). Typically, we consider $X$ to be a subshift. The winning shift $W(X)$ of a subshift $X$ is defined as the set of choice sequences for which $A$ has a winning strategy when the target set is the language of $X$. The winning shift $W(X)$ mirrors some properties of $X$. For instance, $W(X)$ and $X$ have the same entropy. Virtually nothing is known about the structure of the winning shifts of subshifts common in combinatorics on words. In this paper, we study the winning shifts of subshifts generated by marked uniform substitutions, and show that these winning shifts, viewed as subshifts, also have a substitutive structure. Particularly, we give an explicit description of the winning shift for the generalized Thue-Morse substitutions. It is known that $W(X)$ and $X$ have the same factor complexity. As an example application, we exploit this connection to give a simple derivation of the first difference and factor complexity functions of subshifts generated by marked substitutions. We describe these functions in particular detail for the generalized Thue-Morse substitutions., Comment: Extended version of a paper presented at RuFiDiM IV

Given a surjective mapping $f : E \to F$ between Banach spaces, we investigate the existence of a subspace $G$ of $E$, with the same density character as $F$, such that the restriction of $f$ to $G$ remains surjective. We obtain a positive answer whenever $f$ is continuous and uniformly open. In the smooth case, we deduce a positive answer when $f$ is a $C^1$-smooth surjection whose set of critical values is countable. Finally we show that, when $f$ takes values in the Euclidean space $\mathbb R^n$, in order to obtain this result it is not sufficient to assume that the set of critical values of $f$ has zero-measure., Comment: 11 pages, a small mistake in the proof of Theorem 8 has been fixed

The geometry of the set of restrictions of rank-one tensors to some of their coordinates is studied. This gives insight into the problem of rank-one completion of partial tensors. Particular emphasis is put on the semialgebraic nature of the problem, which arises for real tensors with constraints on the parameters. The algebraic boundary of the completable region is described for tensors parametrized by probability distributions and where the number of observed entries equals the number of parameters. If the observations are on the diagonal of a tensor of format $d\times\dots\times d$, the complete semialgebraic description of the completable region is found., Comment: 24 pages, 3 figures; v2: Final version, accepted in SIAM Journal of Applied Algebra and Geometry (SIAGA)

We give an explicit description of fundamental domains associated to the $p$-adic uniformisation of families of Shimura curves of discriminant $Dp$ and level $N\geq 1$, for which the one-sided ideal class number $h(D,N)$ is $1$. The obtained results generalise those in \cite[Ch. IX]{Gerritzen_vanderPut1980} for Shimura curves of discriminant $2p$ and level $N=1$. The method we present here enables us to find Mumford curves covering Shimura curves, together with a free system of generators for the associated Schottky groups, $p$-adic good fundamental domains and their stable reduction-graphs. This is based on a detailed study of the modular arithmetic of an Eichler order of level $N$ inside the definite quaternion algebra of discriminant $D$, for which we generalise classical results of Hurwitz \cite{Hurwitz1896}. As an application, we prove general formulas for the reduction-graphs with lengths at $p$ of the considered families of Shimura curves., Comment: article, 39 pages, 3 tables, 8 figures. This is second version of the paper, which will appear on Trans. of the AMS. Conjecture in the Introduction becomes a theorem in this version

Markerless tracking of hands and fingers is a promising enabler for human-computer interaction. However, adoption has been limited because of tracking inaccuracies, incomplete coverage of motions, low framerate, complex camera setups, and high computational requirements. In this paper, we present a fast method for accurately tracking rapid and complex articulations of the hand using a single depth camera. Our algorithm uses a novel detection-guided optimization strategy that increases the robustness and speed of pose estimation. In the detection step, a randomized decision forest classifies pixels into parts of the hand. In the optimization step, a novel objective function combines the detected part labels and a Gaussian mixture representation of the depth to estimate a pose that best fits the depth. Our approach needs comparably less computational resources which makes it extremely fast (50 fps without GPU support). The approach also supports varying static, or moving, camera-to-scene arrangements. We show the benefits of our method by evaluating on public datasets and comparing against previous work., Comment: 9 pages, Accepted version of paper published at CVPR 2015

A word w is rich if it has |w|+1 many distinct palindromic factors, including the empty word. A word is square-free if it does not have a factor uu, where u is a non-empty word. Pelantov\'a and Starosta (Discrete Math. 313 (2013)) proved that every infinite rich word contains a square. We will give another proof for that result. Pelantov\'a and Starosta denoted by r(n) the length of a longest rich square-free word on an alphabet of size n. The exact value of r(n) was left as an open question. We will give an upper and a lower bound for r(n), and make a conjecture that our lower bound is exact. We will also generalize the notion of repetition threshold for a limited class of infinite words. The repetition thresholds for episturmian and rich words are left as an open question., Comment: 25 pages

This paper presents the generalized optimal sub-pattern assignment (GOSPA) metric on the space of finite sets of targets. Compared to the well-established optimal sub-pattern assignment (OSPA) metric, GOSPA is unnormalized as a function of the cardinality and it penalizes cardinality errors differently, which enables us to express it as an optimisation over assignments instead of permutations. An important consequence of this is that GOSPA allows us to penalize localization errors for detected targets and the errors due to missed and false targets, as indicated by traditional multiple target tracking (MTT) performance measures, in a sound manner. In addition, we extend the GOSPA metric to the space of random finite sets, which is important to evaluate MTT algorithms via simulations in a rigorous way., Comment: The paper received the Jean Pierre Le Cadre best paper award at the 20th International Conference on Information Fusion, July 2017. A Matlab implementation of the proposed GOSPA metric is available in https://github.com/abusajana/GOSPA Also visit https://youtu.be/M79GTTytvCM for a 15-min presentation about the paper

We consider the normalized $p$-Poisson problem $$-��^N_p u=f \qquad \text{in}\quad ��.$$ The normalized $p$-Laplacian $��_p^{N}u:=|D u|^{2-p}��_p u$ is in non-divergence form and arises for example from stochastic games. We prove $C^{1,��}_{loc}$ regularity with nearly optimal $��$ for viscosity solutions of this problem. In the case $f\in L^{\infty}\cap C$ and $p>1$ we use methods both from viscosity and weak theory, whereas in the case $f\in L^q\cap C$, $q>\max(n,\frac p2,2)$, and $p>2$ we rely on the tools of nonlinear potential theory., 41 pages

The norm in classical Sobolev spaces can be expressed as a difference quotient. This expression can be used to generalize the space to the fractional smoothness case. Because the difference quotient is based on shifting the function, it cannot be used in generalized Orlicz spaces. In its place, we introduce a smoothed difference quotient and show that it can be used to characterize the generalized Orlicz-Sobolev space. Our results are new even in Orlicz spaces and variable exponent spaces., Comment: 18 pages

In this paper, we develop {finite-time horizon} causal filters using the nonanticipative rate distortion theory. We apply the {developed} theory to {design optimal filters for} time-varying multidimensional Gauss-Markov processes, subject to a mean square error fidelity constraint. We show that such filters are equivalent to the design of an optimal \texttt{\{encoder, channel, decoder\}}, which ensures that the error satisfies {a} fidelity constraint. Moreover, we derive a universal lower bound on the mean square error of any estimator of time-varying multidimensional Gauss-Markov processes in terms of conditional mutual information. Unlike classical Kalman filters, the filter developed is characterized by a reverse-waterfilling algorithm, which ensures {that} the fidelity constraint is satisfied. The theoretical results are demonstrated via illustrative examples., Comment: 35 pages, 6 figures, submitted for publication in SIAM Journal on Control and Optimization (SICON)

We show that if $X$ is a uniformly perfect complete metric space satisfying the finite doubling property, then there exists a fully supported measure with lower regularity dimension as close to the lower dimension of $X$ as we wish. Furthermore, we show that, under the condensation open set condition, the lower dimension of an inhomogeneous self-similar set $E_C$ coincides with the lower dimension of the condensation set $C$, while the Assouad dimension of $E_C$ is the maximum of the Assouad dimensions of the corresponding self-similar set $E$ and the condensation set $C$. If the Assouad dimension of $C$ is strictly smaller than the Assouad dimension of $E$, then the upper regularity dimension of any measure supported on $E_C$ is strictly larger than the Assouad dimension of $E_C$. Surprisingly, the corresponding statement for the lower regularity dimension fails.

We study the Goldbach problem for primes represented by the polynomial $x^2+y^2+1$. The set of such primes is sparse in the set of all primes, but the infinitude of such primes was established by Linnik. We prove that almost all even integers $n$ satisfying certain necessary local conditions are representable as the sum of two primes of the form $x^2+y^2+1$. This improves a result of Matom\"aki, which tells that almost all even $n$ satisfying a local condition are the sum of one prime of the form $x^2+y^2+1$ and one generic prime. We also solve the analogous ternary Goldbach problem, stating that every large odd $n$ is the sum of three primes represented by our polynomial. As a byproduct of the proof, we show that the primes of the form $x^2+y^2+1$ contain infinitely many three term arithmetic progressions, and that the numbers $\alpha p \pmod 1$ with $\alpha$ irrational and $p$ running through primes of the form $x^2+y^2+1$, are distributed rather uniformly., Comment: 49 pages; Referee comments incorporated; To appear in Mathematika

In preprocessing tensor-valued data, e.g. images and videos, a common procedure is to vectorize the observations and subject the resulting vectors to one of the many methods used for independent component analysis (ICA). However, the tensor structure of the original data is lost in the vectorization and, as a more suitable alternative, we propose the matrix- and tensor fourth order blind identification (MFOBI and TFOBI). In these tensorial extensions of the classic fourth order blind identification (FOBI) we assume a Kronecker structure for the mixing and perform FOBI simultaneously on each direction of the observed tensors. We discuss the theory and assumptions behind MFOBI and TFOBI and provide two different algorithms and related estimates of the unmixing matrices along with their asymptotic properties. Finally, simulations are used to compare the method's performance with that of classical FOBI for vectorized data and we end with a real data clustering example., Comment: 26 pages, 4 figures

Locally repairable codes (LRCs) are error correcting codes used in distributed data storage. Besides a global level, they enable errors to be corrected locally, reducing the need for communication between storage nodes. There is a close connection between almost affine LRCs and matroid theory which can be utilized to construct good LRCs and derive bounds on their performance. A generalized Singleton bound for linear LRCs with parameters $(n,k,d,r,\delta)$ was given in [N. Prakash et al., "Optimal Linear Codes with a Local-Error-Correction Property", IEEE Int. Symp. Inf. Theory]. In this paper, a LRC achieving this bound is called perfect. Results on the existence and nonexistence of linear perfect $(n,k,d,r,\delta)$-LRCs were given in [W. Song et al., "Optimal locally repairable codes", IEEE J. Sel. Areas Comm.]. Using matroid theory, these existence and nonexistence results were later strengthened in [T. Westerb\"ack et al., "On the Combinatorics of Locally Repairable Codes", Arxiv: 1501.00153], which also provided a general lower bound on the maximal achievable minimum distance $d_{\rm{max}}(n,k,r,\delta)$ that a linear LRC with parameters $(n,k,r,\delta)$ can have. This article expands the class of parameters $(n,k,d,r,\delta)$ for which there exist perfect linear LRCs and improves the lower bound for $d_{\rm{max}}(n,k,r,\delta)$. Further, this bound is proved to be optimal for the class of matroids that is used to derive the existence bounds of linear LRCs., Comment: 5 pages, final version, submitted to the 2016 IEEE International Symposium on Information Theory, presents the results of arXiv:1512.05325 in an abridged form

We discretize the Lagrange multiplier formulation of the obstacle problem by mixed and stabilized finite element methods. A priori and a posteriori error estimates are derived and numerically verified.

Nonlinear analysis of the phase-locked loop (PLL) based circuits is a challenging task, thus in modern engineering literature simplified mathematical models and simulation are widely used for their study. In this work the limitations of numerical approach is discussed and it is shown that, e.g. hidden oscillations may not be found by simulation. Corresponding examples in SPICE and MatLab, which may lead to wrong conclusions concerning the operability of PLL-based circuits, are presented.

In this paper, we derive an ambiguity function (AF) for the transmit beamspace (TB)-based multipleinput multiple-output (MIMO) radar for the case of far-field targets and narrow-band waveforms. The effects of transmit coherent processing gain and waveform diversity are incorporated into the AF definition. To cover all the phase information conveyed by different factors, we introduce the equivalent transmit phase centers. The newly defined AF serves as a generalized AF form for which the phased-array (PA) and traditional MIMO radar AFs are important special cases. We establish relationships among the defined TB-based MIMO radar AF and the existing AF results including the Woodward's AF, the AFs defined for the traditional colocated MIMO radar, and also the PA radar AF, respectively. Moreover, we compare the TB-based MIMO radar AF with the square-summation-form AF definition and identify two limiting cases to bound its 'clear region' in Doppler-delay domain that is free of sidelobes. Corresponding bounds for these two cases are derived, and it is shown that the bound for the worst case is inversely proportional to the number of transmitted waveforms K, whereas the bound for the best case is independent of K. The actual 'clear region' of the TB-based MIMO radar AF depends on the array configuration and is in between of the worst- and best-case bounds. We propose a TB design strategy to reduce the levels of the AF sidelobes, and show in simulations that proper design of the TB matrix leads to reduction of the relative sidelobe levels of the TB-based MIMO radar AF., Comment: 28 pages, 4 figures, Submitted to IEEE Trans. Signal Processing on July 2014

We establish both sufficient and necessary conditions for the validity of the so-called Hardy-Sobolev inequalities on open sets of the Euclidean space. These inequalities form a natural interpolating scale between the (weighted) Sobolev inequalities and the (weighted) Hardy inequalities. The Assouad dimension of the complement of the open set turns out to play an important role in both sufficient and necessary conditions., Comment: 27 pages

Volume of metric balls relates to rate-distortion theory and packing bounds on codes. In this paper, the volume of balls in complex Grassmann manifolds is evaluated for an arbitrary radius. The ball is defined as a set of hyperplanes of a fixed dimension with reference to a center of possibly different dimension, and a generalized chordal distance for unequal dimensional subspaces is used. First, the volume is reduced to one-dimensional integral representation. The overall problem boils down to evaluating a determinant of a matrix of the same size as the subspace dimensionality. Interpreting this determinant as a characteristic function of the Jacobi ensemble, an asymptotic analysis is carried out. The obtained asymptotic volume is moreover refined using moment-matching techniques to provide a tighter approximation in finite-size regimes. Lastly, the pertinence of the derived results is shown by rate-distortion analysis of source coding on Grassmann manifolds., Comment: 26 pages

The exact Lyapunov dimension formula for the Lorenz system has been analytically obtained first due to G.A.Leonov in 2002 under certain restrictions on parameters, permitting classical values. He used the construction technique of special Lyapunov-type functions developed by him in 1991 year. Later it was shown that the consideration of larger class of Lyapunov-type functions permits proving the validity of this formula for all parameters of the system such that all the equilibria of the system are hyperbolically unstable. In the present work it is proved the validity of the formula for Lyapunov dimension for a wider variety of parameters values, which include all parameters satisfying the classical physical limitations. One of the motivation of this work is the possibility of computing a chaotic attractor in the Lorenz system in the case of one unstable and two stable equilibria., Comment: arXiv admin note: text overlap with arXiv:1505.04729

We study the Cauchy problem for a nonlocal heat equation, which is of fractional order both in space and time. We prove four main theorems: (i) a representation formula for classical solutions, (ii) a quantitative decay rate at which the solution tends to the fundamental solution, (iii) optimal L 2 -decay of mild solutions in all dimensions, (iv) L 2 -decay of weak solutions via energy methods. The first result relies on a delicate analysis of the definition of classical solutions. After proving the representation formula we carefully analyze the integral representation to obtain the quantitative decay rates of (ii). Next we use Fourier analysis techniques to obtain the optimal decay rate for mild solutions. Here we encounter the critical dimension phenomenon where the decay rate attains the decay rate of that in a bounded domain for large enough dimensions. Consequently, the decay rate does not anymore improve when the dimension increases. The theory is markedly different from that of the standard caloric functions and this substantially complicates the analysis. Finally, we use energy estimates and a comparison principle to prove a quantitative decay rate for weak solutions defined via a variational formulation. Our main idea is to show that the L 2 -norm is actually a subsolution to a purely time-fractional problem which allows us to use the known theory to obtain the result.

In order to harvest the business potential of device-to-device (D2D) communication, direct communication between devices subscribed to different mobile operators should be supported. This would also support meeting requirements resulting from D2D relevant scenarios, like vehicle-to-vehicle communication. In this paper, we propose to allocate the multi-operator D2D communication over dedicated cellular spectral resources contributed from both operators. Ideally, the operators should negotiate about the amount of spectrum to contribute, without revealing proprietary information to each other and/or to other parties. One possible way to do that is to use the sequence of operators' best responses, i.e., the operators make offers about the amount of spectrum to contribute using a sequential updating procedure until reaching consensus. Besides spectrum allocation, we need a mode selection scheme for the multi-operator D2D users. We use a stochastic geometry framework to capture the impact of mode selection on the distribution of D2D users and assess the performance of the best response iteration algorithm. With the performance metrics considered in the paper, we show that the best response iteration has a unique Nash equilibrium that can be reached from any initial strategy. In general, asymmetric operators would contribute unequal amounts of spectrum for multi-operator D2D communication. Provided that the multi-operator D2D density is not negligible, we show that both operators may experience significant performance gains as compared to the scheme without spectrum sharing., Comment: To appear in IEEE ICC 2015

Some sharp inequalities between the triangular ratio metric and the Cassinian metric are proved in the unit ball., Comment: 15 pages, 2 figures

Let $\Omega$ be a domain in $\mathbf R^d$ and $h(\varphi)=\sum^d_{k,l=1}(\partial_k\varphi, c_{kl}\partial_l\varphi)$ a quadratic form on $L_2(\Omega)$ with domain $C_c^\infty(\Omega)$ where the $c_{kl}$ are real symmetric $L_\infty(\Omega)$-functions with $C(x)=(c_{kl}(x))>0$ for almost all $x\in \Omega$. Further assume there are $a, \delta>0$ such that $a^{-1}d_\Gamma^{\delta}\,I\le C\le a\,d_\Gamma^{\delta}\,I$ for $d_\Gamma\le 1$ where $d_\Gamma$ is the Euclidean distance to the boundary $\Gamma$ of $\Omega$. We assume that $\Gamma$ is Ahlfors $s$-regular and if $s$, the Hausdorff dimension of $\Gamma$, is larger or equal to $d-1$ we also assume a mild uniformity property for $\Omega$ in the neighbourhood of one $z\in\Gamma$. Then we establish that $h$ is Markov unique, i.e. it has a unique Dirichlet form extension, if and only if $\delta\ge 1+(s-(d-1))$. The result applies to forms on Lipschitz domains or on a wide class of domains with $\Gamma$ a self-similar fractal. In particular it applies to the interior or exterior of the von Koch snowflake curve in $\mathbf R^2$ or the complement of a uniformly disconnected set in $\mathbf R^d$., Comment: 25 pages, 2 figures

International audience; We present an algorithm to solve BSDEs with jumps based on Wiener Chaos Expansion and Picard's iterations. This paper extends the results given in Briand-Labart (2014) to the case of BSDEs with jumps. We get a forward scheme where the conditional expectations are easily computed thanks to chaos decomposition formulas. Concerning the error, we derive explicit bounds with respect to the number of chaos, the discretization time step and the number of Monte Carlo simulations. We also present numerical experiments. We obtain very encouraging results in terms of speed and accuracy.

We extend the De Giorgi-Nash-Moser theory to nonlocal, possibly degenerate integro-differential operators., Comment: 26 pages. To appear in Ann. Inst. H. Poincare Anal. Non Lineaire. arXiv admin note: text overlap with arXiv:1405.7842

The set of matrices of given positive semidefinite rank is semialgebraic. In this paper we study the geometry of this set, and in small cases we describe its boundary. For general values of positive semidefinite rank we provide a conjecture for the description of this boundary. Our proof techniques are geometric in nature and rely on nesting spectrahedra between polytopes., Comment: 23 pages

This work considers finding optimal positions for the electrodes within the Bayesian paradigm based on available prior information on the conductivity; the aim is to place the electrodes so that the posterior density of the (discretized) conductivity, i.e., the conditional density of the conductivity given the measurements, is as localized as possible. To make such an approach computationally feasible, the complete electrode forward model of impedance tomography is linearized around the prior expectation of the conductivity, allowing explicit representation for the (approximate) posterior covariance matrix. Two approaches are considered: minimizing the trace or the determinant of the posterior covariance. The introduced optimization algorithm is of the steepest descent type, with the needed gradients computed based on appropriate Fr\'echet derivatives of the complete electrode model. The functionality of the methodology is demonstrated via two-dimensional numerical experiments., Comment: 22 pages, 6 figures

We introduce the quasiminimal subshifts, subshifts having only finitely many subsystems. With $\mathbb{N}$-actions, their theory essentially reduces to the theory of minimal systems, but with $\mathbb{Z}$-actions, the class is much larger. We show many examples of such subshifts, and in particular construct a universal system with only a single proper subsystem, refuting a conjecture of [Delvenne, K\r{u}rka, Blondel, '05]., Comment: 40 pages, 1 figure, submitted to JCSS

This paper considers space-time coding over several independently Rayleigh faded blocks. In particular, we will concentrate on giving upper bounds for the coding gain of lattice space-time codes as the number of blocks grow. This problem was previously considered in the single antenna case by Bayer-Fluckiger et al. in 2006. Crucial to their work was Odlyzko’s bound on the discriminant of an algebraic number field, as this provides an upper bound for the normalized coding gain of number field codes. In the MIMO context natural codes are constructed from division algebras defined over number fields and the coding gain is measured by the discriminant of the corresponding (noncommutative) algebra. In this paper, we will develop analogues of the Odlyzko bounds in this context and show how these bounds limit the normalized coding gain of a very general family of division algebra based space-time codes. These bounds can also be used as benchmarks in practical code design and as tools to analyze asymptotic bounds of performance as the number of independently faded blocks increases.

This paper proposes an analysis of the effects of consensus and preference aggregation on the consistency of pairwise comparisons. We define some boundary properties for the inconsistency of group preferences and investigate their relation with different inconsistency indices. Some results are presented on more general dependencies between properties of inconsistency indices and the satisfaction of boundary properties. In the end, given three boundary properties and nine indices among the most relevant ones, we will be able to present a complete analysis of what indices satisfy what properties and offer a reflection on the interpretation of the inconsistency of group preferences., Comment: 21 pages, 6 figures

We investigate systems of degenerate parabolic equations idealizing reactive solute transport in porous media. Taking advantage of the inherent structure of the system that allows to deduce a scalar Generalized Porous Medium Equation for the sum of the solute concentrations, we show existence of a unique weak solution to the coupled system and derive regularity estimates. We also prove that the system supports solutions propagating with finite speed thus giving rise to free boundaries and interaction of compactly supported initial concentrations of different species.

For any open set $\Omega\subset\mathbb R^n$ and $n\ge 2$, we establish everywhere differentiability of viscosity solutions to the Aronsson equation $$ =0 \quad \rm in\ \ \Omega, $$ where $H$ is given by $$H(x,\,p)==\sum_{i,\,j=1}^na^{ij}(x)p_i p_j,\ x\in\Omega, \ p\in\mathbb R^n, $$ and $A=(a^{ij}(x))\in C^{1,1}(\bar\Omega,\mathbb R^{n\times n})$ is uniformly elliptic. This extends an earlier theorem by Evans and Smart \cite{es11a} on infinity harmonic functions., Comment: 24 pages

We present a new, dynamical way to study powers (that is, repetitions) in Sturmian words based on results from Diophantine approximation theory. As a result, we provide an alternative and shorter proof of a result by Damanik and Lenz characterizing powers in Sturmian words [Powers in Sturmian sequences, Eur. J. Combin. 24 (2003), 377--390]. Further, as a consequence, we obtain a previously known formula for the fractional index of a Sturmian word based on the continued fraction expansion of its slope., Comment: 9 pages, 1 figure

We construct examples to show the sharpness of uniform continuity of quasiconformal mappings onto $s$-John domains. Our examples also give a negative answer to a prediction in [7]., Comment: 16 pages, 7 figures

We develop the canonical Volterra representation for a self-similar Gaussian process by using the Lamperti transformation of the corresponding stationary Gaussian process, where this latter one admits a canonical integral representation under the assumption of pure non-determinism. We apply the representation obtained for the self-similar Gaussian process to derive an expression for Gaussian processes that are equivalent in law to the self-similar Gaussian process in question.

User performance on a wireless network depends on whether a neighboring cochannel interferer applies a single (spatial) stream or a multi stream transmission. This work analyzes the impact of interference rank on a beamforming and orthogonal space-time block coded (OSTBC) user transmission. We generalize existing analytical results on signal-to-interference-plus-noise-ratio (SINR) distribution and outage probability under arbitrary number of unequal power interferers. We show that higher rank interference causes lower outage probability, and can support better outage threshold especially in the case of beamforming., Comment: Revised after first round of reviews

In the setting of a metric space equipped with a doubling measure and supporting a Poincare inequality, and based on results by Bjorn and Shanmugalingam (2007) [7] , we show that functions of bounded variation can be extended from any bounded uniform domain to the whole space. Closely related to extensions is the concept of boundary traces, which have previously been studied by Hakkarainen et al. (2014) [17] . On spaces that satisfy a suitable locality condition for sets of finite perimeter, we establish some basic results for the traces of functions of bounded variation. Our analysis of traces also produces novel pointwise results on the behavior of functions of bounded variation in their jump sets.

In this note we establish the boundedness properties of local maximal operators MG on the fractional Sobolev spaces Ws;p(G) whenever G is an open set in Rn, 0 < s < 1 and 1 < p < 1. As an application, we characterize the fractional (s;p)-Hardy inequality on a bounded open set by a Maz'ya-type testing condition localized to Whitney cubes. pq(G) whenever G is an open set in R n , 0 < s < 1 and 1 < p;q