Optimal symmetric / general matrix multiply

I understand that the low level BLAS library has operations like CSYMM which allow one to do a symmetric complex matrix multiplying a general complex matrix faster than a general on general matrix multiply.

How does one invoke this ability from Matlab in a case like the last line in

t=1:2000;
A = exp(-1.0i * (transpose(t)*t) );
B = rand(2000) + 1.0i * rand(2000);
C = A*B;

I find that the above runs at the same speed regardless of whether A is symmetric or not. Thank you for any ideas!