Download All

Code covered by the BSD License  

Highlights from
3D Crouzeix-Raviart mortar finite element method

image thumbnail
from 3D Crouzeix-Raviart mortar finite element method by Jan Valdman
Implementation of 3D Crouzeix-Raviart mortar finite element

middle_points=face_middlepoint(X,F) 
function middle_points=face_middlepoint(X,F)
middle_points=[sum(reshape(X(F,1),size(F,1),3),2), ...
               sum(reshape(X(F,2),size(F,1),3),2), ...
               sum(reshape(X(F,3),size(F,1),3),2)]/3;



% middle_points=zeros(size(F,2),3);
% for i=1:size(T,1)
%     face1=T(i,[2 3 4])'; middle_points(F(i,1),:)=sum(X(face1,:),1)/3;
%     face2=T(i,[1 3 4])'; middle_points(F(i,2),:)=sum(X(face2,:),1)/3;
%     face3=T(i,[1 2 4])'; middle_points(F(i,3),:)=sum(X(face3,:),1)/3;
%     face4=T(i,[1 2 3])'; middle_points(F(i,4),:)=sum(X(face4,:),1)/3;
% end

Contact us at files@mathworks.com