Interpolazione 3d c++

[repka]

Salve a tutti,

il mio problema e' il seguente: Ho un insieme di punti Xi,Yi,Zi (con i=0...i=n). Possiamo vedere questi punti come
una superficie nello spazio, e la sua rappresentazione non e' importante, puo' essere una matrice o 3 vettori.

Io ho bisogno di una interpolazione 3d (il metodo dell'interpolazione puo' essere spline, polinomiale ecc ecc)
che mi permetta di sapere la Zk (per k!=n) di un dato Xk, Yk (per k!=n).

Qualcuno ha una qualche idea di libreria c++ gia' pronta o di algoritmo di facile implementazione?

Grazie a tutti della cortese attenzione