sub chkopt.for

来自「非线性回归问题SQP解法」· FOR 代码 · 共 19 行

	SUBROUTINE CHKOPT(N,ME,M,NAMAX,A,B,X,MJ,JS,US,OPTM,ZERO)
	DOUBLE PRECISION A(NAMAX,M),B(M),X(N),RES,ZERO,US
	LOGICAL MJ(M),OPTM
	JS=0
	US=-ZERO
	DO 2 I=ME+1,M
		IF (MJ(I))GOTO 2
		RES=-B(I)
		DO 1 J=1,N
			RES=RES+A(J,I)*X(J)
1		CONTINUE
		IF(RES.LT.US)THEN
			US=RES
			JS=I
		ENDIF
2	CONTINUE
	IF(JS.EQ.0) OPTM=.TRUE.
	RETURN
	END