Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

In the following example, inside the minimumX's Module, x is initialized to 1, but the results say otherwise:

solutionQ = Function[{x, d}, IntegerQ[Sqrt[(x^2 - 1)/d]] && x != 1];
minimumX = 
  Compile[{{d, _Integer}}, 
   Module[{x = 1}, While[! solutionQ[x, d], x++]; x]];

minimumX[2]
minimumX[3]
minimumX[2]

(*
3

7

17
*)

If I define minimumX as:

minimumX = 
  Function[{d}, Module[{x = 1}, While[! solutionQ[x, d], x++]; x]];

or

minimumX = 
  Compile[{{d, _Integer}}, 
   Module[{x = 1}, x = 1; While[! solutionQ[x, d], x++]; x]];

It works as I was expecting:

(*
3

2

3
*)

Is this expected behavior, or should I report a bug?

(Mathematica 9.0.1)

share|improve this question
4  
This is definitely not expected and you should report it. The localization of x in the instruction {46, Function[{d}, While[ !solutionQ[x, d], x++]], {x, 2, 0, 1, Module}, 2, 0, 0, 6, 0, 17} is clearly not being respected. The same behavior is observed for With in place of Module, but not for Block, which works as it should. It also works correctly in versions 7 and 5.2; the bug seems to have been introduced in version 8, when major changes were made to the VM. (The compiled code generated by all versions is functionally identical and differs only by opcode numbering.) –  Oleksandr R. May 12 '13 at 3:38
    
For completeness, let me just note that the function can be written more optimally as With[{solutionQ = Function[{x, d}, FractionalPart[Sqrt[(x^2 - 1)/d]] == 0 && x != 1]}, minimumX = Compile[{{d, _Integer}}, Module[{x = 1}, While[! solutionQ[x, d], x++]; x]]]. This is somewhat beside the point as far as the question goes since the bug does not manifest in this case, but I'm sure others would have been moved to comment on it if I hadn't. Indeed, I think it's worth emphasizing to WRI support that this particular non-optimal formulation is necessary to expose the bug. –  Oleksandr R. May 12 '13 at 4:04
1  
I think, what's worth noticing is: if You define minimumX by SetDelayed(:=) not only Set(=) it works like it should. –  Kuba May 12 '13 at 7:32
    
Thank you all. I've just reported the bug. –  P. Fonseca May 12 '13 at 10:19

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.