o *j@sUddlZddlmZmZddlmZddlmZmZm Z ddl Z ddl Z ddl mZddlmZmZmZmZddlmZmZmZmZddlmZdd lmZmZeee j fZ!eee j fZ"e"e#d <gd Z$d d Z%ddZ&ddZ'ddZ(ddZ)ddZ*e(ddZ+ddZ,ddZ-ddZ.d d!Z/dd"d#Z0e( $ $dd%d&Z1e(d'd(Z2e(d)d*Z3e(d+d,Z4d-d.Z5e(d/d0Z6e(dd2d3Z7d4d5Z8d6d7Z9e(dd8d9Z:d:d;Z;dddDdEZ?dFdGZ@dHdIZAdJdKZBdLdMZCe( $ $ddNdOZDe( ddPdQZEdRdSZFdTdUZGdVdWZHe(dXdYZIe(dZd[ZJd\d]ZKd^d_ZLd`daZMddcddZNe( >ddedfZOdgdhZPdidjZQe( $ddkdlZRe(dmdnZSe(dodpZTdqe!dreUdse eUdte jVfdudvZW  $ddwe jVdreUdse eUdxeXdte jVf dydzZY  $dd e!dreUdse eUdxeXdte jVf d{d|ZZ }dd e!dreUdseUd~e[de[dteeUf ddZ\de jVdte jVfddZ]de jVde jVde jVdte jVfddZ^ddZ_dS)N)reducewraps)add)ListMutableMappingOptional) data_utils)batched_gatherone_hottensor_tree_maptree_map) N_EXTRA_MSAN_MSAN_RESN_TPL)residue_constants)FrameRotationprotein)msadeletion_matrixmsa_mask msa_row_mask bert_masktrue_msa msa_chainscCsX|D]%\}}|dr|tj||<q|jtjtjtjfvr)|tj ||<q|S)NZ_mask) itemsendswithtypetorchfloat32dtypeint32Zuint8int8int64)rkvr'p/var/www/html/Deteccion_Ine/venv/lib/python3.10/site-packages/modelscope/models/science/unifold/data/data_ops.pycast_to_64bit_ints#s r)cCstj|djtjd|d<|S)Naatyper!Zseq_maskronesshaper rr'r'r( make_seq_mask.s  r0cCs"tj|djdtjd|d<|S)Ntemplate_aatyperr+ template_maskr,r/r'r'r(make_template_mask4s r3cstfdd}|S)z#Supply all arguments but the first.csfddS)Ncs|gRiSNr'x)argsfkwargsr'r(?sz$curry1..fc..r'r7r9r8r;r(fc=szcurry1..fc)r)r8r=r'r<r(curry1:sr>cCs\|d|d<tj}tj|tjddd|djd}t |d|d|d<|S)z1Correct MSA restype to have the same order as rc.rr+r) longrcZ MAP_HHBLITS_AATYPE_TO_OUR_AATYPErtensorr# unsqueezeexpandr.Zgather)rZnew_order_listZ new_orderr'r'r(correct_msa_restypesDsrFcCst|djdkrtj|ddd|d<d|vr)t|djdkr)|dd|d<dD];}||vrft||jrf||jd}t|trf|dkrft||r[tj||dd||<q+tj||dd ||<q+d D]}||vr~t||jr~||d||<qi|S) z=Remove singleton and repeated dimensions in protein features.r*r?dim resolutionr@r) Z domain_namernum_alignments seq_lengthsequenceZ superfamilyrbetween_segment_residuesZ residue_indextemplate_all_atom_mask)Zaxis)rLrK) lenr.rargmax isinstanceintZ is_tensorZsqueezenp)rr%Z final_dimr'r'r(squeeze_featuresPs" rUcCs|dkrLtj|dj|k}d}d}t||d|k}t|t|d||d|d<tj|dj|k}t|t|d||d|d<|S)z&Replace a portion of the MSA with 'X'.rr*)rTrandomrandr.r logical_andwhereZ ones_like)rZreplace_proportionrZx_idxZgap_idxZ aatype_maskr'r'r(!randomly_replace_msa_with_unknownrs(r]cCs:d}ttjdd|}tt|| | }|S)zGenerate Gumbel Noise of given Shape. This generates samples from Gumbel(0, 1). Args: shape: Shape of noise to return. Returns: Gumbel noise of given shape. ư>rr@)r from_numpyrTrYuniformlog)r.epsilonZ uniform_noiseZgumbelr'r'r( gumbel_noisesrccCst|j}tj||ddS)a;Samples from a probability distribution given by 'logits'. This uses Gumbel-max trick to implement the sampling in an efficient manner. Args: logits: Logarithm of probabilities to sample from, probabilities can be unnormalized. Returns: Sample from logprobs in one-hot form. r?rH)rcr.rrQlogitszr'r'r(gumbel_max_samples rgcCst|j}tj||dddS)aSamples with replacement from a distribution given by 'logits'. This uses Gumbel trick to implement the sampling an efficient manner. For a distribution over k items this samples k times without replacement, so this is effectively sampling a random permutation with probabilities over the permutations derived from the logprobs. Args: logits: Logarithm of probabilities to sample from, probabilities can be unnormalized. Returns: Sample from logprobs in index r?T)rIZ descending)rcr.rargsortrdr'r'r(gumbel_argsort_sample_idxs ricCs4ttj|dd}tjtdg|fddS)Nr@rrH)rr_rTrY permutationcatrC)num_seqshuffledr'r'r(uniform_permutationsrncCsbtj|dddk}tj|tjd}d||<t|jdks!J|dd}|dd}|jddkr:tdgS|dur|ddd}d||<t j |dd\}}|}|dk}||}|dk} d |d ||<d||<|D]} | dkr|| k} | } | dkr|| || | 9<qxt |d }t |d} tj tdg| fddS) Nr?rHrr+g.r@TZ return_counts?r^)rsumrA zeros_liker rPr.rCreshaperTuniquerarirk)rrZhas_msarekeyscountsZ num_has_msaZnum_pairZ num_unpairZ num_chainsr%Zcur_maskZcur_cntrmr'r'r(gumbel_permutations8         rwFc Cs|djd}t||}|st|}n|rd|vr|dnd}t|d|}t||}t||||g\} } tD] } | |vrY|rNt|| d| |d| <t|| d| || <q9|S)zLSample MSA randomly, remaining sequences are stored are stored as `extra_*`.rrrNrextra_)r.minrnrwrsplitMSA_FEATURE_NAMES index_select) rmax_seq keep_extra gumbel_sampleZbiased_msa_by_chainrlnum_selZ index_orderrZsel_seqZ not_sel_seqr%r'r'r( sample_msas0      rcCs(d|vr|ddkrt|dd|}|S)Nis_distillationr@F)r~)r)rr}r'r'r(sample_msa_distillationsrcCs|djd}|djd}|j|}||krSttjj|d|dddd}t|d}tj t dg|fdd}t D]}||vrRt ||d|||<qA|S)Nrrr@FreplacerH) r.Z max_msa_entryrr_rTrYchoicerAsortrkrCr{r|)rconfigrlZseq_lenr}Z keep_indexr%r'r'r(random_delete_msas& rcCsj|djd}t||}ttj|d|}tD]}d||vr2t|d|d||d|<q|S)N extra_msarrx) r.ryrr_rTrYrjr{r|)rZ max_extra_msarlrZselect_indicesr%r'r'r(crop_extra_msas   rcCs$tD] }d||vr|d|=q|S)Nrx)r{)rr%r'r'r(delete_extra_msas   rcCsd|vr |ddkr |S|djd}||jkr|Sttj|tjd|jtj}|j r:t j d|j d}n|j }tt j d||g}|dddftd|}t|d|d}t|d}ttd|d|dtddf}|jdd \} } | | dk} | d} ttdd| dgd} t| ddksJtD]} | |vrtj|| d| d || <q|S) Nrr@rrr+r?rGTro)index)r.Z min_num_msarfloorrCr Zmsa_fraction_per_blocktor"Zrandomize_num_blocksrTrYrandintZ num_blocksr_arangecliprtviewZhstackzerosrArSr{r|)rrrlZ block_num_seqnbZdel_block_startsZ del_blocksZ del_indicescombinedZuniquesrv differenceZ keep_indicesr%r'r'r(block_delete_msa&sL    rrVcCsttd|tdtdgd}t|dd}|ddddddf|}t|dd}|ddddddf|}|j\}}} |j\} } } || |d} ||||ddd} | | } tj| dd  |d <|S) NrXr@rrrrextra_msa_maskrHextra_cluster_assignment) rrkr-rr r.rZ transposerQrA)rgap_agreement_weightweights msa_one_hotZsample_one_hotZextra_msa_one_hot extra_one_hotrlnum_res_Z extra_num_seqab agreementr'r'r(nearest_neighbor_clustersJs"  rcCst|jdkr|jd|jdksJ|j|jdgdt|jddR}||j}|gt|jdd}tj|d||}| |j }|S)Nr@rr@) rPr.rrElistrrZ scatter_add_floatrr!)dataZ segment_idsZ num_segmentsr.rCr'r'r(unsorted_segment_sumes  rcsdjdfdd}d}dd||}||dddddftd d }|tdd 7}||dddddfd <~||d }|d 7}||d<~S)z=Produce profile and deletion_matrix_mean within each cluster.rrcst|dS)Nr)rr5rlrr'r(csumus z summarize_clusters..csumrr^rNrrcluster_profileextra_deletion_matrixrcluster_deletion_mean)r.r )rrmaskZ mask_countsmsa_sumdel_sumr'rr(summarize_clustersqs(  rcCstjdgd|gdgtjd}|d}|d}t|dd}t|d d}|d d d d d f|}|d d d d d f|}||} | | jd | jd | jd } ||jd |jd |jd } | | j} tjjj d| d d} | t d||9} tj | dd} | d7} t d| |}||7}|| d d d d f}|d}|d}t d| ||}||7}|| d d d f}||d<||d<|S)zFAssign each extra MSA sequence to its nearest neighbor in sampled MSA.rprXrVr+rrrrrNrr@rGg@@rHz mr, nr->mnr?z nm, mrc->nrcrrz nm, mc->ncrr) rrCr r rr.Tnn functionalZsoftmaxeinsumrq)batchrrrZ extra_maskrrZmsa_one_hot_maskedZextra_one_hot_maskedt1t2rZcluster_assignmentZ cluster_countrrrrrrr'r'r(nearest_neighbor_clusters_v2sF " rcCsDd|vrtj|djtjd|d<tj|djdtjd|d<|S)z:Mask features are all ones, but will later be zero-padded.rrr+rrr,r/r'r'r( make_msa_masks   rc Cs|jddkr?t|tjd}tjd}tjd}tt|ddgt|jdg|d|d d f|d|d d f}n |j g|jdR}|d urs|jddkrht||d|f|d|f}||fSt | }||fS|S) zCreate pseudo beta features.rGCAZCB.Nr@.N) r.reqrB restype_order atom_orderr\tilerP new_zerosrrr)r*all_atom_positions all_atom_maskZis_glyca_idxZcb_idx pseudo_betapseudo_beta_maskr'r'r(pseudo_beta_fns&    rcCsP|dvsJt||r dnd||d||rdnd\||d<||d<|S) z9Create pseudo-beta (alpha for glycine) position and mask.)rZ template_r1r*rrOrrr)r)rprefixr'r'r(make_pseudo_betas    rcCst|||<|Sr4rrC)rkeyvaluer'r'r(add_constant_fieldsr绽|=cCs`|j}|d}t||d|g}t}|tjdtj|d|d}t||ddS)Nr?ir@) generator) r.rrs GeneratorZ manual_seedrTrYr multinomial)ZprobsrbZdsZ num_classesgenrvr'r'r(shaped_categoricalsrcCs0d|vr|St|dd}tj|dd|d<|S)7Compute the HHblits MSA profile if not already present.hhblits_profilerrrH)r rmean)rrr'r'r(make_hhblits_profiles rcCsLt|dd}|ddddddf}||9}|jdd|jdddS)zCompute the MSA profile.rrrNrrHr)r rq)rZohrr'r'r(make_msa_profile srcCsd|vr|St||d<|S)rr)rr/r'r'r(make_hhblits_profile_v2s rc Csd|vr|S|d}|d}|d}|}|dk}|D]*}||k}t||d} | dkrF||@} |dd| f} | d| |dd|f<q|S)Nnum_sym entity_idsym_idr@r)rtrSrepeat) mask_positionrrrrZunique_entity_idsZfirst_sym_maskZ cur_entity_idZcur_entity_maskZ cur_num_symZ cur_sym_maskZcur_sym_bert_maskr'r'r(share_mask_by_entitys$rc Csntjdgdddgtjd}|j||j|d|jt|dd}ttt dd t t |j D}d |d <d |j|j|j}|dksKJtj jj|||d }|dj } ttjj| |k} | |d M} d|vry| |dM} |rt| |} |rt|d} t| } nt|} t| | |d} | |d 9} | tj|d<|d|d<| |d<|S)z Create data for BERT on raw MSA.g?rWrVr+rrrcSg|]}dqS)rrr'.0rr'r'r( @z#make_masked_msa..r@rp)rrrr^r)rrCr Z uniform_probZ profile_probZ same_probr rrrrangerPr.rrpadr_rTrYrZboolrrargrr\rAr) rrZreplace_fractionrZ share_maskZ random_aaZcategorical_probsZ pad_shapesZ mask_probshrreZbert_msar'r'r(make_masked_msa/sB      rc sXddd}|dur$d|vr|djdn|djd}||kr$||d}d |vr8|d jd}||kr8||d }t|t|t|t|i|D]c\} | d krOqFtj} || } d } t| t| ksqJ| d | d| d| fddt| | D} fddt | D}| tt j |}|rt jj||| <t || | || <qF|S)z;Guess at the MSA and sequence dimension to make fixed size.cSst|||d||S)Nr@max)Zcur_sizeZ multiplierr'r'r( get_pad_sizejsz%make_fixed_size..get_pad_sizeNr*rrr@rrz0Rank mismatch between shape and shape schema for z: z vs cs g|] \}}|dp |qSr4)get)rs1s2) pad_size_mapr'r(rsz#make_fixed_size..cs"g|] \}}d|j|fqSrr.)rip)r&r'r(rs")r)r.rrr rrrrPzip enumeratereverse itertoolschainrrrrrs)r shape_schemaZmsa_cluster_sizeZextra_msa_sizer num_templatesrZ input_num_resZinput_extra_msa_sizer%r.ZschemamsgZpad_sizepaddingr')rr&r(make_fixed_size_sF     . rcCs|d|d<d|vrt|dtjdd}n%tj|dtjd}d|vr;|d}tj|ddd d}d ||<|}t|dd }tj |dd|g}tj |dd|d <|S) $Create and concatenate MSA features.r*rNrr@r+asym_lenr?rHNrprX target_feat) rArrrr rrZcumsumrr rDrk)rZ has_breakr Z entity_endsZ aatype_1hotr r'r'r(make_target_feats" r cCst|dd}t|ddd}t|dddtj}|tj|dd tj|dd g}d |vrMt|d ddtj}||d tj|dd gd |vrlt|d dd|d <t|d ddtj|d<tj|dd |d<|S)r rrrrVrp@@r?rHrrrextra_msa_has_deletionextra_msa_deletion_valuemsa_feat) r rratanrTpirDextendrk)rmsa_1hot has_deletiondeletion_valuerdeletion_mean_valuer'r'r( make_msa_feats@       rcCst|dd}|d}t|ddd}t|ddtjd}t|d ddtjd}||||d |g}tj|d d |d <|S)r rrrrVrprr rrrr?rHr)r rrrrTrZarctanrk)rrrrrrrr'r'r(make_msa_feat_v2s"rcCsz|dd|}|dd|}t|dd}t|ddtj}|dd|}||d<||d<||d<||d <|S) NrrrVrpr rrrr)rrrrTr)rZ num_extra_msarrrrrr'r'r(make_extra_msa_featsrcsfdd|DS)Ncsi|] \}}|vr||qSr'r'rr%r& feature_listr'r( szselect_feat..)r)rrr'rr( select_featsr cCsd|vr|Stjtjtj|djd}tjtjtj|djd}tjtjtj|djd}tjtj tj|djd}|d }|| |d<|| |d<|||d<|||d<|S)z>Construct denser atom positions (14 dimensions instead of 37).atom14_atom_existsr*r!deviceresidx_atom14_to_atom37Zresidx_atom37_to_atom14Zatom37_atom_exists) rrCrBrestype_atom14_to_atom37r$r#restype_atom37_to_atom14restype_atom14_maskr restype_atom37_maskrA)rr%r&r'r(Zprotein_aatyper'r'r(make_atom14_maskssH   r)cCs,tdd|tj}t|}tdd|}|S)NcS t|Sr4r)nr'r'r(r: z&make_atom14_masks_np..cSr*r4)rTarray)tr'r'r(r:r,)r rTndarrayr)r )routr'r'r(make_atom14_masks_npsr1c Cs8|d|d<|d|d<|d|d<|d}|d}|t|d|dt|djddd}|d t|d|d t|djdd d}||d<||d <||d <tjtj|dj |dj d }||d}t d||}||d<t d||}||d<tjtj |dj |dj d } | |d|d<|S)z?Constructs denser atom positions (14 dimensions instead of 37).r*rrr!r$r?NrIZnum_batch_dimsrZatom14_gt_existsZatom14_gt_positionsr"z...rac,...rab->...rbcZatom14_alt_gt_positionsz...ra,...rab->...rbZatom14_alt_gt_existsZatom14_atom_is_ambiguous) rArr rPr.rrCrBrenaming_matricesr!r#rrestype_atom14_is_ambiguous) rZresidx_atom14_maskr$Zresidx_atom14_gt_maskZresidx_atom14_gt_positionsr4Zrenaming_transformZalternative_gt_positionsZalternative_gt_maskr5r'r'r(make_atom14_positions!s^ r6:0yE>cs|d}|d}|d}t|jdd}tjgddtd}gd|ddd ddf<gd |ddd ddf<ttjD]-\}}tj|} t d D]} tj || rjtj | | } | d d||| d ddf<qKq>| g|jddddR} d | d<d | d<| tj | dddd df<tjd d<tfdd} | |}| |}|jgd||jR}t||d|d}t||dt|jddd}tj|dd ddf|dd ddf|ddddf|d}t| |d|d}t||dt|jddd}tj|ddd |}tjd |j|jd}t|gd|dd d R}d|d<d|d<t|d }|t|d}|j gd|ddR}tjd |j|jd}t|gd|ddd d R}tjD]7\} }tjtj | }t!t"tj |d } d |d|| d f<d|d|| d d d f<d|d|| d ddf<qst||d|d}t||d!|d}t|d }|t|d}|#}|#}||d"<||d#<||d$<||d%<||d&<|S)'Nr*rrr?)rXrrrr+)CrNr)rr8Orrr@rXr.r).r.rWcs|Sr4r'r5Z lookuptabler'r(r:sz"atom37_to_frames..rr2r3rG)Z p_neg_x_axisoriginZ p_xy_planeepsrHr").rrr).rrGrG)Zmatrigidgroups_gt_framesrigidgroups_gt_existsZrigidgroups_group_existsZrigidgroups_group_is_ambiguousZrigidgroups_alt_gt_frames)$rPr.rTfullobjectrrBZrestypesZ restype_1to3rchi_angles_maskZchi_angles_atomsr new_tensorrcopyZ vectorizerr r from_3_pointsrryeyer!r#rrZcomposeZresidue_atom_renaming_swapsrrZ restype_3to1rSrqZ to_tensor_4x4)rr?r*rrZ batch_dimsZ"restype_rigidgroup_base_atom_namesrestypeZrestype_letterZresnameZchi_idxnamesZrestype_rigidgroup_masklookupZ"restype_rigidgroup_base_atom37_idxZ!residx_rigidgroup_base_atom37_idxZ base_atom_posZ gt_framesZ group_existsZgt_atoms_existZ gt_existsZrotsZrestype_rigidgroup_is_ambiguousZrestype_rigidgroup_rotsrZresidx_rigidgroup_is_ambiguousZresidx_rigidgroup_ambiguity_rotZ alt_gt_framesZgt_frames_tensorZalt_gt_frames_tensorr'r<r(atom37_to_framesbs       "  rMc Cs8||d}||d}||d}|jddkrL|j}|jg|ddR||d<|jg|ddR||d <|jg|dR||d <|Stj|d d }|g|jdd ddd}tj||dddddddfgd d}|g|jdddd}tj||dddddfgdd}tj|dddddf|dddddfgdd} tj|dddddf|dddddfgdd} tj|dddddf|dddddfgdd} tj|dddfddtj|dddfdd} |dtj|dddfd|jd} tj|dddfd|jd|d}tjtj |j d}|d|ddddf}t ||dt |jdd}t tj}|gd||}||ddf}t ||dt |jddd}tj|d|jd}||}tj| ddddddf| ddddddf| ddddddf|gd d}tj| d| d|d|gdd}tj|ddddf|ddddf|ddddfdd}||ddddf}tj|d|dgdd}ttjt|d|jd d!d}||}||gd"d#t |jddtddf}|tj|df}tj|jg|jdRd$d%|gdd}||d}|d&krtjt|jddt|jddgdd}||d|d|d}||d|d|d}|||d<|||d <|||d <|S)'Nr*rrr?rrGtorsion_angles_sin_cosalt_torsion_angles_sin_costorsion_angles_maskrWrr=r@%r.rHr3r).rG)rIr!).rr#)rVrVrVrVr2rr7)r?).r@T)rIr!Zkeepdims)rprpgrprprprpr4rprr)r.rrclamprkprodr!Z as_tensorrBchi_atom_indicesr#r rPrrEappendrFrrHinvertapplystacksqrtrqZsquaresliceZchi_pi_periodicZnew_onesr-r)rrr*rrZ base_shaperZprev_all_atom_positionsZprev_all_atom_maskZpre_omega_atom_posZ phi_atom_posZ psi_atom_posZpre_omega_maskZphi_maskZpsi_maskrWZ atom_indicesZ chis_atom_posrEZ chis_maskZchi_angle_atoms_maskZtorsions_atom_posrQZtorsion_framesZfourth_atom_rel_posrOdenomZchi_is_ambiguousZmirror_torsion_anglesrPZplaceholder_torsionsr'r'r(atom37_to_torsion_anglessP     $***             r_cCs4|dddddddf|d<|dd|d<|S)NrA.rZtrue_frame_tensorrBr;Z frame_maskr'r/r'r'r(get_backbone_framess  r`cCsP|dj}|ddddddf||d<|ddddf||d<|S)NrrO.rZchi_angles_sin_cosrQZchi_mask)r!r)rr!r'r'r(get_chi_angless rac Csd|vr |djd}nd}|dkr|r0t|tjd|d}tjtjj||ddtjd}n tj t||tjd}| D]H\}}| drz||}Wn4t y}z(t |j|t d |t ||jt d |t d d d | DWYd}~nd}~ww|||<q?|S)Nr2r?rr@Frr+templaterzprotein:zprotein_shape:cSs$i|]\}}dt|vr||jqSr)dirr.rr'r'r(rsz"crop_templates..)r.ryrTrYrrrCrr$rr startswith Exceptionprint __class__)rZ max_templatesZsubsample_templatesrZ template_idxr%r&exr'r'r(crop_templatessD        ricCs@d|vr |djdn|djd}t|||}t|||}|S) crop to size.r*rrr@)r.get_single_crop_idxapply_crop_idx)r crop_sizerseedrcrop_idxr'r'r(crop_to_size_singles   rpc Cstj|ddtj|k}Wdn1swYd|vo'|ddk}|r3t|||d|S|r=t||||}nt|||}t|||S)rjZ multimer_croprNrr@)rmrrn) r numpy_seedrTrYrZrpget_spatial_crop_idxget_contiguous_crop_idxrl) rrmrrnZspatial_crop_probca_ca_thresholdZuse_spatial_croprror'r'r(crop_to_size_multimers(   rvrrm random_seedreturncCsj||kr t|St|ttjd||d}t|||WdS1s.wYdS)Nrr@)rrrrrrSrTrYr)rrmrwZ crop_startr'r'r(rks   $rkr use_multinomialc CsD|stj|ddtjt|}Wdn1swY|}t|}dd|D}t |D]K\}} || } || 8}t || } t | t d||} tj||ddt tjj t | t | dd } Wdn1sswY|| 8}| || <q5t|}|S|t|}ttjj||d }t ||}|S) z"get crop sizes for contiguous cropZmultimer_contiguous_permrqNcSrrr'rr'r'r(rrz-get_crop_sizes_each_chain..rZmultimer_contiguous_crop_sizer@)lowhigh)Zpvals)rrrrTrYrjrPrqrrCrryrrSrr_r)r rmrwryZ shuffle_idxZnum_leftZ num_budget crop_sizesjidxZthis_lenmax_sizeZmin_sizeZthis_crop_sizeZ entity_probsr'r'r(get_crop_sizes_each_chainsB     rc Cs|djd}||krt|Sd|vsJ|d}t||||}g}tjdtjd}tj|dd2t||D]$\} } t j dt | | d} | t|| || | || 7}q8Wdn1sgwYt|S)Nr*rr r+Z"multimer_contiguous_crop_start_idxrqr@)r.rrrrCr$rrrrrTrYrrSrXrk) rrmrwryrr r|Z crop_idxsZ asym_offsetZllZcszZ this_startr'r'r(rts0      rtL@ruinfcCs4tjd}|dd|ddf}|dd|f}|jdddkr*t|||S|d|ddddf}t|} t| |d ||} t | rit j |d d t t j| } Wdn1scwYnt|||S| | } || |<tjd | jd| jd d} | | 7} t| d|}|jS)Nrr.rr?rHr@rasym_idZmultimer_spatial_croprqrrTgMbP?)rBrrrqallrtget_pairwise_distancesget_interface_candidatesranyrrrrSrTrYrrr.r#rrhrvalues)rrmrwrurrZ ca_coordsZca_mask pair_mask ca_distancesinterface_candidatesZ target_resZto_target_distancesZ break_tieretr'r'r(rs.s<      rscoordscCs,|d|d}ttj|dddS)Nr3r=rGr?rH)rDrr\rq)rZ coord_diffr'r'r(rWsrrrrcCs\|d|ddddfk}|d||}tj|dk||k@dd}|jddd}|S) Nr.rprr?rHT)as_tuple)rrrqZnonzero)rrrruZ in_same_asymZcnt_interfacesrr'r'r(r\srcCsXi}|D]#\}}||vrqt||D]\}}|tkr$t|||}q|||<q|Sr4)rrrrr|)rrroZcropped_proteinr%r&rZdim_sizer'r'r(rlls rlr4)FF)rV)r)rr)r7)F)NF)r)`r functoolsrroperatorrtypingrrrnumpyrTrZ unicore.datarZ unicore.utilsr r r r Z(modelscope.models.science.unifold.configr rrrZ&modelscope.models.science.unifold.datarrBZ/modelscope.models.science.unifold.modules.framerrstrr/Z NumpyDictZ TorchDict__annotations__r{r)r0r3r>rFrUr]rcrgrirnrwrrrrrrrrrrrrrrrrrrrrrr rrrr r)r1r6rMr_r`rarirprvrSZTensorrkrrrtrrsrrrlr'r'r'r(s&        "   "      #   3        / 6  % A    -      ) # )