Re: Re: Compatibility woes

*To*: mathgroup at smc.vnet.net*Subject*: [mg76276] Re: [mg76190] Re: Compatibility woes*From*: "Chris Chiasson" <chris at chiasson.name>*Date*: Fri, 18 May 2007 06:35:02 -0400 (EDT)*References*: <f2btsp$iha$1@smc.vnet.net> <f2em96$2pm$1@smc.vnet.net>

If you don't mind every argument to FindRoot being evaluated before FindRoot "sees" them, then you could do: ReleaseHold[Hold[FindRoot][args]] or Block[{FindRoot},FindRoot[args]] etc... On 5/17/07, Mark Adler <madler at alumni.caltech.edu> wrote: > On May 16, 3:26 am, Mark Adler <mad... at alumni.caltech.edu> wrote: > > I doubt that I have plumbed the depths of the > > incompatibilities, since the rate of discovery is not decreasing. > > In this case, I hate being right. The latest 5.2 to 6.0 > incompatibility showed up in FindRoot[]. A relatively complicated > function used in a FindRoot[] worked just dandy in 5.2, but now in 6.0 > I get complaints due to how and when the function is executed. Trying > Release[] or Hold[] around the function results in exactly the same > errors. Providing two starting points to avoid computing a Jacobian > also results in no change. The function itself works fine when > executed on its own. Putting it in a FindRoot[] is what exposes the > errors. Example errors are below. > > I am spending much more time than I like trying to get things to work > again that worked just fine before upgrading. So far, this one is > turning out to be one of the tougher nuts to crack. > > Mark > > > NDSolve::ndinnt: Initial condition 0.0174533 x is not a number or a > rectangular array of numbers. >> > > ReplaceAll::reps: {<<1>>} is neither a list of replacement rules nor > a valid dispatch table, and so cannot be used for replacing. >> > > Set::shape: Lists {vf,\[Gamma]f,sf,hf} and {v,\[Gamma],s,h}<<2>>{v^ > \[Prime][t]=={{-4.28284*10^13 Power[<<2>>] Sin[<<1>>]+0.000709677 > <<2>> <<21>>[<<5>>][<<1>>],<<1>>-<<1>>},<<1>>},<<6>>,<<1>>} are > not > the same shape. >> > > > -- http://chris.chiasson.name/