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

 
 

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

Submitted by:  None
Submitted on:  Fri 13 Dec 2013 10:42:35 PM UTC  
 
Category: Octave Forge PackageSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Segfault, Bus Error, etc.
Status: DuplicateAssigned to: None
Originator Name: Saul MtakulaOriginator Email: -unavailable-
Open/Closed: ClosedRelease: 3.6.4
Operating System: Microsoft Windows

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

(Jump to the original submission 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>
Project Member
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>
Project Administrator
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>
Project Administrator
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

as well as

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>
Project Administrator
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 File(s):
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

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

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 8 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Mon 24 Feb 2014 03:47:24 AM UTCgoffioulStatusNeed Info=>Duplicate
      Open/ClosedOpen=>Closed
    Tue 17 Dec 2013 03:55:03 AM UTCmtmillerItem 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
    Mon 16 Dec 2013 05:46:30 PM UTCmtmillerItem 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.1-cleanup1