1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
/****************************************************************************
**
** Copyright (C) 2013 Digia Plc and/or its subsidiary(-ies).
** Contact: http://www.qt-project.org/legal
**
** This file is part of the test suite of the Qt Toolkit.
**
** $QT_BEGIN_LICENSE:LGPL$
** Commercial License Usage
** Licensees holding valid commercial Qt licenses may use this file in
** accordance with the commercial license agreement provided with the
** Software or, alternatively, in accordance with the terms contained in
** a written agreement between you and Digia. For licensing terms and
** conditions see http://qt.digia.com/licensing. For further information
** use the contact form at http://qt.digia.com/contact-us.
**
** GNU Lesser General Public License Usage
** Alternatively, this file may be used under the terms of the GNU Lesser
** General Public License version 2.1 as published by the Free Software
** Foundation and appearing in the file LICENSE.LGPL included in the
** packaging of this file. Please review the following information to
** ensure the GNU Lesser General Public License version 2.1 requirements
** will be met: http://www.gnu.org/licenses/old-licenses/lgpl-2.1.html.
**
** In addition, as a special exception, Digia gives you certain additional
** rights. These rights are described in the Digia Qt LGPL Exception
** version 1.1, included in the file LGPL_EXCEPTION.txt in this package.
**
** GNU General Public License Usage
** Alternatively, this file may be used under the terms of the GNU
** General Public License version 3.0 as published by the Free Software
** Foundation and appearing in the file LICENSE.GPL included in the
** packaging of this file. Please review the following information to
** ensure the GNU General Public License version 3.0 requirements will be
** met: http://www.gnu.org/copyleft/gpl.html.
**
**
** $QT_END_LICENSE$
**
****************************************************************************/
#include "utils.h"
#include <assert.h>
#include <qglobal.h>
#include "qnum.h"
#define FloatToXFixed(i) (int)((i) * 65536)
#define IntToXFixed(i) ((i) << 16)
static double compute_x_at(XFixed y, XPointFixed p1, XPointFixed p2)
{
double d = XFixedToDouble(p2.x - p1.x);
return
XFixedToDouble(p1.x) + d*XFixedToDouble(y - p1.y)/XFixedToDouble(p2.y - p1.y);
}
double compute_area(XTrapezoid *trap)
{
double x1 = compute_x_at(trap->top, trap->left.p1, trap->left.p2);
double x2 = compute_x_at(trap->top, trap->right.p1, trap->right.p2);
double x3 = compute_x_at(trap->bottom, trap->left.p1, trap->left.p2);
double x4 = compute_x_at(trap->bottom, trap->right.p1, trap->right.p2);
double top = XFixedToDouble(trap->top);
double bottom = XFixedToDouble(trap->bottom);
double h = bottom - top;
double top_base = x2 - x1;
double bottom_base = x4 - x3;
if ((top_base < 0 && bottom_base > 0)
|| (top_base > 0 && bottom_base < 0)) {
double y0 = top_base*h/(top_base - bottom_base) + top;
double area = qAbs(top_base * (y0 - top) / 2.);
area += qAbs(bottom_base * (bottom - y0) /2.);
return area;
}
return 0.5 * h * qAbs(top_base + bottom_base);
}
double compute_area_for_x(const QVector<XTrapezoid> &traps)
{
double area = 0;
for (int i = 0; i < traps.size(); ++i) {
XTrapezoid trap = traps[i];
area += compute_area(&trap);
}
return area;
}
|
