I am using a modified version of the sample Penny Lane VQLS code in order to solve a system of PDE’s. In this case the state |b> will not be U = H1H2H3. I have been able to represent |b> as the result of a sum unitaries, with scalar coefficients. Similar to decomposing the A operator (matrix) in A|x>=|b>. While I can see how to modify the sample code to add each of the gates for U, I am not sure how to apply the coefficients for summing here. Would these be rotations prior (or after?) to applying the CNOT, CY, CX?
Thank you very much for any help/advice