OneScience-Group/DiffDock
053
1o
2 �aOj4 � @ s� d dl mZ d dlmZ d dlmZ d dlZd dlZd dlm Z d dl3m Z m Z
d dlmZ dd lmZ dd4lmZ dd lmZ d Zd
d� Zdd� Zdd� Zd dd�Zd dd�Z d dd�Z d dd �Z dS ) � )� Namespace)�partial)�PathN)� DataParallel)�get_timestep_embedding�5t_to_sigma)�ExponentialMovingAverage� )�AAModel)�CGModel)�6AAOldModel)�moad_esm_embeddings_path�pdbbind_esm_embeddings_path� pdbsidechain_esm_embeddings_path�esm_embeddings_path�esm_embeddings_modelc C sZ t | �} | d }|jddd��}t�|�pi }W d � n1 s!w Y tdi |��S )Nzmodel_parameters.yml�rzutf-8)�encoding� )r �open�yamlZ full_loadr )� model_dirZ config_pathZhandleZconfigr r �+modelscope/diffdock/models/score_wrapper.py�load_model_args s �r c s t � fdd�tD ��S )Nc 3 s � | ] }t � |d �d uV qd S �N)�getattr)�.0�key��7model_argsr r � <genexpr># s � z+model_uses_lm_embeddings.<locals>.<genexpr>)�any�_LM_EMBEDDING_KEYSr r r r �model_uses_lm_embeddings" s r# c C s( z|| v W S t y t| |� Y S w r )� TypeError�hasattr)�args�namer r r �_has_arg& s8 9 �r( c C s. t | t�r | �||�S t| |�rt| |�S |S r )�10isinstance�dict�getr( r )r& r' �defaultr r r �_get_arg- s11 12 1314r- Fc
C s| t t| dd�| jt| dd�d�}t| dd�}|r |s td��d }t| d �d us<t| d15�d us<t| d �d us<t| d �d ur>d
}t| d�d urH| j}|rMt} n|rRt} nt} tdLi d|�d|�d| j �d|�d| j16�d| j �d| j �d| j�dt| dd��d| j
�d| j�d| j�d| j�d | j �d | j�d | j�d | j�d | j�d!t| d!d��d"t| d"d��d#|�d$|�d%t| d%d&��d't| d'd��d(t| d(d��d)t| d)d*��d+tt| d,�t�r�t| j�d* nd*�d-tt| d.�t�r�t| j�d* nd*�d/t| d/d0��d1t| d2d� �d3t| d3d��d4t| d4d��d5t| d5d6��d7t| d8d� �d9t| d9d6��d:t| d:d;��d<t| d<d��d=t| d=d��d>t| d?d&�d&k�d@t| dA��rR| j d;k�p\t| dB��o\| j d;k�dCt| dCd���}17| tu �rtt| dDd �|18dD< n'| tu �r�dED ] } |19� | d � �q{t| d �d u�r�dFnd |20d#< t| dGdH�|21dG< | dLi |22��} |j dIk�r�|�s�t| dJ�dKk�r�t | �} | �!|� | S )MN�embedding_typeZ23sinusoidal�embedding_scalei' )r. Z
embedding_dimr/ � all_atomsFz�The old coarse-grained DiffDock model path is not migrated yet. Use old=True only with all_atoms=True, or migrate old_cg_model.py first.r
r r r Z precomputedr r �device�24no_torsion�timestep_emb_func�num_conv_layersZlig_max_radius�scale_by_sigma�sigma_embed_dimZ
norm_by_sigma�ns�nv�distance_embed_dim�cross_distance_embed_dimZ25batch_norm�dropout�use_second_order_repr�cross_max_distance�dynamic_max_crossZ smooth_edgesZ26odd_parity�lm_embedding_type�confidence_modeZconfidence_dropoutg Zconfidence_no_batchnormZaffinity_predictionZparallelr Znum_confidence_outputs�rmsd_classification_cutoff�atom_num_confidence_outputs� atom_rmsd_classification_cutoffZparallel_aggregators� Zfixed_center_convZnot_fixed_center_convZno_aminoacid_identitiesZinclude_miscellaneous_atomsZsh_lmax� �differentiate_convolutionsZ no_differentiate_convolutions�tp_weights_layers�num_prot_emb_layersr �reduce_pseudoscalars�embed_also_ligand�atom_confidenceZatom_confidence_loss_weight�sidechain_pred�sidechain_loss_weight�backbone_loss_weight�depthwise_convolutionZ crop_beyond) rB rF rG rH rI rJ rK rL rO ZesmZuse_old_atom_encoderT�cudaZdatasetZ torsionalr )"r r- r6 �NotImplementedErrorr r r27 r r* r2 r4 Z28max_radiusr5 r7 r8 r9 r: Z
no_batch_normr; r< r= r> r) �list�lenrA rC r( rM rN �pop�typer �to)
r& r1 r � no_parallelr@ �oldr3 r0 r? Z model_classZ model_kwargsr �modelr r r � get_model5 s� 2930� ��������� �31� � �
������ � ��� � � � � ���"��� %�&� '� (� )�*� +� ,� -� .�/�1�� 4�32633 �"34rZ c C s( t t| d�}t| |||||d�}||fS )N)r& )r rW r@ rX )r �t_to_sigma_complrZ )r r1 rW r@ rX r rY r r r �build_score_model� s �r\ Tc
C s |d u rt �t j�� r
dnd�}t| �}t|||||d�\}} t| �| }35t j|36t �d�d�} t| t �rkd| v rkd| v rk|j37| d |d� d| v rjt |d d �d urjt |�
� |jd38�} | j39| d |d � | �|�
� � n|j40| |d� |�|�}|�� ||| fS ) NrP Zcpu)r r1 rW r@ rX )Z map_locationrY Z optimizer)�strict� ema_weights�ema_rate)Zdecay)r1 )�torchr1 rP Z is_availabler r\ r �loadr) r* Zload_state_dictr r Z41parametersr_ Zcopy_torV �eval)
r Zckptr1 rW r@ rX r] r rY r Zcheckpoint_pathZ42state_dictr^ r r r �load_score_model� s. 43� �4445rc r )FFF)NTFFT) �argparser � functoolsr Zpathlibr r` r Z torch_geometric.nn.data_parallelr Z)onescience.utils.diffdock.diffusion_utilsr r r[ Z onescience.utils.diffdock.utilsr Zaa_modelr46 Zcg_modelr Z old_aa_modelr r" r r# r( r- rZ r\ rc r r r r �<module> s6 4748u49��