U
    •Ê¹j„R  ã                   @   sX  d dl Z d dl mZ ddlmZmZmZmZmZmZm	Z	m
Z
mZ d dlmZmZ ddgZG dd„ deƒZd	d
e
› de› de	› d� e_dee ee ee ee ee ee eee eeeeeeeedœdd„Zee ee ee ee ee ee eeeeeeeeedœdd„Zee ee ee ee ee ee eeeeeeeeedœdd„ZdS )é    N)ÚTensoré   )	Ú	OptimizerÚ_use_grad_for_differentiableÚ
_get_valueÚ_dispatch_sqrtÚ_stack_if_compilingÚ_capturable_docÚ_differentiable_docÚ_foreach_docÚ_default_to_fused_or_foreach)ÚListÚOptionalÚNAdamÚnadamc                       sZ   e Zd Zdddddœeee eed	œ‡ fd
d„Z‡ fdd„Zdd„ Zeddd„ƒZ	‡  Z
S )r   çü©ñÒMb`?©gÍÌÌÌÌÌì?g+‡ÙÎ÷ï?ç:Œ0âŽyE>r   çü©ñÒMbp?FN)ÚforeachÚ
capturableÚdifferentiable)Údecoupled_weight_decayr   r   r   c                   sÜ   d|kst d|› �ƒ‚d|ks,t d|› �ƒ‚d|d   krDdk sXn t d|d › �ƒ‚d|d   krpdk s„n t d|d › �ƒ‚d|ksšt d	|› �ƒ‚d|ks°t d
|› �ƒ‚t||||||||	|
d�	}tƒ  ||¡ d S )Nç        zInvalid learning rate: zInvalid epsilon value: r   ç      ð?z#Invalid beta parameter at index 0: r   z#Invalid beta parameter at index 1: zInvalid weight_decay value: zInvalid momentum_decay value: )	ÚlrÚbetasÚepsÚweight_decayÚmomentum_decayr   r   r   r   )Ú
ValueErrorÚdictÚsuperÚ__init__)ÚselfÚparamsr   r   r   r   r   r   r   r   r   Údefaults©Ú	__class__© úP/var/www/html/bwave/python/venv/lib/python3.8/site-packages/torch/optim/nadam.pyr#   
   s*       ýzNAdam.__init__c                    sÚ   t ƒ  |¡ | jD ]4}| dd ¡ | dd¡ | dd¡ | dd¡ qt| j ¡ ƒ}t|ƒdkort 	|d d ¡}|s˜|D ]}t 
t|d ƒ¡|d< q|t|ƒdko´t 	|d d ¡}|sÖ|D ]}t 
|d ¡|d< q¾d S )	Nr   r   Fr   r   r   ÚstepÚ
mu_product)r"   Ú__setstate__Úparam_groupsÚ
setdefaultÚlistÚstateÚvaluesÚlenÚtorchZ	is_tensorÚtensorÚfloat)r$   r1   ÚgroupZstate_valuesZstep_is_tensorÚsZmu_product_is_tensorr'   r)   r*   r-       s    
zNAdam.__setstate__c           
      C   s  |d D � ]þ}|j d k	r| |¡ |j jr2tdƒ‚| |j ¡ | j| }	t|	ƒdkrÐ|d rptjdtj|j	d�nt 
d¡|	d< |d rštjdtj|j	d�nt 
d	¡|	d
< tj|tjd�|	d< tj|tjd�|	d< | |	d ¡ | |	d ¡ | |	d
 ¡ | |	d ¡ qd S )Nr%   z'NAdam does not support sparse gradientsr   r   r)   )ZdtypeÚdevicer   r+   r   r,   )Zmemory_formatÚexp_avgÚ
exp_avg_sq)ÚgradÚappendZ	is_sparseÚRuntimeErrorr1   r3   r4   Zzerosr6   r9   r5   ZonesZ
zeros_likeZpreserve_format)
r$   r7   Úparams_with_gradÚgradsÚexp_avgsÚexp_avg_sqsÚmu_productsÚstate_stepsÚpr1   r)   r)   r*   Ú_init_group1   s,    


ÿþÿþzNAdam._init_groupc                 C   sÀ   |   ¡  d}|dk	r.t ¡ � |ƒ }W 5 Q R X | jD ]†}g }g }g }g }g }g }	|d \}
}|  |||||||	¡ t||||||	|
||d |d |d |d |d |d |d	 |d
 d� q4|S )z±Performs a single optimization step.

        Args:
            closure (Callable, optional): A closure that reevaluates the model
                and returns the loss.
        Nr   r   r   r   r   r   r   r   r   )
Úbeta1Úbeta2r   r   r   r   r   r   r   r   )Z _cuda_graph_capture_health_checkr4   Zenable_gradr.   rF   r   )r$   ÚclosureZlossr7   r?   r@   rA   rB   rC   rD   rG   rH   r)   r)   r*   r+   Q   s@    

ñz
NAdam.step)r   r   r   r   r   F)N)Ú__name__Ú
__module__Ú__qualname__Úboolr   r#   r-   rF   r   r+   Ú__classcell__r)   r)   r'   r*   r   	   s$         ÿ ý ý a�  Implements NAdam algorithm.

    .. math::
       \begin{aligned}
            &\rule{110mm}{0.4pt}                                                                 \\
            &\textbf{input}      : \gamma_t \text{ (lr)}, \: \beta_1,\beta_2 \text{ (betas)},
                \: \theta_0 \text{ (params)}, \: f(\theta) \text{ (objective)}                   \\
            &\hspace{13mm} \: \lambda \text{ (weight decay)}, \:\psi \text{ (momentum decay)}    \\
            &\hspace{13mm} \: \textit{decoupled\_weight\_decay}                                  \\
            &\textbf{initialize} :  m_0 \leftarrow 0 \text{ ( first moment)},
                v_0 \leftarrow 0 \text{ ( second moment)}                                 \\[-1.ex]
            &\rule{110mm}{0.4pt}                                                                 \\
            &\textbf{for} \: t=1 \: \textbf{to} \: \ldots \: \textbf{do}                         \\
            &\hspace{5mm}g_t           \leftarrow   \nabla_{\theta} f_t (\theta_{t-1})           \\
            &\hspace{5mm} \theta_t \leftarrow \theta_{t-1}                                       \\
            &\hspace{5mm} \textbf{if} \: \lambda \neq 0                                          \\
            &\hspace{10mm}\textbf{if} \: \textit{decoupled\_weight\_decay}                       \\
            &\hspace{15mm} \theta_t \leftarrow \theta_{t-1} - \gamma \lambda \theta_{t-1}                    \\
            &\hspace{10mm}\textbf{else}                                                          \\
            &\hspace{15mm} g_t \leftarrow g_t + \lambda \theta_{t-1}                             \\
            &\hspace{5mm} \mu_t \leftarrow \beta_1 \big(1 - \frac{1}{2}  0.96^{t \psi} \big)     \\
            &\hspace{5mm} \mu_{t+1} \leftarrow \beta_1 \big(1 - \frac{1}{2} 0.96^{(t+1)\psi}\big)\\
            &\hspace{5mm}m_t           \leftarrow   \beta_1 m_{t-1} + (1 - \beta_1) g_t          \\
            &\hspace{5mm}v_t           \leftarrow   \beta_2 v_{t-1} + (1-\beta_2) g^2_t          \\
            &\hspace{5mm}\widehat{m_t} \leftarrow \mu_{t+1} m_t/(1-\prod_{i=1}^{t+1}\mu_i)\\[-1.ex]
            & \hspace{11mm} + (1-\mu_t) g_t /(1-\prod_{i=1}^{t} \mu_{i})                         \\
            &\hspace{5mm}\widehat{v_t} \leftarrow   v_t/\big(1-\beta_2^t \big)                   \\
            &\hspace{5mm}\theta_t \leftarrow \theta_t - \gamma \widehat{m_t}/
                \big(\sqrt{\widehat{v_t}} + \epsilon \big)                                       \\
            &\rule{110mm}{0.4pt}                                                          \\[-1.ex]
            &\bf{return} \:  \theta_t                                                     \\[-1.ex]
            &\rule{110mm}{0.4pt}                                                          \\[-1.ex]
       \end{aligned}

    For further details regarding the algorithm we refer to `Incorporating Nesterov Momentum into Adam`_.
    a  
    Args:
        params (iterable): iterable of parameters to optimize or dicts defining
            parameter groups
        lr (float, optional): learning rate (default: 2e-3)
        betas (Tuple[float, float], optional): coefficients used for computing
            running averages of gradient and its square (default: (0.9, 0.999))
        eps (float, optional): term added to the denominator to improve
            numerical stability (default: 1e-8)
        weight_decay (float, optional): weight decay (L2 penalty) (default: 0)
        momentum_decay (float, optional): momentum momentum_decay (default: 4e-3)
        decoupled_weight_decay (bool, optional): whether to use decoupled weight
            decay as in AdamW to obtain NAdamW (default: False)
        z	
        zÏ

    .. _Incorporating Nesterov Momentum into Adam:
        https://openreview.net/forum?id=OM0jvwB8jIp57ZJjtNEZ
    .. _Decoupled Weight Decay Regularization:
        https://arxiv.org/abs/1711.05101

    F)r%   r@   rA   rB   rC   rD   r   r   r   r   rG   rH   r   r   r   r   c
                C   s¦   t dd„ |D ƒƒstdƒ‚t dd„ |D ƒƒs4tdƒ‚|dkrNt| |	dd�\}}|rdtj ¡ rdtd	ƒ‚|rxtj ¡ sxt}nt}|| ||||||
||||||||	d
� dS )zpFunctional API that performs NAdam algorithm computation.

    See :class:`~torch.optim.NAdam` for details.
    c                 s   s   | ]}t |tjƒV  qd S ©N©Ú
isinstancer4   r   ©Ú.0Útr)   r)   r*   Ú	<genexpr>Ó   s     znadam.<locals>.<genexpr>zPAPI has changed, `state_steps` argument must contain a list of singleton tensorsc                 s   s   | ]}t |tjƒV  qd S rO   rP   rR   r)   r)   r*   rU   Ö   s     zPAPI has changed, `mu_products` argument must contain a list of singleton tensorsNF)Z	use_fusedz6torch.jit.script not supported with foreach optimizers)	rG   rH   r   r   r   r   r   r   r   )Úallr>   r   r4   ZjitZis_scriptingÚ_multi_tensor_nadamÚ_single_tensor_nadam)r%   r@   rA   rB   rC   rD   r   r   r   r   rG   rH   r   r   r   r   Ú_Úfuncr)   r)   r*   r   º   s6    ò)r%   r@   rA   rB   rC   rD   rG   rH   r   r   r   r   r   r   r   c       	         C   s  t | ƒD �]ô\}}|| }|| }|| }|| }|| }tj ¡ st|rt|jrZ|jrZ|jst|jrl|jrl|jsttdƒ‚|d7 }|r†|}nt|ƒ}d||  }|	dkrÈ|rº| d||	  ¡ n|j	||	d�}|ddd||
     }|ddd|d |
     }||9 }| 
|d| ¡ | |¡j||d| d� | |¡ ¡ }|�sH|�r¤| 	|¡}|| }|| d|  d|   }|| | d|   }| ||¡ | ||¡ qt|ƒ| }| |¡ |j||| d|  dt|ƒ  d� |j||| | d|  d� qd S )	NzUIf capturable=True, params, mu_products, and state_steps must be CUDA or XLA tensors.r   r   ©Úalphar   ç      à?ç¸…ëQ¸î?)Úvalue)Ú	enumerater4   Ú_utilsÚis_compilingÚis_cudaZis_xlaÚAssertionErrorr   Zmul_ÚaddZlerp_Zaddcmul_ÚdivÚsqrtZaddcdiv_Zadd_)r%   r@   rA   rB   rC   rD   rG   rH   r   r   r   r   r   r   r   ÚiÚparamr<   r:   r;   r,   Zstep_tr+   Zbias_correction2ÚmuÚmu_nextÚdenomZmu_product_nextr)   r)   r*   rX   õ   s^    ÿÿÿÿÿÿþ

&rX   c       	             s  t | ƒdkrd S |rtdƒ‚tj ¡ sL|rLtdd„ t| ||ƒD ƒƒsLtdƒ‚t | |||||g¡}| 	¡ D �]š\\}}}}}}}t 
|d¡ |	dkr¾|r®t |dˆ|	  ¡ ntj|||	d�}t ||dˆ  ¡ t |ˆ¡ t |||dˆ ¡ t |¡}|�r¨t |ˆ¡}t d|¡}t |d	¡ t 
|d
¡ t |ˆ ¡ t 
|ˆ¡ t d|¡}t |d	¡ t 
|d
¡ t |ˆ ¡ ~t ˆ|¡}t |d
¡ t |¡ t |¡ n:‡fdd„|D ƒ}‡ ‡fdd„|D ƒ}‡ ‡fdd„|D ƒ}t ||¡ t ||¡ t 
||¡ ~|�r®t |d
¡ t |ˆ¡ t |d
¡}t |¡ t ||¡ |}~t ||¡}t |ˆ¡ t |d
¡ t ||¡ |}~t ||¡}t |||¡ t |||¡ qjt‡fdd„t||ƒD ƒƒ}t‡fdd„t||ƒD ƒƒ}t ||||¡ t ||||¡ qjd S )Nr   z#_foreach ops don't support autogradc                 s   s&   | ]\}}}|j o|j o|j V  qd S rO   )rc   )rS   rE   Úmpr+   r)   r)   r*   rU   Y  s   ÿz&_multi_tensor_nadam.<locals>.<genexpr>zNIf capturable=True, params, mu_products, and state_steps must be CUDA tensors.r   r[   r^   g      à¿r   c                    s    g | ]}t d ˆ t|ƒ  ƒ‘qS )r   )r   r   ©rS   r+   )rH   r)   r*   Ú
<listcomp>Œ  s     z'_multi_tensor_nadam.<locals>.<listcomp>c                    s(   g | ] }ˆ d ddt |ƒˆ     ‘qS )r   r]   r^   ©r   rn   ©rG   r   r)   r*   ro   �  s     c                    s,   g | ]$}ˆ d ddt |ƒd ˆ     ‘qS )r   r]   r^   r   rp   rn   rq   r)   r*   ro   Ž  s   ÿc                    s,   g | ]$\}}ˆ d |  d t |ƒ  d ‘qS ©r   éÿÿÿÿrp   )rS   r,   rj   ©r   r)   r*   ro   »  s   ÿc                    s,   g | ]$\}}ˆ | d t |ƒ|   d ‘qS rr   rp   )rS   r,   rk   rt   r)   r*   ro   ½  s   ÿ)r3   rd   r4   ra   rb   rV   Úzipr   Z"_group_tensors_by_device_and_dtyper2   Z_foreach_add_Z_foreach_mul_Z_foreach_addZ_foreach_lerp_Z_foreach_addcmul_Z_foreach_sqrtZ_foreach_mulZ_foreach_powZ_foreach_sub_Z_foreach_neg_Z_foreach_sqrt_Z_foreach_div_Z_foreach_subZ_foreach_addcdiv_r   ) r%   r@   rA   rB   rC   rD   rG   rH   r   r   r   r   r   r   r   Zgrouped_tensorsZgrouped_paramsZgrouped_gradsZgrouped_exp_avgsZgrouped_exp_avg_sqsZgrouped_mu_productsZgrouped_state_stepsrY   Zexp_avg_sq_sqrtÚexponentZmusZmu_nextsZbias_correction_sqrtrl   Zstep_size_gradsZstep_size_expavgÚ	numeratorr)   )rG   rH   r   r   r*   rW   A  s–    
ÿ
þÿ   

ÿ
ÿÿrW   )FNFF)r4   r   Z	optimizerr   r   r   r   r   r	   r
   r   r   Útypingr   r   Ú__all__r   Ú__doc__rM   r6   r   rX   rW   r)   r)   r)   r*   Ú<module>   s†   ,u#óòñÝD    õî;ñLñ