For a given FanF the function computes the subfan F1 of all smooth cones.
Let's consider the fan consisting of the following three dimensional cone and all of its faces:
i1 : C = posHull matrix {{1,-1,0},{1,1,0},{1,1,1}}
o1 = {ambient dimension => 3 }
dimension of lineality space => 0
dimension of the cone => 3
number of facets => 3
number of rays => 3
o1 : Cone
i2 : F = fan C
o2 = {ambient dimension => 3 }
number of generating cones => 1
number of rays => 3
top dimension of the cones => 3
o2 : Fan
This cone is not smooth, therefore also the fan is not. But if we remove the interior and one of the two dimensional faces the resulting subfan is smooth.
i3 : F1 = smoothSubfan F
o3 = {ambient dimension => 3 }
number of generating cones => 2
number of rays => 3
top dimension of the cones => 2
o3 : Fan