Added addrow function - adds a row at the bottom of a matrix

miscmaths.cc

@@ -667,6 +667,21 @@ namespace MISCMATHS {
       return r;
     }
 
+  int addrow(Matrix& m, int ncols)
+  {
+    if (m.Nrows()==0) {
+      Matrix mm(1,ncols);
+      mm=0;
+      m = mm;
+    } else {
+      Matrix mm(m.Nrows()+1,ncols);
+      mm = 0;
+      mm.SubMatrix(1,m.Nrows(),1,ncols) = m;
+      m = mm;
+    }
+    return 0;
+  }
+  
   //------------------------------------------------------------------------//
 
 

miscmaths.h

@@ -134,6 +134,7 @@ namespace MISCMATHS {
 
   void reshape(Matrix& r, const Matrix& m, int nrows, int ncols);
   ReturnMatrix reshape(const Matrix& m, int nrows, int ncols);
+  int addrow(Matrix& m, int ncols);
 
   int construct_rotmat_euler(const ColumnVector& params, int n, Matrix& aff);
   int construct_rotmat_euler(const ColumnVector& params, int n, Matrix& aff,