--- branches/development/src/brains/Thermo.cpp	2012/07/03 18:32:27	1764
+++ branches/development/src/brains/Thermo.cpp	2012/09/13 14:10:11	1798
@@ -54,8 +54,10 @@
 #include "types/FixedChargeAdapter.hpp"
 #include "types/FluctuatingChargeAdapter.hpp"
 #include "types/MultipoleAdapter.hpp"
+#ifdef HAVE_QHULL
 #include "math/ConvexHull.hpp"
 #include "math/AlphaHull.hpp"
+#endif
 
 using namespace std;
 namespace OpenMD {
@@ -370,12 +372,8 @@ namespace OpenMD {
             pCount++;
           }
           
-          MultipoleAdapter ma = MultipoleAdapter(atom->getAtomType());
-          if (ma.isDipole() ) {
-            Vector3d u_i = atom->getElectroFrame().getColumn(2);
-            moment = ma.getDipoleMoment();
-            moment *= debyeToCm;
-            dipoleVector += u_i * moment;
+          if (atom->isDipole()) {
+            dipoleVector += atom->getDipole() * debyeToCm;
           }
         }
       }
@@ -830,9 +828,9 @@ namespace OpenMD {
           data[0] = pos1.x();
           data[1] = pos1.y();
           data[2] = pos1.z();          
-          MPI_Bcast(data, 3, MPI_REALTYPE, proc1, MPI_COMM_WORLD);
+          MPI::COMM_WORLD.Bcast(data, 3, MPI::REALTYPE, proc1);
         } else {
-          MPI_Bcast(data, 3, MPI_REALTYPE, proc1, MPI_COMM_WORLD);
+          MPI::COMM_WORLD.Bcast(data, 3, MPI::REALTYPE, proc1);
           pos1 = Vector3d(data);
         }
 
@@ -841,10 +839,10 @@ namespace OpenMD {
           pos2 = sd2->getPos();
           data[0] = pos2.x();
           data[1] = pos2.y();
-          data[2] = pos2.z();          
-          MPI_Bcast(data, 3, MPI_REALTYPE, proc2, MPI_COMM_WORLD);
+          data[2] = pos2.z();  
+          MPI::COMM_WORLD.Bcast(data, 3, MPI::REALTYPE, proc2);
         } else {
-          MPI_Bcast(data, 3, MPI_REALTYPE, proc2, MPI_COMM_WORLD);
+          MPI::COMM_WORLD.Bcast(data, 3, MPI::REALTYPE, proc2);
           pos2 = Vector3d(data);
         }
 #else
@@ -864,9 +862,9 @@ namespace OpenMD {
 
   RealType Thermo::getHullVolume(){
     Snapshot* snap = info_->getSnapshotManager()->getCurrentSnapshot();
-    
-    if (!snap->hasHullVolume) {
 
+#ifdef HAVE_QHULL    
+    if (!snap->hasHullVolume) {
       Hull* surfaceMesh_;
 
       Globals* simParams = info_->getSimParams();
@@ -902,5 +900,8 @@ namespace OpenMD {
       snap->setHullVolume(surfaceMesh_->getVolume());
     }
     return snap->getHullVolume();
-  }  
+#else
+    return 0.0;
+#endif
+  }
 }