def draw(self):
  stave = [f"#q{i}" for i in range(self.num_qubits)]

  for gate in self.data:
    if gate[0] == "cx":
      mxlen = max([len(line) for line in stave])
      stave[gate[1]] += ("-" * (mxlen - len(stave[gate[1]]))
                         + ("--" if stave[gate[2]][-2] == "R" else "---") + "*")
      stave[gate[2]] += ("-" * (mxlen - len(stave[gate[2]]))
                         + ("--" if stave[gate[2]][-2] == "R" else "---") + "C")

      for i in range(gate[1] + 1, gate[2]):
        stave[i] += "-" * (mxlen - len(stave[i])) + "---|"

    elif gate[0][0] == "r":
      stave[gate[2]] += ("--" if stave[gate[2]][-2] == "R" else "---") + gate[0].upper()

    else:
      stave[gate[1]] += ("--" if stave[gate[1]][-2] == "R" else "---") + gate[0][0].upper()


  # Auto-measure
  # for l in range(len(stave)):
  #   stave[l] += "-" * (max([len(line) for line in stave]) - len(stave[l])) + "-/m"


  mxlen = max([len(line) for line in stave])
  print("=" * (mxlen + 3))
  for qubit in stave:
    print(qubit)
  print("=" * (mxlen + 3))

# define a function that determines a brightness for any given point
# uses a seed that is a list of four numbers
def get_brightness(x,y,qc,seed):
    qc.data.clear() # empty the circuit
    
    # perform rotations whose angles depend on x and y
    qc.rx( x * pi/seed[0], 0 )
    qc.ry( y * pi/seed[1], 0 )

    # calculate probability for outcome 1
    qc.measure(0,0)
    p = simulate(qc,shots=1000,get='counts')['1']/1000
    # return brightness depending on this probability
    # the chosen values here are fairly arbitrary
    if p>0.7:
        if p<0.8:
            return 1
        elif p<0.9:
            return 2
        else:
            return 3
    else:
        return 0

def test_multiqubit():
    qc = QuantumCircuit(7,7)
    qc.h(0)
    qc.cx(0,2)
    qc.cx(2,1)
    qc.h(5)
    qc.cx(5,3)
    qc.cx(3,4)
    qc.cx(3,6)
    ket = simulate(qc,get='statevector')
    check = True
    for string in ['0000000','0000111','1111000','1111111']:
        check = check and round(ket[int(string,2)][0],2)==0.50
    assert( check )
    for j in range(7):
        qc.measure(j,j)
    for get in ['counts','probabilities_dict']:
        counts  = simulate(qc,shots=shots,get='counts')
        check = True
        for string in ['0000000','0000111','1111000','1111111']:
            p = counts[string]/shots
            check = check and round(p,2)==0.25
        assert( check )

def test_multiqubit():
    qc = QuantumCircuit(7,7)
    qc.h(0)
    qc.cx(0,2)
    qc.cx(2,1)
    qc.h(5)
    qc.cx(5,3)
    qc.cx(3,4)
    qc.cx(3,6)
    ket = simulate(qc,get='statevector')
    check = True
    for string in ['0000000','0000111','1111000','1111111']:
        check = check and round(ket[int(string,2)][0],2)==0.50
    assert( check )
    for j in range(7):
        qc.measure(j,j)
    counts  = simulate(qc,shots=shots,get='counts')
    check = True
    for string in ['0000000','0000111','1111000','1111111']:
        p = float(counts[string])/shots
        check = check and round(p,2)==0.25
    assert( check )

import microqiskit as mq
myqc = mq.QuantumCircuit(3,3)
myqc.h(1)
myqc.h(2)
res = mq.simulate(myqc,shots=16,get='counts')
print(res)