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.

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

Fixed identity operator evaluation and extended the accepted input types in compute_rdms #269

Merged

Conversation

alexfleury-sb
Copy link
Contributor

@alexfleury-sb alexfleury-sb commented Jan 18, 2023

There were several bugs preventing users to use the compute rdms function:
1.

{pauli_string_to_of(term): exp_val for term, exp_val in exp_data.items()}

to

{pauli_string_to_of(term): exp_val for term, exp_val in exp_vals.items()}

(wrong dictionary).

List are unhashable, changed output type of the pauli_string_to_of function to a tuple.

if qterm in exp_vals:
    exp_val = exp_vals[qterm] if qterm else 1.

was failing for the () term. I added a conditional statement.

  1. Can now accept dictionaries with string ("IXYZ") or openfermion-like ((1, "X), (2, "Y"), ...) terms.

@alexfleury-sb alexfleury-sb marked this pull request as ready for review January 25, 2023 20:37
@alexfleury-sb
Copy link
Contributor Author

I tested all the possibilities, i.e.:

  • exp_data with openfermion-like terms ((0, "X"), (1, "Y"), ...)
  • exp_data with string terms ("XY...")
  • exp_vals with openfermion-like terms
  • exp_vals with string terms

Also, a missing entry in exp_vals raises an error. A missing entry in exp_data only raises an error when no compatible bases are detected.

@ValentinS4t1qbit
Copy link
Contributor

Two things:

  • Please change the PR title to something a bit better, that briefly describes what the PR did, for all readers.
  • I would like you to add tests addressing the use cases this fix covers. The fact that all tests were passing before this PR and are still passing means that what you did is not reflected in them: we thus need more tests, that illustrate what you fixed / contributed.

@alexfleury-sb alexfleury-sb changed the title Bugs in the compute rdms function (WIP to confirm with exp. data) Fixed identity operator evaluation and extended the accepted input types in compute_rdms Jan 26, 2023
Copy link
Contributor

@ValentinS4t1qbit ValentinS4t1qbit left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you !

@ValentinS4t1qbit ValentinS4t1qbit merged commit 9e8b0f3 into goodchemistryco:develop Jan 27, 2023
@alexfleury-sb alexfleury-sb deleted the bugs_compute_rdms branch January 27, 2023 14:51
ValentinS4t1qbit added a commit that referenced this pull request Feb 9, 2023
* Update develop to catch up with main version 0.3.3 (#249)

* Fix link of fig in qchem & excited states notebook (#250)

* allow single flip index dis for qcc (#247)

* allow single flip index dis for jkmn

* Richardson extrapolation: bug fix + error estimation (#252)

* Bugfix: DMET with QCC (#253)

* iQCC using only Clifford circuits notebook (#254)

* pUCCD ansatz (#251)

* UHF reference (#240)

* uhf implementation with VQESolver functionality
* support for all types of orbital freezing
* add active_spin and uhf attributes to SecondQuantizedDMETFragment
* add spin to adapt_ansatz arguments

* added multi-product, grid_circuits and discrete_clock (#257)

* added multi-product, grid_circuits and discrete_clock

* translation to pennylane (#260)

* tangelo to pennylane format translation

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* bump testing version to 3.8 (#262)

Updating python version to 3.8 in automated tests, as 3.7 is no longer maintained by the Python dev team

* Auto threshold cutoff for small coefficients in LCU (#261)

* check for small value lcu
* changed to keep same vector length but apply no operations

* Openshell DMET (#208)

* Open-shell DMET.
* Fix for get_rdm CCSD.
* Added NAO localization.
* Added LiO2 spin=1 DMET test.
* Added UHF MF for DMET. New get_rdm for VQESolver.

* Save mid-circuit measurement (#256)

* A flag now allows users to save mid-circuit measurements for each shot run.

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>

* Fix for IBMQConnection (#264)

* Updated copyright year to 2023 (#267)

* added draw method to circuit (#266)

* added draw method to circuit (warning about font comes up in linux env)

* Cirq, qulacs, pennylane and projectq operator translation functions (#268)

* Support to bidirectional translation of operators for cirq, qulacs, PennyLane, projectq.
* Fixed a bug in QulacsSimulator.expectation_value_from_prepared_state.
* Important bug in qiskit -> tangelo op translation fixed

* Fixed identity operator evaluation and extended the accepted input types in compute_rdms (#269)

* Compute RDMs function can now hangle both term formats (string and openfermion tuple).
* Unhashable list for of -> conversion to tuples.
* Handling all the inputs in compute_rdms.

* Temporarily disabling OS-DMET (#271)

* Temporarily disabling for code verification purpose, as the implementation seems to return very odd and incorrect results for simple systems.

* Link to new example repo

* Update README.rst

* Revert changes that broke tests (#275)

* reverted changes that broke tests

---------

Co-authored-by: AlexandreF-1qbit <76115575+AlexandreF-1qbit@users.noreply.github.com>
Co-authored-by: James Brown <84878946+JamesB-1qbit@users.noreply.github.com>
Co-authored-by: KrzysztofB-1qbit <86750444+KrzysztofB-1qbit@users.noreply.github.com>
ValentinS4t1qbit added a commit that referenced this pull request Feb 16, 2023
* Update develop to catch up with main version 0.3.3 (#249)

* Fix link of fig in qchem & excited states notebook (#250)

* allow single flip index dis for qcc (#247)

* allow single flip index dis for jkmn

* Richardson extrapolation: bug fix + error estimation (#252)

* Bugfix: DMET with QCC (#253)

* iQCC using only Clifford circuits notebook (#254)

* pUCCD ansatz (#251)

* UHF reference (#240)

* uhf implementation with VQESolver functionality
* support for all types of orbital freezing
* add active_spin and uhf attributes to SecondQuantizedDMETFragment
* add spin to adapt_ansatz arguments

* added multi-product, grid_circuits and discrete_clock (#257)

* added multi-product, grid_circuits and discrete_clock

* translation to pennylane (#260)

* tangelo to pennylane format translation

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* bump testing version to 3.8 (#262)

Updating python version to 3.8 in automated tests, as 3.7 is no longer maintained by the Python dev team

* Auto threshold cutoff for small coefficients in LCU (#261)

* check for small value lcu
* changed to keep same vector length but apply no operations

* Openshell DMET (#208)

* Open-shell DMET.
* Fix for get_rdm CCSD.
* Added NAO localization.
* Added LiO2 spin=1 DMET test.
* Added UHF MF for DMET. New get_rdm for VQESolver.

* Save mid-circuit measurement (#256)

* A flag now allows users to save mid-circuit measurements for each shot run.

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>

* Fix for IBMQConnection (#264)

* Updated copyright year to 2023 (#267)

* added draw method to circuit (#266)

* added draw method to circuit (warning about font comes up in linux env)

* Cirq, qulacs, pennylane and projectq operator translation functions (#268)

* Support to bidirectional translation of operators for cirq, qulacs, PennyLane, projectq.
* Fixed a bug in QulacsSimulator.expectation_value_from_prepared_state.
* Important bug in qiskit -> tangelo op translation fixed

* Fixed identity operator evaluation and extended the accepted input types in compute_rdms (#269)

* Compute RDMs function can now hangle both term formats (string and openfermion tuple).
* Unhashable list for of -> conversion to tuples.
* Handling all the inputs in compute_rdms.

* Temporarily disabling OS-DMET (#271)

* Temporarily disabling for code verification purpose, as the implementation seems to return very odd and incorrect results for simple systems.

* Link to new example repo

* Update README.rst

* Revert changes that broke tests (#275)

* reverted changes that broke tests

* DMET + frozen orbitals for each fragment (#276)

* DMET + frozen orbitals,  added DMET/frozen tests.

* Example folder relocated to new repo Tangelo-Examples (#274)

* removed examples folder, removed notebooks from automated tests, added link higher in README.rst

* Docs update (#272)

* Updated conf.py, tested building docs and fixed comments/docs.
* Removed tutorials.rst with the notebooks.

* Desired measurement result (#263)

* New feature: return frequencies matching desired mid-circuit measurement values.

---------

Co-authored-by: James Brown <james.brown@1qbit.com>

* Bumping Tangelo version number in _version.py

* Update Changelog (v0.3.4)

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>
Co-authored-by: AlexandreF-1qbit <76115575+AlexandreF-1qbit@users.noreply.github.com>
Co-authored-by: James Brown <84878946+JamesB-1qbit@users.noreply.github.com>
Co-authored-by: KrzysztofB-1qbit <86750444+KrzysztofB-1qbit@users.noreply.github.com>
Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: James Brown <james.brown@1qbit.com>
Co-authored-by: GitHub Actions <noreply@github.com>
ValentinS4t1qbit added a commit that referenced this pull request Jun 29, 2023
* Update develop to catch up with main version 0.3.3 (#249)

* Fix link of fig in qchem & excited states notebook (#250)

* allow single flip index dis for qcc (#247)

* allow single flip index dis for jkmn

* Richardson extrapolation: bug fix + error estimation (#252)

* Bugfix: DMET with QCC (#253)

* iQCC using only Clifford circuits notebook (#254)

* pUCCD ansatz (#251)

* UHF reference (#240)

* uhf implementation with VQESolver functionality
* support for all types of orbital freezing
* add active_spin and uhf attributes to SecondQuantizedDMETFragment
* add spin to adapt_ansatz arguments

* added multi-product, grid_circuits and discrete_clock (#257)

* added multi-product, grid_circuits and discrete_clock

* translation to pennylane (#260)

* tangelo to pennylane format translation

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* bump testing version to 3.8 (#262)

Updating python version to 3.8 in automated tests, as 3.7 is no longer maintained by the Python dev team

* Auto threshold cutoff for small coefficients in LCU (#261)

* check for small value lcu
* changed to keep same vector length but apply no operations

* Openshell DMET (#208)

* Open-shell DMET.
* Fix for get_rdm CCSD.
* Added NAO localization.
* Added LiO2 spin=1 DMET test.
* Added UHF MF for DMET. New get_rdm for VQESolver.

* Save mid-circuit measurement (#256)

* A flag now allows users to save mid-circuit measurements for each shot run.

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>

* Fix for IBMQConnection (#264)

* Updated copyright year to 2023 (#267)

* added draw method to circuit (#266)

* added draw method to circuit (warning about font comes up in linux env)

* Cirq, qulacs, pennylane and projectq operator translation functions (#268)

* Support to bidirectional translation of operators for cirq, qulacs, PennyLane, projectq.
* Fixed a bug in QulacsSimulator.expectation_value_from_prepared_state.
* Important bug in qiskit -> tangelo op translation fixed

* Fixed identity operator evaluation and extended the accepted input types in compute_rdms (#269)

* Compute RDMs function can now hangle both term formats (string and openfermion tuple).
* Unhashable list for of -> conversion to tuples.
* Handling all the inputs in compute_rdms.

* Temporarily disabling OS-DMET (#271)

* Temporarily disabling for code verification purpose, as the implementation seems to return very odd and incorrect results for simple systems.

* Link to new example repo

* Update README.rst

* Revert changes that broke tests (#275)

* reverted changes that broke tests

* DMET + frozen orbitals for each fragment (#276)

* DMET + frozen orbitals,  added DMET/frozen tests.

* Example folder relocated to new repo Tangelo-Examples (#274)

* removed examples folder, removed notebooks from automated tests, added link higher in README.rst

* Docs update (#272)

* Updated conf.py, tested building docs and fixed comments/docs.
* Removed tutorials.rst with the notebooks.

* Desired measurement result (#263)

* New feature: return frequencies matching desired mid-circuit measurement values.

---------

Co-authored-by: James Brown <james.brown@1qbit.com>

* fix error for imaginary qubit_op with desired meas result (#278)

* Integrate bug fix into main before release (#278) (#279) (#280)

Co-authored-by: James Brown <james.brown@1qbit.com>
* fix error for imaginary qubit_op with desired meas result (#278)

* Depth function performance increase (#285)

* Depth function takes now less than half a second to compute it from scratch for ~120K depth circuit.

* Fix for (inverse of S and T gate) #287. (#288)

* Translation perf tests (operator, circuit) (#289)

* Performance test for translation layer of linq, for most formats

---------

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>

* deterministic desired_meas_result (#290)

* deterministic desired_meas_result

* Mifnohelper print: remove problem handle (#293)

* Removed problem handle from print. It now only focuses on chemistry information

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>

* Combinatorial mapping of fermionic Hamiltonians (#286)

* Initial implementation. Known performance scaling issues, may attempt to introduce a different QubitOperator in the future.

* Re-enabling open-shell DMET (#291)

* Re-enabling OS-DMET.

* Support for symbolic simulator (simpy) (#292)

Support for symbolic backend based on sympy, to provide algebraic expressions for states / circuits with less than 10 qubits. Notebook / examples to follow later on.

* Arg simulate_options and projective_circuit to VQESolver (#298)

* It is now possible to request desired measurements and other options for the "simulate" method called by the backend
* Projective circuits can be used with VQE algorithm(s). Either passed by the user explicitly or integrated through the ansatz.

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* Supporting symbols as parameters and add nsimplify to results (#300)

* Symbolic expression simplified.

* Simplification of exp. value and fixed a bug in translate op.

* Add ability to get inverse with Symbol parameters.

* Switch nsimplify to evalf.

* Trim trivial qubits from circuit and Hamiltonian (#302)

Remove qubits that are in a deterministic state in quantum circuits, and simplify qubit Hamiltonian accordingly, in order to reduce resource requirements while computing expectation values.

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* Push contributors towards develop and not to main for PRs (#303)

* Update to CONTRIB (main -> develop) and added workflow to help avoid PR to main that are not from develop

* Set the DMET virtual orbital truncation threshold at the user level (#304)

* Add option to turn off virt. orb. truncation.
* Changed the threshold instead of a flag. Occ. set to their abs values.

* Trim bug (#305)

* minor bug fixes + test

* Added typing for option dictionaries in quantum agorithms (#307)

* added typing for option dictionaries, which improves user experience (autocomplete, auto-loading of docstrings for different GUIs) and handle options better overall

* IntegralSolver class: base + support for Psi4 (#297)

* Tests for pyscf, test for psi4, tests for custom IntegralSolver
* temporary: install NO_PYSCF available to users
* new: MP2Solver

This PR will be followed by 1-2 other PRs to decouple pyscf from almost all algorithms and streamline the use of classical chem backends (support for all classical solvers, turning pyscf as an optional dependency for Tangelo, fixing docs and README etc as a result)

* Change multi-controlled CNOT to multi-controlled CX (#308)

* change multi-controlled CNOT to multi-controlled CX

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* Braket connection (#312)

Support for Braket connection, including batch submission of jobs. Tests, docs, refresh of env_var.sh
A PR for the QPU notebook will follow

---------

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>

* Support for UMP2 initial parameters (#310)

* Support for UMP2 initial parameters.

* Checkfile for IntegralSolverPySCF (#311)

* chkfile for IntegralSolverPySCF.

* QubitOperator import from Tangelo, for better encapsulation (#299)

* Use Tangelo's `QubitOperator` instead of OpenFermion when possible.
* Constructor method `from_openfermion` and export method `to_openfermion` available in `QubitOperator` class.

Remaining usage of Openfermion is tied to qubit mappings mostly, and code relying on other Openfermion features.

---------

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>

* FCISolverPsi4 (#309)

* moved chem tests to molecular_computation tests
* Changed FCISolver to class with solver attribute

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* Removed deprecated Linq functions (#316)

* Removed deprecated versions of trasnlation functions and Simulator.
---------

Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>

* CCSD solver psi4 (#313)

* add ccsd_solver test to psi4 testing
* get_rdm currently not supported, suggestions added in comments

* ADAPT: add spin to available options  (#317)

* add spin to available options and docstring

* MP2 psi4 (#315)

* Support for MP2 with psi4 backend.

* get_rdm and get_mp2_parameters currently unsupported, suggestions left in docstrings.

* skip performance tests (Linq) in pytest (#319)

Perrformance tests take time and should not be run systematically

* Typos in docs (#321)

* Added spaces in messages.

* Typos and code format.

* Add period.

* DMET Effective Core Potential fix(related to #306) (#318)

* potential fix for #306, test added

* Change installation (#320)

* remove no_pyscf, attempt windows test
* updated psi4 test

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* DMET-ecp test fix: change optimizer to minimize square (#323)

* improved Error message for MP2 initialization in UCCSD (pySCf currently required) (#322)

* improved Error message for MP2 initialization if pyscf is not found

* Qiskit-related: deprecation update to SparsePauliOp and bug workaround (#325)

* Stim clifford simulator integration (#314)

* direct tableau
* noise and circuit sampler
* decomposition into clifford gates
* clifford decomp tests, and SDAG in cirq
* is_clifford in gate class
* expand clifford decomp to integer values
* small comments, add pip install stim to workflow

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>

* Bumping Tangelo version number in _version.py

* Update CHANGELOG.md

---------

Co-authored-by: Valentin Senicourt <41597680+ValentinS4t1qbit@users.noreply.github.com>
Co-authored-by: AlexandreF-1qbit <76115575+AlexandreF-1qbit@users.noreply.github.com>
Co-authored-by: James Brown <84878946+JamesB-1qbit@users.noreply.github.com>
Co-authored-by: KrzysztofB-1qbit <86750444+KrzysztofB-1qbit@users.noreply.github.com>
Co-authored-by: Valentin Senicourt <valentin.senicourt@1qbit.com>
Co-authored-by: James Brown <james.brown@1qbit.com>
Co-authored-by: ValentinS4t1qbit <valentin@goodchemistry.gmail>
Co-authored-by: elloyd-1qbit <58313607+elloyd-1qbit@users.noreply.github.com>
Co-authored-by: GitHub Actions <noreply@github.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants