Computing weight modular forms of level

For a prime we describe an algorithm for computing the Brandt matrices giving the action of the Hecke operators on the space of modular forms of weight and level . For we define a special Hecke stable subspace of which contains the space of modular forms with CM by the ring of integers of and we describe the calculation of the corresponding Brandt matrices.

     

The convergence of harmonic Ritz values, harmonic Ritz vectors, and refined harmonic Ritz vectors

This paper concerns a harmonic projection method for computing an approximation to an eigenpair of a large matrix . Given a target point and a subspace that contains an approximation to , the harmonic projection method returns an approximation to . Three convergence results are established as the deviation of from approaches zero. First, the harmonic Ritz value converges to if a certain Rayleigh quotient matrix is uniformly nonsingular. Second, the harmonic Ritz vector converges to if the Rayleigh quotient matrix is uniformly nonsingular and remains well separated from the other harmonic Ritz values. Third, better error bounds for the convergence of are derived when converges. However, we show that the harmonic projection method can fail to find the desired eigenvalue --in other words, the method can miss if it is very close to . To this end, we propose to compute the Rayleigh quotient of with respect to and take it as a new approximate eigenvalue. is shown to converge to once tends to , no matter how is close to . Finally, we show that if the Rayleigh quotient matrix is uniformly nonsingular, then the refined harmonic Ritz vector, or more generally the refined eigenvector approximation introduced by the author, converges. We construct examples to illustrate our theory.

     

On the nonexistence of -cycles for the problem

This article generalizes a proof of Steiner for the nonexistence of -cycles for the problem to a proof for the nonexistence of -cycles. A lower bound for the cycle length is derived by approximating the ratio between numbers in a cycle. An upper bound is found by applying a result of Laurent, Mignotte, and Nesterenko on linear forms in logarithms. Finally numerical calculation of convergents of shows that -cycles cannot exist.

     

Finding -strong pseudoprimes

Let be odd primes and . Put

Then we call the kernel, the triple the signature, and the height of , respectively. We call a -number if it is a Carmichael number with each prime factor . If is a -number and a strong pseudoprime to the bases for , we call a -spsp . Since -numbers have probability of error (the upper bound of that for the Rabin-Miller test), they often serve as the exact values or upper bounds of (the smallest strong pseudoprime to all the first prime bases). If we know the exact value of , we will have, for integers , a deterministic efficient primality testing algorithm which is easy to implement.

In this paper, we first describe an algorithm for finding -spsp(2)'s, to a given limit, with heights bounded. There are in total -spsp's with heights . We then give an overview of the 21978 - spsp(2)'s and tabulate of them, which are -spsp's to the first prime bases up to ; three numbers are spsp's to the first 11 prime bases up to 31. No -spsp's to the first prime bases with heights were found. We conjecture that there exist no -spsp's to the first prime bases with heights and so that

which was found by the author in an earlier paper. We give reasons to support the conjecture. The main idea of our method for finding those -spsp's is that we loop on candidates of signatures and kernels with heights bounded, subject those candidates of -spsp's and their prime factors to Miller's tests, and obtain the desired numbers. At last we speed our algorithm for finding larger -spsp's, say up to , with a given signature to more prime bases. Comparisons of effectiveness with Arnault's and our previous methods for finding -strong pseudoprimes to the first several prime bases are given.

     

-class groups of certain extensions of degree

Let be an odd prime number. In this article we study the distribution of -class groups of cyclic number fields of degree , and of cyclic extensions of degree of an imaginary quadratic field whose class number is coprime to . We formulate a heuristic principle predicting the distribution of the -class groups as Galois modules, which is analogous to the Cohen-Lenstra heuristics concerning the prime-to--part of the class group, although in our case we have to fix the number of primes that ramify in the extensions considered. Using results of Gerth we are able to prove part of this conjecture. Furthermore, we present some numerical evidence for the conjecture.

     

Solutions of the congruence

To supplement existing data, solutions of are tabulated for primes with and . For , five new solutions 2^{32}$"> are presented. One of these, for , also satisfies the ``reverse' congruence . An effective procedure for searching for such ``double solutions' is described and applied to the range , . Previous to this, congruences are generally considered for any and fixed prime to see where the smallest prime solution occurs.

     

Improved methods and starting values to solve the matrix equations iteratively

The two matrix iterations are known to converge linearly to a positive definite solution of the matrix equations , respectively, for known choices of and under certain restrictions on . The convergence for previously suggested starting matrices is generally very slow. This paper explores different initial choices of in both iterations that depend on the extreme singular values of and lead to much more rapid convergence. Further, the paper offers a new algorithm for solving the minus sign equation and explores mixed algorithms that use Newton's method in part.

     

Minimality and other properties of the width- nonadjacent form

Let be an integer and let be the set of integers that includes zero and the odd integers with absolute value less than . Every integer can be represented as a finite sum of the form , with , such that of any consecutive 's at most one is nonzero. Such representations are called width- nonadjacent forms (-NAFs). When these representations use the digits and coincide with the well-known nonadjacent forms. Width- nonadjacent forms are useful in efficiently implementing elliptic curve arithmetic for cryptographic applications. We provide some new results on the -NAF. We show that -NAFs have a minimal number of nonzero digits and we also give a new characterization of the -NAF in terms of a (right-to-left) lexicographical ordering. We also generalize a result on -NAFs and show that any base 2 representation of an integer, with digits in , that has a minimal number of nonzero digits is at most one digit longer than its binary representation.

     

Predicting nonlinear pseudorandom number generators

Let be a prime and let and be elements of the finite field of elements. The inversive congruential generator (ICG) is a sequence of pseudorandom numbers defined by the relation . We show that if sufficiently many of the most significant bits of several consecutive values of the ICG are given, one can recover the initial value (even in the case where the coefficients and are not known). We also obtain similar results for the quadratic congruential generator (QCG), , where . This suggests that for cryptographic applications ICG and QCG should be used with great care. Our results are somewhat similar to those known for the linear congruential generator (LCG), , but they apply only to much longer bit strings. We also estimate limits of some heuristic approaches, which still remain much weaker than those known for LCG.

     

On the greatest prime factor of with effective constants

Let denote a prime. In this article we provide the first published lower bounds for the greatest prime factor of exceeding in which the constants are effectively computable. As a result we prove that it is possible to calculate a value such that for every x_0$"> there is a with the greatest prime factor of exceeding . The novelty of our approach is the avoidance of any appeal to Siegel's Theorem on primes in arithmetic progression.

     

Two kinds of strong pseudoprimes up to

Let be an odd composite integer. Write with odd. If either mod or mod for some , then we say that is a strong pseudoprime to base , or spsp() for short. Define to be the smallest strong pseudoprime to all the first prime bases. If we know the exact value of , we will have, for integers , a deterministic efficient primality testing algorithm which is easy to implement. Thanks to Pomerance et al. and Jaeschke, the are known for . Conjectured values of were given by us in our previous papers (Math. Comp. 72 (2003), 2085-2097; 74 (2005), 1009-1024).

The main purpose of this paper is to give exact values of for ; to give a lower bound of : ; and to give reasons and numerical evidence of K2- and -spsp's to support the following conjecture: for any , where (resp. ) is the smallest K2- (resp. -) strong pseudoprime to all the first prime bases. For this purpose we describe procedures for computing and enumerating the two kinds of spsp's to the first 9 prime bases. The entire calculation took about 4000 hours on a PC Pentium IV/1.8GHz. (Recall that a K2-spsp is an spsp of the form: with primes and ; and that a -spsp is an spsp and a Carmichael number of the form: with each prime factor mod .)

     

On the absolute Mahler measure of polynomials having all zeros in a sector. II

Let be an algebraic integer of degree , not or a root of unity, all of whose conjugates are confined to a sector . In the paper On the absolute Mahler measure of polynomials having all zeros in a sector, G. Rhin and C. Smyth compute the greatest lower bound of the absolute Mahler measure ( of , for belonging to nine subintervals of . In this paper, we improve the result to thirteen subintervals of and extend some existing subintervals.

     

Rational double points on supersingular surfaces

We investigate configurations of rational double points with the total Milnor number on supersingular surfaces. The complete list of possible configurations is given. As an application, we also give the complete list of extremal (quasi-)elliptic fibrations on supersingular surfaces.

     

A finite dimensional realization of the mollifier method for compact operator equations

We introduce and analyze a stable procedure for the approximation of where is the least residual norm solution of the minimal norm of the ill-posed equation , with compact operator between Hilbert spaces, and has some smoothness assumption. Our method is based on a finite number of singular values of and some finite rank operators. Our results are in a more general setting than the one considered by Rieder and Schuster (2000) and Nair and Lal (2003) with special reference to the mollifier method, and it is also applicable under fewer smoothness assumptions on .

     

Visible evidence for the Birch and Swinnerton-Dyer conjecture for modular abelian varieties of analytic rank zero

This paper provides evidence for the Birch and Swinnerton-Dyer conjecture for analytic rank  abelian varieties  that are optimal quotients of attached to newforms. We prove theorems about the ratio , develop tools for computing with , and gather data about certain arithmetic invariants of the nearly abelian varieties of level . Over half of these have analytic rank , and for these we compute upper and lower bounds on the conjectural order of  . We find that there are at least such for which the Birch and Swinnerton-Dyer conjecture implies that is divisible by an odd prime, and we prove for of these that the odd part of the conjectural order of really divides by constructing nontrivial elements of using visibility theory. We also give other evidence for the conjecture. The appendix, by Cremona and Mazur, fills in some gaps in the theoretical discussion in their paper on visibility of Shafarevich-Tate groups of elliptic curves.

     

Elliptic curves with nonsplit mod representations

We calculate explicitly the -invariants of the elliptic curves corresponding to rational points on the modular curve by giving an expression defined over of the -function in terms of the function field generators and of the elliptic curve . As a result we exhibit infinitely many elliptic curves over with nonsplit mod representations.

     

Real zeros of Dedekind zeta functions of real quadratic fields

Let be a primitive, real and even Dirichlet character with conductor , and let be a positive real number. An old result of H. Davenport is that the cycle sums are all positive at and this has the immediate important consequence of the positivity of . We extend Davenport's idea to show that in fact for , 0$"> for all with so that one can deduce the positivity of by the nonnegativity of a finite sum for any . A simple algorithm then allows us to prove numerically that has no positive real zero for a conductor up to 200,000, extending the previous record of 986 due to Rosser more than 50 years ago. We also derive various estimates explicit in of the as well as the shifted cycle sums considered previously by Leu and Li for . These explicit estimates are all rather tight and may have independent interests.

     

Polyharmonic splines on grids and their limits

Radial Basis Functions (RBF) have found a wide area of applications. We consider the case of polyharmonic RBF (called sometimes polyharmonic splines) where the data are on special grids of the form having practical importance. The main purpose of the paper is to consider the behavior of the polyharmonic interpolation splines on such grids for the limiting process 0.$"> For a large class of data functions defined on it turns out that there exists a limit function This limit function is shown to be a polyspline of order on strips. By the theory of polysplines we know that the function is smooth up to order everywhere (in particular, they are smooth on the hyperplanes , which includes existence of the normal derivatives up to order while the RBF interpolants are smooth only up to the order The last fact has important consequences for the data smoothing practice.

     

Pointwise error estimates of the local discontinuous Galerkin method for a second order elliptic problem

In this paper we derive some pointwise error estimates for the local discontinuous Galerkin (LDG) method for solving second-order elliptic problems in (). Our results show that the pointwise errors of both the vector and scalar approximations of the LDG method are of the same order as those obtained in the norm except for a logarithmic factor when the piecewise linear functions are used in the finite element spaces. Moreover, due to the weighted norms in the bounds, these pointwise error estimates indicate that when at least piecewise quadratic polynomials are used in the finite element spaces, the errors at any point depend very weakly on the true solution and its derivatives in the regions far away from . These localized error estimates are similar to those obtained for the standard conforming finite element method.