Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue ngcd #318

Merged
merged 40 commits into from
Mar 23, 2021
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
40 commits
Select commit Hold shift + click to select a range
657112d
WIP
jverzani Feb 5, 2021
2f5080e
WIP
jverzani Feb 6, 2021
154052f
WIP
jverzani Feb 7, 2021
c52febd
move symbol into type
jverzani Feb 9, 2021
409d49e
clean up
jverzani Feb 9, 2021
2f91ef3
change what iteration means for a polynomial
jverzani Feb 9, 2021
5a84d9f
adjust integrate to return polynomial of same type, even if NaN
jverzani Feb 10, 2021
a2fc233
clean up offset vector tests
jverzani Feb 11, 2021
5a47053
drop OffsetArrays dependency; clean up scalar operations for Immutabl…
jverzani Feb 14, 2021
7a350ca
cleanup
jverzani Feb 14, 2021
166522b
change error into ArgumentError
jverzani Feb 14, 2021
b398f03
replace isnothing; standardize check for same variable
jverzani Feb 15, 2021
1426b38
registerN, bug
jverzani Feb 15, 2021
b3ca9c0
update documentation
jverzani Feb 15, 2021
fc3014b
minor edits
jverzani Feb 17, 2021
206de82
cleanup, docfix
jverzani Feb 18, 2021
5d52918
clean up integration
jverzani Feb 19, 2021
a031b8e
adjust doc string
jverzani Feb 19, 2021
93c15ea
WIP
jverzani Feb 20, 2021
264189d
WIP
jverzani Feb 22, 2021
b9324b6
WIP
jverzani Feb 22, 2021
5349be1
work on promotion for + and * operations
jverzani Feb 22, 2021
c0d688c
edits
jverzani Feb 23, 2021
38110e0
laurent polynomial depwarns removed
jverzani Feb 23, 2021
dfe1d73
run doctests
jverzani Feb 23, 2021
dff6629
WIP
jverzani Feb 24, 2021
65a7c92
fix basis symbol handling
jverzani Feb 24, 2021
e57784d
WIP
jverzani Feb 25, 2021
7470f25
WIP
jverzani Feb 25, 2021
316a1ff
adjust evaluation so that implemented evalpoly instead of call syntax…
jverzani Feb 25, 2021
3cfb7f8
add example to extending; specialize evalpoly rather than call syntax
jverzani Feb 25, 2021
226ac83
Merge branch 'v2.0.0-add-examples' into v2.0.0
jverzani Feb 25, 2021
456bc09
refactor truncate! chop!
jverzani Feb 26, 2021
c5e155e
oops
jverzani Feb 26, 2021
1e89476
move things into common (truncate, chop, zero, one,variable, basis)
jverzani Mar 2, 2021
8d30e68
merge changes
jverzani Mar 23, 2021
3450f5e
close issue #316
jverzani Mar 23, 2021
9baddaa
remove changes from a different branch
jverzani Mar 23, 2021
3cd3e77
fix issue with ngcd
jverzani Mar 23, 2021
62dea56
version bump
jverzani Mar 23, 2021
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
registerN, bug
  • Loading branch information
jverzani committed Feb 15, 2021
commit 1426b38e353886279d5494f1600088dedae0c2d9
63 changes: 8 additions & 55 deletions src/abstract.jl
Original file line number Diff line number Diff line change
Expand Up @@ -3,7 +3,7 @@ export AbstractPolynomial
const SymbolLike = Union{AbstractString,Char,Symbol}

"""
AbstractPolynomial{T}
AbstractPolynomial{T,X}

An abstract container for various polynomials.

Expand All @@ -27,7 +27,7 @@ Given a polynomial with `name`, creates some common convenience constructors and

# Example
```julia
struct MyPolynomial{T} <: AbstractPolynomial{T} end
struct MyPolynomial{T,X} <: AbstractPolynomial{T,X} end

Polynomials.@register MyPolynomial
```
Expand Down Expand Up @@ -72,68 +72,21 @@ macro registerN(name, params...)
quote
Base.convert(::Type{P}, q::Q) where {$(αs...),T, P<:$poly{$(αs...),T}, Q <: $poly{$(αs...),T}} = q
Base.convert(::Type{$poly{$(αs...)}}, q::Q) where {$(αs...),T, Q <: $poly{$(αs...),T}} = q
Base.promote(p::P, q::Q) where {$(αs...),T, P <:$poly{$(αs...),T}, Q <: $poly{$(αs...),T}} = p,q
Base.promote_rule(::Type{<:$poly{$(αs...),T}}, ::Type{<:$poly{$(αs...),S}}) where {$(αs...),T,S} =
$poly{$(αs...),promote_type(T, S)}
Base.promote_rule(::Type{<:$poly{$(αs...),T}}, ::Type{S}) where {$(αs...),T,S<:Number} =
$poly{$(αs...),promote_type(T,S)}
Base.promote(p::P, q::Q) where {$(αs...),T, X, P <:$poly{$(αs...),T,X}, Q <: $poly{$(αs...),T,X}} = p,q
Base.promote_rule(::Type{<:$poly{$(αs...),T,X}}, ::Type{<:$poly{$(αs...),S,X}}) where {$(αs...),T,S,X} =
$poly{$(αs...),promote_type(T, S),X}
Base.promote_rule(::Type{<:$poly{$(αs...),T,X}}, ::Type{S}) where {$(αs...),T,X,S<:Number} =
$poly{$(αs...),promote_type(T,S),X}

function $poly{$(αs...),T}(x::AbstractVector{S}, var::SymbolLike = :x) where {$(αs...),T,S}
$poly{$(αs...),T}(T.(x), Symbol(var))
end
$poly{$(αs...)}(coeffs::AbstractVector{T}, var::SymbolLike=:x) where {$(αs...),T} =
$poly{$(αs...),T}(coeffs, Symbol(var))
$poly{$(αs...),T,Symbol(var)}(coeffs)
$poly{$(αs...),T}(n::Number, var::SymbolLike = :x) where {$(αs...),T} = n*one($poly{$(αs...),T}, Symbol(var))
$poly{$(αs...)}(n::Number, var::SymbolLike = :x) where {$(αs...)} = n*one($poly{$(αs...)}, Symbol(var))
$poly{$(αs...),T}(var::SymbolLike=:x) where {$(αs...), T} = variable($poly{$(αs...),T}, Symbol(var))
$poly{$(αs...)}(var::SymbolLike=:x) where {$(αs...)} = variable($poly{$(αs...)}, Symbol(var))
end
end


# deprecated. If desired, replace with @registerN type parameters... macro
# Macros to register POLY{α, T} and POLY{α, β, T}
macro register1(name)
@warn "@register1 is deprecated use @registerN"
poly = esc(name)
quote
Base.convert(::Type{P}, p::P) where {P<:$poly} = p
Base.promote(p::P, q::Q) where {α,T, P <:$poly{α,T}, Q <: $poly{α,T}} = p,q
Base.promote_rule(::Type{<:$poly{α,T}}, ::Type{<:$poly{α,S}}) where {α,T,S} =
$poly{α,promote_type(T, S)}
Base.promote_rule(::Type{<:$poly{α,T}}, ::Type{S}) where {α,T,S<:Number} =
$poly{α,promote_type(T,S)}
function $poly{α,T}(x::AbstractVector{S}, var::SymbolLike = :x) where {α,T,S}
$poly{α,T}(T.(x), Symbol(var))
end
$poly{α}(coeffs::AbstractVector{T}, var::SymbolLike=:x) where {α,T} =
$poly{α,T}(coeffs, Symbol(var))
$poly{α,T}(n::Number, var::SymbolLike = :x) where {α,T} = n*one($poly{α,T}, Symbol(var))
$poly{α}(n::Number, var::SymbolLike = :x) where {α} = n*one($poly{α}, Symbol(var))
$poly{α,T}(var::SymbolLike=:x) where {α, T} = variable($poly{α,T}, Symbol(var))
$poly{α}(var::SymbolLike=:x) where {α} = variable($poly{α}, Symbol(var))
end
end


# Macro to register POLY{α, β, T}
macro register2(name)
@warn "@register2 is deprecated use @registerN"
poly = esc(name)
quote
Base.convert(::Type{P}, p::P) where {P<:$poly} = p
Base.promote(p::P, q::Q) where {α,β,T, P <:$poly{α,β,T}, Q <: $poly{α,β,T}} = p,q
Base.promote_rule(::Type{<:$poly{α,β,T}}, ::Type{<:$poly{α,β,S}}) where {α,β,T,S} =
$poly{α,β,promote_type(T, S)}
Base.promote_rule(::Type{<:$poly{α,β,T}}, ::Type{S}) where {α,β,T,S<:Number} =
$poly{α,β,promote_type(T, S)}
$poly{α,β}(coeffs::AbstractVector{T}, var::SymbolLike = :x) where {α,β,T} =
$poly{α,β,T}(coeffs, Symbol(var))
$poly{α,β,T}(x::AbstractVector{S}, var::SymbolLike = :x) where {α,β,T,S<:Number} = $poly{α,β,T}(T.(x), var)
$poly{α,β,T}(n::Number, var::SymbolLike = :x) where {α,β,T} = n*one($poly{α,β,T}, Symbol(var))
$poly{α,β}(n::Number, va::SymbolLiker = :x) where {α,β} = n*one($poly{α,β}, Symbol(var))
$poly{α,β,T}(var::SymbolLike=:x) where {α,β, T} = variable($poly{α,β,T}, Symbol(var))
$poly{α,β}(var::SymbolLike=:x) where {α,β} = variable($poly{α,β}, Symbol(var))
end
end

2 changes: 1 addition & 1 deletion src/polynomials/ImmutablePolynomial.jl
Original file line number Diff line number Diff line change
Expand Up @@ -66,7 +66,7 @@ function ImmutablePolynomial{T,X}(coeffs::AbstractVector{S}) where {T,X,S}
@warn "ignoring the axis offset of the coefficient vector"
end
N = findlast(!iszero, coeffs)
N = nothing && return ImmutablePolynomial{R,X,0}(())
N == nothing && return ImmutablePolynomial{R,X,0}(())
N′ = N + 1 - firstindex(coeffs)
cs = NTuple{N′,T}(coeffs[i] for i ∈ firstindex(coeffs):N)
ImmutablePolynomial{T, X, N′}(cs)
Expand Down