     % Number of NNZ values for each Ke using simmetry
sets.tdofs = sets.nnod * sets.dxn;              % Number of total DOFs in the mesh
%% Element stiffness matrix computation on CPU

tic;
Ke_hf = eStiff_vsa(Mesh, MP, sets);        % Computation of Ke for K
times = toc;
fprintf('Elapsed time for computing Ke for K on serial CPU: %f\n',times);
%% Element stiffness matrix computation on CPU (symmetry)

tic;
Ke_hs  = eStiff_vssa(Mesh, MP, sets);      % Computation of Ke for tril(K)
times = toc;
fprintf('Elap