#import plotting
import subprocess
import math as mth
import cmath as cmt
import numpy as npy
import scipy.linalg as sla

alpha = 1 / npy.sqrt(2)
beta  = 1 / npy.sqrt(2)
state = npy.array([alpha, beta])

dt = 0.067
iterations = 41

H = npy.array([[2,-1j],[1j,3]])


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) / 2)
    phi = cmt.phase(temp)

    return npy.array([phi, theta])

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

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

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

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