1
// icse.init.bas : initialisation et tests pour icse
2
//**************************************************
5
// creation y0 (etat initial)
7
if exists('y0')==0, y0=ones(1,ny); end;
8
// creation de b,fy et fu et passage en vecteur ligne
9
if exists('b')==0, b=ones(1,ny); end;
10
if exists('fy')==0, fy1=ones(1,ny*ny);
11
else fy1=matrix(fy,1,ny*ny) ; end;
12
if exists('fu')==0, fu1=ones(1,ny*(nuc+nuv));
13
else fu1=matrix(fu,1,ny*(nuc+nuv)); end;
19
if xx+yy>2, error('itu et/ou dtu n'est pas un vecteur ligne'), end;
20
u=[]; if nuc>0, u=uc; end;
21
if nuv>0, u=[u,uv]; end;
22
nu=nuc+nuv*(nti+ntf+1);
23
if size(u)<>[1,nu],error('dimensions du controle incompatibles'),end;
25
// initialisation du common icsez
26
[nitv,nrtv,ndtv]=fort('icse0',nu,1,'i',t0,2,'d',tf,3,'d',dti,4,'d',..
27
dtf,5,'d',ermx,6,'d',iu,7,'i',nuc,8,'i',nuv,9,'i',ilin,10,'i',nti,..
28
11,'i',ntf,12,'i',ny,13,'i',nea,14,'i',itmx,15,'i',nex,16,'i',nob,..
29
17,'i',ntob,18,'i',ntobi,19,'i',nitu,20,'i',ndtu,21,'i','sort',..
30
[1,1],22,'i',[1,1],23,'i',[1,1],24,'i');