Weak coupling of solutions of first-order least-squares method

A theoretical analysis of a first-order least-squares finite element method for second-order self-adjoint elliptic problems is presented. We investigate the coupling effect of the approximate solutions for the primary function and for the flux . We prove that the accuracy of the approximate solution for the primary function is weakly affected by the flux . That is, the bound for is dependent on , but only through the best approximation for multiplied by a factor of meshsize . Similarly, we provide that the bound for is dependent on , but only through the best approximation for multiplied by a factor of the meshsize . This weak coupling is not true for the non-selfadjoint case. We provide the numerical experiment supporting the theorems in this paper.

     

Computations of Eisenstein series on Fuchsian groups

We present numerical investigations of the value distribution and distribution of Fourier coefficients of the Eisenstein series on arithmetic and non-arithmetic Fuchsian groups. Our numerics indicate a Gaussian limit value distribution for a real-valued rotation of as , and also, on non-arithmetic groups, a complex Gaussian limit distribution for when near and , at least if we allow at some rate. Furthermore, on non-arithmetic groups and for fixed with near , our numerics indicate a Gaussian limit distribution for the appropriately normalized Fourier coefficients.

     

Computational estimation of the constant

The paper describes a computational estimation of the constant characterizing the bounds of . It is known that as

with , while the truth of the Riemann hypothesis would also imply that . In the range , two sets of estimates of are computed, one for increasingly small minima and another for increasingly large maxima of . As increases, the estimates in the first set rapidly fall below and gradually reach values slightly below , while the estimates in the second set rapidly exceed and gradually reach values slightly above . The obtained numerical results are discussed and compared to the implications of recent theoretical work of Granville and Soundararajan.

     

Short effective intervals containing primes in arithmetic progressions and the seven cubes problem

For any and any non-exceptional modulus , we prove that, for large enough ( ), the interval contains a prime in any of the arithmetic progressions modulo . We apply this result to establish that every integer larger than is a sum of seven cubes.

     

Sign changes in sums of the Liouville function

The Liouville function is the completely multiplicative function whose value is at each prime. We develop some algorithms for computing the sum , and use these methods to determine the smallest positive integer where . This answers a question originating in some work of Turán, who linked the behavior of to questions about the Riemann zeta function. We also study the problem of evaluating Pólya's sum , and we determine some new local extrema for this function, including some new positive values.

     

The hyperdeterminant and triangulations of the 4-cube

The hyperdeterminant of format is a polynomial of degree in unknowns which has terms. We compute the Newton polytope of this polynomial and the secondary polytope of the -cube. The regular triangulations of the -cube are classified into -equivalence classes, one for each vertex of the Newton polytope. The -cube has coarsest regular subdivisions, one for each facet of the secondary polytope, but only of them come from the hyperdeterminant.

     

Equal moments division of a set

Let be the minimal positive integer , for which there exists a splitting of the set into  subsets, , , ..., , whose first moments are equal. Similarly, let be the maximal positive integer , such that there exists a splitting of into subsets whose first moments are equal. For , these functions were investigated by several authors, and the values of and have been found for and , respectively. In this paper, we deal with the problem for any prime . We demonstrate our methods by finding for any and for .

     

On the irreducibility of Hecke polynomials

Let be the characteristic polynomial of the th Hecke operator acting on the space of cusp forms of weight for the full modular group. We record a simple criterion which can be used to check the irreducibility of the polynomials . Using this criterion with some machine computation, we show that if there exists such that is irreducible and has the full symmetric group as Galois group, then the same is true of for each prime .

     

Links and cubic 3-polytopes

It is well known that a prime link diagram corresponds to a signed plane graph without cut vertices (Kauffman, 1989). In this paper, we present a new relation between prime links and cubic 3-polytopes. Let be the set of links such that each has a diagram whose corresponding signed plane graph is the graph of a cubic 3-polytope. We show that all nontrivial prime links, except -torus links and -pretzel links, can be obtained from by using some operation of untwining. Furthermore, we define the generalized cubic 3-polytope chains and then show that any nontrivial link can be obtained from by some untwining operations, where is the set of links corresponding to generalized cubic 3-polytope chains. These results are used to simplify the computation of the Kauffman brackets of links so that the computing can be done in a unified way for many infinite families of links.

     

Rationality problem of three-dimensional purely monomial group actions: the last case

A -automorphism of the rational function field is called purely monomial if sends every variable to a monic Laurent monomial in the variables . Let be a finite subgroup of purely monomial -automorphisms of . The rationality problem of the -action is the problem of whether the -fixed field is -rational, i.e., purely transcendental over , or not. In 1994, M. Hajja and M. Kang gave a positive answer for the rationality problem of the three-dimensional purely monomial group actions except one case. We show that the remaining case is also affirmative.

     

New expansions of numerical eigenvalues for by nonconforming elements

The paper explores new expansions of the eigenvalues for in with Dirichlet boundary conditions by the bilinear element (denoted ) and three nonconforming elements, the rotated bilinear element (denoted ), the extension of (denoted ) and Wilson's elements. The expansions indicate that and provide upper bounds of the eigenvalues, and that and Wilson's elements provide lower bounds of the eigenvalues. By extrapolation, the convergence rate can be obtained, where is the maximal boundary length of uniform rectangles. Numerical experiments are carried out to verify the theoretical analysis made.

     

Parity-regular Steinhaus graphs

Steinhaus graphs on vertices are certain simple graphs in bijective correspondence with binary -sequences of length . A conjecture of Dymacek in 1979 states that the only nontrivial regular Steinhaus graphs are those corresponding to the periodic binary sequences of any length . By an exhaustive search the conjecture was known to hold up to 25 vertices. We report here that it remains true up to 117 vertices. This is achieved by considering the weaker notion of parity-regular Steinhaus graphs, where all vertex degrees have the same parity. We show that these graphs can be parametrized by an -vector space of dimension approximately and thus constitute an efficiently describable domain where true regular Steinhaus graphs can be searched by computer.

     

The generalized triangular decomposition

Given a complex matrix , we consider the decomposition , where is upper triangular and and have orthonormal columns. Special instances of this decomposition include the singular value decomposition (SVD) and the Schur decomposition where is an upper triangular matrix with the eigenvalues of on the diagonal. We show that any diagonal for can be achieved that satisfies Weyl's multiplicative majorization conditions:

where is the rank of , is the -th largest singular value of , and is the -th largest (in magnitude) diagonal element of . Given a vector which satisfies Weyl's conditions, we call the decomposition , where is upper triangular with prescribed diagonal , the generalized triangular decomposition (GTD). A direct (nonrecursive) algorithm is developed for computing the GTD. This algorithm starts with the SVD and applies a series of permutations and Givens rotations to obtain the GTD. The numerical stability of the GTD update step is established. The GTD can be used to optimize the power utilization of a communication channel, while taking into account quality of service requirements for subchannels. Another application of the GTD is to inverse eigenvalue problems where the goal is to construct matrices with prescribed eigenvalues and singular values.

     

Explicit values of multi-dimensional Kloosterman sums for prime powers, II

For any integer fix , and let denote the group of reduced residues modulo . Let , a power of a prime . The hyper-Kloosterman sums of dimension are defined for by

where denotes the multiplicative inverse of modulo .

Salie evaluated in the classical setting for even , and for odd with . Later, Smith provided formulas that simplified the computation of in these cases for . Recently, Cochrane, Liu and Zheng computed upper bounds for in the general case , stopping short of their explicit evaluation. Here I complete the computation they initiated to obtain explicit values for the Kloosterman sums for , relying on basic properties of some simple specialized exponential sums. The treatment here is more elementary than the author's previous determination of these Kloosterman sums using character theory and -adic methods. At the least, it provides an alternative, independent evaluation of the Kloosterman sums.

     

A counterexample concerning the -projector onto linear spline spaces

For the -orthogonal projection onto spaces of linear splines over simplicial partitions in polyhedral domains in , , we show that in contrast to the one-dimensional case, where independently of the nature of the partition, in higher dimensions the -norm of cannot be bounded uniformly with respect to the partition. This fact is folklore among specialists in finite element methods and approximation theory but seemingly has never been formally proved.

     

Fast algorithms for computing isogenies between elliptic curves

We survey algorithms for computing isogenies between elliptic curves defined over a field of characteristic either 0 or a large prime. We introduce a new algorithm that computes an isogeny of degree ( different from the characteristic) in time quasi-linear with respect to . This is based in particular on fast algorithms for power series expansion of the Weierstrass -function and related functions.

     

On the distinctness of modular reductions of maximal length sequences modulo odd prime powers

We discuss the distinctness problem of the reductions modulo of maximal length sequences modulo powers of an odd prime , where the integer has a prime factor different from . For any two different maximal length sequences generated by the same polynomial, we prove that their reductions modulo are distinct. In other words, the reduction modulo of a maximal length sequence is proved to contain all the information of the original sequence.

     

Monotonicity of some functions involving the gamma and psi functions

Let , where is Euler's gamma function. We determine conditions for the numbers so that the function is strongly completely monotonic on . Through this result we obtain some inequalities involving the ratio of gamma functions and provide some applications in the context of trigonometric sum estimation. We also give several other examples of strongly completely monotonic functions defined in terms of and functions. Some limiting and particular cases are also considered.

     

The Euler implicit/explicit scheme for the 2D time-dependent Navier-Stokes equations with smooth or non-smooth initial data

This paper considers the stability and convergence results for the Euler implicit/explicit scheme applied to the spatially discretized two-dimensional (2D) time-dependent Navier-Stokes equations. A Galerkin finite element spatial discretization is assumed, and the temporal treatment is implicit/explict scheme, which is implicit for the linear terms and explicit for the nonlinear term. Here the stability condition depends on the smoothness of the initial data , i.e., the time step condition is in the case of , in the case of and in the case of for mesh size and some positive constant . We provide the -stability of the scheme under the stability condition with and obtain the optimal error estimate of the numerical velocity and the optimal error estimate of the numerical pressure under the stability condition with .

     

Comparison theorems of Kolmogorov type and exact values of -widths on Hardy-Sobolev classes

Let be a strip in complex plane. denotes those -periodic, real-valued functions on which are analytic in the strip and satisfy the condition , . Osipenko and Wilderotter obtained the exact values of the Kolmogorov, linear, Gel'fand, and information -widths of in , , and 2-widths of in , , .

In this paper we continue their work. Firstly, we establish a comparison theorem of Kolmogorov type on , from which we get an inequality of Landau-Kolmogorov type. Secondly, we apply these results to determine the exact values of the Gel'fand -width of in , . Finally, we calculate the exact values of Kolmogorov -width, linear -width, and information -width of in , , .