bug fixed

Author: aghaeifar

src/file_utils.cpp

@@ -313,7 +313,7 @@ bool file_utils::read_phantom(std::string phantom_filename, std::vector<float> &
         if (product(dims) != mask.size())
         {
             BOOST_LOG_TRIVIAL(warning) << "Hint in reading mask: " << phantom_filename;
-            BOOST_LOG_TRIVIAL(warning) << "Mask size does not match the expected size: " << product(dims) << " vs " << mask.size() << ". Aborting...!";
+            BOOST_LOG_TRIVIAL(warning) << "Mask size does not match the expected size: " << product(dims) << " vs " << mask.size() << ". resize it...!";
             mask.resize(product(dims));
         }
         read_h5(phantom_filename, mask.data(), "mask", "uint8_t");

src/kernels.cu

@@ -136,7 +136,7 @@ void sim(const simulation_parameters *param, const float *pFieldMap, const uint8
             ind = sub2ind(xyz_new[0]*param->scale2grid[0], xyz_new[1]*param->scale2grid[1], xyz_new[2]*param->scale2grid[2], param->fieldmap_size[0], param->fieldmap_size[1], param->fieldmap_size[2]);
             if(ind >= param->matrix_length || ind < 0)
             {
-                printf("Error:spin=%d, ind=%" PRId64 ", %d,\n\tscale2grid=(%.10f, %.10f, %.10f)\n\txyz=(%.10f, %.10f, %.10f)\n\txyz_old=(%.10f, %.10f, %.10f)\n\tscale2grid*xyz(%.10f, %.10f, %.10f)\n\tPhantomSize=(%d, %d, %d)\n\tFoV=(%.10f, %.10f, %.10f)\n\tError?(%d, %d, %d)\n",spin_no, ind, current_timepoint, param->scale2grid[0], param->scale2grid[1], param->scale2grid[2], xyz_new[0], xyz_new[1], xyz_new[2], xyz_old[0], xyz_old[1], xyz_old[2], param->scale2grid[0]*xyz_new[0], param->scale2grid[1]*xyz_new[1], param->scale2grid[2]*xyz_new[2], param->fieldmap_size[0], param->fieldmap_size[1], param->fieldmap_size[2], param->fov[0], param->fov[1], param->fov[2], xyz_new[0] >= param->fov[0], xyz_new[1] >= param->fov[1], xyz_new[2] >= param->fov[2]);
+                printf("Error:spin=%d, ind=%" PRId64 ", %d,\n\tscale2grid=(%.10f, %.10f, %.10f)\n\txyz=(%.10f, %.10f, %.10f)\n\txyz_old=(%.10f, %.10f, %.10f)\n\tscale2grid*xyz(%.10f, %.10f, %.10f)\n\tPhantomSize=(%" PRId64 ", %" PRId64 ", %" PRId64 ")\n\tFoV=(%.10f, %.10f, %.10f)\n\tError?(%d, %d, %d)\n",spin_no, ind, current_timepoint, param->scale2grid[0], param->scale2grid[1], param->scale2grid[2], xyz_new[0], xyz_new[1], xyz_new[2], xyz_old[0], xyz_old[1], xyz_old[2], param->scale2grid[0]*xyz_new[0], param->scale2grid[1]*xyz_new[1], param->scale2grid[2]*xyz_new[2], param->fieldmap_size[0], param->fieldmap_size[1], param->fieldmap_size[2], param->fov[0], param->fov[1], param->fov[2], xyz_new[0] >= param->fov[0], xyz_new[1] >= param->fov[1], xyz_new[2] >= param->fov[2]);
                 return;
             }
             // ------ accumulate phase ------

src/simulation_parameters.h

@@ -49,7 +49,7 @@ typedef struct simulation_parameters
     int32_t timestep_us, TR_us, TE_us[MAX_TE], RF_us[MAX_RF], dephasing_us[MAX_DEPHASE], gradient_us[MAX_GRADIENT];
     int32_t n_dummy_scan;
     uint32_t n_spins, n_timepoints, n_fieldmaps, n_TE, n_RF, n_dephasing, n_gradient, n_fov_scale, n_tissue_type;
-    uint32_t fieldmap_size[3], seed, max_iterations;
+    size_t fieldmap_size[3], seed, max_iterations;
     int64_t matrix_length;
     bool enDebug, enCrossFOV, enRecordTrajectory, enProfiling;
     bool fieldmap_exist, mask_exist, no_gpu;
@@ -72,6 +72,7 @@ typedef struct simulation_parameters
         fieldmap_exist(true),
         mask_exist(true),
         enProfiling(false),
+        matrix_length(0),
         no_gpu(false)
     {
         memset(fieldmap_size,   0, 3*sizeof(fieldmap_size[0]));