[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