}
}
}
-
chimera = new Pintail(fastaFileNames[s], templatefile, filter, processors, maskfile, consfile, quanfile, window, increment, outputDir);
if (outputDir == "") { outputDir = m->hasPath(fastaFileNames[s]); }//if user entered a file with a path then preserve it
#if you are a mac user use the following line
TARGET_ARCH += -arch x86_64
+ #if you using cygwin to build Windows the following line
+ #CXX = x86_64-w64-mingw32-g++
+ #CC = x86_64-w64-mingw32-g++
+ #TARGET_ARCH += -m64
+
#if you are a linux user use the following line
#CXXFLAGS += -mtune=native -march=native -m64
}
}
/***********************************************************/
-double MothurFisher::hyper(double n11){
+double MothurFisher::myhyper(double n11){
try {
- return(hyper0(n11,0,0,0));
+ double hyper0Result = hyper0(n11,0,0,0);
+ return hyper0Result;
}catch(exception& e) {
- m->errorOut(e, "MothurFisher", "hyper");
+ m->errorOut(e, "MothurFisher", "myhyper");
exit(1);
}
}
}
prob=hyper0(n11,n1_,n_1,n);
sleft=0;
- p=hyper(min);
+ p=myhyper(min);
for(i=min+1; p<0.99999999*prob; i++)
{
sleft += p;
- p=hyper(i);
+ p=myhyper(i);
}
i--;
if(p<1.00000001*prob) sleft += p;
else i--;
sright=0;
- p=hyper(max);
+ p=myhyper(max);
for(j=max-1; p<0.99999999*prob; j--)
{
sright += p;
- p=hyper(j);
+ p=myhyper(j);
}
j++;
if(p<1.00000001*prob) sright += p;
return prob;
}catch(exception& e) {
- m->errorOut(e, "MothurFisher", "hyper0");
+ m->errorOut(e, "MothurFisher", "exact");
exit(1);
}
}
double lnfact(double);
double lnbico(double, double);
double hyper_323(double, double, double, double);
- double hyper(double);
+ double myhyper(double);
double hyper0(double, double, double, double);
double exact(double, double, double, double);
};
string command = "mv " + oldName + " " + newName;
system(command.c_str());
#else
- m->mothurRemove(newName);
+ mothurRemove(newName);
int renameOk = rename(oldName.c_str(), newName.c_str());
#endif
return 0;
\r
#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)\r
#else\r
-#ifdef BIT_VERSION\r
-#include <io.h>\r
-int fseeko(FILE *stream, off_t offset, int whence)\r
- {\r
- off_t FilePos = _fseeki64(stream, offset, whence);\r
- return (FilePos == -1L) ? -1 : 0;\r
- }\r
-#define ftello(fm) (off_t) _ftelli64(fm)\r
-#else \r
+//#ifdef BIT_VERSION\r
+//#include <io.h>\r
+//int fseeko(FILE *stream, off_t offset, int whence)\r
+//// {\r
+// off_t FilePos = _fseeki64(stream, offset, whence);\r
+// return (FilePos == -1L) ? -1 : 0;\r
+// }\r
+//#define ftello(fm) (off_t) _ftelli64(fm)\r
+//#else \r
int fseeko(FILE *stream, off_t offset, int whence)\r
{\r
off_t FilePos = fseek(stream, offset, whence);\r
return (FilePos == -1L) ? -1 : 0;\r
}\r
#define ftello(fm) (off_t) ftell(fm)\r
-#endif\r
+//#endif\r
#endif\r
\r
void LogStdioFileState(FILE *f)\r
off_t GetStdioFileSize(FILE *f)\r
{\r
off_t CurrentPos = GetStdioFilePos(f);\r
- int Ok = fseeko(f, 0, SEEK_END);\r
+ off_t zeroPos = 0;\r
+ int Ok = fseeko(f, zeroPos, SEEK_END);\r
if (Ok < 0)\r
Die("fseek in GetFileSize");\r
\r
string noOutliers, outliers;
if ((!filter) && (seqMask == "")) {
- noOutliers = templateFileName + "pintail.quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.quan";
}else if ((!filter) && (seqMask != "")) {
- noOutliers =templateFileName + "pintail.masked.quan";
+ noOutliers =m->getRootName(m->getSimpleName(templateFileName)) + "pintail.masked.quan";
}else if ((filter) && (seqMask != "")) {
- noOutliers = templateFileName + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "masked.quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "masked.quan";
}else if ((filter) && (seqMask == "")) {
- noOutliers = templateFileName + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "quan";
}
decalc->removeObviousOutliers(quantilesMembers, templateSeqs.size());