Remove phase encoding from magnetization_to_signal

NLeSC-COMPAS · Apr 16, 2024 · 5d6f20a · 5d6f20a
1 parent b8066a6
commit 5d6f20a
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/simulate/phase_encoding_kernels.cuh b/src/simulate/phase_encoding_kernels.cuh
@@ -16,10 +16,10 @@ __global__ void phase_encoding(
     auto voxel = index_t(blockIdx.x * blockDim.x + threadIdx.x);
 
     if (readout < echos.size(0) && voxel < echos.size(1)) {
+        auto k_y = trajectory.k_start[readout].imag();
         auto y = parameters.get(voxel).y;
-        auto k = trajectory.k_start[readout].imag();
 
-        ph_en_echos[readout][voxel] = echos[readout][voxel] * exp(cfloat(0.0f, y * k));
+        ph_en_echos[readout][voxel] = echos[readout][voxel] * exp(cfloat(0.0f, y * k_y));
     }
 }
 

diff --git a/src/trajectories/cartesian_view.cuh b/src/trajectories/cartesian_view.cuh
@@ -19,7 +19,7 @@ cfloat rewind(cfloat m, float R2, float delta_t, TissueVoxel p) {
 // apply gradient prephaser (i.e. phase encoding + readout prephaser for Cartesian)
 COMPAS_DEVICE
 cfloat prephaser(cfloat m, float k_x, float k_y, float x, float y) {
-    return m * exp(cfloat(0, k_x * x + k_y * y));
+    return m * exp(cfloat(0, k_x * x /*+ k_y * y*/));
 }
 
 struct CartesianTrajectoryView {