Skip to content

FIx photonics and noisy demos #742

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
rmoyard merged 2 commits into from
Apr 17, 2023
Merged

FIx photonics and noisy demos #742

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
ed11420
FIx two demos
rmoyard Apr 17, 2023
a744526
small change
rmoyard Apr 17, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion demonstrations/tutorial_noisy_circuits.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -231,7 +231,7 @@ def depolarizing_circuit(p):
# ensure that the trainable parameters give rise to a valid channel parameter, i.e., a number
# between 0 and 1.
#
ev = np.tensor(0.7781, requires_grad=False) # observed expectation value
ev = 0.7781 # observed expectation value

def sigmoid(x):
return 1/(1+np.exp(-x))
Expand Down
31 changes: 19 additions & 12 deletions demonstrations/tutorial_photonics.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -168,8 +168,8 @@ def vacuum_measure_p():


# Sample measurements in phase space
x_sample = vacuum_measure_x().numpy()
p_sample = vacuum_measure_p().numpy()
x_sample = vacuum_measure_x()
p_sample = vacuum_measure_p()

# Import some libraries for a nicer plot
from scipy.stats import gaussian_kde
Expand Down Expand Up @@ -287,8 +287,8 @@ def measure_coherent_p(alpha, phi):


# Choose alpha and phi and sample 1000 measurements
x_sample_coherent = measure_coherent_x(3, np.pi / 3).numpy()
p_sample_coherent = measure_coherent_p(3, np.pi / 3).numpy()
x_sample_coherent = measure_coherent_x(3, np.pi / 3)
p_sample_coherent = measure_coherent_p(3, np.pi / 3)

# Plot as before
xp = vstack([x_sample_coherent, p_sample_coherent])
Expand Down Expand Up @@ -514,8 +514,8 @@ def measure_squeezed_p(r):


# Choose alpha and phi and sample 1000 measurements
x_sample_squeezed = measure_squeezed_x(0.4).numpy()
p_sample_squeezed = measure_squeezed_p(0.4).numpy()
x_sample_squeezed = measure_squeezed_x(0.4)
p_sample_squeezed = measure_squeezed_p(0.4)

# Plot as before
xp = vstack([x_sample_squeezed, p_sample_squeezed])
Expand Down Expand Up @@ -617,22 +617,29 @@ def measurement(a, phi):


@qml.qnode(dev_exact2)
def measurement2(a, theta, alpha, phi):
def measurement2_0(a, theta, alpha, phi):
qml.Displacement(a, theta, wires = 0) # We choose the initial to be a displaced vacuum
qml.CoherentState(alpha, phi, wires = 1) # Prepare coherent as second qumode
qml.Beamsplitter(np.pi / 4, 0, wires=[0, 1]) # Interfere both states
return qml.expval(qml.NumberOperator(0)), qml.expval(qml.NumberOperator(1)) # Read out N
return qml.expval(qml.NumberOperator(0)) # Read out N

@qml.qnode(dev_exact2)
def measurement2_1(a, theta, alpha, phi):
qml.Displacement(a, theta, wires = 0) # We choose the initial to be a displaced vacuum
qml.CoherentState(alpha, phi, wires = 1) # Prepare coherent as second qumode
qml.Beamsplitter(np.pi / 4, 0, wires=[0, 1]) # Interfere both states
return qml.expval(qml.NumberOperator(1)) # Read out N


print(
"Expectation value of x-quadrature after displacement: {}\n".format(measurement(3, 0).numpy())
"Expectation value of x-quadrature after displacement: {}\n".format(measurement(3, 0))
)
print("Expected current in each detector:")
print("Detector 1: {}".format(measurement2(3, 0, 1, 0)[0].numpy()))
print("Detector 2: {}".format(measurement2(3, 0, 1, 0)[1].numpy()))
print("Detector 1: {}".format(measurement2_0(3, 0, 1, 0)))
print("Detector 2: {}".format(measurement2_1(3, 0, 1, 0)))
print(
"Difference between currents: {}".format(
measurement2(3, 0, 1, 0)[1].numpy() - measurement2(3, 0, 1, 0)[0].numpy()
measurement2_1(3, 0, 1, 0) - measurement2_0(3, 0, 1, 0)
)
)

Expand Down