Euclidean -1.0 0.0 1.0 0.0 0 0 1 0 14 1 2 0 0.5 6 1 -1 -1 -1 def calc( arg1, arg2 ): # Construct the hyperbolic circle centered in arg1 through arg2 a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by <= 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) num2 = bsq - asq - den*bx dby2 = num2/by dy = by - dby2/2. d = Coordinate(ax,dy) r2 = sqrt( den*den + dby2*dby2 ) return Circle(d,r2/2.) 0.01 def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa) def calc( arg1, arg2 ): # Construct the hyperbolic segment (arc) connecting the two points a=arg1.coordinate(); b=arg2.coordinate(); ax=a.x ay=a.y if ay < 0: return InvalidObject() bx=b.x by=b.y if by < 0: return InvalidObject() asq=ax*ax + ay*ay bsq=bx*bx + by*by den = 2*(bx - ax) if den == 0: return Segment(a,b) x = (bsq - asq)/den c = Coordinate(x,0) rsq = asq + x*x - 2*x*ax r = sqrt (rsq) if rsq > 1000*asq: return Segment(a,b) alfa = atan2(ay, ax-x) beta = atan2(by, bx-x) return Arc(c,r,alfa,beta-alfa)