Re: [OctDev] odepkg

2008-04-23 Thread Thomas Treichl
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 (a

Re: [OctDev] odepkg

2008-04-23 Thread Marco Caliari
> 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). >

Re: [OctDev] odepkg

2008-04-23 Thread Thomas Treichl
Marco Caliari schrieb: >>> other solvers): at line 336, it should be >>> >>> vdelta = abs(y5 - y4); > >> So number of evaluations and steps increased dramatically. I can see >> the same for ode45..ode78. Can you send an example where I can see >> results becoming better if I use your modificatio

Re: [OctDev] [Fwd: Re: Add a logdet function?]

2008-04-23 Thread Fredrik Lingvall
David Bateman wrote: > >> Yes that would be preferable but I don't think Octave (3.0.0) keeps >> track of if the matrix is symmetric (pos def), at least a test I did >> with for example, >> >> > > A = diag(1:4,1) + 10*eye(5) + diag(1:4,-1); > matrix_type (A) > ans = Positive Definite > b =

Re: [OctDev] readxml

2008-04-23 Thread David Bateman
Tomas Kubin wrote: > Hello, > > function xmlread.oct (miscellaneous-1.0.5) does not read an attribute > encoding in XML declaration > > > > however, I think, it is valid attribute. > Is it my mistake or a bug? > > Regards > Tomas > This xml read function is extremely limited, and is basically

Re: [OctDev] Release of Octave-Forge for Octave 3.0.1

2008-04-23 Thread David Bateman
Tatsuro MATSUOKA wrote: > Hi David > > --- David Bateman <[EMAIL PROTECTED]> wrote: > > >> Err can you point me to it.. I seem to have missed it. >> >> > At that time you were too busy. > > http://sourceforge.net/mailarchive/forum.php?thread_name=20080322120239.GE4753%40heima.gjk.dk&forum

[OctDev] readxml

2008-04-23 Thread Tomas Kubin
Hello, function xmlread.oct (miscellaneous-1.0.5) does not read an attribute encoding in XML declaration however, I think, it is valid attribute. Is it my mistake or a bug? Regards Tomas - This SF.net email is sponsored

Re: [OctDev] [Fwd: Re: Add a logdet function?]

2008-04-23 Thread David Bateman
Fredrik Lingvall wrote: > >> I don't know the numerical details of the different methods (e.g. >> chol2logdet vs. qr2logdet). In a similar way, I don't know if I should >> invert a matrix using SVD or QR. I really love that when I compute the >> inverse of a matrix using 'inv', Octave choo

Re: [OctDev] odepkg

2008-04-23 Thread Marco Caliari
>> other solvers): at line 336, it should be >> >> vdelta = abs(y5 - y4); > So number of evaluations and steps increased dramatically. I can see the same > for ode45..ode78. Can you send an example where I can see results becoming > better if I use your modification? vdelta is an error estimat