示例#1
0
void indexOnDisk()
{
    std::cout << "Building index for " << opt::readsFile << " on disk\n";
    BWTDiskParameters parameters;
    parameters.inFile = opt::readsFile;
    parameters.outPrefix = opt::prefix;
    parameters.bwtExtension = BWT_EXT;
    parameters.saiExtension = SAI_EXT;
    parameters.numReadsPerBatch = opt::numReadsPerBatch;
    parameters.numThreads = opt::numThreads;
    parameters.storageLevel = opt::gapArrayStorage;
    parameters.bBuildReverse = false;
    parameters.bUseBCR = (opt::algorithm == "bcr");
		
	if(opt::bBuildForward)
	{
		buildBWTDisk(parameters);
	}
		
	if(opt::bBuildReverse)
    {
        parameters.bwtExtension = RBWT_EXT;
        parameters.saiExtension = RSAI_EXT;
        parameters.bBuildReverse = true;
        buildBWTDisk(parameters);
    }
}
示例#2
0
文件: index.cpp 项目: avilella/sga
void indexOnDisk()
{
    std::cout << "Building index for " << opt::readsFile << " on disk\n";
    buildBWTDisk(opt::readsFile, opt::prefix, BWT_EXT, SAI_EXT, false, opt::numThreads, opt::numReadsPerBatch, opt::gapArrayStorage);
    
    if(opt::bBuildReverse)
        buildBWTDisk(opt::readsFile, opt::prefix, RBWT_EXT, RSAI_EXT, true, opt::numThreads, opt::numReadsPerBatch, opt::gapArrayStorage);
}