from scipy.stats import kruskal
import pandas as pd
import numpy as np

#if p value less than alpha reject
qiskit_data = np.full(100000, 0)
print(qiskit_data)

cirq_data = np.full(100000, 0)

qs_data = np.full(100000, 0)

kruskal(qiskit_data, cirq_data, qs_data)


qiskit_data = np.concatenate((np.full(41192, 8), np.full(1555, 9), np.full(7096, 10), 
                        np.full(226, 11), np.full(41044, 12), np.full(1637, 13),
                        np.full(6961, 14), np.full(289, 15)))

print(qiskit_data)

cirq_data = np.concatenate((np.full(41016, 8), np.full(1163, 9), np.full(7092, 10),
                      np.full(301, 11), np.full(41024, 12), np.full(1635, 13),
                      np.full(7017, 14), np.full(282, 15))) 

print(cirq_data)

qs_data = np.concatenate((np.full(40973, 8), np.full(1706, 9), np.full(7123, 10),
                      np.full(314, 11), np.full(40899, 12), np.full(1621, 13),
                      np.full(7076, 14), np.full(288, 15)))

print(qs_data)

kruskal(qiskit_data, cirq_data, qs_data)

qiskit_data = np.concatenate((np.full(1676, 8), np.full(40877, 9), np.full(285, 10), 
                        np.full(7103, 11), np.full(1610, 12), np.full(41135, 13),
                        np.full(256, 14), np.full(7058, 15)))

print(qiskit_data)

cirq_data = np.concatenate((np.full(1639, 8), np.full(40954, 9), np.full(285, 10),
                      np.full(6906, 11), np.full(1582, 12), np.full(41332, 13),
                      np.full(265, 14), np.full(7037, 15))) 

print(cirq_data)

qs_data = np.concatenate((np.full(1656, 8), np.full(41277, 9), np.full(251, 10),
                      np.full(7103, 11), np.full(1608, 12), np.full(40752, 13),
                      np.full(271, 14), np.full(7082, 15)))

print(qs_data)

kruskal(qiskit_data, cirq_data, qs_data)

qiskit_data = np.concatenate((np.full(41064, 8), np.full(1645, 9), np.full(7072, 10), 
                        np.full(261, 11), np.full(40850, 12), np.full(1612, 13),
                        np.full(7223, 14), np.full(273, 15)))

print(qiskit_data)

cirq_data = np.concatenate((np.full(41036, 8), np.full(1582, 9), np.full(7094, 10),
                      np.full(250, 11), np.full(41178, 12), np.full(1621, 13),
                      np.full(6967, 14), np.full(272, 15))) 


print(cirq_data)

qs_data = np.concatenate((np.full(41120, 8), np.full(1643, 9), np.full(6953, 10),
                      np.full(251, 11), np.full(41011, 12), np.full(1631, 13),
                      np.full(7131, 14), np.full(260, 15)))

print(qs_data)

kruskal(qiskit_data, cirq_data, qs_data)

