Calculer les facettes d'une surface via Numpy

Bonjour,
avant de réinventer la roue, je voulais savoir si quelqu'un connaissait, ou bien aurait, un script qui étant donnée une surface ayant une équation cartésienne z = f(x ; y) calculerait des facettes triangulaires relativement à un maillage en (x ;y) (c'est pour représenter la surface).

Autrement dit, on donne le maillage et l'équation, et le script calcule les sommets tout en revoyant des triplets correspondant aux triangles.

Ce serait un joli cadeau de Noël... :lol: