File: of.sci

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (35 lines) | stat: -rw-r--r-- 1,209 bytes parent folder | download | duplicates (4) 
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
30
31
32
33
34
35
 
 function [P,Q]=of(A,B,C)
// Copyright INRIA
// Generated by lmitool on Thu Feb 09 11:22:01 MET 1995
   
   Mbound = 1e3;
   abstol = 1e-10;
   nu = 10;
   maxiters = 100;
   reltol = 1e-10;
   options=[Mbound,abstol,nu,maxiters,reltol];
    
 ///////////DEFINE INITIAL GUESS AND PRELIMINARY CALCULATIONS BELOW
 P_init=zeros(A)                                                   
 Q_init=zeros(A)                                                   
 NB=kernel(B');                                                    
 NC=kernel(C);                                                     
 I=eye(A);                                                         
 /////////// 
  
 XLIST0=list(P_init,Q_init)
 XLIST=lmisolver(XLIST0,of_eval,options)
 [P,Q]=XLIST(:)
  
  
  
 /////////////////EVALUATION FUNCTION////////////////////////////
  
 function [LME,LMI,OBJ]=of_eval(XLIST)
 [P,Q]=XLIST(:)
  
 /////////////////DEFINE LME, LMI and OBJ BELOW
 LME=list(P-P',Q-Q')                                               
 LMI= list(-NB'*(A*Q+Q*A'+Q)*NB,-NC'*(A'*P+P*A+P)*NC,[P I; I Q]);  
 OBJ=trace(P+Q);