18 #ifndef __itkDiscreteMeanCurvatureQuadEdgeMeshFilter_h
19 #define __itkDiscreteMeanCurvatureQuadEdgeMeshFilter_h
35 template<
class TInputMesh,
class TOutputMesh=TInputMesh >
71 #ifdef ITK_USE_CONCEPT_CHECKING
101 if ( qe != qe->GetOnext() )
116 qe_it2 = qe_it->GetOnext();
117 q0 = output->GetPoint( qe_it->GetDestination() );
118 q1 = output->GetPoint( qe_it2->GetDestination() );
120 temp_coeff = coefficent(output, qe_it);
121 Laplace += temp_coeff * ( iP - q0 );
123 temp_area = this->ComputeMixedArea(qe_it, qe_it2);
126 face_normal = TriangleType::ComputeNormal(q0, iP, q1);
127 normal += face_normal;
131 while ( qe_it != qe );
139 if ( normal.GetSquaredNorm() > 0. )
142 Laplace *= 0.25 / area;
143 oH = Laplace * normal;
158 void operator=(
const Self &);