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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
subroutine satur (n,x,binf,bsup,d,ttmin,ttsup,topt,tg,td,
& tmi,icoi,icos,iproj)
c
c subroutine calculant ,dans un intervalle donne, un pas proche
c de tmi saturant une contrainte
c topt:pas calculer (=0 s'il n'existe pas un tel pas (s)
c ttmin,ttsup:bornes de l'intervalle dans lequel doit
c etre topt (e)
c tmi:pas au voisinnage duquel on calcul topt (e)
c iproj:indicateur de projection (e)
c =0:on cherche un pas saturant une contrainte dans
c l'intervalle ttmin,ttsup
c =1:on cherche un pas dans l'intervalle tg,td et on
c le ramene dans l'intervalle ttmin,ttsup
c par projection
c icos:indice de la variable saturee par la borne superieure
c icoi:indice de la variable saturee par la borne inferieure
c inf:indicateur pour l initialisation de icoi et icos
c =0:la borne superieure est atteinte
c =1:la borne superieure est atteinte
c =2:le pas est obtenu par projection sur ttmin ttsup
c
implicit double precision(a-h,o-z)
integer iproj
dimension x(n),binf(n),bsup(n),d(n)
c
icoi=0
icos=0
ep=tmi
c
c boucle
do 70 i=1,n
inf=0
c calcul du pas saturant la ieme contrainte:tb
if (d(i)) 61,70,62
61 tb=(binf(i)-x(i))/d(i)
inf=1
go to 63
62 tb=(bsup(i)-x(i))/d(i)
63 if ((tb.gt.ttsup).or.(tb.lt.ttmin))then
c projection de tb sur l intervalle ttmin,ttsup
if ((iproj.eq.0).or.(tb.lt.tg).or.(tb.gt.td)) go to 70
tb=max(tb,ttmin)
tb=min(tb,ttsup)
inf=2
end if
c recherche du pas le plus proche de tmi
e=abs(tb-tmi)
if (e.ge.ep) go to 70
topt=tb
ep=e
c mise a jour de icoi,icos
icoi=0
icos=0
if (inf.eq.0) icos=i
if (inf.eq.1) icoi=i
70 continue
return
end
|
