-
Notifications
You must be signed in to change notification settings - Fork 56
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
Hermite improve #93
Hermite improve #93
Conversation
Codecov Report
@@ Coverage Diff @@
## master #93 +/- ##
==========================================
+ Coverage 95.75% 96.19% +0.44%
==========================================
Files 11 11
Lines 800 894 +94
==========================================
+ Hits 766 860 +94
Misses 34 34
Continue to review full report at Codecov.
|
Fun fact. The very first implementation of |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
💯
@@ -21,6 +21,7 @@ All probability amplitudes of a pure Gaussian state | |||
All matrix elements of a general Gaussian state :func:`thewalrus.quantum.density_matrix()` | |||
A particular probability amplitude of pure Gaussian state :func:`thewalrus.quantum.pure_state_amplitude()` | |||
A particular matrix element of a general Gaussian state :func:`thewalrus.quantum.density_matrix_element()` | |||
The Fock representation of a Gaussian unitary :func:`thewalrus.quantum.fock_tensor()` | |||
================================================================================ ============================= |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
💯
for (int ii = 0; ii < dim; ii++) { | ||
H[nextCoordinate] = H[nextCoordinate] + R(k, ii) * y(ii, 0); | ||
} | ||
H[nextCoordinate] = H[nextCoordinate] + y(k, 0); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Was + R(k, ii) * y(ii, 0)
the original matrix multiplication that is now redundant?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
yes!
@@ -776,7 +776,7 @@ TEST(LoopHafnianEigenComplex, EvenRandom) { | |||
|
|||
} | |||
|
|||
// Check loop hafnian with eignevalues for complex matrices with odd dimensions. | |||
// Check loop hafnian with eigenvalues for complex matrices with odd dimensions. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
lol good catch
Improving the speed of the hermite polynomials and also changing their definition
Context:
Description of the Change:
Benefits:
Possible Drawbacks:
The docs need to be updated to reflect the new definition
Related GitHub Issues: