
    M/Phz                        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j        d
            ej        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    r> ej'         ej(        d e)e                    e            ed ee                      ed            ed           g Z e	j*        eeeeed
          \  ZZ Z+Z,Z-Z%Z&e%dk    rO ej.                      ej'         ej(        d e)e                    e            ed ee                     dS dS )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                 >   t          | j        t          | z            t          | j        t          | z            z  }t          | z  |t          z  | z  z
  t	          j        t          | z            z  }t                              t	          j        |                     |S )N)r   TABr   normdataappend)vrqress      W/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/sandbox/mle.pyRr      sq    	QS1c!#acll	"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      s:    !AaDDF1Hc!#acll*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      sd    !A$$q&qs2a55)))%1ac*:*::;CAaCLL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      sQ   
     % % % % % % % % % % & & & & & & & & & & % % % % % %             BIj GAJ	E!AJJ AFKNN o   1S{!LLL >$dJ
Hgq==CL1SSYY''---	E
"11Q44((( n    n   
 !CrSTUUU 2$'Q;;CJLLLCL1SSYY''---	E
!!!A$$''''' ;r   