BuildMeshP2.m

Go to the documentation of this file.

function [qf,mef]=BuildMeshP2(Th,e2q,e2t,e2n)
% Numerotation globale
qf=[Th.q (Th.q(:,e2q(:,1))+Th.q(:,e2q(:,2)))/2];

ned=length(e2q);

% Code vectorisé
me1=zeros(3,Th.nme);
%me1(1:3,:)=Th.me;
%me1=Th.me;
L=[1:ned]+Th.nq;
me1(3*(e2t(:,1)-1)+e2n(:,1))=L;
I=3*(e2t(:,2)-1)+e2n(:,2);
K=find(I~=-3);
me1(I(K))=L(K);

mef([1,3,6,2,5,4],:)=[Th.me;me1];

% Code non vectorisé
% mef=zeros(3,Th.nme);
% for l=1:ned
%   mef(e2n(l,1),e2t(l,1))=Th.nq+l;
%   if (e2t(l,2) ~= 0 )
%     mef(e2n(l,2),e2t(l,2))=Th.nq+l;
%   end
% end