Store global, not local column indices in each worker process.
Use print function from __future__; move statements out of per query loop.
Add parameter to_n_p_values to BH_array to implment an approximate Benjamini-Hochberg method.
Get the number of p values from each worker process in process row 0 only; print times from both process 0 and process 1.