liste1

 LISTE1 Assembles lists of internal, boundary, interface nodes (similar to liste)

 [lbor,lint,lintint,lgamma]=liste1(ifro);

 It computes lists of internal, boundary, interface nodes.
 All lists are referred to global ordering on Omega

 Input : ifro = column array of length noe=nov(npdx*npdy,ne):
            if (x_i,y_i) is internal to Omega then ifro(i)=0,
            if (x_i,y_i) is on \partial\Omega then ifro(i)=1,

 Output:
 lbor: list of boundary nodes
 lint: list of internal nodes

0001 function [lbor,lint,lintint,lgamma]=liste1(ifro);
0002 % LISTE1 Assembles lists of internal, boundary, interface nodes (similar to liste)
0003 %
0004 % [lbor,lint,lintint,lgamma]=liste1(ifro);
0005 %
0006 % It computes lists of internal, boundary, interface nodes.
0007 % All lists are referred to global ordering on Omega
0008 %
0009 % Input : ifro = column array of length noe=nov(npdx*npdy,ne):
0010 %            if (x_i,y_i) is internal to Omega then ifro(i)=0,
0011 %            if (x_i,y_i) is on \partial\Omega then ifro(i)=1,
0012 %
0013 % Output:
0014 % lbor: list of boundary nodes
0015 % lint: list of internal nodes
0016 %
0017 
0018 %   Written by Paola Gervasio
0019 %   $Date: 2007/04/01$
0020 
0021 noe=length(ifro);
0022 
0023 
0024 lint=[];lbor=[];lgamma=[]; lintint=[];
0025 for i=1:noe
0026 if(ifro(i)==1)
0027 lbor=[lbor;i];
0028 else
0029 lint=[lint;i];
0030 end
0031 end
0032 
0033 for i=1:length(lint)
0034 if(ifro(i)==-1)
0035 lgamma=[lgamma;i];
0036 else
0037 lintint=[lintint;i];
0038 end
0039 end
0040 
0041 return