Octave and Freemat
Sebastien Loisel
loisel at temple.edu
Tue Mar 4 16:09:17 CST 2008
Dear David,
Thank you for your quick reply.
That is an error, as the code on octave-forge is not part of Octave. The
> copyright strings were recently updated in octave-forge and I believe
>
OK, I see that ORTH
http://velveeta.che.wisc.edu/cgi-bin/cvsweb.cgi/~checkout~/octave/scripts/linear-algebra/orth.m?rev=HEAD&content-type=text/plain
is part of Octave, while ODE45
http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/extra/ode/inst/ode45.m?revision=HEAD&content-type=text/plainis
not. I'll just leave the copyrights alone. For those files that I lift
out of Octave (like ORTH), I will put a note that the file was copied into
Freemat, and that Freemat's not part of Octave, and otherwise leave the
copyright alone. Does that sound like the right idea?
> roots is an octave core function and not from Octave forge, so send the
> proposed patches to roots to maintainers at octave.org. As for octave-forge
The correct implementation of roots.m is
function z = roots(p)
if(any(isnan(p) | isinf(p)))
error('Input to ROOTS must not contain NaN or Inf.');
end
while(any(isinf(p./p(1))))
p=p(2:end);
end
p = vec(p);
n = numel(p)-1;
A = diag(ones(n-1,1),-1);
A(1,:) = -p(2:n+1)./p(1);
z = eig(A);
changes, those are easier to deal with as the acceptance process of
> patches is easy.. Just send to the octave-dev at lists.sourceforge.org list
> and in generally you then commit the code yorself..
OK, I'll see what can be done for that.
--
Sébastien Loisel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://www.cae.wisc.edu/pipermail/octave-maintainers/attachments/20080304/bb47d18e/attachment.html
More information about the Octave-maintainers
mailing list