Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- *
- * INITIAL OPERATIONS WITH SIMPLE BOUNDS
- *
- IF (KBF.GT.0) THEN
- DO 2 I = 1,NF
- IF ((IX(I).EQ.3.OR.IX(I).EQ.4) .AND. XU(I).LE.XL(I)) THEN
- XU(I) = XL(I)
- IX(I) = 5
- ELSE IF (IX(I).EQ.5 .OR. IX(I).EQ.6) THEN
- XL(I) = X(I)
- XU(I) = X(I)
- IX(I) = 5
- END IF
- 2 CONTINUE
- CALL PCBS04(NF,X,IX,XL,XU,EPS9,KBF)
- CALL PYADC0(NF,N,X,IX,XL,XU,INEW)
- END IF
- CALL OBJ(NF,X,F)
- NFV=NFV+1
- CALL DOBJ(NF,X,GF)
- NFG=NFG+1
- LD=KD
- 11020 CONTINUE
- CALL PYTRCG(NF,NF,IX,GF,UMAX,GMAX,KBF,IOLD)
- CALL MXVCOP(NF,GF,GN)
- IF (ABS(IPRNT).GT.1)
- & WRITE (6,'(1X,''NIT='',I5,2X,''NFV='',I5,2X,''NFG='',I5,2X,
- & ''F='', G16.9,2X,''G='',E10.3)') NIT,NFV,NFG,F,GMAX
- CALL PYFUT1(NF,F,FO,UMAX,GMAX,DMAX,TOLX,TOLF,TOLB,TOLG,KD,NIT,KIT,
- & MIT,NFV,MFV,NFG,MFG,NTESX,MTESX,NTESF,MTESF,ITES,IRES1,IRES2,
- & IREST,ITERS,ITERM)
- IF (ITERM.NE.0) GO TO 11080
- IF (KBF.GT.0) THEN
- CALL PYRMC0(NF,N,IX,GN,EPS8,UMAX,GMAX,RMAX,IOLD,IREST)
- IF (UMAX.GT.EPS8*GMAX) IREST=MAX(IREST,1)
- END IF
- CALL MXVCOP(NF,X,XO)
- 11040 CONTINUE
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement