Code covered by the BSD License  

Highlights from
3D Crouzeix-Raviart mortar finite element method

3.0 | 1 rating Rate this file 29 Downloads (last 30 days) File Size: 32.8 KB File ID: #22430
image thumbnail

3D Crouzeix-Raviart mortar finite element method



13 Dec 2008 (Updated )

Implementation of 3D Crouzeix-Raviart mortar finite element

| Watch this File

File Information

This is a software generating results for the paper
Leszek Marcinkowski, Talal Rahman and Jan Valdman, A 3D Crouzeix-Raviart mortar finite element. Computing 86, No. 4, 313-330 (2009)

A link to the paper can be found at the author web page located at

Please cite the paper when you find the code useful.

Description; This solves the Poisson problem in the domain Omega= (0,1) x (0,1) x (-1,1) assuming the volume force equal to 3 pi*pi* sin(pi*x) sin(pi*y) sin(pi*z) and the zero Dirichlet conditions. The discrete solution is computed on two equal subdomains, Omega1= (0,1) x (0,1) x (0,1), Omega2= (0,1) x (0,1) x (-1,0) using 3D Crouzeix-Raviart mortar finite element method.The subdomains are triangulated using different regular triangulation in tetrahedrons. Therefore, there are no matching grids across the intersection plane z=0. To "connect" both subdomains solutions, a new interpolation operator is implemented and compared with a classical one.
To test the functionality of the algorithm, run "start.m".

It takes about 2 minutes on my PC to compute a solution on a mesh with 57216 edges. The major obstacle for the faster implementation is the use of matlab function polyxpoly works fine but is too general - it might be replaced by a faster function, which computes an intersection of two triangles in 2D.

Required Products MATLAB Compiler
MATLAB release MATLAB 7.5 (R2007b)
MATLAB Search Path
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
21 Aug 2011 Qwuasi Moore  
16 Dec 2008

improved functionality, getFaces function

04 May 2009


16 Nov 2009

full article citation added, title picture changed

07 Mar 2015

Paper information completed

Contact us