I /think/ this could be vectorized equivalently as:
```
idx_points_nm = np.array([dcos * fp[:, 0] - dsin * fp[:, 1] + dt,
                          dsin * fp[:, 0] + dcos * fp[:, 1],
                          fp[:, 2]], order='F', dtype='>f8').T
```
Should be faster.