Verify results match those of numpy