--- ray/src/common/face.c 2003/02/22 02:07:22 2.6 +++ ray/src/common/face.c 2003/03/12 17:26:58 2.8 @@ -1,66 +1,11 @@ #ifndef lint -static const char RCSid[] = "$Id: face.c,v 2.6 2003/02/22 02:07:22 greg Exp $"; +static const char RCSid[] = "$Id: face.c,v 2.8 2003/03/12 17:26:58 greg Exp $"; #endif /* * face.c - routines dealing with polygonal faces. */ -/* ==================================================================== - * The Radiance Software License, Version 1.0 - * - * Copyright (c) 1990 - 2002 The Regents of the University of California, - * through Lawrence Berkeley National Laboratory. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * - * 1. Redistributions of source code must retain the above copyright - * notice, this list of conditions and the following disclaimer. - * - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in - * the documentation and/or other materials provided with the - * distribution. - * - * 3. The end-user documentation included with the redistribution, - * if any, must include the following acknowledgment: - * "This product includes Radiance software - * (http://radsite.lbl.gov/) - * developed by the Lawrence Berkeley National Laboratory - * (http://www.lbl.gov/)." - * Alternately, this acknowledgment may appear in the software itself, - * if and wherever such third-party acknowledgments normally appear. - * - * 4. The names "Radiance," "Lawrence Berkeley National Laboratory" - * and "The Regents of the University of California" must - * not be used to endorse or promote products derived from this - * software without prior written permission. For written - * permission, please contact radiance@radsite.lbl.gov. - * - * 5. Products derived from this software may not be called "Radiance", - * nor may "Radiance" appear in their name, without prior written - * permission of Lawrence Berkeley National Laboratory. - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR - * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, - * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT - * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF - * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND - * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, - * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT - * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF - * SUCH DAMAGE. - * ==================================================================== - * - * This software consists of voluntary contributions made by many - * individuals on behalf of Lawrence Berkeley National Laboratory. For more - * information on Lawrence Berkeley National Laboratory, please see - * . - */ +#include "copyright.h" #include "standard.h" @@ -174,6 +119,7 @@ FACE *f; { int ncross, n; double x, y; + int tst; register int xi, yi; register FLOAT *p0, *p1; @@ -187,13 +133,15 @@ FACE *f; ncross = 0; /* positive x axis cross test */ while (n--) { - if ((p0[yi] > y) ^ (p1[yi] > y)) - if (p0[xi] > x && p1[xi] > x) + if ((p0[yi] > y) ^ (p1[yi] > y)) { + tst = (p0[xi] > x) + (p1[xi] > x); + if (tst == 2) ncross++; - else if (p0[xi] > x || p1[xi] > x) + else if (tst) ncross += (p1[yi] > p0[yi]) ^ ((p0[yi]-y)*(p1[xi]-x) > (p0[xi]-x)*(p1[yi]-y)); + } p0 = p1; p1 += 3; }