3D Crouzeix-Raviart mortar finite element method: f(point) - File Exchange

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

f(point)
function VolumeForce = f(point) x=point(:,1); y=point(:,2); z=point(:,3); VolumeForce =3pi^2sin(pix).sin(piy).sin(pi*z); %VolumeForce =ones(size(VolumeForce));