Do you have explicit blas calls? That won't work. Instead you should just use the more high level functions (either * or the Ax_mul_Bx versions) so that Julias multiple dispatch can do its thing and run the generic versions (julia implementations) of the functions when it is called with "graident numbers".
Easiest would be if you could post a small runnable example of something that models what you actually want to do.
