Replace OpenMP reduction on arrays with atomic updates, since it segfaults on MacOSX.
Make create_fitsheader_for still handle scalar as 1-element 1D array. Update tests.