Code covered by the BSD License
-
M_CR=stima_massCR_3D(elements...
Crouzeix-Raviart mass matrix computation
-
[M_CR Z]=stima_gradgradCR_3D(...
-
[b_volumeforces, volumes]=rhs...
GaussPoints2=GaussPoints1;
-
[element2edges, edge2nodes]=g...
function: [element2edges, edge2nodes]=edge_numbering(elements)
-
[element2faces, face2nodes]=g...
function: [element2faces, face2nodes]=getFaces(elements)
-
[elements4,T1_Plane_index]=ro...
elements4_OnIntersection=elements4(T1_Plane_index,:);
-
[matrix,I]=deleterepeatedrows...
function: [element2edges, edge2nodes]=edge_numbering(elements)
-
[x_center, y_center, area]=ce...
-
entryrows=entryInWhichRows(A)
function: entryrows=entryInWhichRows(A)
-
f(point)
VolumeForce =ones(size(VolumeForce));
-
middle_points=face_middlepoin...
middle_points=zeros(size(F,2),3);
-
show_mesh(elements,coordinate...
-
show_middlepoints(face_number...
-
show_points(pointsX,pointsY,p...
-
show_vertices(vertices,coordi...
-
tetrarefine3(XYZ,Elm,ERef);
author: Ales Janka, ales.janka@unifr.ch, http://perso.unifr.ch/ales.janka
-
value=exact_solution(point)
onesvector=ones(size(x));
-
picture_paper.m
-
start.m
-
View all files
from
3D Crouzeix-Raviart mortar finite element method
by Jan Valdman
Implementation of 3D Crouzeix-Raviart mortar finite element
|
| value=exact_solution(point) |
function value=exact_solution(point)
x=point(:,1);
y=point(:,2);
z=point(:,3);
%onesvector=ones(size(x));
%x_shifted=onesvector-x;
%y_shifted=onesvector-y;
%value=(x.*x_shifted).*(y.*y_shifted);
value=sin(pi*x).*sin(pi*y).*sin(pi*z);
|
|
Contact us at files@mathworks.com