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.