Note

You can download this notebook here.

Boson Sampling

import numpy as np
import piquasso as pq


with pq.Program() as program:
    pq.Q(all) | pq.StateVector([0, 1, 1, 0])

    pq.Q(0, 1) | pq.Beamsplitter(theta=np.pi / 3, phi=np.pi / 4)
    pq.Q(2, 3) | pq.Beamsplitter(theta=np.pi / 2, phi=np.pi / 3)

    pq.Q() | pq.ParticleNumberMeasurement()

simulator = pq.SamplingSimulator(d=4)

result = simulator.execute(program, shots=5)

print(result.samples)

[(1, 0, 0, 1), (1, 0, 0, 1), (1, 0, 0, 1), (1, 0, 0, 1), (1, 0, 0, 1)]