- int start = time(NULL);
-
- chimera = new Pintail(fastafile, templatefile, filter, processors, maskfile, consfile, quanfile, window, increment, outputDir);
-
- //set user options
- if (maskfile == "default") { m->mothurOut("I am using the default 236627 EU009184.1 Shigella dysenteriae str. FBD013."); m->mothurOutEndLine(); }
-
-
- string outputFileName, accnosFileName;
- if (maskfile != "") {
- outputFileName = outputDir + getRootName(getSimpleName(fastafile)) + maskfile + ".pintail.chimeras";
- accnosFileName = outputDir + getRootName(getSimpleName(fastafile)) + maskfile + ".pintail.accnos";
- }else {
- outputFileName = outputDir + getRootName(getSimpleName(fastafile)) + "pintail.chimeras";
- accnosFileName = outputDir + getRootName(getSimpleName(fastafile)) + "pintail.accnos";
- }
- bool hasAccnos = true;
-
- if (m->control_pressed) { delete chimera; return 0; }
-
- if (chimera->getUnaligned()) {
- m->mothurOut("Your template sequences are different lengths, please correct."); m->mothurOutEndLine();
- delete chimera;
- return 0;
- }
- templateSeqsLength = chimera->getLength();
-
- #ifdef USE_MPI
- int pid, end, numSeqsPerProcessor;
- int tag = 2001;
- vector<long> MPIPos;
- MPIWroteAccnos = false;
-
- MPI_Status status;
- MPI_Comm_rank(MPI_COMM_WORLD, &pid); //find out who we are
- MPI_Comm_size(MPI_COMM_WORLD, &processors);
-
- MPI_File inMPI;
- MPI_File outMPI;
- MPI_File outMPIAccnos;
-
- int outMode=MPI_MODE_CREATE|MPI_MODE_WRONLY;
- int inMode=MPI_MODE_RDONLY;
-
- //char* outFilename = new char[outputFileName.length()];
- //memcpy(outFilename, outputFileName.c_str(), outputFileName.length());
-
- char outFilename[1024];
- strcpy(outFilename, outputFileName.c_str());
-
- //char* outAccnosFilename = new char[accnosFileName.length()];
- //memcpy(outAccnosFilename, accnosFileName.c_str(), accnosFileName.length());
-
- char outAccnosFilename[1024];
- strcpy(outAccnosFilename, accnosFileName.c_str());