- n += itk::CrossProduct(
- this->GetPoint( ( *eIt )->GetDestination( ) ) - p0,
- this->GetPoint( ( *nIt )->GetDestination( ) ) - p0
- );
+ VectorType pe = this->GetPoint( ( *eIt )->GetDestination( ) ) - p0;
+ VectorType pn = this->GetPoint( ( *nIt )->GetDestination( ) ) - p0;
+ n[ 0 ] += ( pe[ 1 ] * pn[ 2 ] ) - ( pe[ 2 ] * pn[ 1 ] );
+ n[ 1 ] += ( pe[ 2 ] * pn[ 0 ] ) - ( pe[ 0 ] * pn[ 2 ] );
+ n[ 2 ] += ( pe[ 0 ] * pn[ 1 ] ) - ( pe[ 1 ] * pn[ 0 ] );