o i@s@ddlZddlZdddZddZd d Zd d ZdddZdS)NMbP?Mb`?PcCsdd}dtt|dd|tt|dd|d}tt|ddd|}||||}||||fS)NcsfddS)Ncs,dd|ddS)N)logsqrt)sigmabeta_dbeta_min4/mnt/c/Users/fbmor/ComfyUI/comfy/k_diffusion/deis.pys,z)edm2t....r r r r rrszedm2t..rrg?)nprtorchtensorcpuclonedetach) edm_steps epsilon_s sigma_min sigma_max vp_sigma_inv vp_beta_d vp_beta_mint_stepsr r redm2t s H$2rcCsFd}t|jdD]}||krq ||||||||9}q |S)Nrr)rangeshape)prev_tjtauspolykr r rcal_polys "r&cCs"td|d||||S)Nr)rexp)beta_0beta_1tr r r t2alpha_fn!s"r,c Cstjdd@|}|}|}t |dt|||}|}||j }Wdn1s:wYWdn1sIwYd|t |d|}|S)NF)modeTr'r) rinference_moder enable_gradrequires_grad_r,rsumbackwardgradr)r)r*r#alpha log_alphad_log_alpha_dtau integrandr r rcal_intergrand&s     r8'tabc s|dkrvt|\}}}g}tt|dd|ddD]W\\}}td|} | dkr3|gqt|||} |||} |fddt| D} g} t||| }t| D]}t | || }| t ||| qY|| q|S|dkrcdd }d d }g}tt|dd|ddD]\\}}t|} | d kr|gq|fd dt| dD} | dkr|| dd|| ddd|| d}||dd| d|}||g} nv| dkr|| d| d|||}||| d||| d}||| d||| d}|||g} nG| dkr]|| d| d| d|||}||| d| d||| d}||| d| d||| d}||| d| d||| d}||||g} || q|S)a Get the coefficient list for DEIS sampling. Args: t_steps: A pytorch tensor. The time steps for sampling. max_order: A `int`. Maximum order of the solver. 1 <= max_order <= 4 N: A `int`. Use how many points to perform the numerical integration when deis_mode=='tab'. deis_mode: A `str`. Select between 'tab' and 'rhoab'. Type of DEIS. Returns: A pytorch tensor. A batch of generated samples or sampling trajectories if return_inters=True. r:Nrcg|]}|qSr r .0r%ir r Lz'get_deis_coeff_list..rhoabcSsT|d|dd|d|d||d||||}|||||S)Nrr )abstartendccoeffr r rget_def_intergral_2Vs@z0get_deis_coeff_list..get_def_intergral_2cSs|d|dd|d|d|||d|d|d||||||d|||||}|||||||S)NrDrr )rErFrIrGrHdrJr r rget_def_intergral_3[s 2*z0get_deis_coeff_list..get_def_intergral_3rcr<r r r=r?r rrAfrBrrD) r enumeratezipminappendrlinspacerr8r&r1)r max_orderN deis_moder)r*Ct_curt_nextorderr#dtaur! coeff_tempr7r"r$rKrN coeff_cur coeff_prev1 coeff_prev2 coeff_prev3r r?rget_deis_coeff_list6sV *      % *  0      ra)rrr)r9r:)rnumpyrrr&r,r8rar r r rs