Advertisement
Guest User

Untitled

a guest
Oct 5th, 2017
83
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. *
  2. *     INITIAL OPERATIONS WITH SIMPLE BOUNDS
  3. *
  4.       IF (KBF.GT.0) THEN
  5.       DO 2 I = 1,NF
  6.       IF ((IX(I).EQ.3.OR.IX(I).EQ.4) .AND. XU(I).LE.XL(I)) THEN
  7.       XU(I) = XL(I)
  8.       IX(I) = 5
  9.       ELSE IF (IX(I).EQ.5 .OR. IX(I).EQ.6) THEN
  10.       XL(I) = X(I)
  11.       XU(I) = X(I)
  12.       IX(I) = 5
  13.       END IF
  14.     2 CONTINUE
  15.       CALL PCBS04(NF,X,IX,XL,XU,EPS9,KBF)
  16.       CALL PYADC0(NF,N,X,IX,XL,XU,INEW)
  17.       END IF
  18.       CALL OBJ(NF,X,F)
  19.       NFV=NFV+1
  20.       CALL DOBJ(NF,X,GF)
  21.       NFG=NFG+1
  22.       LD=KD
  23. 11020 CONTINUE
  24.       CALL PYTRCG(NF,NF,IX,GF,UMAX,GMAX,KBF,IOLD)
  25.       CALL MXVCOP(NF,GF,GN)
  26.       IF (ABS(IPRNT).GT.1)
  27.      & WRITE (6,'(1X,''NIT='',I5,2X,''NFV='',I5,2X,''NFG='',I5,2X,
  28.     & ''F='', G16.9,2X,''G='',E10.3)') NIT,NFV,NFG,F,GMAX
  29.       CALL PYFUT1(NF,F,FO,UMAX,GMAX,DMAX,TOLX,TOLF,TOLB,TOLG,KD,NIT,KIT,
  30.      & MIT,NFV,MFV,NFG,MFG,NTESX,MTESX,NTESF,MTESF,ITES,IRES1,IRES2,
  31.      & IREST,ITERS,ITERM)
  32.       IF (ITERM.NE.0) GO TO 11080
  33.       IF (KBF.GT.0) THEN
  34.       CALL PYRMC0(NF,N,IX,GN,EPS8,UMAX,GMAX,RMAX,IOLD,IREST)
  35.       IF (UMAX.GT.EPS8*GMAX) IREST=MAX(IREST,1)
  36.       END IF
  37.       CALL MXVCOP(NF,X,XO)
  38. 11040 CONTINUE
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement