summaryrefslogtreecommitdiff
path: root/single_dequantify.py
blob: 3cf7d3afb375d595495729cbd6670a20497295cb (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
import sys
import subprocess
import math as mth
import cmath as cmt
import numpy as npy
import scipy.linalg as sla

alpha = 1
beta  = 1
dt = 0.1
iterations = 61

H = npy.array([[1,0],[0,2]])

if len(sys.argv) == 3:
    alpha = complex(sys.argv[1])
    beta = complex(sys.argv[2])

norm  = npy.linalg.norm([alpha, beta])
state = npy.array([alpha / norm, beta / norm])

def time_evolution(state, dt = dt):
    return npy.dot(state, sla.expm(-1j * dt * H))

def bloch_map(state):
    temp  = npy.conj(state[0] / state[1])
    theta = 2 * npy.arctan(abs(temp))
    phi = cmt.phase(temp)

    return (phi, theta)

def sphere2cart(phi, theta):
        return [
             mth.sin(theta) * mth.cos(phi),
             mth.sin(theta) * mth.sin(phi),
             -1 * mth.cos(theta)
        ]

historie = npy.array([bloch_map(state)])
f = open("data", "w")

for i in range(iterations + 1):
    historie = npy.vstack([historie,bloch_map(state)])
    (phi, theta) = bloch_map(state)
    coords = sphere2cart(phi, theta)
    f.write(f"{coords[0]}; {coords[1]}; {coords[2]}; {(i + 1) / (iterations + 1)}\n")
    state = time_evolution(state)

f.close()

subprocess.run(["gnuplot", "gnuplot.plt"])