# PlanarAngle issues?

Posted 16 days ago
594 Views
|
13 Replies
|
1 Total Likes
|
 I am observing a very strange behavior from PlanarAngle: The result of PlanarAngle depends on the previous history. PlanerAngle assumes for the first few calls that its arguments are complex. After that, it doesn't. Fortunately, however, all results are correct. Attachments:
13 Replies
Sort By:
Posted 16 days ago
 I am also using Windows 10 and Mathematica 12.3.1. And I can not reproduce the issue. No matter how many times PlanarAngle[{0,1} -> {{-Sin ... is evaluated the output always contains Conjugate Table[PlanarAngle[{0, 1} -> {{-Sin[\[Phi]], Cos[\[Phi]]}, {Cos[\[Phi]], Sin[\[Phi]]}}], 50] 
Posted 15 days ago
 Hmm, this is another miracle.For me, your code executed in a freshly loaded kernel looks like:
Posted 15 days ago
 Werner, so far I had only tried the code on Mathematica installed on my computer. The so called Mathematica|Desktop. And there are no problems.But after your recent reply I tried it on the Wolfram Cloud. And I get exactly the same strange results that you have described!
Posted 15 days ago
 I did not try that within the cloud but on the desktop only. I use a Wolfram|One license.
Posted 15 days ago
 After testing a bit more on another W10 computer I must correct myself. There the change from complex to real also appears, but later than it does on the Cloud. Typically around the 14th evaluation. Table[{i, PlanarAngle[{0, 1} -> {{-Sin[\[Phi]], Cos[\[Phi]]}, {Cos[\[Phi]], Sin[\[Phi]]}}]}, {i, 50}] 
Posted 15 days ago
 Very very strange. Does anybody have any idea what happens here?I think repeatability and reproducibility shoul be one of the most basic features of programming languages. PlanarAngle seems to show a kind of randomness instead.
Posted 15 days ago
 This is clearly a bug and should be reported here. To work around the issue until it is fixed, I recommend calling ClearSystemCache[] before PlanarAngle.
Posted 14 days ago
 @Jason: This makes results at least reproducible, just like Remove, ClearAll do (see above). But it is a very expensive operation since it clears all cached intermediate results so far. We now have: Of course it doesn't prevent PlanarAngle from changing its behavior.I reported this as a bug to Wolfram Support.
Posted 14 days ago
 A somewhat rough workaround could be to use VectorAngle instead: myPlanarAngle[p_, q1_, q2_] := VectorAngle[q1 - p, q2 - p] Table[{i, myPlanarAngle[{0, 1}, {-Sin[\[Phi]], Cos[\[Phi]]}, {Cos[\[Phi]], Sin[\[Phi]]}]}, {i, 50}] 
Posted 14 days ago
 @Hans: Yes, VectorAngle[u,v] remains stable. But VectorAngle is very different from PlanarAngle since it does not honor the direction of rotation. It gives the smaller of the two angles between u and v. Hence it is useless for calculating rotation angles. You must enhance it with quite complicated case diffentiations according to the relative positions of p, q1, q2. This is why PlanarAngle exists. See: p = {0, 0}; q1 = {1, 0}; Print[{q2 = #, VectorAngle[q1 - p, q2 - p], PlanarAngle[p -> {q1, q2}, "Counterclockwise"]}] & /@ {{-1, 1}, {-1, -1}}; {{-1,1},(3 \[Pi])/4,(3 \[Pi])/4} {{-1,-1},(3 \[Pi])/4,(5 \[Pi])/4} 
Posted 13 days ago
 Meanwhile I found out that the symbolic result of PlanarAngle is wrong. The correct symbolic result of PlanarAngle would have to be a Piecewise expression that depends on the relative positions of the points p, q1, q2.See the following and attached Notebook "210717b PlanarAngle Issues 2.nb". Attachments:
 Your findings are touched upon in the PlanarAngle documentation