product of a matrix by the inverse of a symmetric one

[Francisco J Molina]

I am trying to multiply A by (B^(-1)) where B^(-1) is the inverse of B and B
is symmetric.
I think that the best way to do this is to use

gsl_linalg_LU_decomp

to get an LU decomposition of B and then use gsl_blas_dtrsm several times.
Am I right?