Marco Caliari schrieb: >> Then yes, you're absolutely right. I've never investigated complex >> problems before, so I expect there must be more changes that need to >> be done in the m-files to handle these kind of problems. >> >>> I think I will have some time to further investigate the codes (at >>> least, the m-files). >> >> Great, thanks, then I'll wait with that single modification to the >> m-file solvers until I either get a patch or modified m-files or >> whatever from you for all the other modifications. > > So, after the change ' -> .' everywhere, I modified this block > > if (~vnormcontrol) > --> vdelta = abs(y5 - y4); > --> vtau = max (vodeoptions.RelTol * abs(vu.'), vodeoptions.AbsTol); > else > vdelta = norm (y5 - y4, Inf); > vtau = max (vodeoptions.RelTol * max (norm (vu.', Inf), 1.0), ... > vodeoptions.AbsTol); > end > > with the absolute values where pointed. I have to say that I don't > understand > > max (norm (vu.', Inf), 1.0) > > Why the maximum between the norm and 1.0? > I modified in a similar way ode23, ode54 and ode78 and tested the > equation y'(t)=sqrt(i)*y(t) and a more complicated Schroedinger equation > and they work without any problem. > The change in vdelta makes the codes to behave like with "NormControl" = > "on", since the infinity-norm is used. I think the 2-norm should be used in > this case.
Can you share your Schroedinger-example or must I write one for myself? Can you also please send me one of your modified solver files then I hope that I won't forget anything. I'll compare that line-by-line... Thomas ------------------------------------------------------------------------- This SF.net email is sponsored by the 2008 JavaOne(SM) Conference Don't miss this year's exciting event. There's still time to save $100. Use priority code J8TL2D2. http://ad.doubleclick.net/clk;198757673;13503038;p?http://java.sun.com/javaone _______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev
