returned False when he expected True. Let’s take a look at what is going on here.<\/p>\n
The SameQ command (===) and the Equal command (==) are subtly different.\u00a0 Equal (==) tests to see if two expressions are mathematically equivalent<\/strong> whereas SameQ tests to see if they are structurally equivalent<\/strong>.\u00a0 In general, SameQ is the stricter of the two tests and is often too strict for our needs.<\/p>\n A very simple example<\/strong><\/p>\n The simplest example might be to compare 2 and 2. (That is 2 with a period at the end)<\/p>\n Mathematically these are the same since they both represent the number two. However, from a Mathematica structural point of view, they are different since one has infinite precision whereas the other has double precision. So, we have<\/p>\n A more complicated example<\/strong><\/p>\n Mathematically the above is always true but the structure of the expressions on the LHS and RHS are different under normal evaluation rules.<\/p>\n So, what do I mean by ‘structurally equivalent?’.\u00a0 To be honest, I am slightly woolly on this myself but in my head I consider two expressions to be structurally equivalent if their TreeForms are the same<\/p>\n The important thing to remember here is that the argument to TreeForm is evaluated by Mathematica before it gets passed to TreeForm. So, the above images show us that under normal evaluation rules, these two expressions have different structural forms.\u00a0 In the next example we’ll see how this detail can matter<\/p>\n Exp[I Pi] and all that Let’s look at<\/p>\n Would that be True or False based on what we have seen so far?\u00a0 It is certainly mathematically true but remember that a triple equals requires both sides to be structurally equivalent and at first sight it appears that this isn’t the case here. You might expect Exp[I Pi] to have a very different TreeForm from -1 and so you’d expect the above to evaluate to False. However if you do<\/p>\n then you get a single box containing -1. This occurs because Mathematica evaluates Exp[I Pi] to -1 before passing to TreeForm. This evaluation also occurs when you do Exp[I Pi] === -1 so what you are really testing is -1===-1 which is obviously True.<\/p>\n Back to the original example If you use === then it returns False and this is ‘obviously’ the case because if you evaluate the TreeForms of both expressions then you’ll see that they get stored with different structural forms after being evaluated.<\/p>\n Since we are interested in mathematical equivalence, however, we use ==<\/p>\n What happens this time, however, is that Mathematica returns the original expression unevaluated<\/strong>. This is because, under the standard set of transformations that Mathematica applies when it evaluates ==, it can’t be sure one way or the other and so it gives up.<\/p>\n To proceed we wrap the whole thing in FullSimplify which essentially says to Mathematica ‘Use everything you’ve got to try and figure this out – no matter how long it takes’\u00a0 I.e. you allow it to use more transforms.<\/p>\n Finally, we have the result that the user expected.<\/p>\n","protected":false},"excerpt":{"rendered":" A Mathematica user contacted me recently and asked why (1\/16 Cos[a] (Cos[b] + 15 Cos[3 b])) === (1\/8 Cos[b] (-7 + 15 Cos[2b])Cos[a]) returned False when he expected True. Let’s take a look at what is going on here. The SameQ command (===) and the Equal command (==) are subtly different.\u00a0 Equal (==) tests to […]<\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[8,7],"tags":[],"class_list":["post-3217","post","type-post","status-publish","format-standard","hentry","category-mathematica","category-programming"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p3swhs-PT","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/posts\/3217","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=3217"}],"version-history":[{"count":23,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/posts\/3217\/revisions"}],"predecessor-version":[{"id":3238,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=\/wp\/v2\/posts\/3217\/revisions\/3238"}],"wp:attachment":[{"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=3217"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=3217"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/walkingrandomly.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=3217"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}2==2. gives True\r\n2===2. gives False<\/pre>\n
Sqrt[2] + Sqrt[3] == Sqrt[5 + 2 Sqrt[6]] which gives True\r\nSqrt[2] + Sqrt[3] === Sqrt[5 + 2 Sqrt[6]] which gives False\r\n<\/pre>\n
TreeForm[Sqrt[2] + Sqrt[3] ]\r\n<\/pre>\n
![\"TreeForm[Sqrt[2]](\"https:\/\/www.walkingrandomly.com\/images\/mathematica8\/treeform_lhs.png\")
<\/p>\nTreeForm[Sqrt[5 + 2 Sqrt[6]] ]\r\n<\/pre>\n
![\"TreeForm[Sqrt[5](\"https:\/\/www.walkingrandomly.com\/images\/mathematica8\/treeform_rhs.png\")
<\/p>\n
\n<\/strong><\/p>\nExp[I Pi] === -1\r\n<\/pre>\n
TreeForm[Exp[I Pi]]\r\n<\/pre>\n
\n<\/strong><\/p>\n(1\/16 Cos[a] (Cos[b] + 15 Cos[3 b])) === (1\/8 Cos[b] (-7 + 15 Cos[2b])Cos[a])\r\n<\/pre>\n
TreeForm[(1\/16 Cos[a] (Cos[b] + 15 Cos[3 b]))]\r\n<\/pre>\n
![\"TreeForm[(1\/16](\"https:\/\/www.walkingrandomly.com\/images\/mathematica8\/treeform1_lhs.png\")
<\/p>\nTreeForm[(1\/8 Cos[b] (-7 + 15 Cos[2 b]) Cos[a])]\r\n<\/pre>\n
![\"TreeForm[(1\/8](\"https:\/\/www.walkingrandomly.com\/images\/mathematica8\/treeform1_rhs.png\")
<\/p>\n(1\/16 Cos[a] (Cos[b] + 15 Cos[3 b])) == (1\/8 Cos[b] (-7 + 15 Cos[2 b]) Cos[a])\r\n<\/pre>\n
FullSimplify[(1\/16 Cos[a] (Cos[b] + 15 Cos[3 b])) == (1\/8 Cos[b] (-7 + 15 Cos[2 b]) Cos[a])]\r\n\r\nTrue\r\n<\/pre>\n