bugGNU Octave - Bugs: bug #40892, control package: sigma function...

 
 

bug #40892: control package: sigma function crashes Octave when given a state-space model

Submitter:  None
Submitted:  Fri 13 Dec 2013 10:42:35 PM UTC
   
 
Category:  Octave Package Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Segfault, Bus Error, etc.
Status:  Duplicate Assigned to:  None
Originator Name:  Saul Mtakula Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * 3.6.4
Operating System:  * Microsoft Windows Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Jump to the original submission

Mon 24 Feb 2014 03:47:24 AM UTC, comment #11: 

This is a duplicate of bug #41192. The problem goes away when upgrading OpenBLAS (the other bug contains a link to an updated OpenBLAS DLL).

Michael Goffioul <goffioul>
Fri 31 Jan 2014 05:29:09 PM UTC, comment #10: 


What I found is that the problems are with the Microsoft Visual Studio version of Octave 3.6.4  (octave-3.6.4-vs2010-setup.exe date 2013-05-29). The Windows MinGW version Octave 3.6.4 for  works just fine once I installed it.
In fact the bug reports #41192 and #40892 are a related to the Visual Studio Version.
To summarize the Visual Studio Version will cause crashes for matrix inversion, when there is too much unsuppressed output (if there is no semicolon after a large matrix for e.g), with the control toolbox sigma function if the system in state space form (presumably because there could be some intermediate inversion operation).
In my opinion the Visual Studio Version of Octave should be withdrawn; there must not be that many users with it otherwise there would be a huge outcry.

Anonymous
Tue 17 Dec 2013 04:03:22 AM UTC, comment #9: 

I am not sure which BLAS library is used. This is for the windows Octave MINGW binary, control toolbox 2.4.2. All the issues have been on 64 bit windows; could that be relevant.

Anonymous
Tue 17 Dec 2013 03:58:22 AM UTC, comment #8: 

This may or may not be relevant, but do you know which BLAS library your Octave installation is using? There have been reports of certain operations crashing when using one or the other of the BLAS library choices.

Mike Miller <mtmiller>
Group Member
Tue 17 Dec 2013 03:55:03 AM UTC, comment #7: 

Updated again to reflect that it does indeed crash Octave and is not about errors or compatibility.

Mike Miller <mtmiller>
Group Member
Tue 17 Dec 2013 12:03:33 AM UTC, comment #6: 

Hi Lucas
I did not mean to  appear rude in my last response so apologies. I am new to Octave.
I notice there is a some kind of "buffer" error before the crash but I can't see the complete listing. I have run into lots of "buffer" error issues and "panicks" so it could be related to that.

Anonymous
Mon 16 Dec 2013 08:58:56 PM UTC, comment #5: 

Well, then the control package should not crash :-) "zero" and "prescale" (which is deactivated by G.scaled = true) are the only oct-files called by sigma.

Can anyone confirm this problem on windows?

Lukas Reichlin <paramaniac>
Mon 16 Dec 2013 08:13:29 PM UTC, comment #4: 

By crash I don't mean an error message but  the whole Octave application just crashes and exits; one has to restart Octave after.
I am on Octave 3.6.4 but the control version is  2.4.2 which was the default available from the windows install.
I also did zero(G), pole(G) and G.scaled=true. If there was a problem with matrix condition numbers the function would not work for the transfer function version the plant G.

Anonymous
Mon 16 Dec 2013 06:04:45 PM UTC, comment #3: 

It works for me on OS X 10.8, Octave 3.6.4 and control 2.6.1 as expected.

Which version of control are you using? Can you try


zero (G)
pole (G)


as well as


G.scaled = true


before calling sigma?

Lukas Reichlin <paramaniac>
Mon 16 Dec 2013 05:46:30 PM UTC, comment #2: 

Can you confirm when you say the function crashes, you mean it returns an error message and does not return the expected result?

According to the online help for Matlab's sigma function, it appears that it should indeed handle different system model types.

Adding the control package maintainer to the CC list.

Mike Miller <mtmiller>
Group Member
Mon 16 Dec 2013 04:03:33 AM UTC, comment #1: 

It turns out that if the system is converted to transfer function format the sigma function works without crashing.

Anonymous
Fri 13 Dec 2013 10:42:35 PM UTC, original submission:  

Below is an example of a sensitivity function G formed using the script attached. The sigma function crashes no matter what form the function call is. Is there a workaround? Scilab/Matlab have no issue generating the plots.


clear
pkg load control
a=[-0.1164, -236.8, 153.3, 8.581, 38.94,  -74.04;
   13.28, 63.75, -47.21, -1.251, -9.194, 24.45;
    21.19,  212.5, -146.4, -6.629,-32.04, 72.11;
   1.764, -29.33,  22.07, 0.7402, 5.977,-9.715;
   -20.14,  219.7, -148.9,-5.633,-47.92, 72.86;
    -0.01975, 161.9, -107,-5.119, -30.03, 51.12];
 
  b=[-0.2017   -3.976
   0.4577    -3.77
   -0.3924   -2.023
   -0.4243   0.7527
   -0.1368   -1.959
    0.1489    1.218];
 
  c=[ -0.3336   0.08691  0.004088   -0.0117   0.09318    -0.656
     -1.293    0.4011     1.066   -0.2659    0.1607   -0.7863];
 
  d=[1   0;0   1];
  G=ss(a,b,c,d);
 
  [y,w]=sigma(G)

Anonymous

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by paramaniac (Posted a comment)
  • -email is unavailable- added by mtmiller (Posted a comment)
  • -email is unavailable- added by mtmiller
  •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 8 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2014-02-24 goffioul StatusNeed Info Duplicate
        Open/ClosedOpen Closed
    2013-12-17 mtmiller Item GroupMatlab Compatibility Segfault, Bus Error, etc.
        Summarycontrol package: sigma function does not accept a state-space model control package: sigma function crashes Octave when given a state-space model
    2013-12-16 mtmiller Item GroupSegfault, Bus Error, etc. Matlab Compatibility
        StatusNone Need Info
        SummaryCall to sigma function (control toolbox) causes a crash control package: sigma function does not accept a state-space model
        Carbon-Copy- Added paramaniac

    Back to the top

    Powered by Savane 3.13-758e.
    Corresponding source code