flux_optimization.m

%geometry properties (can not be avoided)
[element2edges, edge2nodes]=getEdges(elements);    
nodes2edge=sparse(edge2nodes(:,1),edge2nodes(:,2),1:size(edge2nodes,1),size(coordinates,1),size(coordinates,1));
nodes2edge=symetrizeMatrix(nodes2edge);
edge2elements=entryInWhichRows(element2edges); 

%
[DIVDIV_RT0 MASS_RT0]=stima_divdivRT0(elements,coordinates,areas,element2edges,edge2elements);