potential undefined variable bug in cohprods fixed

ralna · Jan 19, 2024 · 5534b9d · 5534b9d
1 parent efbfd8f
commit 5534b9d
Show file tree

Hide file tree

Showing 4 changed files with 7 additions and 5 deletions.
diff --git a/include/cutest_modules.h b/include/cutest_modules.h
@@ -20,7 +20,7 @@
 #define CUTEST_KINDS_precision CUTEST_KINDS_half_64
 #define CUTEST_precision CUTEST_half_64
 #define CUTEST_PROBLEM_precision CUTEST_PROBLEM_half_64
-#define CUTEST_iNTERFACE_precision CUTEST_INTERFACE_half_64
+#define CUTEST_INTERFACE_precision CUTEST_INTERFACE_half_64
 #define CUTEST_LQP_precision CUTEST_LQP_half_64
 #else
 #define CUTEST_KINDS_precision CUTEST_KINDS_half
@@ -62,7 +62,7 @@
 #define CUTEST_KINDS_precision CUTEST_KINDS_double_64
 #define CUTEST_precision CUTEST_double_64
 #define CUTEST_PROBLEM_precision CUTEST_PROBLEM_double_64
-#define CUTEST_iNTERFACE_precision CUTEST_INTERFACE_double_64
+#define CUTEST_INTERFACE_precision CUTEST_INTERFACE_double_64
 #define CUTEST_LQP_precision CUTEST_LQP_double_64
 #else
 #define CUTEST_KINDS_precision CUTEST_KINDS_double

diff --git a/src/test/test_main.F90 b/src/test/test_main.F90
@@ -77,7 +77,7 @@ PROGRAM CUTEST_test_main
       ELSE
 !       only_print_small = .TRUE.
         only_print_small = n <= 5 .AND. m <= 5
-!       only_print_small = n <= 100 .AND. m <= 100
+        only_print_small = n <= 100 .AND. m <= 100
       END IF
 
       IF ( m == 0 ) THEN

diff --git a/src/tools/cohprods.f90 b/src/tools/cohprods.f90
@@ -261,7 +261,9 @@ SUBROUTINE CUTEST_cohprods_threadsafe( data, work, status, n, goth, X,   &
 !  set the indices for the nonzeros for the objective Hessian product
 
        CALL CUTEST_cohprodsp( status, nnzohp, lohp, IND )
-
+       data%nnzohp = nnzohp
+     ELSE
+       nnzohp = data%nnzohp
      END IF
 
 !  reset IUSED and W_ws to zero

diff --git a/src/tools/cutest.F90 b/src/tools/cutest.F90
@@ -70,7 +70,7 @@ MODULE CUTEST_precision
 
       TYPE, PUBLIC :: CUTEST_data_type
         INTEGER ( KIND = ip_ ) :: n, ng, ng1, nel, nel1, ntotel, nvrels
-        INTEGER ( KIND = ip_ ) :: nnza, ngpvlu, nepvlu, ngng, out
+        INTEGER ( KIND = ip_ ) :: nnza, ngpvlu, nepvlu, ngng, nnzohp, out
         INTEGER ( KIND = ip_ ) :: nvargp, nvar2, nnonnz, nbprod, ntotin
         INTEGER ( KIND = ip_ ) :: lo, ch, lwork, la, lb, nobjgr, lu, ltypee
         INTEGER ( KIND = ip_ ) :: ltypeg, lintre, lft, lepvlu, lstep, lstgp