sub rtest.for

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

	SUBROUTINE RTEST(M,ME,NJ,V,RV,JSET,T,JK,JP,ZERO)
	DOUBLE PRECISION V(M),RV(M),T,ZERO,RATIO,LARGE
	INTEGER JSET(M)
	DATA LARGE/1.0D+20/
	JK=0
	T=LARGE
	DO 1 I=ME+1,NJ
		L=JSET(I)
		IF (RV(I).LT.ZERO)GOTO 1
		RATIO=V(L)/RV(I)
		IF (RATIO.LT.T)THEN
			T=RATIO
			JK=L
			JP=I
	ENDIF
1	CONTINUE
	END