Indices of points/edges/triangles in set of subdomains
c=pdesdp(p,e,t) [i,c]=pdesdp(p,e,t) c=pdesdp(p,e,t,sdl) [i,c]=pdesdp(p,e,t,sdl) i=pdesdt(t) i=pdesdt(t,sdl) i=pdesde(e) i=pdesde(e,sdl)
[i,c]=pdesdp(p,e,t,sdl) given mesh data p, e, and t and a list of subdomain numbers sdl, the function returns all points belonging to those subdomains. A point can belong to several subdomains, and the points belonging to the domains in sdl are divided into two disjoint sets. i contains indices of the points that wholly belong to the subdomains listed in sdl, and c lists points that also belongs to the other subdomains.
If sdl is not given, a list of all subdomains is assumed.