saveMatrix

1 files changed, 22 insertions, 3 deletions

diff --git a/lib/Numeric/GSL/gsl-aux.c b/lib/Numeric/GSL/gsl-aux.c
index b3b524d..d129aeb 100644
--- a/lib/Numeric/GSL/gsl-aux.c
+++ b/lib/Numeric/GSL/gsl-aux.c
@@ -351,7 +351,7 @@ int polySolve(KRVEC(a), CVEC(z)) {
 }
 
 int vector_fscanf(char*filename, RVEC(a)) {
-    DEBUGMSG("gsl_matrix_fscanf");
+    DEBUGMSG("gsl_vector_fscanf");
     DVVIEW(a);
     FILE * f = fopen(filename,"r");
     CHECK(!f,BAD_FILE);
@@ -373,7 +373,7 @@ int vector_fprintf(char*filename, char*fmt, RVEC(a)) {
 }
 
 int vector_fread(char*filename, RVEC(a)) {
-    DEBUGMSG("gsl_matrix_fscanf");
+    DEBUGMSG("gsl_vector_fread");
     DVVIEW(a);
     FILE * f = fopen(filename,"r");
     CHECK(!f,BAD_FILE);
@@ -384,7 +384,7 @@ int vector_fread(char*filename, RVEC(a)) {
 }
 
 int vector_fwrite(char*filename, RVEC(a)) {
-    DEBUGMSG("gsl_vector_fprintf");
+    DEBUGMSG("gsl_vector_fwrite");
     DVVIEW(a);
     FILE * f = fopen(filename,"w");
     CHECK(!f,BAD_FILE);
@@ -394,6 +394,25 @@ int vector_fwrite(char*filename, RVEC(a)) {
     OK
 }
 
+int matrix_fprintf(char*filename, char*fmt, int ro, RMAT(m)) {
+    DEBUGMSG("matrix_fprintf");
+    FILE * f = fopen(filename,"w");
+    CHECK(!f,BAD_FILE);
+    int i,j,sr,sc;
+    if (ro==1) { sr = mc; sc = 1;} else { sr = 1; sc = mr;}
+    #define AT(M,r,c) (M##p[(r)*sr+(c)*sc])
+    for (i=0; i<mr; i++) {
+        for (j=0; j<mc-1; j++) {
+            fprintf(f,fmt,AT(m,i,j));
+            fprintf(f," ");
+        }
+        fprintf(f,fmt,AT(m,i,j));
+        fprintf(f,"\n");
+    }
+    fclose (f);
+    OK
+}
+
 //---------------------------------------------------------------
 
 typedef double Trawfun(int, double*);