From 7762722f264adc402ea3b0f21923b18f072253ba Mon Sep 17 00:00:00 2001
From: xiangzhou
Date: Mon, 22 Sep 2014 11:06:02 -0400
Subject: version 0.95alpha
---
src/main.cpp | 86 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 86 insertions(+)
create mode 100644 src/main.cpp
(limited to 'src/main.cpp')
diff --git a/src/main.cpp b/src/main.cpp
new file mode 100644
index 0000000..e1fb336
--- /dev/null
+++ b/src/main.cpp
@@ -0,0 +1,86 @@
+/*
+ Genome-wide Efficient Mixed Model Association (GEMMA)
+ Copyright (C) 2011 Xiang Zhou
+
+ This program is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see .
+*/
+
+#include
+#include
+#include
+#include
+#include
+
+#include "param.h"
+
+#ifdef FORCE_FLOAT
+#include "gemma_float.h"
+#else
+#include "gemma.h"
+#endif
+
+using namespace std;
+
+
+
+int main(int argc, char * argv[])
+{
+ GEMMA cGemma;
+ PARAM cPar;
+
+ if (argc <= 1) {
+ cGemma.PrintHeader();
+ return EXIT_SUCCESS;
+ }
+ if (argc==2 && argv[1][0] == '-' && argv[1][1] == 'h') {
+ cGemma.PrintHelp(0);
+ return EXIT_SUCCESS;
+ }
+ if (argc==3 && argv[1][0] == '-' && argv[1][1] == 'h') {
+ string str;
+ str.assign(argv[2]);
+ cGemma.PrintHelp(atoi(str.c_str()));
+ return EXIT_SUCCESS;
+ }
+ if (argc==2 && argv[1][0] == '-' && argv[1][1] == 'l') {
+ cGemma.PrintLicense();
+ return EXIT_SUCCESS;
+ }
+
+ cGemma.Assign(argc, argv, cPar);
+
+ ifstream check_dir((cPar.path_out).c_str());
+ if (!check_dir) {
+ mkdir((cPar.path_out).c_str(), S_IRWXU|S_IRGRP|S_IROTH);
+ }
+
+ if (cPar.error==true) {return EXIT_FAILURE;}
+
+ if (cPar.mode_silence) {stringstream ss; cout.rdbuf (ss.rdbuf());}
+
+ cPar.CheckParam();
+
+ if (cPar.error==true) {return EXIT_FAILURE;}
+
+ cGemma.BatchRun(cPar);
+
+ if (cPar.error==true) {return EXIT_FAILURE;}
+
+ cGemma.WriteLog(argc, argv, cPar);
+
+ return EXIT_SUCCESS;
+}
+
+
+
--
cgit v1.2.3