1ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov 2ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov//---------------------------------------------------------------------------- 3ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Anti-Grain Geometry - Version 2.3 4ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Copyright (C) 2002-2005 Maxim Shemanarev (http://www.antigrain.com) 5ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// 6ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Permission to copy, use, modify, sell and distribute this software 7ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// is granted provided this copyright notice appears in all copies. 8ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// This software is provided "as is" without express or implied 9ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// warranty, and with no claim as to its suitability for any purpose. 10ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// 11ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov//---------------------------------------------------------------------------- 12ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Contact: mcseem@antigrain.com 13ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// mcseemagg@yahoo.com 14ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// http://www.antigrain.com 15ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov//---------------------------------------------------------------------------- 16ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Bessel function (besj) was adapted for use in AGG library by Andy Wilk 17ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov// Contact: castor.vulgaris@gmail.com 18ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov//---------------------------------------------------------------------------- 19ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov#ifndef AGG_MATH_INCLUDED 20ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov#define AGG_MATH_INCLUDED 21ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov#include "agg_basics.h" 22ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganovnamespace agg 23ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov{ 24ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganovconst FX_FLOAT intersection_epsilon = 1.0e-30f; 25ee451cb395940862dad63c85adfe8f2fd55e864cSvet GanovAGG_INLINE FX_FLOAT calc_point_location(FX_FLOAT x1, FX_FLOAT y1, 26ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT x2, FX_FLOAT y2, 27ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT x, FX_FLOAT y) 28ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov{ 29ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return FXSYS_Mul(x - x2, y2 - y1) - FXSYS_Mul(y - y2, x2 - x1); 30ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov} 31ee451cb395940862dad63c85adfe8f2fd55e864cSvet GanovAGG_INLINE FX_FLOAT calc_distance(FX_FLOAT x1, FX_FLOAT y1, FX_FLOAT x2, FX_FLOAT y2) 32ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov{ 33ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT dx = x2 - x1; 34ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT dy = y2 - y1; 35ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return FXSYS_sqrt2(dx, dy); 36ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov} 37ee451cb395940862dad63c85adfe8f2fd55e864cSvet GanovAGG_INLINE FX_FLOAT calc_line_point_distance(FX_FLOAT x1, FX_FLOAT y1, 38ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT x2, FX_FLOAT y2, 39ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT x, FX_FLOAT y) 40ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov{ 41ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT dx = x2 - x1; 42ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT dy = y2 - y1; 43ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT d = FXSYS_sqrt2(dx, dy); 44ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov if(d < intersection_epsilon) { 45ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return calc_distance(x1, y1, x, y); 46ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov } 47ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return FXSYS_MulDiv(x - x2, dy, d) - FXSYS_MulDiv(y - y2, dx, d); 48ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov} 49ee451cb395940862dad63c85adfe8f2fd55e864cSvet GanovAGG_INLINE bool calc_intersection(FX_FLOAT ax, FX_FLOAT ay, FX_FLOAT bx, FX_FLOAT by, 50ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT cx, FX_FLOAT cy, FX_FLOAT dx, FX_FLOAT dy, 51ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT* x, FX_FLOAT* y) 52ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov{ 53ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT num = FXSYS_Mul(ay - cy, dx - cx) - FXSYS_Mul(ax - cx, dy - cy); 54ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov FX_FLOAT den = FXSYS_Mul(bx - ax, dy - cy) - FXSYS_Mul(by - ay, dx - cx); 55ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov if (FXSYS_fabs(den) < intersection_epsilon) { 56ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return false; 57ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov } 58ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov *x = ax + FXSYS_MulDiv(bx - ax, num, den); 59ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov *y = ay + FXSYS_MulDiv(by - ay, num, den); 60ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov return true; 61ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov} 62ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov} 63ee451cb395940862dad63c85adfe8f2fd55e864cSvet Ganov#endif 64