h);
fprintf('\tCPU speedup (CPU vs CPU, K vs tril(K)): %f\n',time/time_h);
%% Creation of global stiffness matrix on GPU

tic;
elementsGPU = gpuArray(Mesh.elements');	% Transfer to GPU memory
nodesGPU    = gpuArray(Mesh.nodes');  	% Transfer to GPU memory
K_d = StiffMa_vps(elementsGPU, nodesGPU, MP, sets);      % Assembly on GPU (symmetric)
time_d = toc;
wait(d);
fprintf('Elapsed time for building tril(K) on parallel GPU: %f\n',time_d);
fprintf('\tGPU speedup (K vs tril(K)): %f\n',time/time_d);
fprintf('\tG