keflavich / agpy Goto Github PK

/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy
/core/fromnumeric.py:2641: VisibleDeprecationWarning: `rank` is deprecated; use the `ndim` 
attribute or function instead. To find the rank of a matrix see `numpy.linalg.matrix_rank`.
  VisibleDeprecationWarning)

from agpy import gaussfitter
mpfit2, psf_fit2 = gaussfitter.gaussfit(testimg,returnfitimage=True)

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-2-2f88faf6f9b5> in <module>()
     40 plt.imshow(testimg,interpolation="None")
     41 #tmpcs=reduction.EvalDirtyImage2(testimg,verbose=True)
---> 42 mpfit2, psf_fit2 = gaussfitter.gaussfit(testimg,returnfitimage=True)

/Users/jgirard/anaconda/lib/python2.7/site-packages/agpy/gaussfitter.pyc in gaussfit(data, err, params, autoderiv, return_all, circle, fixed, limitedmin, limitedmax, usemoment, minpars, maxpars, rotate, vheight, quiet, returnmp, returnfitimage, **kwargs)
    241 #        p, cov, infodict, errmsg, success = optimize.leastsq(errorfunction,\
    242 #                params, full_output=1)
--> 243         mp = mpfit(mpfitfun(data,err),parinfo=parinfo,quiet=quiet)
    244 
    245 

/Users/jgirard/anaconda/lib/python2.7/site-packages/agpy/mpfit/mpfit.py in __init__(self, fcn, xall, functkw, parinfo, ftol, xtol, gtol, damp, maxiter, factor, nprint, iterfunct, iterkw, nocovar, rescale, autoderivative, quiet, diag, epsfcn, debug)
   1165                 catch_msg = 'calculating LM parameter (MPFIT_)'
   1166                 [fjac, par, wa1, wa2] = self.lmpar(fjac, ipvt, diag, qtf,
-> 1167                                                      delta, wa1, wa2, par=par)
   1168                 # Store the direction p and x+p. Calculate the norm of p
   1169                 wa1 = -wa1

/Users/jgirard/anaconda/lib/python2.7/site-packages/agpy/mpfit/mpfit.py in lmpar(self, r, ipvt, diag, qtb, delta, x, sdiag, par)
   2153             temp = numpy.sqrt(par)
   2154             wa1 = temp * diag
-> 2155             [r, x, sdiag] = self.qrsolv(r, ipvt, wa1, qtb, sdiag)
   2156             wa2 = diag*x
   2157             dxnorm = self.enorm(wa2)

/Users/jgirard/anaconda/lib/python2.7/site-packages/agpy/mpfit/mpfit.py in qrsolv(self, r, ipvt, diag, qtb, sdiag)
   1974         #x.flags.writeable=True
   1975         # END JGMOD
-> 1976         x[ipvt] = wa
   1977         return (r, x, sdiag)
   1978 

ValueError: assignment destination is read-only