/****************************************************************************

**

** Copyright (C) 2009 Nokia Corporation and/or its subsidiary(-ies).

** All rights reserved.

** Contact: Nokia Corporation (qt-info@nokia.com)

**

** This file is part of the test suite of the Qt Toolkit.

**

** $QT_BEGIN_LICENSE:LGPL$

** No Commercial Usage

** This file contains pre-release code and may not be distributed.

** You may use this file in accordance with the terms and conditions

** contained in the Technology Preview License Agreement accompanying

** this package.

**

** 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, Nokia gives you certain additional

** rights.  These rights are described in the Nokia Qt LGPL Exception

** version 1.1, included in the file LGPL_EXCEPTION.txt in this package.

**

** If you have questions regarding the use of this file, please contact

** Nokia at qt-info@nokia.com.

**

**

**

**

**

**

**

**

** $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;

}