File: smga.sci

package info (click to toggle)
scilab 2.4-1
  • links: PTS
  • area: non-free
  • in suites: potato, slink
  • size: 55,196 kB
  • ctags: 38,019
  • sloc: ansic: 231,970; fortran: 148,976; tcl: 7,099; makefile: 4,585; sh: 2,978; csh: 154; cpp: 101; asm: 39; sed: 5
file content (29 lines) | stat: -rw-r--r-- 676 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
 
function [a,b1,b2,c1,c2,d11,d12,d21,d22]=smga(slp,r)
//Matrix extractions (for use with h_inf,lft and other);
// Utility function
// Copyright INRIA
slp1=slp(1);
if slp1(1)=='lss' then
  p=r(2),r=r(1),
  [s1,s2,t]=size(slp);
  [a,b,c,d]=slp(2:5),
  b1=b(:,1:s2-p),b2=b(:,s2-p+1:s2),
  c1=c(1:s1-r,:),c2=c(s1-r+1:s1,:),
  d11=coeff(d(1:s1-r,1:s2-p),0),
  d12=coeff(d(1:s1-r,s2-p+1:s2),0),
  d21=coeff(d(s1-r+1:s1,1:s2-p),0),
  d22=coeff(d(s1-r+1:s1,s2-p+1:s2),0)
end
if slp1(1)=='r' then
  [nl,nk]=size(slp);
  k1=1:nl-r(1);
  k2=nl-r(1)+1:nl;
  m1=1:nk-r(2);
  m2=nk-r(2)+1:nk;
  a=slp(k1,m1);  //p11
  b1=slp(k1,m2); //p12
  b2=slp(k2,m1); //p21
  c1=slp(k2,m2); //p22
end