sub pval.for

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

	SUBROUTINE PVAL(ME,MI,M,OFUN,CFUN,PPRA,PFUN)
	DOUBLE PRECISION OFUN,CFUN(M),PPRA,PFUN
	PFUN=0.0D0
	DO 1 I=1,ME
		PFUN=PFUN+DABS(CFUN(I))
1	CONTINUE
      DO 2 I=1,MI
		PFUN=PFUN+DABS(DMIN1(0.0D0,CFUN(ME+I)))
2	CONTINUE
	PFUN=OFUN+PPRA*PFUN
	RETURN
	END