[Mathematica] Newbie with Compile[]

[guerom00]

Hello all :)

I'm just beginning trying to use compiled version of several simple functions. And right at the beginning, I stumble onto a “Numerical error” error message which puzzle me…

Here is my simple code :

11/18/11 12:51:07 In[18]:= \[Epsilon]IntrinsicSiliciumC = 
  Compile[{{\[Omega], _Real}, {L, _Real}},
   Module[{\[Omega]0, \[Epsilon]0, \[Epsilon]Inf, c = 299792458},
    \[Omega]0 = (66 10^14)/(c (10^9) ) L;
    \[Epsilon]0 = 11.87;
    \[Epsilon]Inf = 1.035;
    \[Epsilon]Inf + ((\[Epsilon]0 - \[Epsilon]Inf) \[Omega]0^2)/(\
\[Omega]^2 + \[Omega]0^2)
    ]
   , CompilationTarget -> "C"];

11/18/11 13:05:18 In[23]:= \[Epsilon]IntrinsicSiliciumC[1., 1.]

11/18/11 13:05:19 During evaluation of In[23]:= CompiledFunction::cfne: Numerical error encountered; proceeding with uncompiled evaluation. >>

11/18/11 13:05:19 Out[23]= 1.0402

Basically, how to debug such error message (which is quite vague…)

Thanks in advance :)

 

[Adrmay]

</code>The error message you are getting is a bit vague, but it is telling you that your compiled code is not able to compute the result for that particular set of inputs. The most likely cause is that the inputs you are giving it are outside of the range of values that the compiled code is able to handle. To debug this, you'll need to figure out what range of inputs your compiled code is able to handle. This can involve testing different inputs and seeing which ones work and which ones don't. You can also look at the code itself to get an idea of what range of values it can handle.Once you have figured out the input range, you can modify your code so that it will produce a valid output for the inputs you are trying to use. This could involve adding conditions to check for values outside of the allowed range, or using alternative methods to compute the result. Hope this helps!