bug

[jguterl]

running "fit(Bernstein, t, ydata, N)"

bug related to commit 972e991 line 157 in common.jl

ERROR: MethodError: objects of type Symbol are not callable

Some of the types have been truncated in the stacktrace for improved reading. To emit complete information
in the stack trace, evaluate `TruncatedStacktraces.VERBOSE[] = true` and re-run the code.

Stacktrace:
  [1] _broadcast_getindex_evalf
    @ ./broadcast.jl:670 [inlined]
  [2] _broadcast_getindex
    @ ./broadcast.jl:643 [inlined]
  [3] getindex
    @ ./broadcast.jl:597 [inlined]
  [4] copy
    @ ./broadcast.jl:899 [inlined]
  [5] materialize
    @ ./broadcast.jl:860 [inlined]
  [6] (SpecialPolynomials.Bernstein{Float64, :x})(x::Vector{Float64}, var::Polynomials.Var{:x}) (repeats 2 times)
    @ SpecialPolynomials ~/.julia/packages/Polynomials/agWXR/src/abstract.jl:140
  [7] _fit(P::Type{SpecialPolynomials.Bernstein}, x::Vector{Float64}, y::Vector{Float64}, deg::Int64; weights::Nothing, var::Symbol)
    @ Polynomials ~/.julia/packages/Polynomials/agWXR/src/common.jl:157
  [8] #fit#17
    @ ~/.julia/packages/Polynomials/agWXR/src/common.jl:118 [inlined]
  [9] fit
    @ ~/.julia/packages/Polynomials/agWXR/src/common.jl:112 [inlined]
 [10] PlasmaFacingSurfaces.RZBernsteinFit(xdata::Vector{Float64}, ydata::Vector{Float64}, N::Int64)
    @ PlasmaFacingSurfaces ~/.julia/dev/PlasmaFacingSurfaces/src/fit.jl:22
 [11] PlasmaFacingSurfaces.RZBernsteinFit()
    @ PlasmaFacingSurfaces ~/.julia/dev/PlasmaFacingSurfaces/src/fit.jl:16
 [12] PlasmaFacingSurfaces.MainChamberWall(id::Vector{Symbol})
    @ PlasmaFacingSurfaces ~/.julia/dev/PlasmaFacingSurfaces/src/main_chamber_walls.jl:12
 [13] PlasmaFacingSurfaces.MainChamberWalls(eqt::IMASDD.equilibrium__time_slice{Float64})
    @ PlasmaFacingSurfaces ~/.julia/dev/PlasmaFacingSurfaces/src/main_chamber_walls.jl:14
 [14] PlasmaFacingSurfaces.PFSDesign(eqt::IMASDD.equilibrium__time_slice{Float64})
    @ PlasmaFacingSurfaces ~/.julia/dev/PlasmaFacingSurfaces/src/PlasmaFacingSurfaces.jl:38
 [15] top-level scope
    @ ./REPL[6]:41

[jverzani]

Thx for the report. I’ll have a look.

[jverzani]

It is a missing method in SpecialPolynomials. I'll work out a fix. In the interim, you could fit to the Polynomial type and call convert, but that is a bit of a hacky solution.

[jverzani]

Closed by #478

[jguterl]

thanks!

…

On Thu, Mar 9, 2023 at 3:36 PM john verzani ***@***.***> wrote: Closed by #478 <#478> — Reply to this email directly, view it on GitHub <#477 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEESMZCGTCNEDEPRZNB6WUDW3JSOHANCNFSM6AAAAAAVUSXG5U> . You are receiving this because you authored the thread.Message ID: ***@***.***>