WIP: Nested Taylors

.travis.yml

@@ -5,13 +5,16 @@ os:
   - osx
 
 julia:
-    - 0.5
     - 0.6
     - nightly
 
 notifications:
     email: false
 
+matrix:
+  allow_failures:
+  - julia: nightly
+
 script:
  - if [[ -a .git/shallow ]]; then git fetch --unshallow; fi
  - julia --check-bounds=yes -e 'Pkg.clone(pwd()); Pkg.build("TaylorSeries"); Pkg.test("TaylorSeries"; coverage=true)'

README.md

@@ -3,8 +3,8 @@
 A [Julia](http://julialang.org) package for Taylor expansions in one or more
 independent variables.
 
-[![TaylorSeries](http://pkg.julialang.org/badges/TaylorSeries_0.5.svg)](http://pkg.julialang.org/?pkg=TaylorSeries)
 [![TaylorSeries](http://pkg.julialang.org/badges/TaylorSeries_0.6.svg)](http://pkg.julialang.org/?pkg=TaylorSeries)
+[![TaylorSeries](http://pkg.julialang.org/badges/TaylorSeries_0.7.svg)](http://pkg.julialang.org/?pkg=TaylorSeries)
 [![Coverage Status](https://coveralls.io/repos/JuliaDiff/TaylorSeries.jl/badge.svg?branch=master)](https://coveralls.io/github/JuliaDiff/TaylorSeries.jl?branch=master)
 
 [![](https://img.shields.io/badge/docs-stable-blue.svg)](http://www.juliadiff.org/TaylorSeries.jl/stable)

REQUIRE

@@ -1,2 +1,2 @@
-julia 0.5
-Compat 0.17.0
+julia 0.6
+Compat 0.54.0

docs/REQUIRE

@@ -0,0 +1 @@
+Documenter 0.16.1 0.16.1+

docs/make.jl

@@ -16,7 +16,7 @@ makedocs(
 deploydocs(
     repo   = "github.com/JuliaDiff/TaylorSeries.jl.git",
     target = "build",
-    julia = "0.5",
+    julia = "0.6",
     osname = "linux",
     deps   = nothing,
     make   = nothing

docs/src/api.md

@@ -6,6 +6,7 @@
 CurrentModule = TaylorSeries
 ```
 
+## Module
 ```@docs
 TaylorSeries
 ```
@@ -17,20 +18,23 @@ Taylor1
 HomogeneousPolynomial
 TaylorN
 AbstractSeries
-ParamsTaylorN
 ```
 
 ## Functions and methods
 
 ```@docs
-Taylor1([::Type{Float64}], [order::Int64=1])
-HomogeneousPolynomial{T<:Number}(::Type{T}, ::Int)
-TaylorN{T<:Number}(::Type{T}, nv::Int; [order::Int=get_order()])
+Taylor1(::Type{T}, ::Int64=1) where {T<:Number}
+HomogeneousPolynomial(::Type{T}, ::Int) where {T<:Number}
+TaylorN(::Type{T}, ::Int; ::Int=get_order()) where {T<:Number}
 set_variables
+get_variables
 show_params_TaylorN
-get_coeff
+show_monomials
+getcoeff
 evaluate
 evaluate!
+taylor_expand
+update!
 derivative
 integrate
 gradient
@@ -39,27 +43,33 @@ jacobian!
 hessian
 hessian!
 inverse
+abs
+norm
+isapprox
+isfinite
+displayBigO
 ```
 
 ## Internals
 
 ```@docs
+ParamsTaylorN
+_InternalMutFuncs
 generate_tables
 generate_index_vectors
 in_base
 make_inverse_dict
-order_posTb
 resize_coeffs1!
 resize_coeffsHP!
-zero_korder
 constant_term
 mul!
-mul!(c::HomogeneousPolynomial, a::HomogeneousPolynomial, b::HomogeneousPolynomial)
+mul!(::HomogeneousPolynomial, ::HomogeneousPolynomial, ::HomogeneousPolynomial)
+mul!(::Vector{Taylor1{T}}, ::Union{Matrix{T},SparseMatrixCSC{T}},::Vector{Taylor1{T}}) where {T<:Number}
 div!
 pow!
 square
 sqr!
-sqr!(c::HomogeneousPolynomial, a::HomogeneousPolynomial)
+sqr!(::HomogeneousPolynomial{T}, ::HomogeneousPolynomial{T}) where {T<:NumberNotSeriesN}
 sqrt!
 exp!
 log!
@@ -70,7 +80,12 @@ acos!
 atan!
 sinhcosh!
 tanh!
-A_mul_B!
+derivative!
+_internalmutfunc_call
+_dict_unary_ops
+_dict_binary_calls
+_dict_unary_calls
+_dict_binary_ops
 ```
 
 ## Index

docs/src/examples.md

@@ -20,7 +20,7 @@ The first example shows that the four-square identity holds:
 \end{eqnarray}
 ```
 which was originally proved by Euler. The code can also be found in
-[this test](../../test/identities_Euler.jl) of the package.
+[this test](https://github.com/JuliaDiff/TaylorSeries.jl/blob/master/test/identities_Euler.jl) of the package.
 
 First, we reset the maximum degree of the polynomial to 4, since the RHS
 of the equation has *a priori* terms of fourth order, and define the 8
@@ -112,7 +112,7 @@ We note that the above functions use expansions in `Int128`. This is actually
 required, since some coefficients are larger than `typemax(Int)`:
 
 ```@repl fateman
-get_coeff(f2, [1,6,7,20]) # coefficient of x y^6 z^7 w^{20}
+getcoeff(f2, (1,6,7,20)) # coefficient of x y^6 z^7 w^{20}
 ans > typemax(Int)
 length(f2)
 sum(TaylorSeries.size_table)
@@ -127,7 +127,8 @@ monomials in 4 variables.
 ### Bechmarks
 
 The functions described above have been compared against Mathematica v11.1.
-The relevant files used for benchmarking can be found [here](../../perf/).
+The relevant files used for benchmarking can be found
+[here](https://github.com/JuliaDiff/TaylorSeries.jl/tree/master/perf).
 Running on a MacPro with Intel-Xeon processors 2.7GHz, we obtain that
 Mathematica requires on average (5 runs) 3.075957 seconds for the computation,
 while for `fateman1` and `fateman2` above we obtain 2.811391 and 1.490256,