DGEMM error in quadratic calculation