--- branches/development/src/parallel/ForceMatrixDecomposition.cpp	2011/07/09 15:05:59	1588
+++ branches/development/src/parallel/ForceMatrixDecomposition.cpp	2011/07/10 16:05:34	1589
@@ -112,10 +112,16 @@ namespace OpenMD {
     AtomCommIntRow->gather(idents, identsRow);
     AtomCommIntColumn->gather(idents, identsCol);
     
-    vector<int>::iterator it;
-    for (it = AtomLocalToGlobal.begin(); it != AtomLocalToGlobal.end(); ++it) {
-      cerr << "my AtomLocalToGlobal = " << (*it) << "\n";
-    }
+    // allocate memory for the parallel objects
+    AtomRowToGlobal.resize(nAtomsInRow_);
+    AtomColToGlobal.resize(nAtomsInCol_);
+    cgRowToGlobal.resize(nGroupsInRow_);
+    cgColToGlobal.resize(nGroupsInCol_);
+    massFactorsRow.resize(nAtomsInRow_);
+    massFactorsCol.resize(nAtomsInCol_);
+    pot_row.resize(nAtomsInRow_);
+    pot_col.resize(nAtomsInCol_);
+
     AtomCommIntRow->gather(AtomLocalToGlobal, AtomRowToGlobal);
     AtomCommIntColumn->gather(AtomLocalToGlobal, AtomColToGlobal);