I don't remember exactly from where I based my code, but this http://www.math.hkbu.edu.hk/parallel/pgi/doc/pgC++_lib/stdlibug/inp_5394.htm give some explanations that could help. Maybe I have added an additional level :(