Coding the principal character formula for affine Kac-Moody lie algebras

In this paper, an algorithm for computing the principal character for affine Lie algebras is discussed and presented. The principal characters discovered using this program are given and/or proven. Results include level 2 and 3 character formulas in and the sole existence of the Rogers-Ramanujan products in , , , , , , .     

Computation of Stark-Tamagawa units

Let be a totally real number field and let denote an odd prime number. We design an algorithm which computes strong numerical evidence for the validity of the ``Equivariant Tamagawa Number Conjecture' for the -equivariant motive , where is a cyclic extension of degree and group . This conjecture is a very deep refinement of the classical analytic class number formula. In the course of the algorithm, we compute a set of special units which must be considered as a generalization of the (conjecturally existing) Stark units associated to first order vanishing Dirichlet -functions.

     

On the total number of prime factors of an odd perfect number

We say is perfect if , where denotes the sum of the positive divisors of . No odd perfect numbers are known, but it is well known that if such a number exists, it must have prime factorization of the form , where , , ..., are distinct primes and . We prove that if or for all , , then . We also prove as our main result that , where . This improves a result of Sayers given in 1986.

     

Linear algebra algorithms for divisors on an algebraic curve

We use an embedding of the symmetric th power of any algebraic curve of genus into a Grassmannian space to give algorithms for working with divisors on , using only linear algebra in vector spaces of dimension , and matrices of size . When the base field is finite, or if has a rational point over , these give algorithms for working on the Jacobian of that require field operations, arising from the Gaussian elimination. Our point of view is strongly geometric, and our representation of points on the Jacobian is fairly simple to deal with; in particular, none of our algorithms involves arithmetic with polynomials. We note that our algorithms have the same asymptotic complexity for general curves as the more algebraic algorithms in Florian Hess' 1999 Ph.D. thesis, which works with function fields as extensions of . However, for special classes of curves, Hess' algorithms are asymptotically more efficient than ours, generalizing other known efficient algorithms for special classes of curves, such as hyperelliptic curves (Cantor 1987), superelliptic curves (Galbraith, Paulus, and Smart 2002), and curves (Harasawa and Suzuki 2000); in all those cases, one can attain a complexity of .

     

Error bounds for Gauss-Turán quadrature formulae of analytic functions

We study the kernels of the remainder term of Gauss-Turán quadrature formulas

for classes of analytic functions on elliptical contours with foci at , when the weight is one of the special Jacobi weights ; ; , ; , . We investigate the location on the contour where the modulus of the kernel attains its maximum value. Some numerical examples are included.

     

Finite element superconvergence on Shishkin mesh for 2-D convection-diffusion problems

In this work, the bilinear finite element method on a Shishkin mesh for convection-diffusion problems is analyzed in the two-dimensional setting. A superconvergence rate in a discrete -weighted energy norm is established under certain regularity assumptions. This convergence rate is uniformly valid with respect to the singular perturbation parameter . Numerical tests indicate that the rate is sharp for the boundary layer terms. As a by-product, an -uniform convergence of the same order is obtained for the -norm. Furthermore, under the same regularity assumption, an -uniform convergence of order in the norm is proved for some mesh points in the boundary layer region.

     

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.

     

Class numbers of some abelian extensions of rational function fields

Let be a monic irreducible polynomial. In this paper we generalize the determinant formula for of Bae and Kang and the formula for of Jung and Ahn to any subfields of the cyclotomic function field By using these formulas, we calculate the class numbers of all subfields of when and are small.

     

Finding strong pseudoprimes to several bases. II

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 . Upper bounds for were first given by Jaeschke, and those for were then sharpened by the first author in his previous paper (Math. Comp. 70 (2001), 863-872).

In this paper, we first follow the first author's previous work to use biquadratic residue characters and cubic residue characters as main tools to tabulate all strong pseudoprimes (spsp's) to the first five or six prime bases, which have the form with odd primes and ; then we tabulate all Carmichael numbers , to the first six prime bases up to 13, which have the form with each prime factor . There are in total 36 such Carmichael numbers, 12 numbers of which are also spsp's to base 17; 5 numbers are spsp's to bases 17 and 19; one number is an spsp to the first 11 prime bases up to 31. As a result the upper bounds for and are lowered from 20- and 22-decimal-digit numbers to a 19-decimal-digit number:

We conjecture that

and give reasons to support this conjecture. The main idea for finding these Carmichael numbers is that we loop on the largest prime factor and propose necessary conditions on to be a strong pseudoprime to the first prime bases. Comparisons of effectiveness with Arnault's, Bleichenbacher's, Jaeschke's, and Pinch's methods for finding (Carmichael) numbers with three prime factors, which are strong pseudoprimes to the first several prime bases, are given.

     

Using the theory of cyclotomy to factor cyclotomic polynomials over finite fields

We examine the problem of factoring the th cyclotomic polynomial, over , and distinct primes. Given the traces of the roots of we construct the coefficients of in time . We demonstrate a deterministic algorithm for factoring in time when has precisely two irreducible factors. Finally, we present a deterministic algorithm for computing the sum of the irreducible factors of in time .

     

A computational approach for solving

We present a computational approach for finding all integral solutions of the equation for even values of . By reducing this problem to that of finding integral solutions of a certain class of quartic equations closely related to the Pell equations, we are able to apply the powerful computational machinery related to quadratic number fields. Using our approach, we determine all integral solutions for assuming the Generalized Riemann Hypothesis, and for unconditionally.

     

Reducing the construction cost of the component-by-component construction of good lattice rules

The construction of randomly shifted rank- lattice rules, where the number of points is a prime number, has recently been developed by Sloan, Kuo and Joe for integration of functions in weighted Sobolev spaces and was extended by Kuo and Joe and by Dick to composite numbers. To construct -dimensional rules, the shifts were generated randomly and the generating vectors were constructed component-by-component at a cost of operations. Here we consider the situation where is the product of two distinct prime numbers and . We still generate the shifts randomly but we modify the algorithm so that the cost of constructing the, now two, generating vectors component-by-component is only operations. This reduction in cost allows, in practice, construction of rules with millions of points. The rules constructed again achieve a worst-case strong tractability error bound, with a rate of convergence for 0$">.

     

Some new kinds of pseudoprimes

We define some new kinds of pseudoprimes to several bases, which generalize strong pseudoprimes. We call them Sylow -pseudoprimes and elementary Abelian -pseudoprimes. It turns out that every which is a strong pseudoprime to bases 2, 3 and 5, is not a Sylow -pseudoprime to two of these bases for an appropriate prime

We also give examples of strong pseudoprimes to many bases which are not Sylow -pseudoprimes to two bases only, where or

     

New quadratic polynomials with high densities of prime values

Hardy and Littlewood's Conjecture F implies that the asymptotic density of prime values of the polynomials , is related to the discriminant of via a quantity The larger is, the higher the asymptotic density of prime values for any quadratic polynomial of discriminant . A technique of Bach allows one to estimate accurately for any , given the class number of the imaginary quadratic order with discriminant , and for any 0$"> given the class number and regulator of the real quadratic order with discriminant . The Manitoba Scalable Sieve Unit (MSSU) has shown us how to rapidly generate many discriminants for which is potentially large, and new methods for evaluating class numbers and regulators of quadratic orders allow us to compute accurate estimates of efficiently, even for values of with as many as decimal digits. Using these methods, we were able to find a number of discriminants for which, under the assumption of the Extended Riemann Hypothesis, is larger than any previously known examples.

     

Upper bounds for the prime divisors of Wendt's determinant

Let be an even integer, . The resultant of the polynomials and is known as Wendt's determinant of order . We prove that among the prime divisors of only those which divide or can be larger than , where and is the th Lucas number, except when and . Using this estimate we derive criteria for the nonsolvability of Fermat's congruence.

     

A priori error estimates for Galerkin approximations to porous medium and fast diffusion equations

Galerkin approximations to solutions of a Cauchy-Dirichlet problem governed by the generalized porous medium equation

on bounded convex domains are considered. The range of the parameter includes the fast diffusion case . Using an Euler finite difference approximation in time, the semi-discrete solution is shown to converge to the exact solution in norm with an error controlled by for and for . For the fully discrete problem, a global convergence rate of in norm is shown for the range . For , a rate of is shown in norm.

     

Analysis of PSLQ, an integer relation finding algorithm

Let be either the real, complex, or quaternion number system and let be the corresponding integers. Let be a vector in . The vector has an integer relation if there exists a vector , , such that . In this paper we define the parameterized integer relation construction algorithm PSLQ, where the parameter can be freely chosen in a certain interval. Beginning with an arbitrary vector , iterations of PSLQ will produce lower bounds on the norm of any possible relation for . Thus PSLQ can be used to prove that there are no relations for of norm less than a given size. Let be the smallest norm of any relation for . For the real and complex case and each fixed parameter in a certain interval, we prove that PSLQ constructs a relation in less than iterations.

     

Maximum-norm estimates for resolvents of elliptic finite element operators

Let be a convex domain with smooth boundary in . It has been shown recently that the semigroup generated by the discrete Laplacian for quasi-uniform families of piecewise linear finite element spaces on is analytic with respect to the maximum-norm, uniformly in the mesh-width. This implies a resolvent estimate of standard form in the maximum-norm outside some sector in the right halfplane, and conversely. Here we show directly that such a resolvent estimate holds outside any sector around the positive real axis, with arbitrarily small angle. This is useful in the study of fully discrete approximations based on -stable rational functions, with small.

     

On the step-by-step construction of quasi-Monte Carlo integration rules that achieve strong tractability error bounds in weighted Sobolev spaces

We develop and justify an algorithm for the construction of quasi-Monte Carlo (QMC) rules for integration in weighted Sobolev spaces; the rules so constructed are shifted rank-1 lattice rules. The parameters characterising the shifted lattice rule are found ``component-by-component': the ()-th component of the generator vector and the shift are obtained by successive -dimensional searches, with the previous components kept unchanged. The rules constructed in this way are shown to achieve a strong tractability error bound in weighted Sobolev spaces. A search for -point rules with prime and all dimensions 1 to requires a total cost of operations. This may be reduced to operations at the expense of storage. Numerical values of parameters and worst-case errors are given for dimensions up to 40 and up to a few thousand. The worst-case errors for these rules are found to be much smaller than the theoretical bounds.