sub qpout.for

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

	SUBROUTINE QPOUT(N,M,NAMAX,MAMAX,G,C,A,B,V,OSCL,SCL)
	DOUBLE PRECISION G(NAMAX,N),A(NAMAX,M),C(N),B(M),SCL(MAMAX),V(M),
	*                 OSCL,ZERO
	DO 1 J=1,N
		C(J)=C(J)*OSCL
		DO 1 I=1,N
			G(I,J)=G(I,J)*OSCL
1		CONTINUE
	DO 2 I=1,M
		B(I)=B(I)*SCL(I)
		DO 2 J=1,N
			A(J,I)=A(J,I)*SCL(I)
2	CONTINUE
	DO 3 I=1,M
	V(I)=V(I)*OSCL/SCL(I)
3	CONTINUE
	END