cog_elements_alu3d_hexa.m

function cog = cog_elements_alu3d_hexa(M);
%function cog = cog_elements_alu3d_hexa(M);
%
% function computing the cog of the elements of a hexaeder alu3d mesh
  
% Bernard Haasdonk 16.3.2006

linvind = M.elements(:)+1; % matlab/C-offset
X = sum(reshape(M.vertices(1,linvind),8,M.num_elements))/8;
Y = sum(reshape(M.vertices(2,linvind),8,M.num_elements))/8;
Z = sum(reshape(M.vertices(3,linvind),8,M.num_elements))/8;
cog = [X;Y;Z];


% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate