anj1 / algebraicnumbers.jl Goto Github PK

i = sqrt(AlgebraicNumber(-1))
# ≈0.0 + 1.0im
uno  = AlgebraicNumber(1)
# ≈1.0 + 0.0im
two = AlgebraicNumber(2)
# ≈2.0 + 0.0im
ipu = 1+i

[ipu == two, two.coeff == ipu.coeff, two.prec == ipu.prec, two.apprx == ipu.apprx]'
#  true  true  true  false

(i+1).coeff'
# -2  1
AlgebraicNumbers.composed_sum(uno.coeff, i.coeff)'
# -2  1

julia> using Pkg
julia> Pkg.add("AlgebraicNumber")
ERROR: The following package names could not be resolved:
 * AlgebraicNumber (not found in project, manifest or registry)
Please specify by known `name=uuid`.

julia> convert(Int, AlgebraicNumber(2))
-2

using LinearAlgebra
using AlgebraicNumbers


A = AlgebraicNumber.([
    1 2
    3 4
])

λ = AlgebraicNumber(2)

B = A - λ*I # error here

ERROR: LoadError: MethodError: promote_rule(::Type{Bool}, ::Type{AlgebraicNumber{BigInt, BigFloat}}) is ambiguous.

Candidates:
  promote_rule(x::Type{T}, y::Type{AlgebraicNumber{S, F}}) where {T<:Integer, S, F}
    @ AlgebraicNumbers ~/.julia/packages/AlgebraicNumbers/H0atE/src/promote.jl:14
  promote_rule(::Type{Bool}, ::Type{T}) where T<:Number
    @ Base bool.jl:4

Possible fix, define
  promote_rule(::Type{Bool}, ::Type{AlgebraicNumber{S, F}}) where {S, F}

Stacktrace:
 [1] promote_type(#unused#::Type{AlgebraicNumber{BigInt, BigFloat}}, #unused#::Type{Bool})
   @ Base ./promotion.jl:307
 [2] _promote
   @ ./promotion.jl:357 [inlined]
 [3] promote
   @ ./promotion.jl:381 [inlined]
 [4] *(x::AlgebraicNumber{BigInt, BigFloat}, y::Bool)
   @ Base ./promotion.jl:411
 [5] *(x::AlgebraicNumber{BigInt, BigFloat}, J::UniformScaling{Bool})
   @ LinearAlgebra ~/julia-1.9.3/share/julia/stdlib/v1.9/LinearAlgebra/src/uniformscaling.jl:263
 [6] top-level scope
   @ ~/MSU/Code/Quadrics/bug.jl:12
 [7] include(fname::String)
   @ Base.MainInclude ./client.jl:478
 [8] top-level scope
   @ REPL[14]:1

using LinearAlgebra
using NiceNumbers

A = NiceNumber.([
    1 2
    3 4
])

λ = NiceNumber(2)

B = A - λ*I

2×2 Matrix{NiceNumber}:
 -1  2
  3  2

using LinearAlgebra
using AlgebraicNumbers


M = AlgebraicNumber.([
    1 2
    3 4
])

b = AlgebraicNumber.([5; 6])

println(typeof(M))
println(typeof(b))

x = M * b

println(typeof(x))

using LinearAlgebra
using AlgebraicNumbers


M = Rational.([
    1 2
    3 4
])

b = Rational.([5; 6])

println(typeof(M))
println(typeof(b))

x = M * b

println(typeof(x))