LAPACK/ScaLAPACK Development

Hi,

I have a scalar function known by its value for various points. I thought I could use DGELSD to find a "good" polynomial interpolation using less polynoms than points, but it seems not to work for that... am I wrong or is there something else for that kind of purposes in LAPACK ?

thanks