Skip to content
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

Added Pauli noise and standard-basis projectors, collapse evolution and other improvements #26

Merged
merged 19 commits into from
Jun 4, 2024
Merged

Added Pauli noise and standard-basis projectors, collapse evolution and other improvements #26

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
19 commits
Select commit Hold shift + click to select a range
87b0cd8
Added standard basis projectors support and other improvements
haggaila May 15, 2024
46b02a9
lint fix
haggaila May 15, 2024
a4a47cd
reformatted tabs
haggaila May 15, 2024
857c219
reformat
haggaila May 15, 2024
914929f
Updated solver version
haggaila May 15, 2024
e04264f
Fixed redundant opening of empty output files
haggaila May 16, 2024
8838a19
Version change
haggaila May 16, 2024
b24b349
Collapse Parameter defined
haggaila May 16, 2024
f2c4916
reformat
haggaila May 16, 2024
309e67b
Fix some Python error messages and complete collapse loop
haggaila May 19, 2024
29ee2a1
Fixed relaxed parameter tests
haggaila May 19, 2024
4604896
Various improvements to gates handling, logging, and output
haggaila May 26, 2024
e0e2fd1
reformatting and solver version increase
haggaila May 26, 2024
6a9c864
Added placeholders for a bit flip Lindbladian term
haggaila May 26, 2024
42b63e7
Added placeholders for a bit flip Lindbladian term
haggaila May 26, 2024
b775a0a
implementation of Sx_Sx and Sy_Sy (to be tested)
gregoire-misguich May 29, 2024
eb49efc
Finished implementation of bit flip and bit-phase flip (Pauli) jump o…
haggaila Jun 4, 2024
fe33619
Added test for g_3, g_4
haggaila Jun 4, 2024
e17c7e3
Reformatting
haggaila Jun 4, 2024
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
82 changes: 41 additions & 41 deletions lindbladmpo/LindbladMPOSolver.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -246,14 +246,14 @@ def build(self, parameters: Optional[dict] = None):
if i_op != n_indices - 1:
file.write(",")
file.write("\n")
elif key == "custom_observables":
elif key == "custom_observables" or key == "collapse":
file.write(key + " = ")
observables: list = parameters[key]
n_observables = len(observables)
for i_obs, (obs_def, obs_components) in enumerate(observables):
file.write(obs_def[0])
# file.write(" ")
# file.write(obs_def[1])
file.write(" ")
file.write(obs_def[1])
file.write(":")
n_components = len(obs_components)
for i_component, obs_component in enumerate(obs_components):
Expand Down Expand Up @@ -533,7 +533,7 @@ def verify_parameters(
if not LindbladMPOSolver.is_float(parameters[key]):
check_msg += "Error 140: " + key + " is not a float\n"
continue
if key != "t_init" and parameters[key] <= 0:
if key == "tau" and parameters[key] <= 0:
check_msg += "Error 150: " + key + " must be larger than 0\n"
continue
if key == "t_init" and parameters[key] > parameters["t_final"]:
Expand All @@ -552,7 +552,11 @@ def verify_parameters(
+ " should be equal to or larger than 0 (integer)\n"
)
continue
elif key == "output_step" or key == "force_rho_hermitian_step":
elif (
key == "output_step"
or key == "force_rho_hermitian_step"
or key == "force_rho_hermitian_gates"
):
if not LindbladMPOSolver._is_int(parameters[key]):
check_msg += "Error 180: " + key + " should be an integer\n"
continue
Expand All @@ -570,6 +574,8 @@ def verify_parameters(
or (key == "g_0")
or (key == "g_1")
or (key == "g_2")
or (key == "g_3")
or (key == "g_4")
):
if LindbladMPOSolver.is_float(parameters[key]):
continue
Expand Down Expand Up @@ -716,7 +722,9 @@ def verify_parameters(
"list/np.array) in the size of number_of_qubits^2 of floats\n"
)
continue
elif key == "apply_gates" or key == "custom_observables":
elif (
key == "apply_gates" or key == "custom_observables" or key == "collapse"
):
if (
not isinstance(parameters[key], tuple)
and not isinstance(parameters[key], list)
Expand Down Expand Up @@ -759,7 +767,8 @@ def verify_parameters(
" (time, gate name, qubit, [qubit])\n"
)
continue
else: # Hence key == "custom_observables"
else: # Hence key == "custom_observables" or key == "collapse"
b_is_collapse = key == "collapse"
for g_tuple in custom_list:
tuple_len = len(g_tuple)
if (
Expand All @@ -775,15 +784,31 @@ def verify_parameters(
)
continue
obs_type = g_tuple[0][1]
if not isinstance(g_tuple[0][0], str) or (
obs_type != "g" and obs_type != "p"
):
if not isinstance(g_tuple[0][0], str):
check_msg += (
"Error 342: each member of the first element of"
+ key
+ " must be a tuple of the form"
" (obs_name, obs_type), with obs_type being either 'g' or"
" 'p' to indicate a gate-based observable or a Pauli expansion\n"
" (obs_name, obs_type)\n"
)
continue
if b_is_collapse:
if obs_type != "o":
check_msg += (
"Error 342: each member of the first element of"
+ key
+ " must be a tuple of the form"
" (obs_name, obs_type), with obs_type being 'o' to indicate"
" a 1Q operator expansion\n"
)
continue
elif obs_type != "g" and obs_type != "o":
check_msg += (
"Error 342: each member of the first element of"
+ key
+ " must be a tuple of the form (obs_name, obs_type),"
" with obs_type being either 'g' or 'o' to indicate"
" a gate-based observable or a 1Q operator expansion\n"
)
continue
for o_tuple in g_tuple[1]:
Expand All @@ -799,15 +824,15 @@ def verify_parameters(
" (gate_name, q0, q1, ...)\n"
)
continue
if obs_type == "p" and (
if obs_type == "o" and (
not isinstance(o_tuple[0], str)
or not LindbladMPOSolver.is_float(o_tuple[1])
or not LindbladMPOSolver._is_int(o_tuple[1])
):
check_msg += (
"Error 343: each member of Pauli-based component of"
"Error 343: each member of an operator-based component of"
+ key
+ " must be a tuple of the form"
" (pauli-string, weight)\n"
" (operator name, qubit)\n"
)
continue

Expand Down Expand Up @@ -1196,29 +1221,4 @@ def verify_parameters(
check_msg += "Error: unknown parameter key passed: " + key + "\n"
# End of: "for key in dict.keys(parameters)"

# More cross-parameter checks:
if ("t_final" in parameters) and ("tau" in parameters):
if (LindbladMPOSolver.is_float(parameters["tau"])) and (
LindbladMPOSolver.is_float(parameters["t_final"])
):
if (parameters["tau"] > 0) and (parameters["t_final"] > 0):
if parameters["tau"] > parameters["t_final"] - parameters.get(
"t_init", 0.0
):
check_msg += (
"Error 640: t_final (total time) is smaller than tau (time step "
"for time evolution)\n "
)
elif "output_step" in parameters:
if LindbladMPOSolver._is_int(parameters["output_step"]):
if parameters["output_step"] > 0:
if (
parameters["output_step"] * parameters["tau"]
> parameters["t_final"]
):
check_msg += (
"Error 650: output_step multiplied by tau is larger "
"than t_final (output_step in units of tau, times "
"tau is larger than the simulation time)\n "
)
return check_msg
Loading
Loading