nov_p1_1d

 NOV_P1_1D  P1-FEM Local to global map

 [nov_p1,nov_p1_g]=nov_p1_1d(nov);

 INPUT: nov = map mesh_PN in Omega_ie -----> global mesh in Omega

 OUTPUT:   nov_p1 = map: mesh_P1 -----> mesh P_N in Omega_ie
           nov_p1_g = map: mesh_P1  -----> global mesh in Omega

0001 function [nov_p1,nov_p1_g]=nov_p1_1d(nov);
0002 % NOV_P1_1D  P1-FEM Local to global map
0003 %
0004 % [nov_p1,nov_p1_g]=nov_p1_1d(nov);
0005 %
0006 % INPUT: nov = map mesh_PN in Omega_ie -----> global mesh in Omega
0007 %
0008 % OUTPUT:   nov_p1 = map: mesh_P1 -----> mesh P_N in Omega_ie
0009 %           nov_p1_g = map: mesh_P1  -----> global mesh in Omega
0010 
0011 [npdx,ne]=size(nov);
0012 noe=nov(npdx,ne);
0013 nx=npdx-1;
0014 nov_p1=zeros(2,nx);
0015 nov_p1_g=zeros(2,nx);
0016 for ie=1:nx
0017     nov_p1(1:2,ie)=[ie,ie+1];
0018     nov_p1_g(1:2,ie)=nov([ie,ie+1]);
0019 end
0020 return