function [ST,FT]=faceTronquee(S,F,k)
% [ST,FT]=FACETRONQUEE(S,F,k) coupe un polyèdre régulier donné par S(ommets et F(aces
% K est le coefficient de découpage des arêtes.

nbs=F(1,1)-1;	% nombre de sommets des faces 
[nbf,co]=size(F);	% nombre de faces

ST=zeros(2*nbs*nbf,3);	% matrice des sommets
FT=zeros(nbf,2*nbs+2);	% matrice des faces

for i=1:nbf
	j1=F(i,2:nbs+1);
	j2=F(i,3:nbs+2);
	d=2*nbs*(i-1);
	ST(d+1:2:d+2*nbs,:)=S(j1,:)+k*(S(j2,:)-S(j1,:));
	ST(d+2:2:d+2*nbs,:)=S(j1,:)+(1-k)*(S(j2,:)-S(j1,:));
	FT(i,:)=[2*nbs+1 d+1:d+2*nbs d+1];
end

