DUf3ddlZddlZddlZddlZddlZddlZGddZGddZGddZ Gdd Z dS) Nc`eZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zd ZdZdS) TestRemoteNhttp://raw.githubusercontent.com/dpryan79/pyBigWig/master/pyBigWigTest/test.bwc@tj|j}|J|S)N)pyBigWigopenfnameselfbws N/var/www/html/software/conda/lib/python3.11/site-packages/pyBigWigTest/test.pydoOpenzTestRemote.doOpen s! ]4: & &~~~ ctj|j5}|dddksJ ddddS#1swxYwYdS)Nc Q?110)rrr chromsr s r doOpenWithzTestRemote.doOpenWiths ]4: & & E"99;; "C"CCCCCC E E E E E E E E E E E E E E E E E EsAA  A c|dddksJ|ddksJ|dJdS)Nrrrrc)rr s r doChromszTestRemote.doChromssZyy{{IY??????yy~~****yy~~%%%%%rc J|ddddddddksJdS) Nir)maxValsumDataminValversion sumSquarednLevels nBasesCovered)headerr s r doHeaderzTestRemote.doHeadersSyy{{saTUehuvJMNNNNNNNNrc|ddddgksJ|dddddgksJ|ddd dd d d gksJ|dtjdtjd dd d d gksJ|ddgksJdS)NrrgUUU?max)type@33?crr-nBins`ff??gq\?)statsnpint64r s r doStatszTestRemote.doStatssxxQ""'9&:::::xxQx..3F2GGGGGxxBsax88=NPS>>>>    %rcrdgdz}gd}gd}gd}tjd}|j}|t j|d}|d d g||d d |d d |d d |d d  ||d d|d d|d d|d d |d|dd|ddd|d|dd|ddd|dd|dddd|dd|dddd ||d d |d d |d d |d d  dd ksJn#t$rYnwxYw |d|dd|ddddd ksJn#t$rYnwxYw |d|dd|ddddd ksJn#t$rYnwxYw |dd|dddddd ksJn#t$rYnwxYw|dgd z|d d |d d |d d  |tj t |d  }|dksJtj|dS)z Test all three modes of storing entries. Also test to ensure that we get error messages when doing something silly This is a modified version of the writing example from libBigWig r) rrB}r0riXiqii iR)x~)?i@b@9@rqrrrsrtrurv4rvrwrxrvFrKrMr@B2`rr+rP )r<span i)r<rstepir3r|rb ef104f198c6ce8310acc149d0377fc16N)rRrSrTrUrrrWrY RuntimeErrorhashlibmd5read hexdigestrZr[) r rrarQr<r\r]r hs r doWrite2zTestRemote.doWrite2Us qKKK++++5999  ]5# & & nn5666 fQqSk6!A#;T!A#YvaPQc{ SSS fQqSk6!A#;T!A#YvaPQc{ SSS c6!A#;vac{ DDD c6!B$<qt 2 FFF c3vbe}2B GGG c3vbe}2B GGG  MM&1+vac{ac6RSTURU;M W W Wa44444    D   MM#vac{6!A#;RM H H Ha44444    D   MM#vac{6!A#;RM H H Ha44444    D   MM#s6"R%=rM K K Ka44444    D  seAgvac{ac6!A#; OOO  KUD))..00 1 1 ; ; = =44444 %sI8AF:: GG 1G== H  H 1I I  I *I<< J J ctjd}|j}|t j|d}|ddg|tjt |d  }|dksJt j|}| dd d ksJ| d dksJ| d d dgksJ|d d dd ddgksJ|tj|dS)NFrKrMryr{r 361c600e5badf0b45d819552a7822937rzr})rr|rrr)r2)rRrSrTrUrrrWrrrrrrEr<r5rZr[)r r\r]r rs r doWriteEmptyzTestRemote.doWriteEmptysa+5999  ]5# & & nn5666   KUD))..00 1 1 ; ; = =44444]5 ! !yy{{G'::::::||C  D((((yya))R////xxQqx11dD\AAAA   %rctjd}|j}|t j|d}|gdtjdgdzdgdzzdgdzzd gdzz}tjgd tj }tjgd tj }tjtj d tj }| |||||dt|||D}t j|}|dddddksJt!gdD]\} } t!|| D]b\} } | d|d| z| zksJ| d|d| z| zksJtj| d|d| z| zsJc|t'j|dS)NFrKrM))chr1rB)chr2rC)chr3r0)chr4rrrrr)r (2<FPZdtype)rkr-7AKU_rPc g|] \}}}|||f Sr).0xyzs r z+TestRemote.doWriteNumpy..s$DDDgaAAq DDDrrBrCr0r)rrrrrr)rRrSrTrUrrrWr6arrayr7random random_samplefloat64rYzipr enumeraterEiscloserZr[) r r\r]r rrarQvalues0valsidx1chromidx2tups r doWriteNumpyzTestRemote.doWriteNumpysN+5999  ]5# & & QQQRRR6(Q,&A51 DxRS|STT999JJJx777rxHHH(2922155RZHHH ff4 @@@  DDVT7)C)CDDD ]5 ! !yy{{sCcRRRRRR$%E%E%EFF E EKD%&r||E':':;; E E c1vD4!888881va$ho!66666z#a&'!d(T/*BCCDDDD E   %rcP|}|||jds||||||||||| || | | | |dS)Nhttp)rrr startswithr)r8rIr=rFrerrrrrUr s r testAllzTestRemote.testAlls [[]] bz$$V,, MM"    LL    JJrNNN MM"      R LL    OO    MMOOO              rN)__name__ __module__ __qualname__r rrrr)r8r=rFrIrerrrrrrr rrs \E EEE''' OOO666||| mmm MMM###J999v.4rrceZdZdZdS) TestLocalct}tjtjdz|_|dS)Nz/pyBigWigTest/test.bw)rrZpathdirnamer__file__r r)r blahs r testFoozTestLocal.testFoos:||W__X%677:QQ  rN)rrrrrrr rrs#rrceZdZdZdS) TestBigBedcztjtjdz}tj|}|J|dksJ|dksJd}|}t|tr| d}||ksJ| ddd}gd }||ksJ| dtjdtjd}||ksJ|dS) Nz/pyBigWigTest/test.bigBedrratable RnaElements "BED6 + 3 scores for RNA Elements data " ( string chrom; "Reference sequence chromosome or scaffold" uint chromStart; "Start position in chromosome" uint chromEnd; "End position in chromosome" string name; "Name of item" uint score; "Normalized score from 0-1000" char[1] strand; "+ or - or . for unknown" float level; "Expression level such as RPKM or FPKM. Set to -1 for no data." float signif; "Statistical significance such as IDR. Set to -1 for no data." uint score2; "Additional measurement/count e.g. number of reads. Set to 0 for no data." ) ASCIIrii))ii(z61035 130 - 0.026 0.42 404)i7͘iQΘz61047 136 - 0.029 0.42 404)iΘisz61048 630 - 5.420 0.00 2672399)rZrrrrrisBigWigisBigBedSQL isinstancebytesdecodeentriesr6r7rU)r r bbroutputoexpecteds r testBigBedzTestBigBed.testBigBeds+ 1225PP ]5 ! !~~~{{}}!!!!{{}}!!!!  fe $ $ ,]]7++F}}}} JJvhx 0 0CCCH}}}} JJvbhx00(1C1C D DH}}}}  rN)rrrrrrr rrs#rrceZdZdZdZdS) TestNumpyc, ddl}tjdkrdSddltjdd dgdjdgdz}jgdj }jgd j }jjdj } |||| jgd j }jgd j }jjdj } |||| jgdj }jjdj } d|jd|jgdj }jjdj } d|jd|jjdj } djdjdjd| jjdj } djdjdjd|  tjd J fd} | d|| d|| d|| d|| d| | d|  dddd} t| tsJ dddd} t| jsJ |jddS)Nrz/tmp/delete.bwrM)rirNrr) rrrrrrrrrrr) rkr#rrrrrrrP) rBnrlrC) isrh) r0rirjriiii"r)rr<) ,i6i@iJiTi^ihiri|iir)rrr<rcg}d||dzD]}||dj|}jj||z dksJdS)NrrBrgh㈵>)rErXrallabs)startv2vtr r6s r compyz"TestNumpy.testNumpy..compy!sA\\#ueck::  1 A26&"&R..4/00 1 1 1 1 1rrBr0rFnumpyT)rZrrrrWrr7rrrrYrUr<rlistndarrayr[)r rZrrarQrvalues1values2values3values4values5rfoor r6s @@r testNumpyzTestNumpy.testNumpysN >Q  1 ]+S 1 1 k]Q ///3%"*%%AAARRRrx???rxPPP"(2922266bjIII ff4 @@@LLLTVT\]]]rxJJJRTRZ[[["(2922266bjIII ff4 @@@LLLTVT\]]]"(2922266bjIII c6 G DDDLLLTVT\]]]"(2922266bjIII c6 G DDD"(2922266bjIII c828C==xrx{{!U\ ]]]"(2922266bjIII c828C==xrx{{!U\ ]]]  ]+ , ,~~~ 2 2 2 2 2 2 a c7 c7 c7 c7 c7iiQ5i11#t$$%%%iiQ4i00#rz**+++   "#####rc tjdkrdSddl}d}tj|d}|j|dddd|j|ddd|jd sJ|j|dddd d d|j|dddd d |jd sJdS) NrrrrrTrr) equal_nanmeanrk) rrrallcloser<rfloat32r5r)r r6r r s r testNumpyValueszTestNumpy.testNumpyValues8s >Q  1` ]5# & &r{ IIc1bI - - BHRYYsAr**"* = = =      r{ HHS!R$H 7 7 BHRXXc1b&!44BJ G G G        rN)rrrrrrrr rrs3>$>$>$@     rr) rrRrZsysrrr6rrrrrrr rs GGGGGGGGR @S S S S S S S S S S r