CGAL users discussion list

Hello all, I am using CGAL for mesh editing and I am using TAUCS packaged in 
CGAL. 

I find that it is very slow to stuff the symmetric matrix 
Taucs_symmetric_solver_traits<T>::Matrix with set_coef(int i, int j, T val),
like this:
for (i=0; i<dim; i++) {
    for (j=0; j<dim; j++) {
       M.set_coef(int i, int j, my_matrix[i][j]);
    }
}
I know this may be because there is some optimizing work during the stuffing 
but it indeed dramatically slow down my program. 

Please tell me how to tackle this.

Thanks in advance!