diff --git a/probtrackxOptions.h b/probtrackxOptions.h
index 80efc0ab3f106a625ca4e86b6eee180af42e3690..f38519e1f01394b2ba918024f1ad0f1f04695867 100644
--- a/probtrackxOptions.h
+++ b/probtrackxOptions.h
@@ -67,6 +67,7 @@ class probtrackxOptions {
   Option<bool> modeuler;
   Option<int> rseed;
   Option<bool> seedcountastext;
+  FmribOption<bool> splitmatrix2;
 
   void parse_command_line(int argc, char** argv,Log& logger);
   void modecheck();
@@ -203,12 +204,15 @@ class probtrackxOptions {
   modeuler(string("--modeuler"), false, 
 	   string("Use modified euler streamlining"), 
 	   false, no_argument),
-   rseed(string("--rseed"), 12345,
-	 string("Random seed"),
-	 false, requires_argument), 
-   seedcountastext(string("--seedcountastext"), false,
-	 string("Output seed-to-target counts as a text file (useful when seeding from a mesh)"),
-	 false, no_argument), 
+  rseed(string("--rseed"), 12345,
+	string("Random seed"),
+	false, requires_argument), 
+  seedcountastext(string("--seedcountastext"), false,
+		  string("Output seed-to-target counts as a text file (useful when seeding from a mesh)"),
+		  false, no_argument), 
+  splitmatrix2(string("--splitmatrix2"), false,
+		  string("split matrix 2 (in case it is too big)"),
+		  false, no_argument), 
    options("probtrackx","probtrackx -s <basename> -m <maskname> -x <seedfile> -o <output> --targetmasks=<textfile>\n probtrackx --help\n")
    {
      
@@ -253,6 +257,7 @@ class probtrackxOptions {
        options.add(modeuler);
        options.add(rseed);
        options.add(seedcountastext);
+       options.add(splitmatrix2);
      }
      catch(X_OptionError& e) {
        options.usage();
diff --git a/ptx_meshmask.cc b/ptx_meshmask.cc
index c7ed837929287e03197355754bbeb3ca2849dc8f..8bd9850717d83a835fe071261348e139cc2e4b79 100644
--- a/ptx_meshmask.cc
+++ b/ptx_meshmask.cc
@@ -19,9 +19,11 @@ void meshmask()
   probtrackxOptions& opts =probtrackxOptions::getInstance();
 
   // load seed mesh
+  cout<<"loading mesh files"<<endl;
   Mesh mseeds;
   mseeds.load(opts.meshfile.value());
   mseeds.load_fs_label(opts.seedfile.value());
+  cout<<"mesh files loaded"<<endl;
 
   // internally create seed mask in voxel space
   volume<float> seeds;