Commit f9ce9dc1 authored by Alexis  CRISCUOLO's avatar Alexis CRISCUOLO
Browse files

1.0

parent fad01857
......@@ -254,12 +254,12 @@ public class AlienRemover {
try { bloom = new LHBF(filename, true); k = bloom.k(); } catch ( ClassNotFoundException e ) {}
catch ( OutOfMemoryError e ) { System.err.println("not enough RAM: " + (alienFile.length()>>19) + "Mb required (java option -Xmx)"); System.exit(1); }
k_1 = k; --k_1;
msk = 1L; msk <<= BSIZE * k; --msk; //## NOTE: °°°1111111111111111 == 2*k 1s
del = msk & (msk << BSIZE); //## NOTE: °°°1111111111111100 == msk-3
drc = msk >>> BSIZE; //## NOTE: °°°0011111111111111 == 2*(k-1) 1s
rcA = BT; rcA <<= BSIZE * k_1; //## NOTE: °°°1100000000000000
rcC = BG; rcC <<= BSIZE * k_1; //## NOTE: °°°1000000000000000
rcG = BC; rcG <<= BSIZE * k_1; //## NOTE: °°°0100000000000000
msk = 1L; msk <<= BSIZE * k; --msk; //## NOTE: ooo1111111111111111 == 2*k 1s
del = msk & (msk << BSIZE); //## NOTE: ooo1111111111111100 == msk-3
drc = msk >>> BSIZE; //## NOTE: ooo0011111111111111 == 2*(k-1) 1s
rcA = BT; rcA <<= BSIZE * k_1; //## NOTE: ooo1100000000000000
rcC = BG; rcC <<= BSIZE * k_1; //## NOTE: ooo1000000000000000
rcG = BC; rcG <<= BSIZE * k_1; //## NOTE: ooo0100000000000000
rcT = BA; rcT <<= BSIZE * k_1; //## NOTE: rcT == 0, because BA == 0
//## counters
cnt_r_in1 = cnt_r_in2 = cnt_b_in1 = cnt_b_in2 = 0;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment