From: Michael Orlitzky Date: Sun, 4 Nov 2018 06:13:16 +0000 (-0500) Subject: orthogonal_polynomials.py: use xrange everywhere. X-Git-Url: http://gitweb.michael.orlitzky.com/?a=commitdiff_plain;ds=sidebyside;h=5f0f517a80ae4f3e12bf4d3a84ea3d182cb96e70;p=sage.d.git orthogonal_polynomials.py: use xrange everywhere. --- diff --git a/mjo/orthogonal_polynomials.py b/mjo/orthogonal_polynomials.py index 7e094fc..fd5fd75 100644 --- a/mjo/orthogonal_polynomials.py +++ b/mjo/orthogonal_polynomials.py @@ -82,7 +82,7 @@ def legendre_p(n, x, a = -1, b = 1): ... for b in basis]) ... sage: f = sin(x) - sage: legendre_basis = [ legendre_p(k, x, a, b) for k in range(0,4) ] + sage: legendre_basis = [ legendre_p(k, x, a, b) for k in xrange(4) ] sage: proj = project(legendre_basis, f) sage: proj.simplify_trig() 5/2*(7*(pi^2 - 15)*x^3 - 3*(pi^4 - 21*pi^2)*x)/pi^6 @@ -92,7 +92,7 @@ def legendre_p(n, x, a = -1, b = 1): We should agree with Maxima for all `n`:: sage: eq = lambda k: bool(legendre_p(k,x) == legendre_P(k,x)) - sage: all([eq(k) for k in range(0,20) ]) # long time + sage: all([eq(k) for k in xrange(20) ]) # long time True We can evaluate the result of the zeroth polynomial:: @@ -196,6 +196,6 @@ def legendre_p(n, x, a = -1, b = 1): # From Abramowitz & Stegun, (22.3.2) with alpha = beta = 0. # Also massaged to support finite field elements. - P = sum([ c(m)*g(m) for m in range(0,n+1) ])/(2**n) + P = sum([ c(m)*g(m) for m in xrange(n+1) ])/(2**n) return P