Code covered by the BSD License  

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

5.0

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

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

by

 

04 May 2009 (Updated )

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

| Watch this File

File Information
Description

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
http://sites.google.com/site/janvaldman/publications

Please cite the paper when using the code.

Acknowledgements

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.

Updates
04 May 2009

description

05 May 2009

description + screenshot

11 Jun 2009

description

14 Jun 2009

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

22 Jul 2009

paper citation added

Contact us