Adding polynomial features means that results are different when adding constant features to the dataset

[bejaeger]

test_constant_feature_handling failed with the default 2_6 which is probably not the desired behavior. Need to move the constant feature removal before polynomial features.