Make parameter controlling which "mass axpy" function to use configurable

[pelesh]

This constant (if (k < 200)) would need to become a parameter at some point (but this is a future issue).

Originally posted by @kswirydo in #110 (comment)