-
Notifications
You must be signed in to change notification settings - Fork 16
Open
Description
I'm trying to solve a system of eight equations involving 8 unknowns. I don't think it's too complicated but I got the following "too many right parenthese" error:
julia> phy_in_alphas = RDAB.solve(eqns, tuple(Symbol.(:phy_, 1:nphy)...) ) # use of "phi" hangs !!!
ERROR: Reduce:
solve({-(+(phy_1,phy_4))=-(+(alpha_1,alpha_2,alpha_3,alpha_4)),-(-(-(*(phy_1,
phy_4),*(phy_2,phy_3)),phy_5),phy_8)=+(*(alpha_1,alpha_2),*(alpha_1,alpha_3),*(
alpha_1,alpha_4),*(alpha_2,alpha_3),*(alpha_2,alpha_4),*(alpha_3,alpha_4)),+(-(-
(*(phy_1,phy_8),*(phy_2,phy_7)),*(phy_3,phy_6)),*(phy_4,phy_5))=-(+(*(alpha_1,
alpha_2,alpha_3),*(alpha_1,alpha_2,alpha_4),*(alpha_1,alpha_3,alpha_4),*(alpha_2
,alpha_3,alpha_4))),-(*(phy_5,phy_8),*(phy_6,phy_7))=*(alpha_1,alpha_2,alpha_3,
alpha_4),$$$);
at line 1
***** Too many right parentheses
Stacktrace:
[1] ReduceCheck at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/Reduce.jl:60 [inlined]
[2] read(::Reduce.PSL) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/Reduce.jl:95
[3] readsp at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/Reduce.jl:99 [inlined]
[4] rcall(::Reduce.RExpr; on::Array{Symbol,1}, off::Array{Symbol,1}) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/rexpr.jl:473
[5] rcall(::Reduce.RExpr) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/rexpr.jl:438
[6] |> at ./operators.jl:823 [inlined]
[7] parse_args(::String, ::Reduce.RExpr; be::Int64) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/parser.jl:343
[8] r_args(::String, ::SubString{String}, ::Array{Any,1}; be::Int64) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/parser.jl:383
[9] parse_args(::String, ::Reduce.RExpr, ::Reduce.RExpr; be::Int64) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/parser.jl:331
[10] #solve#207 at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/parser.jl:437 [inlined]
[11] solve at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/parser.jl:437 [inlined]
[12] solve(::NTuple{8,Expr}, ::NTuple{8,Symbol}) at /Users/Thomas/.julia/packages/Reduce/HgVTl/src/args.jl:230
[13] top-level scope at REPL[25]:1
julia> for i in 1:length(eqns)
println(i, " : ", eqns[i])
end
1 : -((phy_1 + phy_4)) == -((alpha_1 + alpha_2 + alpha_3 + alpha_4))
2 : ((phy_1 * phy_4 - phy_2 * phy_3) - phy_5) - phy_8 == alpha_1 * alpha_2 + alpha_1 * alpha_3 + alpha_1 * alpha_4 + alpha_2 * alpha_3 + alpha_2 * alpha_4 + alpha_3 * alpha_4
3 : ((phy_1 * phy_8 - phy_2 * phy_7) - phy_3 * phy_6) + phy_4 * phy_5 == -((alpha_1 * alpha_2 * alpha_3 + alpha_1 * alpha_2 * alpha_4 + alpha_1 * alpha_3 * alpha_4 + alpha_2 * alpha_3 * alpha_4))
4 : phy_5 * phy_8 - phy_6 * phy_7 == alpha_1 * alpha_2 * alpha_3 * alpha_4
5 : :phy_3 == :alpha_5
6 : :phy_6 == :alpha_6
7 : :phy_7 == :alpha_7
8 : :phy_8 == :alpha_8
is there any way to prevent the error by setting something? thanks.
Metadata
Metadata
Assignees
Labels
No labels