r of DOFs per node for the vector problem
sets.edof = sets.dxn * sets.nxe;                % Number of DOFs per element
sets.sz = (sets.edof * (sets.edof + 1) )/2;     % Number of NNZ values for each Ke using simmetry
sets.tdofs = sets.nnod * sets.dxn;              % Number of total DOFs in the mesh
%% GPU Settings

d = gpuDevice;                          % Selects the GPU device 
sets.tbs      = d.MaxThreadsPerBlock;   % Thread Block Size
sets.numSMs   = d.MultiprocessorCount;  % Number of multiprocessors o