With odepkg 0.8.0 and Octave 3.6.1, I am getting a segmentation fault in
my Debian system with the following code (taken from the test suite in
odepkg_octsolver_mebdfdae.cc):
###################################################################
function [ydot] = fpol (vt, vy, varargin) %# The Van der Pol
ydot = [vy(2); (1 - vy(1)^2) * vy(2) - vy(1)];
end
function [vout] = fout (vt, vy, vflag, varargin)
if (regexp (char (vflag), 'init') == 1)
if (size (vt) != [2, 1] && size (vt) != [1, 2])
error ('"fout" step "init"');
end
elseif (isempty (vflag))
if (size (vt) ~= [1, 1]) error ('"fout" step "calc"'); end
vout = false;
elseif (regexp (char (vflag), 'done') == 1)
if (size (vt) ~= [1, 1]) error ('"fout" step "done"'); end
else error ('"fout" invalid vflag');
end
end
vopt = odeset ('OutputFcn', @fout, 'OutputSel', 1, 'Refine', 5);
vsol = odebda (@fpol, [0 2], [2 0], vopt);
###################################################################
I get :
warning: Option "RelTol" not set, new value 1.0e-06 is used
warning: Option "AbsTol" not set, new value 1.0e-06 is used
warning: Option "InitialStep" not set, new value 1.0e-06 is used
warning: Option "MaxStep" not set, new value 1.6e-01 is used
warning: Option "MaxOrder" not set, new value 3 is used
panic: Segmentation fault -- stopping myself...
attempting to save variables to `octave-core'...
save to `octave-core' complete
Segmentation fault
It seems to be triggered by the presence of 'OutputFcn' in teh cal to
odeset. Is this a known problem?
Rafael
------------------------------------------------------------------------------
Try before you buy = See our experts in action!
The most comprehensive online learning library for Microsoft developers
is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3, MVC3,
Metro Style Apps, more. Free future releases when you subscribe now!
http://p.sf.net/sfu/learndevnow-dev2
_______________________________________________
Octave-dev mailing list
Octave-dev@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/octave-dev
