[Scipy-svn] r2093 - trunk/Lib/optimize/cobyla
scipy-svn at scipy.org
scipy-svn at scipy.org
Thu Jul 13 04:50:54 EDT 2006
Author: oliphant
Date: 2006-07-13 03:50:51 -0500 (Thu, 13 Jul 2006)
New Revision: 2093
Modified:
trunk/Lib/optimize/cobyla/cobyla2.f
Log:
more printing
Modified: trunk/Lib/optimize/cobyla/cobyla2.f
===================================================================
--- trunk/Lib/optimize/cobyla/cobyla2.f 2006-07-13 08:50:16 UTC (rev 2092)
+++ trunk/Lib/optimize/cobyla/cobyla2.f 2006-07-13 08:50:51 UTC (rev 2093)
@@ -367,9 +367,19 @@
ISDIRN=IVMC+MP
IDXNEW=ISDIRN+N
IVMD=IDXNEW+N
+ IF (IPRINT .EQ. 3) THEN
+ print *, 'BEFORE trstlp:'
+ PRINT *, ' **DX = ', (DX(I),I=1,N)
+ PRINT *, ' **IACT = ', (IACT(I),I=1,M+1)
+ PRINT *, 'M,N,RHO,IFULL =', M, N, RHO, IFULL
+ PRINT *, ' **CON = ', (CON(I),I=1,M)
+ PRINT *, ' **A = ', ((A(I,K),I=1,N),K=1,MP)
+ PRINT *, ' **W = ', (W(I),I=1,ITOTAL)
+ END IF
CALL TRSTLP (N,M,A,CON,RHO,DX,IFULL,IACT,W(IZ),W(IZDOTA),
1 W(IVMC),W(ISDIRN),W(IDXNEW),W(IVMD))
IF (IPRINT .EQ. 3) THEN
+ print *, 'AFTER trstlp:'
PRINT *, ' **DX = ', (DX(I),I=1,N)
PRINT *, ' **IACT = ', (IACT(I),I=1,M+1)
PRINT *, 'M,N,RHO,IFULL =', M, N, RHO, IFULL
More information about the Scipy-svn
mailing list