18 #ifndef __itkDiscreteGaussianCurvatureQuadEdgeMeshFilter_h
19 #define __itkDiscreteGaussianCurvatureQuadEdgeMeshFilter_h
22 #include "vnl/vnl_math.h"
37 template<
class TInputMesh,
class TOutputMesh=TInputMesh >
71 #ifdef ITK_USE_CONCEPT_CHECKING
102 qe_it2 = qe_it->GetOnext();
103 q0 = output->GetPoint( qe_it->GetDestination() );
104 q1 = output->GetPoint( qe_it2->GetDestination() );
108 TriangleType::ComputeAngle(q0, iP, q1) );
109 area += this->ComputeMixedArea(qe_it, qe_it2);
112 while ( qe_it != qe );
123 void operator=(
const Self &);