n the device
sets.WarpSize = d.SIMDWidth;            % The warp size in threads
%% Creation of global stiffness matrix on CPU (serial)

tic;
K_f = StiffMa_vs(Mesh, MP, sets);	% Assembly on CPU: K
time = toc;
fprintf('\nElapsed time for building K on serial CPU: %f\n',time);
%% Creation of global stiffness matrix on CPU (serial) taking advantage of symmetry

tic;
K_s = StiffMa_vss(Mesh, MP, sets);  % Assembly on CPU (tril(K))
time_h = toc;
fprintf('Elapsed time for building tril(K) on serial CPU: %f\n',time_