1 too many bvals, `grad * grad == (bvec * bval) * (bvec * bval)`. You could do `grad * bvec` but that doesn't seem any better to me.