# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1219845595 14400
# Node ID e3e3d12364b01734846c43745bbc63128a93ef28
# Parent  c476c2214a9014bea5529ff05603a3944ad81751
make glpk return NA values for unfeasible problems

diff --git a/src/ChangeLog b/src/ChangeLog
--- a/src/ChangeLog
+++ b/src/ChangeLog
@@ -1,3 +1,8 @@
+2008-08-27  Jaroslav Hajek  <highegg@gmail.com>
+
+	* DLD-FUNCTIONS/__glpk__.cc (F__glpk__): Initialize output vectors 
+	by NA value.
+
 2008-08-26 Ben Abbott <bpabott@mac.com>
 
 	* mappers.cc: Adjust tolerance for single precesion tests.
diff --git a/src/DLD-FUNCTIONS/__glpk__.cc b/src/DLD-FUNCTIONS/__glpk__.cc
--- a/src/DLD-FUNCTIONS/__glpk__.cc
+++ b/src/DLD-FUNCTIONS/__glpk__.cc
@@ -801,11 +801,11 @@
   OCTAVE_GLPK_GET_REAL_PARAM ("tolobj", 9);
 
   //-- Assign pointers to the output parameters
-  ColumnVector xmin (mrowsc);
-  ColumnVector fmin (1);
+  ColumnVector xmin (mrowsc, octave_NA);
+  ColumnVector fmin (1, octave_NA);
   ColumnVector status (1);
-  ColumnVector lambda (mrowsA);
-  ColumnVector redcosts (mrowsc);
+  ColumnVector lambda (mrowsA, octave_NA);
+  ColumnVector redcosts (mrowsc, octave_NA);
   ColumnVector time (1);
   ColumnVector mem (1);