Computing rectangle enclosures

The rectangle enclosure problem is the problem of determining the subset of n iso-oriented planar rectangles that enclose a query rectangle Q. In this paper, we use a three layered data structure which is a combination of Range and Priority search trees and answers both the static and dynamic cases of the problem. Both the cases use O(n> log² n) space. For the static case, the query time is O(log² n log log n + K). The dynamic case is supported in O(log³ n + K) query time using O(log³ n) amortized time per update. K denotes the size of the answer. For the d-dimensional space the results are analogous. The query time is O(log^2d-2 n log log n + K) for the static case and O(log^2d-1 n + K) for the dynamic case. The space used is O(n> log^2d-2 n) and the amortized time for an update is O(log^2d-1 n). The existing bounds given for a class of problems which includes the present one, are O(log^2d n + K) query time, O(log^2d n) time for an insertion and O(log^2d-1 n) time for a deletion.     

High order difference schemes for the system of two space second order nonlinear hyperbolic equations with variable coefficients

In this paper, we develop implicit difference schemes of O(k⁴ + k²h² + h⁴), where k > 0, h > 0 are grid sizes in time and space coordinates, respectively, for solving the system of two space dimensional second order nonlinear hyperbolic partial differential equations with variable coefficients having mixed derivatives subject to appropriate initial and boundary conditions. The proposed difference method for the scalar equation is applied for the solution of wave equation in polar coordinates to obtain three level conditionally stable ADI method of O(k⁴ + k²h² + h⁴). Some physical nonlinear problems are provided to demonstrate the accuracy of the implementation.     

Crank–Nicolson least-squares Galerkin procedures for parabolic integro-differential equations

Two Crank–Nicolson least-squares Galerkin finite element schemes are formulated to solve parabolic integro-differential equations. The advantage of this method is that it is not subject to the LBB condition. The convergence analysis shows that the methods yield the approximate solutions with optimal accuracy in H(div; Ω) × H¹(Ω) and (L²(Ω))² × L²(Ω), respectively. Moreover, the two methods both get the approximate solutions with second-order accuracy in time increment.     

Dynamic motion planning in low obstacle density environments

A fundamental task for an autonomous robot is to plan its own motions. Exact approaches to the solution of this motion planning problem suffer from high worst-case running times. The weak and realistic low obstacle density (L.O.D.) assumption results in linear complexity in the number of obstacles of the free space (Van der Stappen et al., 1997). In this paper we address the dynamic version of the motion planning problem in which a robot moves among polygonal obstacles which move along polylines. The obstacles are assumed to move along constant complexity polylines, and to respect the low density property at any given time. We will show that in this situation a cell decomposition of the free space of size O(n²(n) log² n) can be computed in O(n²(n) log² n) time. The dynamic motion planning problem is then solved in O(n²(n) log³ n) time. We also show that these results are close to optimal.     

Diagonal flips in pseudo-triangulations on closed surfaces

Negami has already shown that there is a natural number N(F²) for any closed surface F² such that two triangulations on F² with n vertices can be transformed into each other by a sequence of diagonal flips if nN(F²). We investigate the same theorem for pseudo-triangulations with or without loops, estimating the length of a sequence of diagonal flips. Our arguments will be applied to simple triangulations to obtain a linear upper bound for N(F²) with respect to the genus of F².     

New constructions of divisible designs

A construction is given for a (p^2a(p+1),p²,p^2a+1(p+1),p^2a+1,p^2a(p+1)) (p a prime) divisible difference set in the group H×Z²_pa+1 where H is any abelian group of order p+1. This can be used to generate a symmetric semi-regular divisible design; this is a new set of parameters for λ₁≠0, and those are fairly rare. We also give a construction for a (p^a−1+p^a−2+…+p+2,p^a+2, p^a(p^a+p^a−1+…+p+1), p^a(p^a−1+…+p+1), p^a−1(p^a+…+p²+2)) divisible difference set in the group H×Z_p2×Z^a_p. This is another new set of parameters, and it corresponds to a symmetric regular divisible design. For p=2, these parameters have λ₁=λ₂, and this corresponds to the parameters for the ordinary Menon difference sets.     

Harnack Differential Inequalities for the Parabolic Equation ut=LF(u) on Riemannian Manifolds and Applications

In this paper, let(M~n, g) be an n-dimensional complete Riemannian manifold with the mdimensional Bakry–mery Ricci curvature bounded below. By using the maximum principle, we first prove a Li–Yau type Harnack differential inequality for positive solutions to the parabolic equation u_t= LF(u)=ΔF(u)-f·F(u),on compact Riemannian manifolds Mn, where F∈C~2(0, ∞), F0 and f is a C~2-smooth function defined on M~n. As application, the Harnack differential inequalities for fast diffusion type equation and porous media type equation are derived. On the other hand, we derive a local Hamilton type gradient estimate for positive solutions of the degenerate parabolic equation on complete Riemannian manifolds. As application, related local Hamilton type gradient estimate and Harnack inequality for fast dfiffusion type equation are established. Our results generalize some known results.     

A theory of electromagnetic and gravitational fields

Careful calculations using classical field theory show that if a macroscopic ball with uniform surface charge (say, a billiard ball with 1E6 excess electrons) is released near the surface of the earth, it will almost instantaneously accelerate to relativistic speed and blow a hole in the ground. This absurd prediction is just the macroscopic version of the self-force problem for charged particles [1]. Furthermore, if one attempts to develop from electromagnetism a parallel theory for gravitational [2], the result is the same, self-acceleration.

The basis of the new theory is a measure of energy density for any wave equation [3–5]. Given any solution of any four-vector wave equation in spacetime (for example, the potentials (c^-1φA)=(A⁰,A¹,A²,A³) in electromagnetism), one can form the 16th first order partial derivatives of the vector components, with respect to the time and space variables (ct,x) = (x⁰, x¹, x², x³). The sum of the squares of the 16 terms is a natural energy function [6, p. 283] (satisfying a conservation law . Such energy functions are routinely utilized by mathematicians as Lyapunov functions in the theory of stability of waves with boundary conditions. A Lagrangian using this sum leads to a new energy tensor for electromagnetic and gravitational fields, an alternative to that in [7].     

Numerical approximation of the product of the square root of a matrix with a vector

Given an n×n symmetric positive definite matrix A and a vector , two numerical methods for approximating are developed, analyzed, and computationally tested. The first method applies a Newton iteration to a specific nonlinear system to approximate while the second method applies a step-control method to numerically solve a specific initial-value problem to approximate . Assuming that A is first reduced to tridiagonal form, the first method requires O(n²) operations per iteration while the second method requires O(n) operations per iteration. In contrast, numerical methods that first approximate A^1/2 and then compute generally require O(n³) operations per iteration.     

Analytic solutions of an iterative functional differential equation

This paper is concerned with a functional differential equation x^′(z)=1/x(az+bx(z)), where a, b are two complex numbers. By constructing a convergent power series solution y(z) of a auxiliary equation of the form b²y^′(z)=(y(μ²z)−ay(μz))(μy^′(μz)−ay^′(z)), analytic solutions of the form for the original differential equation are obtained.     

The Second Neighbourhood for Quasi-transitive Oriented Graphs

In 2006, Sullivan stated the conjectures:(1) every oriented graph has a vertex x such that d~(++)(x) ≥ d~-(x);(2) every oriented graph has a vertex x such that d~(++)(x) + d~+(x) ≥ 2 d~-(x);(3) every oriented graph has a vertex x such that d~(++)(x) + d~+(x) ≥ 2 · min{d~+(x), d~-(x)}. A vertex x in D satisfying Conjecture(i) is called a Sullivan-i vertex, i = 1, 2, 3. A digraph D is called quasi-transitive if for every pair xy, yz of arcs between distinct vertices x, y, z, xz or zx("or" is inclusive here) is in D. In this paper, we prove that the conjectures hold for quasi-transitive oriented graphs, which is a superclass of tournaments and transitive acyclic digraphs. Furthermore, we show that a quasi-transitive oriented graph with no vertex of in-degree zero has at least three Sullivan-1 vertices and a quasi-transitive oriented graph has at least three Sullivan-3 vertices unless it belongs to an exceptional class of quasitransitive oriented graphs. For Sullivan-2 vertices, we show that an extended tournament, a subclass of quasi-transitive oriented graphs and a superclass of tournaments, has at least two Sullivan-2 vertices unless it belongs to an exceptional class of extended tournaments.     

Numerical method for solving 3D inverse scattering problems

Let q(x) L²(D), D R³ is a bounded domain, q = 0 outside D, q is real-valued. Assume that A(\Gj;\t';,\Gj;,k) A(\Gj;\t';,\Gj), the scattering amplitude, is known for all \Gj;|t',\Gj; S², S² is the unit sphere, an d a fixed k \r>0. These data determine q(x) uniquely and a numerical method is given for computing q(x).     

Asymptotic and Partial Asymptotic Hankel Operators on $$H^2(\mathbb{D}^n)$$ H 2 ( D n ) )

In this paper, we generalize the concept of asymptotic Hankel operators on H²(D) to the Hardy space H²(Dⁿ) (over polydisk) in terms of asymptotic Hankel and partial asymptotic Hankel operators and investigate some properties in case of its weak and strong convergence. Meanwhile, we introduce i^th-partial Hankel operators on H²(Dⁿ) and obtain a characterization of its compactness for n > 1. Our main results include the containment of Toeplitz algebra in the collection of all strong partial asymptotic Hankel operators on H²(Dⁿ). It is also shown that a Toeplitz operator with symbol φ is asymptotic Hankel if and only if φ is holomorphic function in L^∞(Tⁿ).     

A stochastic maximum principle for processes driven by fractional Brownian motion

We prove a stochastic maximum principle for controlled processes X(t)=X^(u)(t) of the form

dX(t)=b(t,X(t),u(t)) dt+σ(t,X(t),u(t)) dB^(H)(t),

where B^(H)(t) is m-dimensional fractional Brownian motion with Hurst parameter . As an application we solve a problem about minimal variance hedging in an incomplete market driven by fractional Brownian motion.     

Least-squares Galerkin procedures for parabolic integro-differential equations

Two least-squares Galerkin finite element schemes are formulated to solve parabolic integro-differential equations. The advantage of this method is that it is not subject to the LBB condition. The convergence analysis shows that the least-squares mixed element schemes yield the approximate solution with optimal accuracy in H(div;Ω)×H¹(Ω) and (L²(Ω))²×L²(Ω), respectively.     

On the fourth power moment of Fourier coefficients of cusp form

Let a(n)be the Fourier coefficients of a holomorphic cusp form of weightκ=2n≥12 for the full modular group and A(x)=∑_(n≤x)a(n).In this paper,we establish an asymptotic formula of the fourth power moment of A(x)and prove that ∫T1A~4(x)dx=3/(64κπ~4)s_4;2()T~(2κ)+O(T~(2κ-δ_4+ε))with δ_4=1/8,which improves the previous result.     

Phase Retrieval of Real-valued Functions in Sobolev Space

The Sobolev space H^s(R^d) with s > d/2 contains many important functions such as the bandlimited or rational ones. In this paper we propose a sequence of measurement functions {ϕ_j,k^γ} ⊆ H^-s(R^d) to the phase retrieval problem for the real-valued functions in H^s(R^d). We prove that any real-valued function f ∈ H^s(R^d) can be determined, up to a global sign, by the phaseless measurements {|<f, ϕ_j,k^γ>|}. It is known that phase retrieval is unstable in infinite dimensional spaces with respect to perturbations of the measurement functions. We examine a special type of perturbations that ensures the stability for the phase-retrieval problem for all the real-valued functions in H^s(R^d) ∩ C¹(R^d), and prove that our iterated reconstruction procedure guarantees uniform convergence for any function f ∈ H^s(R^d) ∩ C¹(R^d) whose Fourier transform f is L¹-integrable. Moreover, numerical simulations are conducted to test the efficiency of the reconstruction algorithm.     

Minimal enclosing parallelepiped in 3D

We investigate the problem of finding a minimal volume parallelepiped enclosing a given set of n three-dimensional points. We give two mathematical properties of these parallelepipeds, from which we derive two algorithms of theoretical complexity O(n⁶). Experiments show that in practice our quickest algorithm runs in O(n²) (at least for n10⁵). We also present our application in structural biology.     

On parallel rectilinear obstacle- avoiding paths

We give improved space and processor complexities for the problem of computing, in parallel, a data structure that supports queries about shortest rectilinear obstacle-avoiding paths in the plane, where the obstacles are disjoint rectangles. That is, a query specifies any source and destination in the plane, and the data structure enables efficient processing of the query. We now can build the data structure with O(n²/log n) CREW PRAM processors, as opposed to the previous O(n²), and with O(n²) space, as opposed to the previous O(n²(log n)²). The time complexity remains unchanged, at O((log n)²). As before, the data structure we compute enables a query to be processed in O(log n) time, by one processor for obtaining a path length, or by O(k/log n) processors for retrieving a shortest path itself, where k is the number of segments on that path. The new ideas that made our improvement possible include a new partitioning scheme of the recursion tree, which is used to schedule the computations performed on that tree. Since a number of other related shortest paths problems are solved using this technique as a subroutine our improvement translates into a similar improvement in the complexities of these problems as well.     

Self-maps of <Emphasis Type="Italic">S</Emphasis><Superscript>2</Superscript> homotopic to a smooth map with a single <Emphasis Type="Italic">n</Emphasis>-periodic point

We show for which (d,n) ∈ Z×N there exists a smooth self-map f:S² → S² so that deg(f)=d and Fix(fⁿ) is a point.