
    DUf                     Z    d dl ZddlmZ d Zd Zd ed           ed          dfd	ZdS )
    N   )numutilsc                     |                      t                    } t          |          }t          j        || d d         z
  d          dz  }t          j        || dd          z
  d          dz  }d||z
  z  S )Nr   r      g      ?)astypefloatnpmaximum)distbinsscaffold_lengthouter_areasinner_areass       f/var/www/html/software/conda/lib/python3.11/site-packages/cooltools/sandbox/pairs_scaling_functions.py_contact_areasr      sv    u%%HO,,O*_x}<a@@AEK*_x|;Q??1DK++,,    c                 "   ||k    r|\  }}t          | ||z
            }no|\  }}|\  }}	||k    r||}}|	|}	}t          | |	|z
            t          | ||z
            z
  t          | |	|z
            z
  }||k     r|t          | ||z
            z  }|S )N)r   )
r   region1region2startendareasstart1end1start2end2s
             r   contact_areasr      s    '
sxu55V#VFFt$D8TF]33Xv778Xtd{334 	
 &==^Hftm<<<ELr   g      $@g    cA2   c                    dd l m} ||}t          j        |||          }t	          |||          }| | d         |d         k    | d         |d         k     z  | d         |d         k    z  | d         |d         k     z           } | d         | d         z
  j        }	t          |	|j                  r,|                    |	|	|k    |	|k     z           |          \  }
}n*t          j        |	|	|k    |	|k     z           |          \  }
}||
|fS )Nr   )Npos1r   pos2)bins)

dask.arrayarrayr   logbinsr   values
isinstanceArray	histogramr
   )dfr   r   dmindmaxn_binsdar   r   distsobs_s               r   compute_scalingr3   (   s1   df555H(GW55E		Fwqz	!f:
"	$f:#	% f:
"	$
B Z"V*$,E%"" VeUd]ut|$DEHUUQQeUd]ut|$DEHUUUQS%r   )numpyr
   libr   r   r   intr3    r   r   <module>r8      s}             - - -  * *.CCHH33s88TV            r   