tal DOFs in the mesh
%% Index computation on CPU (serial)

tic;
[iKhf, jKhf] = Index_ssa(Mesh.elements, sets);	% Row/column indices of K
times = toc;
fprintf('Elapsed time for computing row/column indices of K on serial CPU: %f\n',times);
%% Index computation on CPU (vectorized)

tic;
[iKhfv, jKhfv] = Index_sa(Mesh.elements', sets);	% Row/column indices of K
time1 = toc;
fprintf('Elapsed time for computing row/column indices of K on CPU (vectorized): %f\n',time1);
fprintf('Speedup (loop vs vectorized): %f\n