Numerics of classical elliptic functions, elliptic integrals and modular forms
|
|
- Job Johns
- 5 years ago
- Views:
Transcription
1 1 / 49 Numerics of classical elliptic functions, elliptic integrals and modular forms Fredrik Johansson LFANT, Inria Bordeaux & Institut de Mathématiques de Bordeaux KMPB Conference: Elliptic Integrals, Elliptic Functions and Modular Forms in Quantum Field Theory DESY, Zeuthen, Germany 24 October 2017
2 2 / 49 Introduction Elliptic functions F(z + ω 1 m + ω 2 n) = F(z), m, n Z Can assume ω 1 = 1 and ω 2 = τ H = {τ C : Im(τ) > 0} Elliptic integrals R(x, P(x))dx; inverses of elliptic functions Modular forms/functions on H F( aτ+b cτ+d ) = (cτ + d)k F(τ) for ( a b c d ) PSL2 (Z) Related to elliptic functions with fixed z and varying lattice parameter ω 2 /ω 1 = τ H Jacobi theta functions (quasi-elliptic functions) Used to construct elliptic and modular functions
3 3 / 49 Numerical evaluation Lots of existing literature, software (Pari/GP, Sage, Maple, Mathematica, Matlab, Maxima, GSL, NAG,... ). This talk will mostly review standard techniques (and many techniques will be omitted). My goal: general purpose methods with Rigorous error bounds Arbitrary precision Complex variables Implementations in the C library Arb (
4 4 / 49 Why arbitrary precision? Applications: Mitigating roundoff error for lengthy calculations Surviving cancellation between exponentially large terms High order numerical differentiation, extrapolation Computing discrete data (integer coefficients) Integer relation searches (LLL/PSLQ) Heuristic equality testing Also: Can increase precision if error bounds are too pessimistic Most interesting range: digits. (Millions, billions...?)
5 Ball/interval arithmetic 5 / 49 A real number in Arb is represented by a rigorous enclosure as a high-precision midpoint and a low-precision radius: [ ± ] Complex numbers: [m 1 ± r 1 ] + [m 2 ± r 2 ]i. Key points: Error bounds are propagated automatically As cheap as arbitrary-precision floating-point To compute f (x) = k=0 N 1 k=0 rigorously, only need analysis to bound k=n Dependencies between variables may lead to inflated enclosures. Useful technique is to compute f ([m ± r]) as [f (m) ± s] where s = r sup x m r f (x).
6 6 / 49 Reliable numerical evaluation Example: sin(π ) IEEE 754 double precision result: e-16 Adaptive numerical evaluation with Arb: 164 bits: [± ] 128 bits: [ ± ] 192 bits: [ ± ] Can be used to implement reliable floating-point functions, even if you don t use interval arithmetic externally: Float input Arb function Increase precision Output midpoint yes Accurate enough? no
7 Elliptic and modular functions in Arb 7 / 49 PSL 2 (Z) transformations and argument reduction Jacobi theta functions θ 1 (z, τ),..., θ 4 (z, τ) Arbitrary z-derivatives of Jacobi theta functions Weierstrass elliptic functions (n) (z, τ), 1 (z, τ), ζ(z, τ), σ(z, τ) Modular forms and functions: j(τ), η(τ), (τ), λ(τ), G 2k (τ) Legendre complete elliptic integrals K (m), E(m), Π(n, m) Incomplete elliptic integrals F(φ, m), E(φ, m), Π(n, φ, m) Carlson incomplete elliptic integrals R F, R J, R C, R D, R G Possible future projects: The suite of Jacobi elliptic functions and integrals Asymptotic complexity improvements
8 8 / 49 An application: Hilbert class polynomials For D < 0 congruent to 0 or 1 mod 4, H D (x) = ( ( b + )) D x j Z[x] 2a (a,b,c) where (a, b, c) is taken over all the primitive reduced binary quadratic forms ax 2 + bxy + cy 2 with b 2 4ac = D. Example: H 31 = x x x Algorithms: modular, complex analytic D Degree Bits Pari/GP classpoly CM Arb s 0.8 s 0.4 s 0.14 s s 8 s 29 s 17 s s 82 s 436 s 274 s
9 Some visualizations 9 / 49 The Weierstrass zeta-function ζ( i, τ) as the lattice parameter τ varies over [ 0.25, 0.25] + [0, 0.15]i.
10 Some visualizations 10 / 49 The Weierstrass elliptic functions ζ(z, i) (left) and σ(z, i) (right) as z varies over [ π, π], [ π, π]i.
11 Some visualizations 11 / 49 The function j(τ) on the complex interval [ 2, 2] + [0, 1]i. The function η(τ) on the complex interval [0, 24] + [0, 1]i.
12 Some visualizations 12 / 49 Plot of j(τ) on [ 13, ] + [0, ]i. Plot of η(τ) on [ 2, ] + [0, ]i.
13 Approaches to computing special functions 13 / 49 Numerical integration (integral representations, ODEs) Functional equations (argument reduction) Series expansions Root-finding methods (for inverse functions) Precomputed approximants (not applicable here)
14 14 / 49 Brute force: numerical integration For analytic integrands, there are good algorithms that easily permit achieving 100s or 1000s of digits of accuracy: Gaussian quadrature Clenshaw-Curtis method (Chebyshev series) Trapezoidal rule (for periodic functions) Double exponential (tanh-sinh) method Taylor series methods (also for ODEs) Pros: Simple, general, flexible approach Can deform path of integration as needed Cons: Usually slower than dedicated methods Possible convergence problems (oscillation, singularities) Error analysis may be complicated for improper integrals
15 Poisson and the trapezoidal rule (historical remark) 15 / 49 In 1827, Poisson considered the example of the perimeter of an ellipse with axis lengths 1/π and 0.6/π: I = 1 2π 2π sin 2 (θ)dθ = 2 E(0.36) = π Poisson used the trapezoidal approximation I I N = 4 N N /4 k= sin 2 (2πk/N ). With N = 16 (four points!), he computed I and proved that the error is < In fact I N I = O(3 N ). See Trefethen & Weideman, The exponentially convergent trapezoidal rule, 2014.
16 A model problem: computing exp(x) 16 / 49 Standard two-step numerical recipe for special functions: (not all functions fit this pattern, but surprisingly many do!) 1. Argument reduction exp(x) = exp(x n log(2)) 2 n 2. Series expansion exp(x) = [ exp(x/2 R ) ] 2 R exp(x) = 1 + x + x2 2! + x3 3! +... Step (1) ensures rapid convergence and good numerical stability in step (2).
17 17 / 49 Reducing complexity for p-bit precision Principles: Balance argument reduction and series order optimally Exploit special (e.g. hypergeometric) structure of series How to compute exp(x) for x 1 with an error of ? Only reduction: apply x x/2 reduction 1000 times Only series evaluation: use 170 terms (170! > ) Better: apply 1000 = 32 reductions and use 32 terms This trick reduces the arithmetic complexity from p to p 0.5 (time complexity from p 2+ε to p 1.5+ε ). With a more complex scheme, the arithmetic complexity can be reduced to O(log 2 p) (time complexity p 1+ε ).
18 Evaluating polynomials using rectangular splitting (Paterson and Stockmeyer 1973; Smith 1989) N i=0 xi in O(N ) cheap steps + O(N 1/2 ) expensive steps ( + x + x 2 + x 3 ) + ( + x + x 2 + x 3 ) x 4 + ( + x + x 2 + x 3 ) x 8 + ( + x + x 2 + x 3 ) x 12 This does not genuinely reduce the asymptotic complexity, but can be a huge improvement (100 times faster) in practice. 18 / 49
19 19 / 49 Elliptic functions Elliptic integrals Argument reduction Move to standard domain (periodicity, modular transformations) Move parameters close together (various formulas) Series expansions Theta function q-series Multivariate hypergeometric series (Appell, Lauricella... ) Special cases Modular forms & functions, theta constants Complete elliptic integrals, ordinary hypergeometric series (Gauss 2 F 1 )
20 Modular forms and functions 20 / 49 A modular form of weight k is a holomorphic function on H = {τ : τ C, Im(τ) > 0} satisfying ( ) aτ + b F = (cτ + d) k F(τ) cτ + d for any integers a, b, c, d with ad bc = 1. A modular function is meromorphic and has weight k = 0. Since F(τ) = F(τ + 1), the function has a Fourier series (or Laurent series/q-expansion) F(τ) = n= m c n e 2iπnτ = n= m c n q n, q = e 2πiτ, q < 1
21 Some useful functions and their q-expansions 21 / 49 Dedekind ( eta ) function η aτ+b cτ+d = ε(a, b, c, d) cτ + dη(τ) η(τ) = e πiτ/12 n= ( 1)n q (3n2 n)/2 The j-invariant ( ) j aτ+b cτ+d = j(τ) j(τ) = 1 q q q2 + j(τ) = 32(θ θ8 3 + θ8 4 )3 /(θ 2 θ 3 θ 4 ) 8 Theta constants (q = e πiτ ) (θ 2, θ 3, θ 4 ) = n= (q (n+1/2)2, q n2, ( 1) n q n2) Due to sparseness, we only need N = O( p) terms for p-bit accuracy (so the evaluation takes p 1.5+ε time).
22 22 / 49 Argument reduction for modular forms PSL 2 (Z) is generated by ( ) and ( By repeated use of τ τ + 1 or τ 1/τ, we can move τ to the fundamental domain { τ H : z 1, Re(z) 1 2}. In the fundamental domain, q exp( π 3) = , which gives rapid convergence of the q-expansion. )
23 Practical considerations 23 / 49 Instead of applying F(τ + 1) = F(τ) or F( 1/τ) = τ k F(τ) step by step, build transformation matrix g = ( ) a b c d and apply to F in one step. This improves numerical stability g can usually be computed cheaply using machine floats If computing F via theta constants, apply transformation for F instead of the individual theta constants.
24 Fast computation of eta and theta function q-series 24 / 49 Consider N n=0 qn2. More generally, q P(n), P Z[x] of degree 2. Naively: 2N multiplications. Enge, Hart & J, Short addition sequences for theta functions, 2016: Optimized addition sequence for P(0), P(1),... (2 speedup) Rectangular splitting: choose splitting parameter m so that P has few distinct residues mod m (logarithmic speedup, in practice another 2 speedup) Schost & Nogneng, On the evaluation of some sparse polynomials, 2017: N 1/2+ε method (p 1.25+ε time complexity) using FFT Faster for p > in practice
25 Jacobi theta functions 25 / 49 Series expansion: θ 3 (z, τ) = q n2 w 2n, n= and similarly for θ 1, θ 2, θ 4. q = e πiτ, w = e πiz The terms eventually decay rapidly (there can be an initial hump if w is large). Error bound via geometric series. For z-derivatives, we compute the object θ(z + x, τ) C[[x]] (as a vector of coefficients) in one step. θ(z+x, τ) = θ(z, τ)+θ (z, τ)x+...+ θ(r 1) (z, τ) x r 1 +O(x r ) C[[x]] (r 1)!
26 Argument reduction for Jacobi theta functions 26 / 49 Two reductions are necessary: Move τ to τ in the fundamental domain (this operation transforms z z, introduces some prefactors, and permutes the theta functions) Reduce z modulo τ using quasiperiodicity General formulas for the transformation τ τ = aτ+b cτ+d are given in (Rademacher, 1973): z = θ n (z, τ) = exp(πir/4) A B θ S (z, τ ) z i cτ + d, A = cτ + d, ( B = exp z 2 ) πic cτ + d R, S are integers depending on n and (a, b, c, d). The argument reduction also applies to θ(z + x, τ) C[[x]].
27 Elliptic functions 27 / 49 The Weierstrass elliptic function (z, τ) = (z + 1, τ) = (z + τ, τ) (z, τ) = 1 z 2 + [ ] 1 (z + m + nτ) 2 1 (m + nτ) 2 n 2 +m 2 0 is computed via Jacobi theta functions as (z, τ) = π 2 θ2(0, 2 τ)θ3(0, 2 τ) θ2 4 (z, τ) π2 [ θ1 2 θ 4 (z, τ) 3 3 (0, τ) + θ3(0, 4 τ) ] Similarly σ(z, τ), ζ(z, τ) and (k) (z, τ) using z-derivatives of theta functions. With argument reduction for both z and τ already implemented for theta functions, reduction for is unnecessary (but can improve numerical stability).
28 Some timings For d decimal digits (z = 5 + 7i, τ = 7 + i/ 11): Function d = 10 d = 10 2 d = 10 3 d = 10 4 d = 10 5 exp(z) log(z) η(τ) j(τ) (θ i (0, τ)) 4 i= (θ i (z, τ)) 4 i= (z, τ) (, ) ζ(z, τ) σ(z, τ) / 49
29 Elliptic integrals 29 / 49 Any elliptic integral R(x, P(x))dx can be written in terms of a small basis set. The Legendre forms are used by tradition. Complete elliptic integrals: K (m) = π/2 dt 0 = 1 1 m sin 2 t 0 E(m) = π/2 0 Π(n, m) = π/2 0 1 m sin 2 t dt = 1 0 (1 n sin 2 t) dt ( 1 t 2 )( 1 mt 2 ) dt = 1 1 m sin 2 t 0 Incomplete integrals: F(φ, m) = φ dt 0 = sin φ 1 m sin 2 t 0 E(φ, m) = φ 0 Π(n, φ, m) = φ 0 1 m sin 2 t dt = sin φ 0 (1 n sin 2 t) 1 mt 2 1 t 2 dt (1 nt 2 ) dt ( 1 t 2 )( 1 mt 2 ) dt = sin φ 1 m sin 2 t 0 1 mt 2 1 t 2 dt (1 nt 2 ) dt 1 t 2 1 mt 2 dt 1 t 2 1 mt 2
30 30 / 49 Complete elliptic integrals and 2 F 1 The Gauss hypergeometric function is defined for z < 1 by 2F 1 (a, b, c, z) = k=0 (a) k (b) k (c) k z k k!, (x) k = x(x + 1) (x + k 1) and elsewhere by analytic continuation. The 2 F 1 function can be computed efficiently for any z C. K (m) = 1 2 π 2F 1 ( 1 2, 1 2, 1, m) E(m) = 1 2 π 2F 1 ( 1 2, 1 2, 1, m) This works, but it s not the best way!
31 Complete elliptic integrals and the AGM 31 / 49 The AGM of x, y is the common limit of the sequences a n+1 = a n + b n, b n+1 = a n b n 2 with a 0 = x, b 0 = y. As a functional equation: ( x + y M(x, y) = M 2, ) xy Each step doubles the number of digits in M(x, y) x y convergence in O(log p) operations (p 1+ε time complexity). K (m) = π 2M(1, 1 m), E(m) = (1 m)(2mk (m) + K (m))
32 Numerical aspects of the AGM 32 / 49 Argument reduction vs series expansion: O(1) terms only. Slightly better than reducing all the way to a n b n < 2 p : π 4K (z 2 ) = 1 2 z2 8 5z z z O(z10 ) Complex variables: simplify to M(z) = M(1, z) using M(x, y) = xm(1, y/x). Some case distinctions for correct square root branches in AGM iteration. Derivatives: can use finite (central) difference for M (z) (better method possible using elliptic integrals), higher derivatives using recurrence relations.
33 Incomplete elliptic integrals Incomplete elliptic integrals are multivariate hypergeometric functions. In terms of the Appell F 1 function F 1 (a, b 1, b 2 ; c; x, y) = where x, y < 1, we have F(z, m) = Problems: z 0 m,n=0 (a) m+n (b 1 ) m (b 2 ) n (c) m+n m! n! x m y n dt 1 m sin 2 t = sin(z) F 1( 1 2, 1 2, 1 2, 3 2 ; sin2 z, m sin 2 z) How to reduce arguments so that x, y 1? How to perform analytic continuation and obtain consistent branch cuts for complex variables? 33 / 49
34 Branch cuts of Legendre incomplete elliptic integrals 34 / 49
35 Branch cuts of F(z, m) with respect to z / 49
36 Branch cuts of F(z, m) with respect to z (continued) 36 / 49
37 Branch cuts of F(z, m) with respect to z (continued) 37 / 49
38 Branch cuts of F(z, m) with respect to z (continued) 38 / 49
39 Branch cuts of F(z, m) with respect to z (continued) 39 / 49
40 Branch cuts of F(z, m) with respect to m 40 / 49 Conclusion: the Legendre forms are not nice as building blocks.
41 Carlson s symmetric forms 41 / 49 In the 1960s, Bille C. Carlson suggested an alternative basis set for incomplete elliptic integrals: R F (x, y, z) = 1 2 R J (x, y, z, p) = 3 2 Advantages: 0 0 dt (t + x)(t + y)(t + z) dt (t + p) (t + x)(t + y)(t + z) R C (x, y) = R F (x, y, y), R D (x, y, z) = R J (x, y, z, z) Symmetry unifies and simplifies transformation laws Symmetry greatly simplifies series expansions The functions have nice complex branch structure Simple universal algorithm for computation
42 Evaluation of Legendre forms 42 / 49 For π 2 Re(z) π 2 : F(z, m) = sin(z) R F (cos 2 (z), 1 m sin 2 (z), 1) Elsewhere, use quasiperiodic extension: F(z + kπ, m) = 2kK (m) + F(z, m), k Z Similarly for E(z, m) and Π(n, z, m). Slight complication to handle (complex) intervals straddling the lines Re(z) = (n )π. Useful for implementations: variants with z πz.
43 Symmetric argument reduction 43 / 49 We have the functional equation ( x + λ R F (x, y, z) = R F 4, y + λ 4, z + λ ) 4 where λ = x y + y z + z x. Each application reduces the distance between x, y, z by a factor 1/4. Algorithm: apply reduction until the distance is ε, then use an order-n series expansion with error term O(ε N ). For p-bit accuracy, need p/(2n ) argument reduction steps. (A similar functional equation exists for R J (x, y, z, p).)
44 Series expansion when arguments are close ( R F (x, y, z) = R 1 1/2 2, 1 2, 1 2, x, y, z) ( R J (x, y, z, p) = R 1 3/2 2, 1 2, 1 2, 1 2, 1 2, x, y, z, p, p) Carlson s R is a multivariate hypergeometric series: R a (b; z) = = M=0 M=0 (a) M ( n j=1 b T M (b 1,..., b n ; 1 z 1,..., 1 z n ) j) M (a) M ( n j=1 b T M j) M z a n ( b 1,..., b n 1 ; 1 z 1 z n,..., 1 z n 1 z n ), T M (b 1,..., b n, w 1,..., w n ) = m m n=m n j=1 (b j ) mj (m j )! wm j j Note that T M Const p(m) max( w 1,..., w n ) M, so we can easily bound the tail by a geometric series. 44 / 49
45 45 / 49 A clever idea by Carlson: symmetric polynomials Using elementary symmetric polynomials E s (w 1,..., w n ), T M ( 1 2, w) = m 1 +2m nm n=m ( 1) M+ j m j ( 1 2) j m j n j=1 E m j j (w) (m j )! We can expand R around the mean of the arguments, taking w j = 1 z j /A where A = 1 n n j=1 z j. Then E 1 = 0, and most of the terms disappear! Carlson suggested expanding to M < N = 8: A 1/2 R F (x, y, z) = 1 E E E E 2E E E E 2 2 E O(ε8 ) Need p/16 argument reduction steps for p-bit accuracy.
46 46 / 49 Rectangular splitting for the R series The exponents of E m 2 2 E m 3 3 appearing in the series for R F are the lattice points m 2, m 3 Z 0 with 2m 2 + 3m 3 < N. m 3 (terms appearing for N = 10) m 2 Compute powers of E 2, use Horner s rule with respect to E 3. Clear denominators so that all coefficients are small integers. O(N 2 ) cheap steps + O(N ) expensive steps For R J, compute powers of E 2, E 3, use Horner for E 4, E 5.
47 47 / 49 Balancing series evaluation and argument reduction Consider R F : p = wanted precision in bits O(ε N ) = error due to truncating the series expansion O(N 2 ) = number of terms in series O(p/N ) = number of argument reduction steps for ε N = 2 p Overall cost O(N 2 + p/n ) is minimized by N p 0.333, giving p arithmetic complexity (p time complexity). Empirically, N 2p 0.4 is optimal (due to rectangular splitting). Speedup over N = 8 at d digits precision: d = 10 d = 10 2 d = 10 3 d = 10 4 d =
48 48 / 49 Some timings We include K (m) (computed by AGM), F(z, m) (computed by R F ) and the inverse Weierstrass elliptic function: 1 (z, τ) = 1 2 z dt (t e1 )(t e 2 )(t e 3 ) = R F (z e 1, z e 2, z e 3 ) Function d = 10 d = 10 2 d = 10 3 d = 10 4 d = 10 5 exp(z) log(z) η(τ) K (m) F(z, m) (z, τ) (z, τ)
49 49 / 49 Quadratic transformations It is possible to construct AGM-like methods (converging in O(log p) steps) for general elliptic integrals and functions. Problems: The overhead may be slightly higher at low precision Correct treatment of complex variables is not obvious Unfortunately, I have not had time to study this topic. However, see the following papers: The elliptic logarithm ( 1 ): John E. Cremona and Thotsaphon Thongjunthug, The complex AGM, periods of elliptic curves over and complex elliptic logarithms, Elliptic and theta functions: Hugo Labrande, Computing Jacobi s θ in quasi-linear time, 2015.
Addition sequences and numerical evaluation of modular forms
Addition sequences and numerical evaluation of modular forms Fredrik Johansson (INRIA Bordeaux) Joint work with Andreas Enge (INRIA Bordeaux) William Hart (TU Kaiserslautern) DK Statusseminar in Strobl,
More informationNumerical Evaluation of Elliptic Functions, Elliptic Integrals and Modular Forms
Numerical Evaluation of Elliptic Functions, Elliptic Integrals and Modular Forms Fredrik Johansson To cite this version: Fredrik Johansson. Numerical Evaluation of Elliptic Functions, Elliptic Integrals
More informationEfficient implementation of the Hardy-Ramanujan-Rademacher formula
Efficient implementation of the Hardy-Ramanujan-Rademacher formula or: Partitions in the quintillions Fredrik Johansson RISC-Linz July 10, 2013 2013 SIAM Annual Meeting San Diego, CA Supported by Austrian
More information1 Introduction. or equivalently f(z) =
Introduction In this unit on elliptic functions, we ll see how two very natural lines of questions interact. The first, as we have met several times in Berndt s book, involves elliptic integrals. In particular,
More informationPartitions in the quintillions or Billions of congruences
Partitions in the quintillions or Billions of congruences Fredrik Johansson November 2011 The partition function p(n) counts the number of ways n can be written as the sum of positive integers without
More informationFast reversion of formal power series
Fast reversion of formal power series Fredrik Johansson LFANT, INRIA Bordeaux RAIM, 2016-06-29, Banyuls-sur-mer 1 / 30 Reversion of power series F = exp(x) 1 = x + x 2 2! + x 3 3! + x 4 G = log(1 + x)
More informationFor COURSE PACK and other PERMISSIONS, refer to entry on previous page. For more information, send to
COPYRIGHT NOTICE: Elias M. Stein and Rami Shakarchi: Complex Analysis is published by Princeton University Press and copyrighted, 2003, by Princeton University Press. All rights reserved. No part of this
More informationNumerical integration in arbitrary-precision ball arithmetic
/ 24 Numerical integration in arbitrary-precision ball arithmetic Fredrik Johansson (LFANT, Bordeaux) Journées FastRelax, Inria Sophia Antipolis 7 June 208 Numerical integration in Arb 2 / 24 New code
More informationFast reversion of power series
Fast reversion of power series Fredrik Johansson November 2011 Overview Fast power series arithmetic Fast composition and reversion (Brent and Kung, 1978) A new algorithm for reversion Implementation results
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationClass polynomials for abelian surfaces
Class polynomials for abelian surfaces Andreas Enge LFANT project-team INRIA Bordeaux Sud-Ouest andreas.enge@inria.fr http://www.math.u-bordeaux.fr/~aenge LFANT seminar 27 January 2015 (joint work with
More informationELLIPTIC FUNCTIONS AND THETA FUNCTIONS
ELLIPTIC FUNCTIONS AND THETA FUNCTIONS LECTURE NOTES FOR NOV.24, 26 Historically, elliptic functions were first discovered by Niels Henrik Abel as inverse functions of elliptic integrals, and their theory
More informationNotations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation. Specialized values
KleinInvariant Notations Traditional name Klein invariant modular function Traditional notation z Mathematica StandardForm notation KleinInvariant z Primary definition 09.50.0.0001.01 ϑ 0, Π z 8 ϑ 3 0,
More informationCongruence Subgroups
Congruence Subgroups Undergraduate Mathematics Society, Columbia University S. M.-C. 24 June 2015 Contents 1 First Properties 1 2 The Modular Group and Elliptic Curves 3 3 Modular Forms for Congruence
More informationRamanujan s first letter to Hardy: 5 + = 1 + e 2π 1 + e 4π 1 +
Ramanujan s first letter to Hardy: e 2π/ + + 1 = 1 + e 2π 2 2 1 + e 4π 1 + e π/ 1 e π = 2 2 1 + 1 + e 2π 1 + Hardy: [These formulas ] defeated me completely. I had never seen anything in the least like
More informationHans Wenzl. 4f(x), 4x 3 + 4ax bx + 4c
MATH 104C NUMBER THEORY: NOTES Hans Wenzl 1. DUPLICATION FORMULA AND POINTS OF ORDER THREE We recall a number of useful formulas. If P i = (x i, y i ) are the points of intersection of a line with the
More informationExact Formulas for Invariants of Hilbert Schemes
Exact Formulas for Invariants of Hilbert Schemes N. Gillman, X. Gonzalez, M. Schoenbauer, advised by John Duncan, Larry Rolen, and Ken Ono September 2, 2018 1 / 36 Paper See the related paper at https://arxiv.org/pdf/1808.04431.pdf
More informationNumerical Integration with the double exponential method
Numerical Integration with the double exponential method Pascal Molin Université Paris 7 December 7 2016 Numerical integration Number-theory context, want fast quadrature high precision (100 1000 digits)
More informationComputation of the error functions erf and erfc in arbitrary precision with correct rounding
Computation of the error functions erf and erfc in arbitrary precision with correct rounding Sylvain Chevillard Arenaire, LIP, ENS-Lyon, France Sylvain.Chevillard@ens-lyon.fr Nathalie Revol INRIA, Arenaire,
More informationChapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn
Chapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide n c: Solve the problem
More informationINDEFINITE THETA FUNCTIONS OF TYPE (n, 1) I: DEFINITIONS AND EXAMPLES
INDEFINITE THETA FUNCTIONS OF TYPE (n, ) I: DEFINITIONS AND EXAMPLES LARRY ROLEN. Classical theta functions Theta functions are classical examples of modular forms which play many roles in number theory
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationOn the zeros of certain modular forms
On the zeros of certain modular forms Masanobu Kaneko Dedicated to Professor Yasutaka Ihara on the occasion of his 60th birthday. The aim of this short note is to list several families of modular forms
More informationarxiv:math/ v1 [math.ca] 6 Sep 1994
NUMERICAL COMPUTATION OF REAL OR COMPLEX arxiv:math/909227v1 [math.ca] 6 Sep 199 ELLIPTIC INTEGRALS B. C. Carlson Ames Laboratory and Department of Mathematics, Iowa State University, Ames, Iowa 50011-020,
More informationProjects on elliptic curves and modular forms
Projects on elliptic curves and modular forms Math 480, Spring 2010 In the following are 11 projects for this course. Some of the projects are rather ambitious and may very well be the topic of a master
More informationELLIPSES AND ELLIPTIC CURVES. M. Ram Murty Queen s University
ELLIPSES AND ELLIPTIC CURVES M. Ram Murty Queen s University Planetary orbits are elliptical What is an ellipse? x 2 a 2 + y2 b 2 = 1 An ellipse has two foci From: gomath.com/geometry/ellipse.php Metric
More informationComplex Variables. Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit.
Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit. 1. The TI-89 calculator says, reasonably enough, that x 1) 1/3 1 ) 3 = 8. lim
More informationThe arithmetic geometric mean (Agm)
The arithmetic geometric mean () Pictures by David Lehavi This is part of exercise 5 question 4 in myc> calculus class Fall 200: a = c a n+ =5a n 2 lim n an = a = a an + bn a n+ = 2 b = b b n+ = a nb n
More informationGREEN FUNCTION, PAINLEVÉ VI EQUATION, AND EISENSTEIN SERIES OF WEIGHT ONE
GREEN FUNCTION, PAINLEVÉ VI EQUATION, AND EISENSTEIN SERIES OF WEIGHT ONE ZHIJIE CHEN, TING-JUNG KUO, CHANG-SHOU LIN, AND CHIN-LUNG WANG ABSTRACT. We study the problem: How many singular points of a solution
More informationOn Newton-type methods with cubic convergence
Journal of Computational and Applied Mathematics 176 (2005) 425 432 www.elsevier.com/locate/cam On Newton-type methods with cubic convergence H.H.H. Homeier a,b, a Science + Computing Ag, IT Services Muenchen,
More informationComputing with precision
1 / 44 Computing with precision Fredrik Johansson Inria Bordeaux X, Mountain View, CA January 24, 2019 Computing with real numbers 2 / 44 How can we represent 3.14159265358979323846264338327950288419716939937510582097
More informationComputing logarithms and other special functions
Computing logarithms and other special functions Mike Giles University of Oxford Mathematical Institute Napier 400 NAIS Symposium April 2, 2014 Mike Giles (Oxford) Computing special functions April 2,
More informationNOTES ON RIEMANN S ZETA FUNCTION. Γ(z) = t z 1 e t dt
NOTES ON RIEMANN S ZETA FUNCTION DRAGAN MILIČIĆ. Gamma function.. Definition of the Gamma function. The integral Γz = t z e t dt is well-defined and defines a holomorphic function in the right half-plane
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationModular Forms, Elliptic Curves, and Modular Curves
1 Modular Forms, Elliptic Curves, and Modular Curves This chapter introduces three central objects of the book. Modular forms are functions on the complex upper half plane. A matrix group called the modular
More informationEllipticTheta2. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation
EllipticTheta Notations Traditional name Jacobi theta function ϑ Traditional notation ϑ z, Mathematica StandardForm notation EllipticTheta, z, Primary definition 09.0.0.0001.01 ϑ z, cos k 1 z ; 1 Specific
More information20 The modular equation
18.783 Elliptic Curves Spring 2015 Lecture #20 04/23/2015 20 The modular equation In the previous lecture we defined modular curves as quotients of the extended upper half plane under the action of a congruence
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationClass invariants for quartic CM-fields
Number Theory Seminar Oxford 2 June 2011 Elliptic curves An elliptic curve E/k (char(k) 2) is a smooth projective curve y 2 = x 3 + ax 2 + bx + c. Q P E is a commutative algebraic group P Q Endomorphisms
More informationConstructing Class invariants
Constructing Class invariants Aristides Kontogeorgis Department of Mathematics University of Athens. Workshop Thales 1-3 July 2015 :Algebraic modeling of topological and computational structures and applications,
More informationMATH5685 Assignment 3
MATH5685 Assignment 3 Due: Wednesday 3 October 1. The open unit disk is denoted D. Q1. Suppose that a n for all n. Show that (1 + a n) converges if and only if a n converges. [Hint: prove that ( N (1 +
More informationOld and new algorithms for computing Bernoulli numbers
Old and new algorithms for computing Bernoulli numbers University of New South Wales 25th September 2012, University of Ballarat Bernoulli numbers Rational numbers B 0, B 1,... defined by: x e x 1 = n
More informationOn the evaluation of modular polynomials
On the evaluation of modular polynomials Andrew V. Sutherland Massachusetts Institute of Technology ANTS X July 10, 2012 http://math.mit.edu:/ drew 1 / 16 Introduction Let l be a prime and let F q be a
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More information2 Write down the range of values of α (real) or β (complex) for which the following integrals converge. (i) e z2 dz where {γ : z = se iα, < s < }
Mathematical Tripos Part II Michaelmas term 2007 Further Complex Methods, Examples sheet Dr S.T.C. Siklos Comments and corrections: e-mail to stcs@cam. Sheet with commentary available for supervisors.
More informationRamanujan s theories of elliptic functions to alternative bases, and beyond.
Ramanuan s theories of elliptic functions to alternative bases, and beyond. Shaun Cooper Massey University, Auckland Askey 80 Conference. December 6, 03. Outline Sporadic sequences Background: classical
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationA Special Module in the Field of Modular Functions
A Special Module in the Field of Modular Functions Research Institute for Symbolic Computation (RISC), Johannes Kepler University, A-4040 Linz, Austria Dedicated to Peter Paule at the Occasion of his 60th
More information20 The modular equation
18.783 Elliptic Curves Lecture #20 Spring 2017 04/26/2017 20 The modular equation In the previous lecture we defined modular curves as quotients of the extended upper half plane under the action of a congruence
More informationIntroduction to modular forms Perspectives in Mathematical Science IV (Part II) Nagoya University (Fall 2018)
Introduction to modular forms Perspectives in Mathematical Science IV (Part II) Nagoya University (Fall 208) Henrik Bachmann (Math. Building Room 457, henrik.bachmann@math.nagoya-u.ac.jp) Lecture notes
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationMicrostates of AdS black holes and supersymmetric localization
Microstates of AdS black holes and supersymmetric localization Seyed Morteza Hosseini Università di Milano-Bicocca IPM, Tehran, May 8-11, 2017 Recent Trends in String Theory and Related Topics in collaboration
More informationCPSC 518 Introduction to Computer Algebra Asymptotically Fast Integer Multiplication
CPSC 518 Introduction to Computer Algebra Asymptotically Fast Integer Multiplication 1 Introduction We have now seen that the Fast Fourier Transform can be applied to perform polynomial multiplication
More informationTHE RAMANUJAN-SERRE DIFFERENTIAL OPERATORS AND CERTAIN ELLIPTIC CURVES. (q = e 2πiτ, τ H : the upper-half plane) ( d 5) q n
THE RAMANUJAN-SERRE DIFFERENTIAL OPERATORS AND CERTAIN ELLIPTIC CURVES MASANOBU KANEKO AND YUICHI SAKAI Abstract. For several congruence subgroups of low levels and their conjugates, we derive differential
More informationDOCTORAL THESIS Extended Abstract
Lodz University of Technology Faculty of Electrical, Electronic, Computer and Control Engineering DOCTORAL THESIS Extended Abstract Dariusz Brzeziński MSc Eng. Problems of Numerical Calculation of Derivatives
More informationClosed Form Solutions
Closed Form Solutions Mark van Hoeij 1 Florida State University Slides of this talk: www.math.fsu.edu/ hoeij/issac2017.pdf 1 Supported by NSF 1618657 1 / 29 What is a closed form solution? Example: Solve
More informationRANKIN-COHEN BRACKETS AND SERRE DERIVATIVES AS POINCARÉ SERIES. φ k M = 1 2
RANKIN-COHEN BRACKETS AND SERRE DERIVATIVES AS POINCARÉ SERIES BRANDON WILLIAMS Abstract. We give expressions for the Serre derivatives of Eisenstein and Poincaré series as well as their Rankin-Cohen brackets
More informationCOMPUTING MODULAR POLYNOMIALS WITH THETA FUNCTIONS
ERASMUS MUNDUS MASTER ALGANT UNIVERSITÀ DEGLI STUDI DI PADOVA UNIVERSITÉ BORDEAUX 1 SCIENCES ET TECHNOLOGIES MASTER THESIS COMPUTING MODULAR POLYNOMIALS WITH THETA FUNCTIONS ILARIA LOVATO ADVISOR: PROF.
More informationThe Arithmetic of Elliptic Curves
The Arithmetic of Elliptic Curves Sungkon Chang The Anne and Sigmund Hudson Mathematics and Computing Luncheon Colloquium Series OUTLINE Elliptic Curves as Diophantine Equations Group Laws and Mordell-Weil
More information(τ) = q (1 q n ) 24. E 4 (τ) = q q q 3 + = (1 q) 240 (1 q 2 ) (1 q 3 ) (1.1)
Automorphic forms on O s+2,2 (R) + and generalized Kac-Moody algebras. Proceedings of the International Congress of Mathematicians, Vol. 1, 2 (Zürich, 1994), 744 752, Birkhäuser, Basel, 1995. Richard E.
More informationPARITY OF THE COEFFICIENTS OF KLEIN S j-function
PARITY OF THE COEFFICIENTS OF KLEIN S j-function CLAUDIA ALFES Abstract. Klein s j-function is one of the most fundamental modular functions in number theory. However, not much is known about the parity
More informationSUMMATION TECHNIQUES
SUMMATION TECHNIQUES MATH 53, SECTION 55 (VIPUL NAIK) Corresponding material in the book: Scattered around, but the most cutting-edge parts are in Sections 2.8 and 2.9. What students should definitely
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationMock and quantum modular forms
Mock and quantum modular forms Amanda Folsom (Amherst College) 1 Ramanujan s mock theta functions 2 Ramanujan s mock theta functions 1887-1920 3 Ramanujan s mock theta functions 1887-1920 4 History S.
More informationINTRODUCTION TO COMPUTATIONAL MATHEMATICS
INTRODUCTION TO COMPUTATIONAL MATHEMATICS Course Notes for CM 271 / AMATH 341 / CS 371 Fall 2007 Instructor: Prof. Justin Wan School of Computer Science University of Waterloo Course notes by Prof. Hans
More informationNUMERICAL MATHEMATICS & COMPUTING 6th Edition
NUMERICAL MATHEMATICS & COMPUTING 6th Edition Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole www.engage.com www.ma.utexas.edu/cna/nmc6 September 1, 2011 2011 1 / 42 1.1 Mathematical
More informationCPSC 518 Introduction to Computer Algebra Schönhage and Strassen s Algorithm for Integer Multiplication
CPSC 518 Introduction to Computer Algebra Schönhage and Strassen s Algorithm for Integer Multiplication March, 2006 1 Introduction We have now seen that the Fast Fourier Transform can be applied to perform
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationElliptic functions. MATH 681: Lecture Notes. Part I. 1 Evaluating ds cos s. 18 June 2014
MATH 68: Lecture Notes 8 June 4 Part I Elliptic functions Evaluating cos s Consider cos s First we can substitute u cos s and arrive at u u p3 u, where p 3 is a third-degree polynomial We can also write
More informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
More informationClosed Form Solutions
Closed Form Solutions Mark van Hoeij 1 Florida State University ISSAC 2017 1 Supported by NSF 1618657 1 / 26 What is a closed form solution? Example: Solve this equation for y = y(x). y = 4 x 3 (1 x) 2
More informationRigorous Polynomial Approximations and Applications
Rigorous Polynomial Approximations and Applications Mioara Joldeș under the supervision of: Nicolas Brisebarre and Jean-Michel Muller École Normale Supérieure de Lyon, Arénaire Team, Laboratoire de l Informatique
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationAN ARITHMETIC FORMULA FOR THE PARTITION FUNCTION
AN ARITHMETIC FORMULA FOR THE PARTITION FUNCTION KATHRIN BRINGMANN AND KEN ONO 1 Introduction and Statement of Results A partition of a non-negative integer n is a non-increasing sequence of positive integers
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More informationA Motivated Introduction to Modular Forms
May 3, 2006 Outline of talk: I. Motivating questions II. Ramanujan s τ function III. Theta Series IV. Congruent Number Problem V. My Research Old Questions... What can you say about the coefficients of
More informationMATH COMPLEX ANALYSIS. Contents
MATH 3964 - OMPLEX ANALYSIS ANDREW TULLOH AND GILES GARDAM ontents 1. ontour Integration and auchy s Theorem 2 1.1. Analytic functions 2 1.2. ontour integration 3 1.3. auchy s theorem and extensions 3
More informationTHE CIRCLE METHOD, THE j FUNCTION, AND PARTITIONS
THE CIRCLE METHOD, THE j FUCTIO, AD PARTITIOS JEREMY BOOHER. Introduction The circle method was first used by Hardy and Ramanujan [] in 98 to investigate the asymptotic growth of the partition function,
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationApplications of Complex Multiplication of Elliptic Curves
Applications of Complex Multiplication of Elliptic Curves MASTER THESIS Candidate: Massimo CHENAL Supervisor: Prof. Jean-Marc COUVEIGNES UNIVERSITÀ DEGLI STUDI DI PADOVA UNIVERSITÉ BORDEAUX 1 Facoltà di
More informationON THE MODULARITY OF CERTAIN FUNCTIONS FROM THE GROMOV-WITTEN THEORY OF ELLIPTIC ORBIFOLDS
ON THE MODULARITY OF CERTAIN FUNCTIONS FROM THE GROMOV-WITTEN THEORY OF ELLIPTIC ORBIFOLDS KATHRIN BRINGMANN, LARRY ROLEN, AND SANDER ZWEGERS Abstract. In this paper, we study modularity of several functions
More informationCompute the behavior of reality even if it is impossible to observe the processes (for example a black hole in astrophysics).
1 Introduction Read sections 1.1, 1.2.1 1.2.4, 1.2.6, 1.3.8, 1.3.9, 1.4. Review questions 1.1 1.6, 1.12 1.21, 1.37. The subject of Scientific Computing is to simulate the reality. Simulation is the representation
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationMath221: HW# 7 solutions
Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More information7 Asymptotics for Meromorphic Functions
Lecture G jacques@ucsd.edu 7 Asymptotics for Meromorphic Functions Hadamard s Theorem gives a broad description of the exponential growth of coefficients in power series, but the notion of exponential
More informationFour-Dimensional GLV Scalar Multiplication
Four-Dimensional GLV Scalar Multiplication ASIACRYPT 2012 Beijing, China Patrick Longa Microsoft Research Francesco Sica Nazarbayev University Elliptic Curve Scalar Multiplication A (Weierstrass) elliptic
More informationElliptical Integrals and Functions
Appendix A Elliptical Integrals and Functions Elliptic integrals and functions are mathematical objects, which nowadays are often omitted in the mathematical curricula of universities. One quite trivial
More informationMA/CSSE 473 Day 05. Factors and Primes Recursive division algorithm
MA/CSSE 473 Day 05 actors and Primes Recursive division algorithm MA/CSSE 473 Day 05 HW 2 due tonight, 3 is due Monday Student Questions Asymptotic Analysis example: summation Review topics I don t plan
More informationComputation of the omega function
Computation of the omega function David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) University of California, Davis, Department of Computer Science 1 / 17 Mordell-Tornheim-Witten
More informationQualifying Exam Complex Analysis (Math 530) January 2019
Qualifying Exam Complex Analysis (Math 53) January 219 1. Let D be a domain. A function f : D C is antiholomorphic if for every z D the limit f(z + h) f(z) lim h h exists. Write f(z) = f(x + iy) = u(x,
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More information1 Theta functions and their modular properties
Week 3 Reading material from the books Polchinski, chapter 7 Ginspargs lectures, chapter 7 Theta functions and their modular properties The theta function is one of the basic functions that appears again
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More information