Simplex.cpp

Simplex.cpp

@@ -15,6 +15,7 @@
 //
 
 #include <iostream>
+#include <iomanip>
 #include <cfloat>
 #include <cmath>
 #include <limits>

Simplex.h

@@ -74,6 +74,8 @@ public:
   void MultiContract();
   /// Update the indicies for best, worst and 2nd to worst points.
   void UpdateRankIndicies();
+  // Prints out the points of the simplex and the associated function values
+  friend ostream& operator<<(ostream& op, const Simplex& smplx); 
 private:
   const MISCMATHS::NonlinCF&                           _cf;
   const NEWMAT::ColumnVector                           _sp;
@@ -88,6 +90,18 @@ private:
   void calculate_reflexion_point(unsigned int ii);
 };
 
+// Global function, friend of Simplex
+ostream& operator<<(ostream& op, const Simplex& smplx)
+{
+  op << "Simplex = " << endl;
+  for (unsigned int i=0; i<smplx._smx.size(); i++) {
+    for (unsigned int j=1; j<smplx._smx.size(); j++) {
+      op << std::setw(10) << std::setprecision(4) << smplx._smx[i](j); 
+    }
+    op << " | " << std::setw(10) << std::setprecision(4) << smplx._fv[i] << endl;
+  }
+}
+
 } // End namespace MISCMATHS
 
 #endif // End #ifndef Simplex_h