CGAL 5.5  dD Geometry Kernel

#include <CGAL/Epick_d.h>
Public Member Functions  
FT  operator() (Weighted_point_d pw, Weighted_point_d qw) 
returns the power product (aka power distance) of the weighted points pw and qw , that is, the squared Euclidean distance between the points minus their weights.  