The sparse data structure appears to be correct (though I did not check to see if it is symmetric), the only problem I notice is that all your values are zero.
By the way, you do know that you are not calling the solve routine, e.g. dss_solve_real, after the dss_factor_real...
