Should we be applying neutralize() before computing mean? I can’t actually figure out how the neutralize() in the notebook applies the 1st equation from the screenshot above.
Also, I think corrected_e_wB may not require mu_ortho to compute the norm in the denominator.
(I’ve tried the above variations, but the results were almost similar to that of above.)
I have just checked to see if any progress has been made. All I can report back here is that the issue ticket has been allocated for investigation. If there is any update, I will report back ASAP.
I want to inform you that finally we have fixed that problem. I want to share with you that in the proposed implementation that we were using to get the “expected values”, we where using in the denominator of 3 expressions the norm of the vectors, when the formulas said norm**2. I confirm that the expected values are -0.7000 and 0.7000. Thanks for your patience and support.
I had the same issue with the neutralize() function where I saw -3.858...e-17 instead of the expected -4.442...e-17. I figured out the issue for that one. Instead of using
np.linalg.norm(g) ** 2
I changed the implementation to use
np.sum(g * g)
And that seems to have “fixed” the mismatch issue. Although, I don’t know why (since the former formula should also work IMO).
I get ±0.23 after equalising; I only get the ±0.700 answer after removing the \mu_{\perp} in the denominator of the expression for e_{w1/2B}^{corrected}.
So question: are the expressions in the assignment for e_{w1/2B}^{corrected} wrong? They certainly don’t line up with the ‘Step 2a’ quoted by jincy-p-janardhanan above…
In you case, I guess the error is in the denominator at the right. Most likely you used \mu in stead of \mu_B. In any cases, the error is in your code to calculate “corrected_e_w1B” and “corrected_e_w2B” if the output is same as Martin’s case.
Thanks for your reply. Unfortunately I could not identify the error in the denominator. From your reply I understand that for all the other values (mu, mu_B, mu_orth, e_w1B, ew2B) the expected output is what is shown in Martin’s post, right?