From c32a97123a3998748b0ac33a33f4ef1d765b5868 Mon Sep 17 00:00:00 2001 From: Unknown Date: Tue, 29 May 2018 17:10:49 -0500 Subject: [PATCH] Define alpha and cleanup tests --- README.md | 6 +----- src/OrthogonalPolynomials.jl | 3 ++- test/runtests.jl | 6 +----- utils.jl | 2 ++ 4 files changed, 6 insertions(+), 11 deletions(-) diff --git a/README.md b/README.md index e549812..80819d0 100644 --- a/README.md +++ b/README.md @@ -1,10 +1,6 @@ # OrthogonalPolynomials -[![Build Status](https://travis-ci.org/miguelraz/OrthogonalPolynomials.jl.svg?branch=master)](https://travis-ci.org/miguelraz/OrthogonalPolynomials.jl) - -[![Coverage Status](https://coveralls.io/repos/miguelraz/OrthogonalPolynomials.jl/badge.svg?branch=master&service=github)](https://coveralls.io/github/miguelraz/OrthogonalPolynomials.jl?branch=master) - -[![codecov.io](http://codecov.io/github/miguelraz/OrthogonalPolynomials.jl/coverage.svg?branch=master)](http://codecov.io/github/miguelraz/OrthogonalPolynomials.jl?branch=master) +[![Build Status](https://travis-ci.org/miguelraz/OrthogonalPolynomials.jl.svg?branch=master)](https://travis-ci.org/miguelraz/OrthogonalPolynomials.jl) [![Coverage Status](https://coveralls.io/repos/miguelraz/OrthogonalPolynomials.jl/badge.svg?branch=master&service=github)](https://coveralls.io/github/miguelraz/OrthogonalPolynomials.jl?branch=master) [![codecov.io](http://codecov.io/github/miguelraz/OrthogonalPolynomials.jl/coverage.svg?branch=master)](http://codecov.io/github/miguelraz/OrthogonalPolynomials.jl?branch=master) [![Build status](https://ci.appveyor.com/api/projects/status/tgs1489v06dlk3u1/branch/master?svg=true)](https://ci.appveyor.com/project/miguelraz/orthogonalpolynomials-jl/branch/master) ## Introduction diff --git a/src/OrthogonalPolynomials.jl b/src/OrthogonalPolynomials.jl index 7c68a14..3408a2f 100644 --- a/src/OrthogonalPolynomials.jl +++ b/src/OrthogonalPolynomials.jl @@ -8,6 +8,7 @@ b(n, m) = n-m+1 c(m, α) = m*(α+m) k(n, α=0) = [-b(n,i)*inv(c(i, α)) * d(n, α) for i in 1:n] f(x) = x +const α = 0 a(x, n, ks = k(n,α), i=0) = i == n ? :(1) : return :(muladd( $(ks[i+1]*f(x)), $((a)(x, n, ks, i+1)) , 1)) @@ -20,6 +21,6 @@ macro a(x,n,ks=k(n,α)) return :($ex) end -export d,b,c,k,f,a +export d,b,c,k,f,a,α end # module diff --git a/test/runtests.jl b/test/runtests.jl index a48dd1f..604f2c2 100644 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -44,9 +44,6 @@ function laguerre_corrected(n, x) l1 end -goal12(x) = @evalpoly x 1.0 -12.0 33.0 -36.666666666666664 20.625 -6.6 1.2833333333333334 -0.15714285714285714 0.012276785714285714 -0.0006062610229276896 1.8187830687830687e-5 -3.0062530062530064e-7 - - # Abramowitz and Stegun pg. 799 LaguerreTable = [ 1 1 1 2 6 24 120 720 5040 40320 362880 3628800 39916800 479001600; @@ -95,9 +92,8 @@ TestLaguerreTable = [ @test all([laguerre_corrected(i-1, testvals[j]) ≈ TestLaguerreTable[i,j] for i in 1:13, j in 2:6]) end @testset "Stubborn Way" begin + const α = 0 @test all(eval.(a(0,i)) .≈ 1 for i in 1:12) @test all(eval.([a(testvals[j],i) for i in 0:12, j in 2:6]) .≈ TestLaguerreTable[:,2:6]) end end - -goal12(x) = @horner x 1.0 -12.0 33.0 -36.666666666666664 20.625 -6.6 1.2833333333333334 -0.15714285714285714 0.012276785714285714 -0.0006062610229276896 1.8187830687830687e-5 -3.0062530062530064e-7 diff --git a/utils.jl b/utils.jl index 8ccc517..fb434ce 100644 --- a/utils.jl +++ b/utils.jl @@ -33,3 +33,5 @@ function T(n,k) 0 end end + +goal12(x) = @evalpoly x 1.0 -12.0 33.0 -36.666666666666664 20.625 -6.6 1.2833333333333334 -0.15714285714285714 0.012276785714285714 -0.0006062610229276896 1.8187830687830687e-5 -3.0062530062530064e-7