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/