diff --git a/libs/tbdefs.hpp b/libs/tbdefs.hpp
index 233c3104..5524195c 100644
--- a/libs/tbdefs.hpp
+++ b/libs/tbdefs.hpp
@@ -21,6 +21,20 @@
 #include <complex>
 #include <time.h>
 
+template<typename T>
+std::ostream& operator<<(std::ostream& os, const std::valarray<T>& v)
+{
+    os << "[";
+    if (v.size() > 0) {
+        os << v[0];
+        for (size_t i = 1; i < v.size(); ++i) {
+            os << ", " << v[i];
+        }
+    }
+    os << "]";
+    return os;
+}
+
 /**************************************
  *          MACROS & DEFINES          *
  **************************************/
diff --git a/make.in b/make.in
index cbb9a62f..4cb3eb06 100644
--- a/make.in
+++ b/make.in
@@ -9,7 +9,7 @@ MAKEDEPEND=g++ -M $(CXXFLAGS)
 # Lapack includes and linking
 
 ILAPACK=
-LLAPACK=-llapack
+LLAPACK=-lopenblas
 
 TBSRC=../libs/
 LDFLAGS=-L$(TBSRC) -ltoolbox $(LLAPACK)
diff --git a/tools/dimproj.cpp b/tools/dimproj.cpp
index bdef5421..cdba4460 100644
--- a/tools/dimproj.cpp
+++ b/tools/dimproj.cpp
@@ -124,7 +124,9 @@ int main(int argc, char**argv)
     else
     {
       csv2floats(fdhd,tfpars);  fhdpars=tfpars;
-      std::cerr<<"high-dim pars"<<tfpars<<"\n";
+      std::cerr << "high-dim pars";
+      toolbox::operator<<(std::cerr, tfpars);
+      std::cerr <<"\n";
       if (tfpars.size()==2)
       {
         opts.tfunH.set_mode(NLDRGamma,tfpars);
@@ -142,7 +144,9 @@ int main(int argc, char**argv)
     else
     {
       csv2floats(fdld,tfpars); fldpars=tfpars;
-      std::cerr<<"lo-dim pars"<<tfpars<<"\n";
+      std::cerr<<"lo-dim pars";
+      toolbox::operator<<(std::cerr, tfpars);
+      std::cerr <<"\n";
       if (tfpars.size()==2)
       {
         opts.tfunL.set_mode(NLDRGamma,tfpars);
diff --git a/tools/dimred.cpp b/tools/dimred.cpp
index d591611d..d86b1f43 100644
--- a/tools/dimred.cpp
+++ b/tools/dimred.cpp
@@ -140,8 +140,10 @@ int main(int argc, char**argv)
     else 
     {
       csv2floats(fdhd,tfpars);  fhdpars=tfpars; 
-      std::cerr<<"high-dim pars"<<tfpars<<"\n";     
-      if (tfpars.size()==2) 
+      std::cerr << "high-dim pars";
+      toolbox::operator<<(std::cerr, tfpars);
+      std::cerr <<"\n";
+      if (tfpars.size()==2)
       {
         iteropts.tfunH.set_mode(NLDRGamma,tfpars);       
       }
@@ -160,8 +162,10 @@ int main(int argc, char**argv)
     else 
     {
       csv2floats(fdld,tfpars); fldpars=tfpars;      
-      std::cerr<<"lo-dim pars"<<tfpars<<"\n";     
-      if (tfpars.size()==2) 
+      std::cerr<<"lo-dim pars";
+      toolbox::operator<<(std::cerr, tfpars);
+      std::cerr <<"\n";
+      if (tfpars.size()==2)
       {
         iteropts.tfunL.set_mode(NLDRGamma,tfpars);  
       }
@@ -200,7 +204,9 @@ int main(int argc, char**argv)
 
     std::cerr<<"hey "<<itermode<<" "<<iteropts.minmode<<"\n";
     std::valarray<double> sat(0.0,2); csv2floats(tempopts,sat);
-    std::cerr<<" simulated annealing ops: "<<sat<<"\n";
+    std::cerr << " simulated annealing ops: ";
+    toolbox::operator<<(std::cerr, sat);
+    std::cerr<<"\n";
     iteropts.saopts.temp_init=sat[0]; iteropts.saopts.temp_final=sat[1];
     iteropts.ptopts.temp_init=sat[0]; iteropts.ptopts.temp_final=sat[1]; 
     iteropts.ptopts.temp_factor=ptfac; iteropts.ptopts.replica=npt;