************************************************************************ C PROGRAM TNECU C C TEST PROGRAM FOR THE SUBROUTINE PNECU C INTEGER NF,IH(1001),JH(8000),IPAR(7),ISPAS,IPRNT,ITERM DOUBLE PRECISION X(1000),RPAR(9),F,GMAX INTEGER NEXT,IERR,I,ITIME,M INTEGER NITER,NFVAL,NGVAL,NSUCC,NITCG COMMON /PROB/ NEXT INTEGER NRES,NDEC,NIN,NIT,NFV,NFG,NFH COMMON /STAT/ NRES,NDEC,NIN,NIT,NFV,NFG,NFH NITER=0 NFVAL=0 NGVAL=0 NSUCC=0 NITCG=0 CALL TYTIM1(ITIME) C C LOOP FOR 22 TEST PROBLEMS C DO 3 NEXT=1,22 C C CHOICE OF INTEGER AND REAL PARAMETERS C DO 1 I=1,7 IPAR(I)=0 1 CONTINUE DO 2 I=1,9 RPAR(I)=0.0D 0 2 CONTINUE ISPAS=2 IPRNT=1 C C PROBLEM DIMENSION C NF=1000 C C INITIATION OF X AND CHOICE OF RPAR(1) AND RPAR(6) C CALL TIUS14(NF,M,X,IH,JH,RPAR(6),RPAR(1),NEXT,IERR) IF (IERR.NE.0) GO TO 3 IF (RPAR(6).EQ.0.0D0) IPAR(4)=1 RPAR(1)=0.0D0 IF (NEXT.EQ.1) RPAR(1)=4.0D0 C C SOLUTION C CALL PNECU(NF,M,X,IH,JH,IPAR,RPAR,F,GMAX,ISPAS,IPRNT,ITERM) NITER=NITER+NIT NFVAL=NFVAL+NFV NGVAL=NGVAL+NFG NITCG=NITCG+NIN IF (ITERM.GT.0.AND.ITERM.LT.9) NSUCC=NSUCC+1 3 CONTINUE WRITE(6,10) NITER,NFVAL,NGVAL,NITCG,NSUCC 10 FORMAT(' NITER =',I5,3X,' NFVAL =',I5,3X,' NGVAL =',I5,3X, & ' NITCG =',I5,3X,' NSUCC =',I5) CALL TYTIM2(ITIME) STOP END C C USER SUPPLIED SUBROUTINE (CALCULATION OF FF) C SUBROUTINE OBJ(NF,X,FF) INTEGER NF DOUBLE PRECISION X(*),FF INTEGER NEXT COMMON /PROB/ NEXT C C FUNCTION EVALUATION C CALL TFFU14(NF,X,FF,NEXT) RETURN END C C USER SUPPLIED SUBROUTINE (CALCULATION OF GF) C SUBROUTINE DOBJ(NF,X,GF) INTEGER NF DOUBLE PRECISION X(*),GF(*) INTEGER NEXT COMMON /PROB/ NEXT C C GRADIENT EVALUATION C CALL TFGU14(NF,X,GF,NEXT) RETURN END C C EMPTY SUBROUTINES C SUBROUTINE FUN(NF,KA,X,FA) INTEGER NF,KA DOUBLE PRECISION X(*),FA KA=NF FA=X(1) RETURN END SUBROUTINE DFUN(NF,KA,X,GA) INTEGER NF,KA DOUBLE PRECISION X(*),GA(*) KA=NF GA(1)=X(1) RETURN END ************************************************************************ C PROGRAM TNECS C C TEST PROGRAM FOR THE SUBROUTINE PNECS C INTEGER NF,IX(1000),IH(1001),JH(8000),IPAR(7),ISPAS,IPRNT,ITERM DOUBLE PRECISION X(1000),XL(1000),XU(1000),RPAR(9),F,GMAX INTEGER NEXT,IERR,I,ITIME,M INTEGER NITER,NFVAL,NGVAL,NSUCC,NITCG COMMON /PROB/ NEXT INTEGER NRES,NDEC,NIN,NIT,NFV,NFG,NFH COMMON /STAT/ NRES,NDEC,NIN,NIT,NFV,NFG,NFH NITER=0 NFVAL=0 NGVAL=0 NSUCC=0 NITCG=0 CALL TYTIM1(ITIME) C C LOOP FOR 22 TEST PROBLEMS C DO 3 NEXT=1,22 C C CHOICE OF INTEGER AND REAL PARAMETERS C DO 1 I=1,7 IPAR(I)=0 1 CONTINUE DO 2 I=1,9 RPAR(I)=0.0D 0 2 CONTINUE ISPAS=2 IPRNT=1 C C PROBLEM DIMENSION C NF=1000 C C INITIATION OF X AND CHOICE OF RPAR(1) AND RPAR(6) C CALL TIUS14(NF,M,X,IH,JH,RPAR(6),RPAR(1),NEXT,IERR) IF (IERR.NE.0) GO TO 3 IF (RPAR(6).EQ.0.0D 0) IPAR(4)=1 CALL MXVINS(NF,3,IX) CALL MXVSET(NF,-1.0D 0,XL) CALL MXVSET(NF, 1.0D 0,XU) RPAR(1)=0.0D 0 IF (NEXT.EQ. 2) RPAR(1)=5.0D 1 IF (NEXT.EQ.12) RPAR(1)=1.0D 0 IF (NEXT.EQ.18) RPAR(1)=1.0D 1 C C SOLUTION C CALL PNECS(NF,M,X,IX,XL,XU,IH,JH,IPAR,RPAR,F,GMAX,ISPAS,IPRNT, & ITERM) NITER=NITER+NIT NFVAL=NFVAL+NFV NGVAL=NGVAL+NFG NITCG=NITCG+NIN IF (ITERM.GT.0.AND.ITERM.LT.9) NSUCC=NSUCC+1 3 CONTINUE WRITE(6,10) NITER,NFVAL,NGVAL,NITCG,NSUCC 10 FORMAT(' NITER =',I5,3X,' NFVAL =',I5,3X,' NGVAL =',I5,3X, & ' NITCG =',I5,3X,' NSUCC =',I5) CALL TYTIM2(ITIME) STOP END C C USER SUPPLIED SUBROUTINE (CALCULATION OF FF) C SUBROUTINE OBJ(NF,X,FF) INTEGER NF DOUBLE PRECISION X(*),FF INTEGER NEXT COMMON /PROB/ NEXT C C FUNCTION EVALUATION C CALL TFFU14(NF,X,FF,NEXT) RETURN END C C USER SUPPLIED SUBROUTINE (CALCULATION OF GF) C SUBROUTINE DOBJ(NF,X,GF) INTEGER NF DOUBLE PRECISION X(*),GF(*) INTEGER NEXT COMMON /PROB/ NEXT C C GRADIENT EVALUATION C CALL TFGU14(NF,X,GF,NEXT) RETURN END C C EMPTY SUBROUTINES C SUBROUTINE FUN(NF,KA,X,FA) INTEGER NF,KA DOUBLE PRECISION X(*),FA KA=NF FA=X(1) RETURN END SUBROUTINE DFUN(NF,KA,X,GA) INTEGER NF,KA DOUBLE PRECISION X(*),GA(*) KA=NF GA(1)=X(1) RETURN END 