convexengineering / gpfit

Fit posynomials to data
http://gpfit.readthedocs.io/en/latest/
MIT License
10 stars 7 forks source link

index error #67

Closed mayork closed 7 years ago

mayork commented 7 years ago

@bqpd I'm doing some D8 fits (the file is naca_cl0_fits.py the Tail Fits folder at commit https://github.com/hoburg/d8/commit/851b1d427caace8369a0f2dd1dbcbf4a917c6438 on D8 master. GPfit is throwing the following error and neither I nor @1ozturkbe know what it is..

---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
/Users/mayork/Documents/GpGit/d8/Tail Fits/naca_cl0_fits.py in <module>()
    100     X, Y = fit_setup(NACA, Re) # call fit(X, Y, 4, "SMA") to get fit
    101     F, A = plot_fits(NACA, Re)
--> 102     make_fit(NACA, Re)
    103     F.savefig("tail_fits/taildragpolar.pdf",
    104               bbox_inches="tight")

/Users/mayork/Documents/GpGit/d8/Tail Fits/naca_cl0_fits.py in make_fit(naca_range, re_range)
     65     print np.size(x)
     66     print np.size(y)
---> 67     fit(x, y, 3)
     68 
     69 def plot_fits(naca_range, re_range):

/Users/mayork/Documents/GpGit/gpfit/gpfit/fit.pyc in fit(xdata, ydata, K, ftype)
     70         w = Variable("w")
     71 
---> 72     params = get_params(ftype, K, xdata, ydata)
     73 
     74     # A: exponent parameters, B: coefficient parameters

/Users/mayork/Documents/GpGit/gpfit/gpfit/fit.pyc in get_params(ftype, K, xdata, ydata)
     24         return r, drdp
     25 
---> 26     ba = ba_init(xdata, ydata.reshape(ydata.size, 1), K).flatten('F')
     27 
     28     if ftype == "ISMA":

/Users/mayork/Documents/GpGit/gpfit/gpfit/ba_init.pyc in ba_init(x, y, K)
     77                       "full rank for local fitting." % (i-iinit, k))
     78         # now create the local fit
---> 79         b[:, k] = lstsq(X[inds.nonzero()], y[inds.nonzero()])[0][:, 0]
     80 
     81     return b

IndexError: index 59 is out of bounds for axis 0 with size 59
mayork commented 7 years ago

figured out this was a pandas error in the xfoil parsing. closing now.

bqpd commented 7 years ago

++