typo fixed

Author: aghaeifar

src/file_utils.cpp

@@ -102,20 +102,23 @@ bool file_utils::read_config(std::string config_filename, simulation_parameters
     }
     filenames["output"] = file_paths;
 
-    // check header of fieldmap matches
-    auto dims = file_utils::get_size_h5(filenames.at("phantom")[0], "mask");
-    if(dims.size() != 3)
+    if(param->n_fieldmaps > 0)
     {
-        BOOST_LOG_TRIVIAL(error) << cf_name << ") " << "mask must be 3D but is " << dims.size() << "D filename:"<< filenames.at("phantom")[0] << ". Aborting...!";
-        return false;
-    }
-    std::vector<float> fov(3, 0);
-    file_utils::read_h5(filenames.at("phantom")[0], fov.data(), "fov", "float");
-    std::copy(dims.begin(), dims.end(), param->fieldmap_size);
-    std::copy(fov.begin(), fov.end(), param->fov);
+        // check header of fieldmap matches
+        auto dims = file_utils::get_size_h5(filenames.at("phantom")[0], "mask");
+        if(dims.size() != 3)
+        {
+            BOOST_LOG_TRIVIAL(error) << cf_name << ") " << "mask must be 3D but is " << dims.size() << "D filename:"<< filenames.at("phantom")[0] << ". Aborting...!";
+            return false;
+        }
+        std::vector<float> fov(3, 0);
+        file_utils::read_h5(filenames.at("phantom")[0], fov.data(), "fov", "float");
+        std::copy(dims.begin(), dims.end(), param->fieldmap_size);
+        std::copy(fov.begin(), fov.end(), param->fov);
 
-    param->fieldmap_exist = file_utils::get_size_h5(filenames.at("phantom")[0], "fieldmap").size() == 3;
-    param->mask_exist     = file_utils::get_size_h5(filenames.at("phantom")[0], "mask").size() == 3;
+        param->fieldmap_exist = file_utils::get_size_h5(filenames.at("phantom")[0], "fieldmap").size() == 3;
+        param->mask_exist     = file_utils::get_size_h5(filenames.at("phantom")[0], "mask").size() == 3;
+    }
 
     // ============== reading section SCAN_PARAMETERS ==============
     param->TR_us = (int32_t)pt.get<float>("SCAN_PARAMETERS.TR", param->TR_us);

src/spinwalk.cu

@@ -178,12 +178,16 @@ bool run(simulation_parameters param, std::map<std::string, std::vector<std::str
             {
                 param_local.fov[i] = fov_scale[sl] * param.fov[i];
                 param_local.scale2grid[i]    = param_local.fieldmap_size[i] / param_local.fov[i];
-                for (int n = 0; n < param.n_tissue_type; i++)
-                    if(param_local.fov[i] < 100 * param.diffusivity[n])
+                for (int n = 0; n < param.n_tissue_type; n++)
+                {
+                    double min_convergence =  2 * param.diffusivity[n] * sqrt(10 * param.n_timepoints); // https://submissions.mirasmart.com/ISMRM2024/Itinerary/PresentationDetail.aspx?evdid=4684
+                    if(param_local.fov[i] < min_convergence )
                     {
-                        BOOST_LOG_TRIVIAL(error) << "FoV (= " << param_local.fov[i] << ") must be 100x larger than sqrt(2*diffusivity*timestep) (= " << param.diffusivity[n] << ")!";
+                        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;
                     }
+                }
             }           
 
             if(param.no_gpu)