Skip to content
Snippets Groups Projects
Commit 7f7e4af0 authored by Pierre Augier's avatar Pierre Augier
Browse files

Correct doc and phase-shift random 1d

parent ac790f47
No related branches found
No related tags found
1 merge request!156Topic/default/phaseshift3d
......@@ -258,9 +258,10 @@
def _get_phase_shift_random(self):
"""Compute two random phase shift terms."""
alpha = uniform(-1, 1)
oper = self.sim.oper
alpha = uniform(-0.5, 0.5)
if alpha < 0:
beta = alpha + 0.5
else:
beta = alpha - 0.5
......@@ -262,14 +263,10 @@
if alpha < 0:
beta = alpha + 0.5
else:
beta = alpha - 0.5
phase_shift_alpha = np.exp(
1j * alpha * self.sim.oper.deltax * self.sim.oper.kx
)
phase_shift_beta = np.exp(
1j * beta * self.sim.oper.deltax * self.sim.oper.kx
)
phase_shift_alpha = np.exp(1j * alpha * oper.deltax * oper.kx)
phase_shift_beta = np.exp(1j * beta * oper.deltax * oper.kx)
return phase_shift_alpha, phase_shift_beta
......@@ -552,7 +549,7 @@
- Approximation 1:
.. math:: \p_t \log S = \sigma + \frac{N_0}{S_0},
.. math:: \p_t \log S = \sigma + \frac{N_{0\alpha}}{S_0},
Integrating from :math:`t` to :math:`t+\dt`, it gives:
......@@ -556,7 +553,7 @@
Integrating from :math:`t` to :math:`t+\dt`, it gives:
.. math:: S_1 = (S_0 + N_0 \dt) e^{\sigma \dt}.
.. math:: S_1 = (S_0 + N_{0\alpha} \dt) e^{\sigma \dt}.
- Approximation 2:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment