From d403f90142d7fda860cb62b597b92705d1036dfd Mon Sep 17 00:00:00 2001
From: Saad Jbabdi <saad@fmrib.ox.ac.uk>
Date: Wed, 8 Apr 2009 08:54:29 +0000
Subject: [PATCH] sorted mesh tracking

---
 probtrackxOptions.h | 17 +++++++++++------
 ptx_meshmask.cc     |  2 ++
 2 files changed, 13 insertions(+), 6 deletions(-)

diff --git a/probtrackxOptions.h b/probtrackxOptions.h
index 80efc0a..f38519e 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 c7ed837..8bd9850 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;
-- 
GitLab