xref: /aoo4110/main/basegfx/source/numeric/ftools.cxx (revision b1cdbd2c)
1*b1cdbd2cSJim Jagielski /**************************************************************
2*b1cdbd2cSJim Jagielski  *
3*b1cdbd2cSJim Jagielski  * Licensed to the Apache Software Foundation (ASF) under one
4*b1cdbd2cSJim Jagielski  * or more contributor license agreements.  See the NOTICE file
5*b1cdbd2cSJim Jagielski  * distributed with this work for additional information
6*b1cdbd2cSJim Jagielski  * regarding copyright ownership.  The ASF licenses this file
7*b1cdbd2cSJim Jagielski  * to you under the Apache License, Version 2.0 (the
8*b1cdbd2cSJim Jagielski  * "License"); you may not use this file except in compliance
9*b1cdbd2cSJim Jagielski  * with the License.  You may obtain a copy of the License at
10*b1cdbd2cSJim Jagielski  *
11*b1cdbd2cSJim Jagielski  *   http://www.apache.org/licenses/LICENSE-2.0
12*b1cdbd2cSJim Jagielski  *
13*b1cdbd2cSJim Jagielski  * Unless required by applicable law or agreed to in writing,
14*b1cdbd2cSJim Jagielski  * software distributed under the License is distributed on an
15*b1cdbd2cSJim Jagielski  * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
16*b1cdbd2cSJim Jagielski  * KIND, either express or implied.  See the License for the
17*b1cdbd2cSJim Jagielski  * specific language governing permissions and limitations
18*b1cdbd2cSJim Jagielski  * under the License.
19*b1cdbd2cSJim Jagielski  *
20*b1cdbd2cSJim Jagielski  *************************************************************/
21*b1cdbd2cSJim Jagielski 
22*b1cdbd2cSJim Jagielski 
23*b1cdbd2cSJim Jagielski 
24*b1cdbd2cSJim Jagielski // MARKER(update_precomp.py): autogen include statement, do not remove
25*b1cdbd2cSJim Jagielski #include "precompiled_basegfx.hxx"
26*b1cdbd2cSJim Jagielski #include <basegfx/numeric/ftools.hxx>
27*b1cdbd2cSJim Jagielski #include <algorithm>
28*b1cdbd2cSJim Jagielski 
29*b1cdbd2cSJim Jagielski namespace basegfx
30*b1cdbd2cSJim Jagielski {
31*b1cdbd2cSJim Jagielski 	// init static member of class fTools
32*b1cdbd2cSJim Jagielski 	double ::basegfx::fTools::mfSmallValue = 0.000000001;
33*b1cdbd2cSJim Jagielski 
snapToNearestMultiple(double v,const double fStep)34*b1cdbd2cSJim Jagielski 	double snapToNearestMultiple(double v, const double fStep)
35*b1cdbd2cSJim Jagielski 	{
36*b1cdbd2cSJim Jagielski 		if(fTools::equalZero(fStep))
37*b1cdbd2cSJim Jagielski 		{
38*b1cdbd2cSJim Jagielski 			// with a zero step, all snaps to 0.0
39*b1cdbd2cSJim Jagielski 			return 0.0;
40*b1cdbd2cSJim Jagielski 		}
41*b1cdbd2cSJim Jagielski 		else
42*b1cdbd2cSJim Jagielski 		{
43*b1cdbd2cSJim Jagielski 			const double fHalfStep(fStep * 0.5);
44*b1cdbd2cSJim Jagielski 			const double fChange(fHalfStep - fmod(v + fHalfStep, fStep));
45*b1cdbd2cSJim Jagielski 
46*b1cdbd2cSJim Jagielski             if(basegfx::fTools::equal(fabs(v), fabs(fChange)))
47*b1cdbd2cSJim Jagielski             {
48*b1cdbd2cSJim Jagielski                 return 0.0;
49*b1cdbd2cSJim Jagielski             }
50*b1cdbd2cSJim Jagielski             else
51*b1cdbd2cSJim Jagielski             {
52*b1cdbd2cSJim Jagielski     			return v + fChange;
53*b1cdbd2cSJim Jagielski             }
54*b1cdbd2cSJim Jagielski 		}
55*b1cdbd2cSJim Jagielski 	}
56*b1cdbd2cSJim Jagielski 
snapToZeroRange(double v,double fWidth)57*b1cdbd2cSJim Jagielski 	double snapToZeroRange(double v, double fWidth)
58*b1cdbd2cSJim Jagielski 	{
59*b1cdbd2cSJim Jagielski 		if(fTools::equalZero(fWidth))
60*b1cdbd2cSJim Jagielski 		{
61*b1cdbd2cSJim Jagielski 			// with no range all snaps to range bound
62*b1cdbd2cSJim Jagielski 			return 0.0;
63*b1cdbd2cSJim Jagielski 		}
64*b1cdbd2cSJim Jagielski 		else
65*b1cdbd2cSJim Jagielski 		{
66*b1cdbd2cSJim Jagielski 			if(v < 0.0 || v > fWidth)
67*b1cdbd2cSJim Jagielski 			{
68*b1cdbd2cSJim Jagielski 				double fRetval(fmod(v, fWidth));
69*b1cdbd2cSJim Jagielski 
70*b1cdbd2cSJim Jagielski 				if(fRetval < 0.0)
71*b1cdbd2cSJim Jagielski 				{
72*b1cdbd2cSJim Jagielski 					fRetval += fWidth;
73*b1cdbd2cSJim Jagielski 				}
74*b1cdbd2cSJim Jagielski 
75*b1cdbd2cSJim Jagielski 				return fRetval;
76*b1cdbd2cSJim Jagielski 			}
77*b1cdbd2cSJim Jagielski 			else
78*b1cdbd2cSJim Jagielski 			{
79*b1cdbd2cSJim Jagielski 				return v;
80*b1cdbd2cSJim Jagielski 			}
81*b1cdbd2cSJim Jagielski 		}
82*b1cdbd2cSJim Jagielski 	}
83*b1cdbd2cSJim Jagielski 
snapToRange(double v,double fLow,double fHigh)84*b1cdbd2cSJim Jagielski 	double snapToRange(double v, double fLow, double fHigh)
85*b1cdbd2cSJim Jagielski 	{
86*b1cdbd2cSJim Jagielski 		if(fTools::equal(fLow, fHigh))
87*b1cdbd2cSJim Jagielski 		{
88*b1cdbd2cSJim Jagielski 			// with no range all snaps to range bound
89*b1cdbd2cSJim Jagielski 			return 0.0;
90*b1cdbd2cSJim Jagielski 		}
91*b1cdbd2cSJim Jagielski 		else
92*b1cdbd2cSJim Jagielski 		{
93*b1cdbd2cSJim Jagielski 			if(fLow > fHigh)
94*b1cdbd2cSJim Jagielski 			{
95*b1cdbd2cSJim Jagielski 				// correct range order. Evtl. assert this (?)
96*b1cdbd2cSJim Jagielski 				std::swap(fLow, fHigh);
97*b1cdbd2cSJim Jagielski 			}
98*b1cdbd2cSJim Jagielski 
99*b1cdbd2cSJim Jagielski 			if(v < fLow || v > fHigh)
100*b1cdbd2cSJim Jagielski 			{
101*b1cdbd2cSJim Jagielski 				return snapToZeroRange(v - fLow, fHigh - fLow) + fLow;
102*b1cdbd2cSJim Jagielski 			}
103*b1cdbd2cSJim Jagielski 			else
104*b1cdbd2cSJim Jagielski 			{
105*b1cdbd2cSJim Jagielski 				return v;
106*b1cdbd2cSJim Jagielski 			}
107*b1cdbd2cSJim Jagielski 		}
108*b1cdbd2cSJim Jagielski 	}
109*b1cdbd2cSJim Jagielski } // end of namespace basegfx
110*b1cdbd2cSJim Jagielski 
111*b1cdbd2cSJim Jagielski // eof
112