Opt/EdgeLengthOpt.m

Go to the documentation of this file.

function lbe=EdgeLengthOpt(be,q)
% function lbe=EdgeLength(be,q)
%   Computation of the lengths of edges in the mesh
%   - Optimized version
%
% Parameters:
%  be: Connectivity array for boundary edges, `2\times\nbe` array.<br/>
% `\be(\il,l)` is the storage index of the
%  `\il`-th  vertex of the `l`-th edge in the array `\q` of vertices coordinates, `\il\in\{1,2\}` and
%       `l\in{\ENS{1}{\nbe}}`.
%  q: Array of vertices coordinates, `2\times\nq` array. <br/>
%  `{\q}(\il,j)` is the
%  `\il`-th coordinate of the `j`-th vertex, `\il\in\{1,2\}` and
%  `j\in\ENS{1}{\nq}`
%
% 
% Return values:
%  lbe: Array of edges lengths,  `1\times\nbe` array. `lbe(j)`  is the length of the `j`-th edge.
% Copyright:
%   See \ref license
  lbe=sqrt(sum((q(:,be(1,:)) - q(:,be(2,:))).^2,1));