Fix minor bugs (#399)

* Fix issues 397 and 398

* Add tests (iss 398)

* Bump patch version

Author: lbenet

Project.toml

@@ -1,6 +1,6 @@
 name = "TaylorSeries"
 uuid = "6aa5eb33-94cf-58f4-a9d0-e4b2c4fc25ea"
-version = "0.21.1"
+version = "0.21.2"
 repo = "https://github.com/JuliaDiff/TaylorSeries.jl.git"
 
 [deps]

src/other_functions.jl

@@ -10,7 +10,7 @@ for T in (:Taylor1, :HomogeneousPolynomial, :TaylorN)
 
     ## real, imag, conj and ctranspose ##
     for f in (:real, :imag, :conj)
-        @eval ($f)(a::$T) = $T($f(a.coeffs))
+        @eval ($f)(a::$T) = $T($f(a.coeffs), get_order(a))
     end
 
     @eval adjoint(a::$T) = conj(a)
@@ -29,7 +29,7 @@ for op in (:mod, :rem)
             function ($op)(a::$T{T}, x::T) where {T<:Real}
                 coeffs = copy(a.coeffs)
                 @inbounds coeffs[1] = ($op)(constant_term(a), x)
-                return $T(coeffs)
+                return $T(coeffs, get_order(a))
             end
 
             function ($op)(a::$T{T}, x::S) where {T<:Real,S<:Real}
@@ -74,7 +74,7 @@ for T in (:Taylor1, :TaylorN)
         function mod2pi(a::$T{T}) where {T<:Real}
             coeffs = copy(a.coeffs)
             @inbounds coeffs[1] = mod2pi( constant_term(a) )
-            return $T(coeffs)
+            return $T(coeffs, get_order(a))
         end
 
         function abs(a::$T{T}) where {T<:Real}

src/parameters.jl

@@ -188,7 +188,7 @@ function set_variables(::Type{R}, names::T; order=get_order(), numvars=-1) where
 
     if length(variable_names) == 1 && numvars ≥ 1
         name = variable_names[1]
-        variable_names = string.("x", TS.subscriptify.(1:numvars))
+        variable_names = string.(name, TS.subscriptify.(1:numvars))
     end
 
     return set_variables(R, variable_names, order=order)

test/manyvariables.jl

@@ -31,6 +31,9 @@ using Test
     @test TS.index_table[2][1] == 3^65
     @test TS.pos_table[2][3^64] == 2
 
+    @test eltype(set_variables("a", order=1, numvars=2))  == TaylorN{Float64}
+    @test TS._params_TaylorN_.variable_names[1] == "a₁"
+    @test TS._params_TaylorN_.variable_symbols[2] == :a₂
     @test eltype(set_variables(Int, "x", numvars=2, order=6)) == TaylorN{Int}
     @test eltype(set_variables("x", numvars=2, order=6)) == TaylorN{Float64}
     @test eltype(set_variables(BigInt, "x y", order=6)) == TaylorN{BigInt}