Something wrong with the dimensions.. simple dgemm problem

[skibbel]

i wanna do Jac=xI^T*dpN

xI_row=dpN_row=9
xI_col=dpN_col=2

in this example
so (9x2)^T*(9x2)=(2x2) Matrix


dgemm(transposed, not_transposed,&xI_col, &dpN_col, &dpN_row, &one_double, xI, &xI_row, dpN, &dpN_row, &zero_double, Jac, &xI_col);

i always have a small error in the result.