home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / ledar34.zip / leda-r-3_4_tar / LEDA-3.4 / src / plane_alg / _float_balaban.c next >
C/C++ Source or Header  |  1996-09-03  |  1KB  |  50 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA 3.4
  4. +
  5. +  _float_balaban.c
  6. +
  7. +  This file is part of the LEDA research version (LEDA-R) that can be 
  8. +  used free of charge in academic research and teaching. Any commercial
  9. +  use of this software requires a license which is distributed by the
  10. +  LEDA Software GmbH, Postfach 151101, 66041 Saarbruecken, FRG
  11. +  (fax +49 681 31104).
  12. +
  13. +  Copyright (c) 1991-1996  by  Max-Planck-Institut fuer Informatik
  14. +  Im Stadtwald, 66123 Saarbruecken, Germany     
  15. +  All rights reserved.
  16. + 
  17. *******************************************************************************/
  18. #define NO_RAT_ALGORITHMS
  19. #include <LEDA/plane_alg.h>
  20.  
  21. int cmp_at_x(const segment& s1, const segment& s2, const point& x)
  22. {   // Precondition : s1.start() <= s1.end()    s2.start() <= s2.end()
  23.  
  24.   if (s1.is_vertical())
  25.      if (s2.is_vertical()) 
  26.         return 0;
  27.      else
  28.         return orientation(s2, x) >  0 ? 1 : -1;
  29.  
  30.   if(s2.is_vertical()) return orientation(s1, x) <= 0 ? 1 : -1;
  31.  
  32.   double s1yd = s2.dx()*(s1.dy()*(x.xcoord() - s1.xcoord1()) + 
  33.              s1.dx()*s1.ycoord1());
  34.  
  35.   double s2yd = s1.dx()*(s2.dy()*(x.xcoord() - s2.xcoord1()) + 
  36.              s2.dx()*s2.ycoord1());
  37.  
  38.   if( s1yd < s2yd ) return -1;
  39.   if( s1yd > s2yd ) return  1;
  40.   return (orientation( s1, x ) > 0) ? cmp_slopes(s1, s2) : -cmp_slopes(s1, s2);
  41. }
  42.  
  43.  
  44. typedef point      POINT;
  45. typedef segment    SEGMENT;
  46.  
  47. #define STAIRCASE   float_stair_case
  48.  
  49. #include "balaban.c"
  50.