gtSetupForwardSimulation: modified dialogues

Signed-off-by: preischig <preischig@xnovotech.com>

5_reconstruction/gtSetupForwardSimulation.m

@@ -86,7 +86,7 @@ else
 end
 
 % get reconstruction parameters
-check = inputwdefault('Use default parameters? [y/n]', 'y');
+check = inputwdefault('Use existing parameters for forward simulation? [y/n]', 'y');
 if ~strcmpi(check,'y')
     list = build_list_v2();
     parameters.fsim = gtModifyStructure(parameters.fsim, list.fsim, 1, 'Forward simulation parameters:');
@@ -94,7 +94,7 @@ if ~strcmpi(check,'y')
 end
 
 % launch forward simulation process to OAR - split to two jobs
-check = inputwdefault('Launch Forward Simulation with OAR? [y/n]', 'y');
+check = inputwdefault('Launch Forward Simulation with OAR (o) or here on the command line (c)? [o/c]', 'o');
 if strcmpi(check,'y')
     disp('Launching several separate OAR jobs');
 
@@ -117,11 +117,27 @@ if strcmpi(check,'y')
     disp(['Using ' num2str(OAR_parameters.njobs) ' OAR jobs']);
 
     pause(0.5);
-    OAR_make('gtForwardSimulate_v2', 1, ngrains(phaseID), OAR_parameters.njobs, ...
-             [parameters.acq.dir ' ' phaseID_str], true, ...
-             'walltime', OAR_parameters.walltime, 'mem', OAR_parameters.mem);
+    
+    if isfield(parameters.labgeo, 'sourcepoint') && ~isempty(parameters.labgeo.sourcepoint)
+        OAR_make('gtForwardSimulatePoly', 1, ngrains(phaseID), OAR_parameters.njobs, ...
+            [parameters.acq.dir ' ' phaseID_str], true, ...
+            'walltime', OAR_parameters.walltime, 'mem', OAR_parameters.mem);
+    else
+        OAR_make('gtForwardSimulate_v2', 1, ngrains(phaseID), OAR_parameters.njobs, ...
+            [parameters.acq.dir ' ' phaseID_str], true, ...
+            'walltime', OAR_parameters.walltime, 'mem', OAR_parameters.mem);
+    end
 else
-    gtForwardSimulate_v2(1, ngrains(phaseID), parameters.acq.dir, phaseID);
+    
+    disp('Launching forward simulation...')
+    disp(' ')
+    
+    if isfield(parameters.labgeo, 'sourcepoint') && ~isempty(parameters.labgeo.sourcepoint)
+        gtForwardSimulatePoly(1, ngrains(phaseID), parameters.acq.dir, phaseID);
+    else
+        gtForwardSimulate_v2(1, ngrains(phaseID), parameters.acq.dir, phaseID);
+    end
+
 end
 
 end % end of function