Code covered by the BSD License  

Highlights from
Computation of the Friedrich's constant for a unit square domain

5.0 | 1 rating Rate this file 7 Downloads (last 30 days) File Size: 5.18 KB File ID: #23991 Version: 1.5
image thumbnail

Computation of the Friedrich's constant for a unit square domain


Jan Valdman (view profile)


04 May 2009 (Updated )

Friedrichs' constant is computed for a unit square geometry and homogeneous Dirichlet conditions.

| Watch this File

File Information

To run the code, type "start_Friedrichs' in the Matlab window.

Mass and stiffness matrices are generated on nested uniform meshes starting from a coarse mesh with 2 elements (level 0) and ending with a finest mesh with 524288 elements (level 9). The Friedrichs' constant is computed as the largest value of a generalized eigenvalue problem. Discrete values are compared with an analytical value known for the unit square domain.

Geometry of interest can be easily changed for instance to L-shape domain by modifying matrices coordinates, elements3, dirichlet.

This is the first part of software generating results for the paper

Jan Valdman - Minimization of Functional Majorant in A Posteriori Error Analysis based on H(div) Multigrid-Preconditioned CG Method, Advances in Numerical Analysis, vol. 2009, Article ID 164519 (2009)

The paper can be downloaded from the author web

Please cite the paper when using the code.


This file inspired Double Porosity Model.

Required Products MATLAB Compiler
MATLAB release MATLAB 7.5 (R2007b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
13 Jun 2009 Nikola Toljic

Works good.

04 May 2009 1.1


05 May 2009 1.2

description + screenshot

11 Jun 2009 1.3


14 Jun 2009 1.4

Faster performance. Almost scalable generation of stiffness and mass matrices.

22 Jul 2009 1.5

paper citation added

Contact us