Hi,

I am facing some issues with HHL implementation here: https://resources.wolframcloud.com/ExampleRepository/resources/HarrowHassidimLloyd-HHL-algorithm/

Whenever I alter the matrix (maintaining hermiticity!), the results start to differ (right after one decimal point) between the simulated quantum and classical.

For example, if I take p*A, where p=1.25, it did not produce the expected outcome.

The same issue persisted when I picked A using Abs[SymmetricMatrix[RandomVariate[GaussianOrthogonalMatrixDistribution]]]

Later, I increased the matrix dimension to 8 (in addition to changing the precision accordingly).

I would like to know if anyone else has faced similar issues and if anyone has any remarks.

Thank you in advance. Arun