1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
|
#define MAXBITCOUNT 210
typedef struct our_bitset{
bitset<MAXBITCOUNT> bo;
} ob;
class bitset_class {
public:
ob binread;
void set_bits();
void conversion(unsigned short int*, int);
};
void bitset_class::set_bits(){
for(int i=0; i<MAXBITCOUNT; i++)
binread.bo.set(i,0);
}
void bitset_class::conversion(unsigned short int* single_read, int seq_length){
int i, j;
bitset<3> temp;
for (i=0; i<seq_length; i++, single_read++){
switch (*single_read){
case 0:
temp.set(0, 0);
temp.set(1, 0);
temp.set(2, 0);
// cout << "0\t" << temp << endl;
break;
case 1:
temp.set(0, 0);
temp.set(1, 0);
temp.set(2, 1);
// cout << "1\t" << temp << endl;
break;
case 2:
temp.set(0, 0);
temp.set(1, 1);
temp.set(2, 0);
// cout << "2\t" << temp << endl;
break;
case 3:
temp.set(0, 0);
temp.set(1, 1);
temp.set(2, 1);
// cout << "3\t" << temp << endl;
break;
case 4:
temp.set(0, 1);
temp.set(1, 0);
temp.set(2, 0);
// cout << "4\t" << temp << endl;
break;
}
for(j=0; j<3; j++){
binread.bo[MAXBITCOUNT-3*(i+1) + j] = temp[j];
}
}
}
int main(){
int i;
bitset_class* bin_reads = new bitset_class[n];
for(i=0; i<n; i++)
bin_reads[i].set_bits();
for(i=0; i<n; i++){
bin_reads[i].conversion(r[i], J);
cout<<i+1<<"th read is"<<endl<<bin_reads[i].binread.bo<<endl<<endl;
}
return 0;
}
|