diff --git a/TP_maillage.m b/TP_maillage.m index 1ce61aa..d9d8b1c 100644 --- a/TP_maillage.m +++ b/TP_maillage.m @@ -337,4 +337,24 @@ view(80, -10); % Sauvegarde des donnees save donnees; +%% Affichage du maillage, au lieu des surfaces +% Calcul des faces du maillage à garder +FACES = [sort(triBis(:, [2 3 4]), 2); sort(triBis(:, [1 3 4]), 2); sort(triBis(:, [1 2 4]), 2); sort(triBis(:, [1 2 3]), 2)]; +FACES = sortrows(FACES); +rep = sum(FACES(1:end - 1, :) == FACES(2:end, :), 2) == 3; +FACES([0; rep] | [rep; 0], :) = []; + +fprintf('Calcul du maillage final termine : %d faces. \n', size(FACES, 1)); + +% Affichage du maillage final +figure(9); +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; + +view(80, -10); \ No newline at end of file diff --git a/TP_maillage_suite.m b/TP_maillage_suite.m deleted file mode 100644 index e694e91..0000000 --- a/TP_maillage_suite.m +++ /dev/null @@ -1,21 +0,0 @@ -load donnees; - -% Calcul des faces du maillage à garder -FACES = [sort(triBis(:, [2 3 4]), 2); sort(triBis(:, [1 3 4]), 2); sort(triBis(:, [1 2 4]), 2); sort(triBis(:, [1 2 3]), 2)]; -FACES = sortrows(FACES); -rep = sum(FACES(1:end - 1, :) == FACES(2:end, :), 2) == 3; -FACES([0; rep] | [rep; 0], :) = []; - -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; - -view(80, -10);