nqrduck/nqr-blochsimulator
1
0
Fork 0
mirror of https://github.com/nqrduck/nqr-blochsimulator.git synced 2024-11-10 12:30:04 +00:00
Code Issues Projects Releases Packages Wiki Activity

Fixed real and imagniray components.

Browse source
This commit is contained in:
jupfi 2024-01-10 12:26:39 +01:00
parent 304b8a4be2
commit b734b1bcaa
1 changed files with 8 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

15
src/nqr_blochsimulator/classes/simulation.py
View file

@ -129,7 +129,7 @@ class Simulation:
# XY-Component # XY-Component
Mtrans = np.squeeze( Mtrans = np.squeeze(
M_sy1[0, :, :] + 1j * M_sy1[1, :, :] M_sy1[1, :, :] + 1j * M_sy1[0, :, :]
) # Indices start at 0 in Python ) # Indices start at 0 in Python
Mtrans_avg = np.mean(Mtrans, axis=0) Mtrans_avg = np.mean(Mtrans, axis=0)
Mtrans_avg = np.delete(Mtrans_avg, -1) # Remove the last element Mtrans_avg = np.delete(Mtrans_avg, -1) # Remove the last element
@ -180,6 +180,7 @@ class Simulation:
B1 = np.tile( B1 = np.tile(
(gadt * (real_pulsepower - 1j * imag_pulsepower) * B1).reshape(-1, 1), Nx (gadt * (real_pulsepower - 1j * imag_pulsepower) * B1).reshape(-1, 1), Nx
) )
K = gadt * xdis * w * self.gradient K = gadt * xdis * w * self.gradient
phi = -np.sqrt(np.abs(B1) ** 2 + K**2) phi = -np.sqrt(np.abs(B1) ** 2 + K**2)
@ -222,25 +223,25 @@ class Simulation:
for n in range(Nu): # time loop for n in range(Nu): # time loop
Mrot = np.zeros((3, Nx)) Mrot = np.zeros((3, Nx))
Mrot[0, :] = ( Mrot[0, :] = (
Bd1.T[:, n] * Mt[0, :] + Bd2.T[:, n] * Mt[1, :] + Bd3.T[:, n] * Mt[2, :] Bd1.conj().T[:, n] * Mt[0, :] + Bd2.conj().T[:, n] * Mt[1, :] + Bd3.conj().T[:, n] * Mt[2, :]
) )
Mrot[1, :] = ( Mrot[1, :] = (
Bd4.T[:, n] * Mt[0, :] + Bd5.T[:, n] * Mt[1, :] + Bd6.T[:, n] * Mt[2, :] Bd4.conj().T[:, n] * Mt[0, :] + Bd5.conj().T[:, n] * Mt[1, :] + Bd6.conj().T[:, n] * Mt[2, :]
) )
Mrot[2, :] = ( Mrot[2, :] = (
Bd7.T[:, n] * Mt[0, :] + Bd8.T[:, n] * Mt[1, :] + Bd9.T[:, n] * Mt[2, :] Bd7.conj().T[:, n] * Mt[0, :] + Bd8.conj().T[:, n] * Mt[1, :] + Bd9.conj().T[:, n] * Mt[2, :]
) )
Mt = np.dot(D, Mrot) + np.tile(b, (Nx, 1)).T Mt = np.dot(D, Mrot) + np.tile(b, (Nx, 1)).T
Mrot[0, :] = ( Mrot[0, :] = (
Bd1.T[:, n] * Mt[0, :] + Bd2.T[:, n] * Mt[1, :] + Bd3.T[:, n] * Mt[2, :] Bd1.conj().T[:, n] * Mt[0, :] + Bd2.conj().T[:, n] * Mt[1, :] + Bd3.conj().T[:, n] * Mt[2, :]
) )
Mrot[1, :] = ( Mrot[1, :] = (
Bd4.T[:, n] * Mt[0, :] + Bd5.T[:, n] * Mt[1, :] + Bd6.T[:, n] * Mt[2, :] Bd4.conj().T[:, n] * Mt[0, :] + Bd5.conj().T[:, n] * Mt[1, :] + Bd6.conj().T[:, n] * Mt[2, :]
) )
Mrot[2, :] = ( Mrot[2, :] = (
Bd7.T[:, n] * Mt[0, :] + Bd8.T[:, n] * Mt[1, :] + Bd9.T[:, n] * Mt[2, :] Bd7.conj().T[:, n] * Mt[0, :] + Bd8.conj().T[:, n] * Mt[1, :] + Bd9.conj().T[:, n] * Mt[2, :]
) )
Mt = Mrot Mt = Mrot