load donnees;

% Calcul des faces du maillage à garder
% nope = [];
% for i = 1:nb_barycentres
%     tetra_i = T.Points(triBis(i, :), :);
%     for j = i+1:nb_barycentres
%         tetra_j = T.Points(triBis(j, :), :);
%         if size(unique([tetra_i; tetra_j], 'rows'), 1) < 6
%             for ki = 1:4
%                 tri_ki = sort(triBis(i,[1:ki - 1 ki + 1:4]));
%                 for kj = 1:4
%                     tri_kj = sort(triBis(j,[1:kj - 1 kj + 1:4]));
%                     if tri_ki == tri_kj
%                         nope = [nope; i ki; j kj];
%                     end
%                 end
%             end
%         end
%     end
%     i
%     size(nope)
% end
FACES = [triBis(:, [2 3 4]); triBis(:, [1 3 4]); triBis(:, [1 2 4]); triBis(:, [1 2 3])];
FACES = sortrows(FACES);
FACESFACES(1:end-1, :) == FACES(2:end, :)

fprintf('Calcul du maillage final termine : %d faces. \n',size(FACES,1));

% Affichage du maillage final
figure;
hold on
for i = 1:size(FACES,1)
   plot3([X(1,FACES(i,1)) X(1,FACES(i,2))],[X(2,FACES(i,1)) X(2,FACES(i,2))],[X(3,FACES(i,1)) X(3,FACES(i,2))],'r');
   plot3([X(1,FACES(i,1)) X(1,FACES(i,3))],[X(2,FACES(i,1)) X(2,FACES(i,3))],[X(3,FACES(i,1)) X(3,FACES(i,3))],'r');
   plot3([X(1,FACES(i,3)) X(1,FACES(i,2))],[X(2,FACES(i,3)) X(2,FACES(i,2))],[X(3,FACES(i,3)) X(3,FACES(i,2))],'r');
end;