ó 1™xTc@ s<ddlmZd„Zd„Zdefd„ƒYZdS(i˙˙˙˙(tdivisionc C sťddlm}|j||ƒ}|j||ƒ}|dksHt‚|ddks^t‚||}|||d|}||d}d||d|||d}|||fS(sFscale factor for (fo-k*fc)**2, only similar to factor for abs(fo-k*fc)i˙˙˙˙(tflexiiiţ˙˙˙i(tscitbx.array_familyRtsumtAssertionError( tfotfcRtk_numtk_dentktk_d_numtk_dtk_d2((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pytcalc_ks "c C sddlm}|||}|j|jƒdƒ}|j|dkdƒ||9}|j|ƒ}|j|ƒ} | dks†t‚|| } |dkr | S||j||ƒ|| | } |dkrÚ| | fS||j||ƒd|| | } | | | fS(Ni˙˙˙˙(Riii(RRtdoubletsizet set_selectedRRtNone( RRR R R Rtdeltastsignstt_numtt_dentttt_dtt_d2((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pytcalc_ts    $  (ttargetcB s5eZdddd„Zd„Zd„Zd„ZRS(g0Žä.˙++cC sě||gjdƒdks!t‚|dk rOddlm}|j|ƒ}nt||ƒ\|_|_|_ t |||j|j|j ƒ\|_ |_ |_ d|_d|_|dk rč|d}||kjƒ}|j|dƒ|j|dƒ|j |||_|jj|dƒ|j|ƒ}|j|ƒ} ||| | || } } } |j | |j | |} |j | |j | |}|j |j ||| }|j| ||ƒ||_|jj|dƒndS(Nii˙˙˙˙(Riyi(iii(tcountRRRRtabsR R R R RRtgtctf_calc_gradientstf_calc_hessianst iselectionRtrealtimagt vec3_double(tOtf_obstf_calct f_calc_abst fca_sq_epsRtfca_sqt isel_zerotatbtaatbbtabthaathbbthab((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pyt__init__#s.! !3     cC s|jS(N(R(R%((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pyt target_work<scC s|jS(N(R(R%((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pytgradients_work?scC s|jS(N(R (R%((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pyt hessians_workBsN(t__name__t __module__RR4R5R6R7(((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pyR!s  N(t __future__RR RtobjectR(((sC/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/xray/targets/r1.pyts