Idom.num2PEstd.tetra: The indicator for two 3D points constituting a dominating set...

In elvanceyhan/pcds: Proximity Catch Digraphs and Their Applications

The indicator for two 3D points constituting a dominating set for Proportional Edge Proximity Catch Digraphs (PE-PCDs) - standard regular tetrahedron case

Description

Returns I({p1,p2} is a dominating set of the PE-PCD) where the vertices of the PE-PCD are the 3D data set Xp in the standard regular tetrahedron T_h=T((0,0,0),(1,0,0),(1/2,\sqrt{3}/2,0),(1/2,\sqrt{3}/6,\sqrt{6}/3)), that is, returns 1 if {p1,p2} is a dominating set of PE-PCD, returns 0 otherwise.

Point, p1, is in the region of vertex rv1 (default is NULL) and point, p2, is in the region of vertex rv2 (default is NULL); vertices (and hence rv1 and rv2) are labeled as 1,2,3,4 in the order they are stacked row-wise in T_h.

PE proximity region is constructed with respect to the tetrahedron T_h with expansion parameter r \ge 1 and vertex regions are based on center of mass CM (equivalent to circumcenter in this case).

ch.data.pnts is for checking whether points p1 and p2 are data points in Xp or not (default is FALSE), so by default this function checks whether the points p1 and p2 would constitute a dominating set if they actually were both in the data set.

See also (\insertCiteceyhan:Phd-thesis,ceyhan:comp-geo-2010;textualpcds).

Usage

Idom.num2PEstd.tetra(
  p1,
  p2,
  Xp,
  r,
  rv1 = NULL,
  rv2 = NULL,
  ch.data.pnts = FALSE
)

Arguments

p1, p2	Two 3D points to be tested for constituting a dominating set of the PE-PCD.
Xp	A set of 3D points which constitutes the vertices of the PE-PCD.
r	A positive real number which serves as the expansion parameter in PE proximity region; must be \ge 1.
rv1, rv2	The indices of the vertices whose regions contains p1 and p2, respectively. They take the vertex labels as 1,2,3,4 as in the row order of the vertices in T_h (default is NULL for both).
ch.data.pnts	A logical argument for checking whether points p1 and p2 are data points in Xp or not (default is FALSE).

Value

I({p1,p2} is a dominating set of the PE-PCD) where the vertices of the PE-PCD are the 3D data set Xp), that is, returns 1 if {p1,p2} is a dominating set of PE-PCD, returns 0 otherwise

Author(s)

Elvan Ceyhan

References

\insertAllCited

See Also

Idom.num2PEtetra, Idom.num2PEtri and Idom.num2PEbasic.tri

Examples

## Not run: 
A<-c(0,0,0); B<-c(1,0,0); C<-c(1/2,sqrt(3)/2,0); D<-c(1/2,sqrt(3)/6,sqrt(6)/3)
tetra<-rbind(A,B,C,D)

n<-5 #try also n<-20
Xp<-runif.std.tetra(n)$g  #try also Xp<-cbind(runif(n),runif(n),runif(n))
r<-1.5

Idom.num2PEstd.tetra(Xp[1,],Xp[2,],Xp,r)

ind.gam2<-vector()
for (i in 1:(n-1))
 for (j in (i+1):n)
 {if (Idom.num2PEstd.tetra(Xp[i,],Xp[j,],Xp,r)==1)
  ind.gam2<-rbind(ind.gam2,c(i,j))}

ind.gam2

#or try
rv1<-rel.vert.tetraCC(Xp[1,],tetra)$rv;rv2<-rel.vert.tetraCC(Xp[2,],tetra)$rv
Idom.num2PEstd.tetra(Xp[1,],Xp[2,],Xp,r,rv1,rv2)

#or try
rv1<-rel.vert.tetraCC(Xp[1,],tetra)$rv;
Idom.num2PEstd.tetra(Xp[1,],Xp[2,],Xp,r,rv1)

#or try
rv2<-rel.vert.tetraCC(Xp[2,],tetra)$rv
Idom.num2PEstd.tetra(Xp[1,],Xp[2,],Xp,r,rv2=rv2)

P1<-c(.1,.1,.1)
P2<-c(.4,.1,.2)
Idom.num2PEstd.tetra(P1,P2,Xp,r)

Idom.num2PEstd.tetra(c(-1,-1,-1),Xp[2,],Xp,r,ch.data.pnts = FALSE)
#gives an error message if ch.data.pnts = TRUE
#since not both points, p1 and p2, are data points in Xp

## End(Not run)

elvanceyhan/pcds documentation built on June 29, 2023, 8:12 a.m.