
    !gz                        d Z ddlZddlmZmZmZ ddlmZmZm	Z	 ddl
mZ ddlmZ d Zd Zd Z ej&                  d	      Zej*                  d   Z eee      Z ej0                  d
        ej2                  e      Z ed       d
Zg Z e	j<                  eeeed
      \  ZZ Z!Z"Z#Z$Z%Z&e%dk(  r9 ejN                   ejP                  d e)e            e        ed ee              ed        ed       g Z e	jT                  eeeeed
      \  ZZ Z+Z,Z-Z%Z&e%dk(  rK ej\                           ejN                   ejP                  d e)e            e        ed ee             yy)zWhat's the origin of this file? It is not ours.
Does not run because of missing mtx files, now included

changes: JP corrections to imports so it runs, comment out print
    N)dotouterrandom)iolinalgoptimize)eyec                 2   t        | j                  t        | z        t        | j                  t        | z        z  }t        | z  |t        z  | z  z
  t	        j
                  t        | z        z  }t        j                  t	        j
                  |             |S )N)r   TABr   normdataappend)vrqress      T/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/sandbox/mle.pyRr      si    	QSS1c!##acl	"BQ3r!tAv:v{{1Q3'
'CKKC !I    c                     dt         | z  t        |       t        z  | z  z
  z  t        | j                  t        | z        z  }|S )z
 Gradient    )r   r   r   r   r   r   results     r   Rpr      s8    !AaDF1Hc!##acl*FMr   c                     dt         t        |       t        z  z
  t        t        | z  t	        |             z
  t        t	        |       t        | z        z
  z  t        | j                  t        | z        z  }|S )z	 Hessian r   )r   r   r   r   r   r   r   r   s     r   Rppr      sX    !A$q&qs2a5))%1ac*::;CAaCLHFMr   znos4.mtx   ztry fmin_bfgs)fprimefull_outputretallzRayleigh quotient BFGSzfmin_bfgs OKztry fmin_ncg)r   fhessr    r!   zRayleigh quotient NCG)/__doc__numpynpr   r   r   scipyr   r   r   scipy.sparser	   speyematplotlib.pyplotpyplotpltr   r   r   mmreadr   shapenr   seedrandv_0printr    r   	fmin_bfgsr   foptgoptHopt
func_calls
grad_callswarnflagallvecssemilogyarangelenfmin_ncgfcallsgcallshcallsfigure r   r   <module>rD      s|  
  % % & & %  BIIj GGAJ	!AJ AFKKN o 1S{!L >$dJ
Hgq=CLL1SY'-	
"1Q4( n  n 
 !CrSTU 2$'Q;CJJLCLL1SY'-	
!!A$' r   