sub qpval.for

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

	SUBROUTINE QPVAL(N,NA,NAMAX,OGRA,B,D1,PPRA,QPOBJ,FEASBL)
	DOUBLE PRECISION OGRA(N),B(NAMAX,N),D1(N),PPRA,QPOBJ
	LOGICAL FEASBL
	QPOBJ=0.0D0
	DO 1 J=1,N
		QPOBJ=QPOBJ+OGRA(J)*D1(J)
		DO 1 I=1,N
			QPOBJ=QPOBJ+5.0D-1*D1(I)*B(I,J)*D1(J)
1	CONTINUE
	IF (FEASBL)RETURN
	DO 2 J=N+1,NA
		QPOBJ=QPOBJ+PPRA*D1(J)
2	CONTINUE
	RETURN
	END