A new lattice Boltzmann model for the laplace equation

In this paper, a new lattice Boltzmann equation which is independent of time is proposed. Based on the new lattice Boltzmann equation, some steady problems can be modeled by the lattice Boltzmann method. In the further study, the Laplace equation is investigated with the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different space scales. The numerical results show that the new method is effective.     

A novel lattice Boltzmann model for the Poisson equation

In this paper, a novel lattice Boltzmann model is proposed to solve the Poisson equation through modifying equilibrium distribution function. Compared with previous models, which can be viewed as the solvers to diffusion equation, the present model is a genuine solver to the Poisson equation, and the transient term derived by previous models is eliminated. Numerical solutions agree well with analytical solutions, which indicates the potential of the present model for solving the Poisson equation.     

A note on lattice chains and Delannoy numbers

Fix nonnegative integers n₁,…,n_d and let L denote the lattice of integer points (a₁,…,a_d)∈Z^d satisfying 0?a_i?n_i for 1?i?d. Let L be partially ordered by the usual dominance ordering. In this paper we offer combinatorial derivations of a number of results concerning chains in L. In particular, the results obtained are established without recourse to generating functions or recurrence relations. We begin with an elementary derivation of the number of chains in L of a given size, from which one can deduce the classical expression for the total number of chains in L. Then we derive a second, alternative, expression for the total number of chains in L when d=2. Setting n₁=n₂ in this expression yields a new proof of a result of Stanley [Enumerative Combinatorics, vol. 2, Cambridge University Press, Cambridge, 1999] relating the total number of chains to the central Delannoy numbers. We also conjecture a generalization of Stanley's result to higher dimensions.     

Valuations on lattice polytopes

Let L be a lattice (that is, a Z-module of finite rank), and let L=P(L) denote the family of convex polytopes with vertices in L; here, convexity refers to the underlying rational vector space V=Q⊗L. In this paper it is shown that any valuation on L satisfies the inclusion-exclusion principle, in the strong sense that appropriate extension properties of the valuation hold. Indeed, the core result is that the class of a lattice polytope in the abstract group L=P(L) for valuations on L can be identified with its characteristic function in V. In fact, the same arguments are shown to apply to P(M), when M is a module of finite rank over an ordered ring, and more generally to appropriate families of (not necessarily bounded) polyhedra.     

Stepsize analysis for descent methods

The convergence rates of descent methods with different stepsize rules are compared. Among the stepsize rules considered are: constant stepsize, exact minimization along a line, Goldstein-Armijo rules, and stepsize equal to that which yields the minimum of certain interpolatory polynomials. One of the major results shown is that the rate of convergence of descent methods with the Goldstein-Armijo stepsize rules can be made as close as desired to the rate of convergence of methods that require exact minimization along a line. Also, a descent algorithm that combines a Goldstein-Armijo stepsize rule with a secant-type step is presented. It is shown that this algorithm has a convergence rate equal to the convergence of descent methods that require exact minimization along a line and that, eventually (i.e., near the minimum), it does not require a search to determine an acceptable stepsize.     

A note on lattice congruence regularization

We give an example that answers in the negative Problem 11 of G. Gr?tzer and E. T. Schmidt, Regular congruence-preserving extensions of lattices. Received November 24, 1999; accepted in final form October 16, 2000.     

On the counting function of the lattice profile of periodic sequences

The lattice profile analyzes the intrinsic structure of pseudorandom number sequences with applications in Monte Carlo methods and cryptology. In this paper, using the discrete Fourier transform for periodic sequences and the relation between the lattice profile and the linear complexity, we give general formulas for the expected value, variance, and counting function of the lattice profile of periodic sequences with fixed period. Moreover, we determine in a more explicit form the expected value, variance, and counting function of the lattice profile of periodic sequences for special values of the period.     

Parallel iterative methods using factorized preconditioning matrices for solving elliptic equations on triangular grids

Parallel analogs of the variants of the incomplete Cholesky-conjugate gradient method and the modified incomplete Cholesky-conjugate gradient method for solving elliptic equations on uniform triangular and unstructured triangular grids on parallel computer systems with the MIMD architecture are considered. The construction of parallel methods is based on the use of various variants of ordering the grid points depending on the decomposition of the computation domain. Results of the theoretic and experimental studies of the convergence rate of these methods are presented. The solution of model problems on a moderate number processors is used to examine the efficiency of the proposed parallel methods.     

Iterative methods based on spline approximations to detect discontinuities from Fourier data

Recently, spline approximations have been proposed for the reconstruction of piecewise smooth functions from Fourier data. That approach makes possible to retrieve the functions from their Fourier coefficients for any given degree of accuracy when the discontinuity points are known. In this paper we present iterative methods based on those spline approximations, for several degrees, to find locations and amplitudes of the jumps of a piecewise smooth function, given its Fourier coefficients. We also present numerical experiments comparing with different previous approaches.     

Bayesian stopping rules for multistart global optimization methods

By far the most efficient methods for global optimization are based on starting a local optimization routine from an appropriate subset of uniformly distributed starting points. As the number of local optima is frequently unknown in advance, it is a crucial problem when to stop the sequence of sampling and searching. By viewing a set of observed minima as a sample from a generalized multinomial distribution whose cells correspond to the local optima of the objective function, we obtain the posterior distribution of the number of local optima and of the relative size of their regions of attraction. This information is used to construct sequential Bayesian stopping rules which find the optimal trade off between reliability and computational effort.     

Existence and construction of shifted lattice rules with an arbitrary number of points and bounded weighted star discrepancy for general decreasing weights

We study the problem of constructing shifted rank-1 lattice rules for the approximation of high-dimensional integrals with a low weighted star discrepancy, for classes of functions having bounded weighted variation, where the weighted variation is defined as the weighted sum of Hardy–Krause variations over all lower-dimensional projections of the integrand. Under general conditions on the weights, we prove the existence of rank-1 lattice rules such that for any δ>0, the general weighted star discrepancy is O(n^−1+δ) for any number of points n>1 (not necessarily prime), any shift of the lattice, general (decreasing) weights, and uniformly in the dimension. We also show that these rules can be constructed by a component-by-component strategy. This implies in particular that a single infinite-dimensional generating vector can be used for integrals in any number of dimensions, and even for infinite-dimensional integrands when they have bounded weighted variation. These same lattices are also good with respect to the worst-case error in weighted Korobov spaces with the same types of general weights. Similar results were already available for various special cases, such as general weights and prime n, or arbitrary n and product weights, but not for the most general combination of n composite, general weights, arbitrary shift, and star discrepancy, considered here. Our results imply tractability or strong tractability of integration for classes of integrands with finite weighted variation when the weights satisfy the conditions we give. These classes are a strict superset of those covered by earlier sufficient tractability conditions.     

A lattice point problem on the regular tree

Huber (1956) [8] considered the following problem on the hyperbolic plane H. Consider a strictly hyperbolic subgroup of automorphisms on H with compact quotient, and choose a conjugacy class in this group. Count the number of vertices inside an increasing ball, which are images of a fixed point x∈H under automorphisms in the chosen conjugacy class, and describe the asymptotic behaviour of this number as the size of the ball goes to infinity. We use a well-known analogy between the hyperbolic plane and the regular tree to solve this problem on the regular tree.     

Enhancing PSO methods for global optimization

The Particle Swarm Optimization (PSO) method is a well-established technique for global optimization. During the past years several variations of the original PSO have been proposed in the relevant literature. Because of the increasing necessity in global optimization methods in almost all fields of science there is a great demand for efficient and fast implementations of relative algorithms. In this work we propose three modifications of the original PSO method in order to increase the speed and its efficiency that can be applied independently in almost every PSO variant. These modifications are: (a) a new stopping rule, (b) a similarity check and (c) a conditional application of some local search method. The proposed were tested using three popular PSO variants and a variety test functions. We have found that the application of these modifications resulted in significant gain in speed and efficiency.     

Periodization strategy may fail in high dimensions

We discuss periodization of smooth functions f of d variables for approximation of multivariate integrals. The benefit of periodization is that we may use lattice rules, which have recently seen significant progress. In particular, we know how to construct effectively a generator of the rank-1 lattice rule with n points whose worst case error enjoys a nearly optimal bound C _d,p n ^−p . Here C _d,p is independent of d or depends at most polynomially on d, and p can be arbitrarily close to the smoothness of functions belonging to a weighted Sobolev space with an appropriate condition on the weights. If F denotes the periodization for f then the error of the lattice rule for a periodized function F is bounded by C _d,p n ^−p ∣∣F∣∣ with the norm of F given in the same Sobolev space. For small or moderate d, the norm of F is not much larger than the norm of f. This means that for small or moderate d, periodization is successful and allows us to use optimal properties of lattice rules also for non-periodic functions. The situation is quite different if d is large since the norm of F can be exponentially larger than the norm of f. This can already be seen for f = 1. Hence, the upper bound of the worst case error of the lattice rule for periodized functions is quite bad for large d. We conjecture not only that this upper bound is bad, but also that all lattice rules fail for large d. That is, if we fix the number of points n and let d go to infinity then the worst case error of any lattice rule is bounded from below by a positive constant independent of n. We present a number of cases suggesting that this conjecture is indeed true, but the most interesting case, when the sum of the weights of the corresponding Sobolev space is bounded in d, remains open.      

A note on the generalized symmetries of lattice hierarchies

We show that the recursion operators of the integrable lattice equations usually considered in the literature can also be used to generate hierarchies of differential-delay equations. All members of these hierarchies of lattice and differential-delay equations commute. It is thus seen that differential-delay hierarchies provide a broader context within which to place lattice hierarchies.     

Cubature Formula and Interpolation on the Cubic Domain

Several cubature formulas on the cubic domains are derived using the dis-crete Fourier analysis associated with lattice tiling, as developed in [10]. The main results consist of a new derivation of the Gaussian type cubature for the product Cheby-shev weight functions and associated interpolation polynomials on [-1,1]2, as well as new results on [-1,1]3. In particular, compact formulas for the fundamental interpo-lation polynomials are derived, based on n3/4 + (n2) nodes of a cubature formula on [-1,1]3.