Tester si un point appartient à un plan

[chellier]

Bonjour,

Pour faire suite à ce fil, j'essaye de tester si un point appartient à un plan.

La routine suivante à l'air de fonctionner :

Code : Tout sélectionner

bool inplan(triple p, path3 surfplane) {
  triple n=normal(surfplane);
  triple[] int=intersectionpoints(p--shift(n)*p,surface(surfplane));
  return (int.length!=0 && abs(int[0]-p)<.0001);
}

Mais n'y a-t'il pas quelque chose qui existe déjà et qui serait plus propre ?

Christophe