little optmised

src/LAPACK/inverse.winxed

@@ -1,7 +1,5 @@
 $load "Double.pbc";
 
-
-
 namespace dgetri_func{
 
 const int PRINT_DEBUG_STUFF = 0;
@@ -12,22 +10,21 @@ function dgetri_exec(var a)
         var lapack = loadlib('liblapack.so');
         if (lapack == null || !lapack)
                 die("Cannot find liblapack");
-        say("In Winxed Main!\n");
 
 	using dgetrf_func.ipiv_size;
         int ipiv_size_v;
         ipiv_size_v=ipiv_size(a);
-        say(ipiv_size_v);
+//      say(ipiv_size_v);
         int ipiv[ipiv_size_v];
         int i;
         for(i=0;i<ipiv_size_v;i++)
                 ipiv[i]=0;
-        say("printing ipiv:\n",ipiv);
-        using dgetrf_func.dgetrf_exec;
-        int info_dgetrf;
-        info_dgetrf=dgetrf_exec(a,ipiv);
+//      say("printing ipiv:\n",ipiv);
+        //using dgetrf_func.dgetrf_exec;
+        //int info_dgetrf;
+        //info_dgetrf=dgetrf_exec(a,ipiv);
 
-        say("Matrix A:");
+        say("Given Matrix:");
 	say(a);
 
 	var dgetri = dlfunc(lapack, "dgetri_", "vppppppp");
@@ -41,22 +38,19 @@ function dgetri_exec(var a)
 	lwork_d=max(1,lwork);
 
 
-	say(n);
-	say(m);
-	say(lda);
-	say("lwork=",lwork);
-	say("lwork_d=",lwork_d);
+//	say(n);
+//	say(m);
+//	say(lda);
+//	say("lwork=",lwork);
+//	say("lwork_d=",lwork_d);
 
 
 	var work=new 'NumMatrix2D';
-	work.resize(lwork_d+1,lwork_d+1);
-	say(work);
-	say("1");
+	work.resize(lwork_d,lwork_d);
 
 
 //calling the function
 	dgetri(n,a,lda,ipiv,work,lwork,info);
-say("1");
 
 	if(info==0)
         {
@@ -85,13 +79,12 @@ function debug(var matrix, string msg, var args [slurpy])
 
 function main[main](var args)
 {
-	say("hello");
 	var pla = loadlib("linalg_group");
 	var a = new 'NumMatrix2D';
         a.initialize_from_args(3, 3,
-                                11.0, 2.0, 3.0,
-                                4.0, 5.0, 6.0,
-                                7.0, 1.0, 2.0);
+                                2.0, -1.0, 1.0,
+                                0.0, 3.0, -3.0,
+                                0.0, 0.0, 2.0);
 	using dgetri_func.dgetri_exec;
 	int info;
 	info=dgetri_exec(a);