This was inconsistent anyway: 
With CBLAS, the function computed a * gf2 + gf1 -> gf1
Without, it computed gf2 + a * gf1 -> gf1