Computation of the Friedrich's constant for a unit square domain: symetrizeMatrix(A) - File Exchange

Code covered by the BSD License

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

[coordinates,elements3,dirich... function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
[element2edges, edge2nodes]=g... function: [element2edges, edge2nodes]=edge_numbering(elements)
[matrix,I]=deleterepeatedrows... function: [element2edges, edge2nodes]=edge_numbering(elements)
show_mesh(elements,coordinate...
symetrizeMatrix(A)
start_Friedrichs.m
View all files

from Computation of the Friedrich's constant for a unit square domain by Jan Valdman
Friedrichs' constant is computed for a unit square geometry and homogeneous Dirichlet conditions.

symetrizeMatrix(A)

function A_sym = symetrizeMatrix(A)
[i,j,k]=find(A);
W=sparse([i; j], [j; i], ones(size(k,1)*2,1));
A_help=sparse([i; j], [j; i], [k; k]);
[i,j,k]=find(A_help);
[i,j,kk]=find(W);
A_sym=sparse(i,j,(kk.^(-1)).*k); %Now Kantennr_sym is a symetric form of Kantennr

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

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