func juliaSinIter(c, seed complex128) float64 {
i := 0
for z := c; cmplx.Abs(z) < 50 && i < maxIter; i++ {
z = cmplx.Sin(z) * seed
}
return float64(maxIter-i) / maxIter
}
// Convenience constructor for a qubit, specified by its spherical coordinates
// on the Bloch sphere.
func NewQubitWithBlochCoords(theta, phi float64) *QReg {
// |psi> = cos(theta/2) + e^{i phi}sin(theta/2)
t := complex(theta/2, 0)
p := complex(phi, 0)
qreg := &QReg{1, []complex128{cmplx.Cos(t),
cmplx.Exp(complex(0, 1)*p) * cmplx.Sin(t)}}
return qreg
}
// Function to compute basis transformation matrix
func BasisTransform(th, phi float64) *[2][2]complex128 {
BTMatrix := new([2][2]complex128)
th_2, phi_2 := 0.5*th, 0.5*phi
csn := cmplx.Cos(complex(th_2, 0.0)) * cmplx.Exp(complex(0.0, -1.0*phi_2))
sn := cmplx.Sin(complex(th_2, 0.0)) * cmplx.Exp(complex(0.0, phi_2))
BTMatrix[0][0] = cmplx.Conj(csn)
BTMatrix[0][1] = cmplx.Conj(sn)
BTMatrix[1][0] = complex(-1.0, 0.0) * sn
BTMatrix[1][1] = csn
return BTMatrix
}