diff --git a/sample/simple3DFilter.cpp b/sample/simple3DFilter.cpp
index de847f9..9073546 100644
--- a/sample/simple3DFilter.cpp
+++ b/sample/simple3DFilter.cpp
@@ -1,3 +1,4 @@
+#include "openmp.hpp"
 #include "omptl/algorithm"
 #include <cassert>
 #include "yorick.hpp"
@@ -53,7 +54,7 @@ void computeInterpolatedField(MyTree *tree1, double boxsize, int Nres, double cx
       {
         double pz = (rz)*boxsize/Nres-cz;
 
-        cout << format("[%d] %d / %d") % omp_get_thread_num() % rz % Nres << endl;
+        cout << format("[%d] %d / %d") % smp_get_thread_id() % rz % Nres << endl;
         for (int ry = 0; ry < Nres; ry++)
           {
             double py = (ry)*boxsize/Nres-cy;