From 9f19c5c2119d89f05291ee446dc8c60dfd1818f4 Mon Sep 17 00:00:00 2001
From: Ali Aghaeifar <ali.aghaeifar@tuebingen.mpg.de>
Date: Wed, 10 Jul 2024 16:14:56 +0200
Subject: [PATCH] record echo time

---
 src/simulation_parameters.h |  5 +++++
 src/spinwalk.cu             | 14 +++++++++++---
 src/version.h               |  2 +-
 3 files changed, 17 insertions(+), 4 deletions(-)

diff --git a/src/simulation_parameters.h b/src/simulation_parameters.h
index b8caf61..0e2b264 100644
--- a/src/simulation_parameters.h
+++ b/src/simulation_parameters.h
@@ -181,6 +181,11 @@ typedef struct simulation_parameters
         if (n_dummy_scan < 0)
             n_dummy_scan = 5.0 * T1_ms[0] / TR_us * 1e3;
 
+        if (n_TE == 1)
+        {
+            BOOST_LOG_TRIVIAL(error) << "Echo time is not set";
+            return false;
+        }        
         if (mask_exist == false)
         {
             BOOST_LOG_TRIVIAL(error) << "Mask with labeled tissues is not provided";
diff --git a/src/spinwalk.cu b/src/spinwalk.cu
index ee15aed..8918141 100644
--- a/src/spinwalk.cu
+++ b/src/spinwalk.cu
@@ -176,9 +176,8 @@ bool run(simulation_parameters param, std::map<std::string, std::vector<std::str
                     double min_convergence = 0.95 * param.diffusivity[n] * sqrt(1.0 * param.TE_us[param.n_TE-2]); // https://submissions.mirasmart.com/ISMRM2024/Itinerary/PresentationDetail.aspx?evdid=4684
                     if(param_local.fov[i] < min_convergence )
                     {
-                        BOOST_LOG_TRIVIAL(error) << "Virtual FoV (= " << param_local.fov[i] << ") is smaller than minimum convergence length (= " << min_convergence << ")!";
-                        BOOST_LOG_TRIVIAL(error) << "Original FoV (= " << param.fov[i] << ") FoV Scale (= " << fov_scale[sl]  << ")!";
-                        return false;
+                        BOOST_LOG_TRIVIAL(warning) << "Virtual FoV (= " << param_local.fov[i] << ") is smaller than minimum convergence length (= " << min_convergence << ")!";
+                        BOOST_LOG_TRIVIAL(warning) << "Original FoV (= " << param.fov[i] << ") FoV Scale (= " << fov_scale[sl]  << ")!";
                     }
                 }
             }           
@@ -253,6 +252,15 @@ bool run(simulation_parameters param, std::map<std::string, std::vector<std::str
 
         dims[0] = fov_scale.size(); dims[1] = 1; dims[2] = 1; dims[3] = 1;
         file_utils::save_h5(f, fov_scale.data(), dims, "scales", "double");
+
+        dims[0] = fov_scale.size(); dims[1] = 1; dims[2] = 1; dims[3] = 1;
+        file_utils::save_h5(f, fov_scale.data(), dims, "scales", "double");
+
+        std::vector<float> TE;
+        for(int i=0; i<param.n_TE; i++) TE.push_back(param.TE_us[i]*param.timestep_us*1e-6); 
+        TE.push_back(param.TR_us*param.timestep_us*1e-6);
+        dims[0] = TE.size(); dims[1] = 1; dims[2] = 1; dims[3] = 1;
+        file_utils::save_h5(f, TE.data(), dims, "TE", "float");
     }
 
 #ifdef __CUDACC__ 
diff --git a/src/version.h b/src/version.h
index 765bf22..4e2362e 100644
--- a/src/version.h
+++ b/src/version.h
@@ -6,7 +6,7 @@
 
 #define SPINWALK_VERSION_MAJOR 1
 #define SPINWALK_VERSION_MINOR 13
-#define SPINWALK_VERSION_PATCH 9
+#define SPINWALK_VERSION_PATCH 10
 
 //---------------------------------------------------------------------------------------------
 //