GetMatPrTh2RefTh.m

Go to the documentation of this file.

function MatProj=GetMatPrTh2RefTh(nq,e2q)
% function MatProj=GetMatPrTh2RefTh(nq,e2q)
%   Build projection matrix from initial mesh to refine mesh
  N=size(e2q,1);
  I=nq+(1:N)'*ones(1,2);
  I=I(:);
  J=double(e2q(:));
  K=ones(N,2)/2;
  K=K(:);
  L=(1:nq)';
  I=[L;I];
  J=[L;J];
  K=[ones(nq,1);K];
  MatProj=sparse(I,J,K,nq+N,nq);
end