1 #include <Sequence/SimDataIO.hpp> 8 void print_problems(
const SimData & d,
11 int main(
int argc,
char ** argv )
18 cerr <<
"Writing in gzip\n";
19 gzFile gzf = gzopen(
"test_zlib_out.gz",
"w");
20 write_SimData_gz(gzf, d);
23 cerr <<
"Reading from gzip\n";
25 gzf = gzopen(
"test_zlib_out.gz",
"r");
26 SimData d2 = read_SimData_gz(gzf);
29 cerr <<
"Writing in binary\n";
31 ofstream obin(
"test_binary_out.bin",ios::binary);
32 write_SimData_binary(obin,d);
35 cerr <<
"Reading from binary\n";
37 ifstream ibin(
"test_binary_out.bin",ios::binary);
38 SimData d3 = read_SimData_binary(ibin);
41 cerr <<
"Writing to binary + gzip\n";
42 gzf = gzopen(
"test_zlib_out.bin.gz",
"bw");
43 write_SimData_gz(gzf, d,
true);
46 cerr <<
"Reading from binary + gzip\n";
48 gzf = gzopen(
"test_zlib_out.bin.gz",
"r");
49 SimData d4 = read_SimData_gz(gzf,
true);
54 cerr <<
"Error: d != d2\n";
59 cerr <<
"Error: d != d3\n";
64 cerr <<
"Error: d != d4\n";
70 void print_problems(
const SimData & d,
73 for(
unsigned i = 0 ; i < d.numsites() ; ++i )
75 if( d.position(i) != d2.position(i) )
77 cerr <<
"Position " << i <<
": " << d.position(i) <<
' ' << d2.position(i) <<
'\n';
79 for(
unsigned i = 0 ; i < d.size() ; ++i )
83 cerr <<
"Haplotype " << i <<
": " << d[i] <<
"\n\n" << d2[i] <<
'\n';
The namespace in which this library resides.
Data from coalescent simulations.