From a801c828849778a53d2ca7891f5c3774ad39f489 Mon Sep 17 00:00:00 2001
From: David Mayerich <mayerich@uh.edu>
Date: Mon, 16 Dec 2019 21:08:35 +0000
Subject: [PATCH] Added code to save a matrix to a file for debugging

---
 src/basic_functions.h | 18 ++++++++++++++++++
 src/main.cpp          |  5 +++--
 2 files changed, 21 insertions(+), 2 deletions(-)

diff --git a/src/basic_functions.h b/src/basic_functions.h
index 6d51168..506d6d0 100644
--- a/src/basic_functions.h
+++ b/src/basic_functions.h
@@ -47,6 +47,24 @@ void mtxMultranspose(T* M3, T* M1, T* M2, size_t r1, size_t c1, size_t r2, size_
 	}		
 }
 
+template<typename T>
+void mtxOutputFile(std::string filename, T* M, size_t rows, size_t cols){
+	std::cout<<"Outputting "<<rows<<" x "<<cols<<" matrix to file: "<<filename<<std::endl;
+	ofstream outfile(filename.c_str());										//open a file for writing
+	//output the matrix
+	for(size_t r = 0; r < rows; r++){
+		for(size_t c = 0; c < cols; c++){
+			outfile<<M[r * cols + c];
+			if(c != cols-1)
+				outfile<<",";
+		}
+		if(r != rows - 1)
+			outfile<<std::endl;
+	}
+	outfile.close();
+	//outfile<<"This is a test."<<std::endl;
+}
+
 		//display within class scatter
 template<typename T>
 void displayS(T* sw, size_t f){
diff --git a/src/main.cpp b/src/main.cpp
index 68dc9e6..b6415cb 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -155,6 +155,7 @@ void gpuComputeEignS( size_t g, size_t fea){
 	
 	///DETERMINANT CURRENTLY REQUIRES OPENCV
 	std::cout<<"ERROR: This code requires a fix to remove an OpenCV dependence."<<std::endl;
+	mtxOutputFile("newSw.csv", nSw, r, r);
 	exit(1);
 	//FIX BY REPLACING THE FOLLOWING THREE LINES OF CODE USING LAPACK
 	//cv::Mat newSw = cv::Mat((int)r, (int)r, CV_32FC1, nSw);		//within scatter of gnome g in the population 
@@ -255,7 +256,7 @@ int main(int argc, char* argv[]){
 	args.add("generations", "number of generationsr", "50");
 	args.add("initial_guess", "initial guess of featues", "");
 	args.add("debug", "display intermediate data for debugging");
-	args.add("binary", "Select features for binary classes", "");
+	args.add("binary", "Calculate features based on class1 vs. all other classes", "");
 	args.add("trim", "this gives wavenumber to use in trim option of siproc which trims all bands from envi file except gagpu selected bands");
 
 	args.parse(argc,argv);                  //parse the command line arguments
@@ -263,7 +264,7 @@ int main(int argc, char* argv[]){
 //Print the help text if set
 	if(args["help"].is_set()){				//display the help text if requested
 		advertisement();
-		std::cout<<std::endl<<"usage: ga-gpu input output --option [A B C ...]"<<std::endl;
+		std::cout<<std::endl<<"usage: ga-gpu input_ENVI output.txt --classes class1.bmp class2.bmp ... --option [A B C ...]"<<std::endl;
 		std::cout<<std::endl<<std::endl;
 		std::cout<<args.str()<<std::endl;
 		exit(1);
--
libgit2 0.21.4