',times/time1);
fprintf('Should be %d!: %d\n', sets.edof^2*sets.nel, sum(iKhf==iKhfv));
fprintf('Should be %d!: %d\n', sets.edof^2*sets.nel, sum(jKhf==jKhfv));
%% Index computation on CPU (symmetry)

tic;
[iKh, jKh] = Index_sssa(Mesh.elements, sets);	% Row/column indices of tril(K)
times = toc;
fprintf('Elapsed time for computing row/column indices of tril(K) on serial CPU: %f\n',times);
%% Credits
% Written by *Francisco Javier Ramírez-Gil*, <http://fjramireg@gmail.com fjramireg@gmail.com>
% 
% Code site: