+# Get the number of canonical faces
+init_nb_f = geompy.GetNbCanonicalFaces(shape_igs)
+res_nb_f = geompy.GetNbCanonicalFaces(shape_cr)
+
+print "Number of canonical faces in the initial shape: %s" % init_nb_f
+print "Number of canonical faces in the result shape: %s" % res_nb_f
+