* Bug #8190 fixed - Optimization: Fixed ICSE demos
[scilab.git] / scilab / modules / optimization / demos / icse / icsez0.f
1 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 c Copyright (C) INRIA
3 c
4 c This file must be used under the terms of the CeCILL.
5 c This source file is licensed as described in the file COPYING, which
6 c you should have received as part of this distribution.  The terms
7 c are also available at
8 c http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
9 c
10       subroutine icsez0(nu,t0,tf,dti,dtf,ermx,iu,nuc,nuv,ilin,nti,ntf,
11      & ny,nea,itmx,nex,nob,ntob,ntobi,nitu,ndtu,nitv,nrtv,ndtv)
12 c
13 c     programme d'initialisation appele par icse.bas
14 c     initialisation des commons icsez icsez0 et nird
15 c
16       implicit double precision (a-h,o-z)
17       dimension iu(5),iu0(5)
18       common/icsez/t00,tf0,dti0,dtf0,ermx0,iu0,nuc0,nuv0,ilin0,nti0,
19      & ntf0,ny0,nea0,itmx0,nex0,nob0,ntob0,ntobi0,nitu0,ndtu0
20       common/nird/nitv0,nrtv0,ndtv0
21 c
22       t00=t0
23       tf0=tf
24       dti0=dti
25       dtf0=dtf
26       ermx0=ermx
27       do 10 i=1,5
28 10    iu0(i)=iu(i)
29       nuc0=nuc
30       nuv0=nuv
31       ilin0=ilin
32       nti0=nti
33       ntf0=ntf
34       ny0=ny
35       nea0=nea
36       itmx0=itmx
37       nex0=nex
38       nob0=nob
39       ntob0=ntob
40       ntobi0=ntobi
41       nitu0=nitu
42       ndtu0=ndtu
43       nitv0=nitu+ny+ntob
44       nrtv0=0
45       ldtvt=ny*(nob+nuc+nuv+nti+ntf+ny+3)+ntob*(nob*nex+nob+1)+2*nu+
46      &ndtu
47       mdtv1=ldtvt+ny*(ny+7)
48       nui=iu(1)*nuc+iu(2)*nuv*(nti+ntf+1)
49       mdtv2=ldtvt+ny*(2*ntob+nu+2*ny+nuc+nuv+4)+max(nuc+nuv,nui)+
50      &nob*ntob+nob
51       ndtv0=max(mdtv1,mdtv2)
52       !ind=0
53       !call icse(ind,nu,zz,zz,zz,zz,zz,zz,zz,zz,zz)
54       nitv=max(1,nitv0)
55       nrtv=max(1,nrtv0)
56       ndtv=max(1,ndtv0)
57       return
58       end