Added code to save a matrix to a file for debugging

David Mayerich
1 parent e9a7821b
Showing 2 changed files with 21 additions and 2 deletions Show diff stats
src/basic_functions.h
src/main.cpp
@@ -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){
@@ -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);