Enabled cheap parallelization.

1 files changed, 2 insertions, 0 deletions

diff --git a/uniform.cpp b/uniform.cpp
index 474399b..08f190a 100644
--- a/uniform.cpp
+++ b/uniform.cpp
@@ -37,6 +37,7 @@ public:
 
   void computeWeights() {
     for (std::vector<Face>& faces : lattices) {
+#pragma omp parallel for
       for (Face& f : faces) {
         Real w = f[0].get().weights.top();
         Real x = f[1].get().weights.top();
@@ -60,6 +61,7 @@ public:
 
   void computeProbabilities() { // destroys *all* weights
     for (auto it = lattices.rbegin(); it != lattices.rend(); it++) {
+#pragma omp parallel for
       for (Face& f : *it) {
         Real p = f[0].get().probability;
         Real q = f[1].get().probability;