Subversion Repositories AndroidProjects

#region --- License ---

/*

Copyright (c) 2006 - 2008 The Open Toolkit library.

Permission is hereby granted, free of charge, to any person obtaining a copy of

this software and associated documentation files (the "Software"), to deal in

the Software without restriction, including without limitation the rights to

use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies

of the Software, and to permit persons to whom the Software is furnished to do

so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all

copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR

IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,

FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE

AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER

LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,

OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE

SOFTWARE.

 */

#endregion

using System;

using System.Runtime.InteropServices;

namespace OpenTK

{

    /// <summary>Represents a 2D vector using two single-precision floating-point numbers.</summary>

    /// <remarks>

    /// The Vector2 structure is suitable for interoperation with unmanaged code requiring two consecutive floats.

    /// </remarks>

    [Serializable]

    [StructLayout(LayoutKind.Sequential)]

    public struct Vector2 : IEquatable<Vector2>

    {

        #region Fields

        /// <summary>

        /// The X component of the Vector2.

        /// </summary>

        public float X;

        /// <summary>

        /// The Y component of the Vector2.

        /// </summary>

        public float Y;

        #endregion

        #region Constructors

        /// <summary>

        /// Constructs a new Vector2.

        /// </summary>

        /// <param name="x">The x coordinate of the net Vector2.</param>

        /// <param name="y">The y coordinate of the net Vector2.</param>

        public Vector2(float x, float y)

        {

            X = x;

            Y = y;

        }

        /// <summary>

        /// Constructs a new Vector2 from the given Vector2.

        /// </summary>

        /// <param name="v">The Vector2 to copy components from.</param>

        [Obsolete]

        public Vector2(Vector2 v)

        {

            X = v.X;

            Y = v.Y;

        }

        /// <summary>

        /// Constructs a new Vector2 from the given Vector3.

        /// </summary>

        /// <param name="v">The Vector3 to copy components from. Z is discarded.</param>

        [Obsolete]

        public Vector2(Vector3 v)

        {

            X = v.X;

            Y = v.Y;

        }

        /// <summary>

        /// Constructs a new Vector2 from the given Vector4.

        /// </summary>

        /// <param name="v">The Vector4 to copy components from. Z and W are discarded.</param>

        [Obsolete]

        public Vector2(Vector4 v)

        {

            X = v.X;

            Y = v.Y;

        }

        #endregion

        #region Public Members

        #region Instance

        #region public void Add()

        /// <summary>Add the Vector passed as parameter to this instance.</summary>

        /// <param name="right">Right operand. This parameter is only read from.</param>

        [Obsolete("Use static Add() method instead.")]

        public void Add(Vector2 right)

        {

            this.X += right.X;

            this.Y += right.Y;

        }

        /// <summary>Add the Vector passed as parameter to this instance.</summary>

        /// <param name="right">Right operand. This parameter is only read from.</param>

        [CLSCompliant(false)]

        [Obsolete("Use static Add() method instead.")]

        public void Add(ref Vector2 right)

        {

            this.X += right.X;

            this.Y += right.Y;

        }

        #endregion public void Add()

        #region public void Sub()

        /// <summary>Subtract the Vector passed as parameter from this instance.</summary>

        /// <param name="right">Right operand. This parameter is only read from.</param>

        [Obsolete("Use static Subtract() method instead.")]

        public void Sub(Vector2 right)

        {

            this.X -= right.X;

            this.Y -= right.Y;

        }

        /// <summary>Subtract the Vector passed as parameter from this instance.</summary>

        /// <param name="right">Right operand. This parameter is only read from.</param>

        [CLSCompliant(false)]

        [Obsolete("Use static Subtract() method instead.")]

        public void Sub(ref Vector2 right)

        {

            this.X -= right.X;

            this.Y -= right.Y;

        }

        #endregion public void Sub()

        #region public void Mult()

        /// <summary>Multiply this instance by a scalar.</summary>

        /// <param name="f">Scalar operand.</param>

        [Obsolete("Use static Multiply() method instead.")]

        public void Mult(float f)

        {

            this.X *= f;

            this.Y *= f;

        }

        #endregion public void Mult()

        #region public void Div()

        /// <summary>Divide this instance by a scalar.</summary>

        /// <param name="f">Scalar operand.</param>

        [Obsolete("Use static Divide() method instead.")]

        public void Div(float f)

        {

            float mult = 1.0f / f;

            this.X *= mult;

            this.Y *= mult;

        }

        #endregion public void Div()

        #region public float Length

        /// <summary>

        /// Gets the length (magnitude) of the vector.

        /// </summary>

        /// <see cref="LengthFast"/>

        /// <seealso cref="LengthSquared"/>

        public float Length

        {

            get

            {

                return (float)System.Math.Sqrt(X * X + Y * Y);

            }

        }

        #endregion

        #region public float LengthFast

        /// <summary>

        /// Gets an approximation of the vector length (magnitude).

        /// </summary>

        /// <remarks>

        /// This property uses an approximation of the square root function to calculate vector magnitude, with

        /// an upper error bound of 0.001.

        /// </remarks>

        /// <see cref="Length"/>

        /// <seealso cref="LengthSquared"/>

        public float LengthFast

        {

            get

            {

                return 1.0f / MathHelper.InverseSqrtFast(X * X + Y * Y);

            }

        }

        #endregion

        #region public float LengthSquared

        /// <summary>

        /// Gets the square of the vector length (magnitude).

        /// </summary>

        /// <remarks>

        /// This property avoids the costly square root operation required by the Length property. This makes it more suitable

        /// for comparisons.

        /// </remarks>

        /// <see cref="Length"/>

        /// <seealso cref="LengthFast"/>

        public float LengthSquared

        {

            get

            {

                return X * X + Y * Y;

            }

        }

        #endregion

        #region public Vector2 PerpendicularRight

        /// <summary>

        /// Gets the perpendicular vector on the right side of this vector.

        /// </summary>

        public Vector2 PerpendicularRight

        {

            get

            {

                return new Vector2(Y, -X);

            }

        }

        #endregion

        #region public Vector2 PerpendicularLeft

        /// <summary>

        /// Gets the perpendicular vector on the left side of this vector.

        /// </summary>

        public Vector2 PerpendicularLeft

        {

            get

            {

                return new Vector2(-Y, X);

            }

        }

        #endregion

        #region public void Normalize()

        /// <summary>

        /// Scales the Vector2 to unit length.

        /// </summary>

        public void Normalize()

        {

            float scale = 1.0f / this.Length;

            X *= scale;

            Y *= scale;

        }

        #endregion

        #region public void NormalizeFast()

        /// <summary>

        /// Scales the Vector2 to approximately unit length.

        /// </summary>

        public void NormalizeFast()

        {

            float scale = MathHelper.InverseSqrtFast(X * X + Y * Y);

            X *= scale;

            Y *= scale;

        }

        #endregion

        #region public void Scale()

        /// <summary>

        /// Scales the current Vector2 by the given amounts.

        /// </summary>

        /// <param name="sx">The scale of the X component.</param>

        /// <param name="sy">The scale of the Y component.</param>

        [Obsolete("Use static Multiply() method instead.")]

        public void Scale(float sx, float sy)

        {

            this.X = X * sx;

            this.Y = Y * sy;

        }

        /// <summary>Scales this instance by the given parameter.</summary>

        /// <param name="scale">The scaling of the individual components.</param>

        [Obsolete("Use static Multiply() method instead.")]

        public void Scale(Vector2 scale)

        {

            this.X *= scale.X;

            this.Y *= scale.Y;

        }

        /// <summary>Scales this instance by the given parameter.</summary>

        /// <param name="scale">The scaling of the individual components.</param>

        [CLSCompliant(false)]

        [Obsolete("Use static Multiply() method instead.")]

        public void Scale(ref Vector2 scale)

        {

            this.X *= scale.X;

            this.Y *= scale.Y;

        }

        #endregion public void Scale()

        #endregion

        #region Static

        #region Fields

        /// <summary>

        /// Defines a unit-length Vector2 that points towards the X-axis.

        /// </summary>

        public static readonly Vector2 UnitX = new Vector2(1, 0);

        /// <summary>

        /// Defines a unit-length Vector2 that points towards the Y-axis.

        /// </summary>

        public static readonly Vector2 UnitY = new Vector2(0, 1);

        /// <summary>

        /// Defines a zero-length Vector2.

        /// </summary>

        public static readonly Vector2 Zero = new Vector2(0, 0);

        /// <summary>

        /// Defines an instance with all components set to 1.

        /// </summary>

        public static readonly Vector2 One = new Vector2(1, 1);

        /// <summary>

        /// Defines the size of the Vector2 struct in bytes.

        /// </summary>

        public static readonly int SizeInBytes = Marshal.SizeOf(new Vector2());

        #endregion

        #region Obsolete

        #region Sub

        /// <summary>

        /// Subtract one Vector from another

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <returns>Result of subtraction</returns>

        [Obsolete("Use static Subtract() method instead.")]

        public static Vector2 Sub(Vector2 a, Vector2 b)

        {

            a.X -= b.X;

            a.Y -= b.Y;

            return a;

        }

        /// <summary>

        /// Subtract one Vector from another

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <param name="result">Result of subtraction</param>

        [Obsolete("Use static Subtract() method instead.")]

        public static void Sub(ref Vector2 a, ref Vector2 b, out Vector2 result)

        {

            result.X = a.X - b.X;

            result.Y = a.Y - b.Y;

        }

        #endregion

        #region Mult

        /// <summary>

        /// Multiply a vector and a scalar

        /// </summary>

        /// <param name="a">Vector operand</param>

        /// <param name="f">Scalar operand</param>

        /// <returns>Result of the multiplication</returns>

        [Obsolete("Use static Multiply() method instead.")]

        public static Vector2 Mult(Vector2 a, float f)

        {

            a.X *= f;

            a.Y *= f;

            return a;

        }

        /// <summary>

        /// Multiply a vector and a scalar

        /// </summary>

        /// <param name="a">Vector operand</param>

        /// <param name="f">Scalar operand</param>

        /// <param name="result">Result of the multiplication</param>

        [Obsolete("Use static Multiply() method instead.")]

        public static void Mult(ref Vector2 a, float f, out Vector2 result)

        {

            result.X = a.X * f;

            result.Y = a.Y * f;

        }

        #endregion

        #region Div

        /// <summary>

        /// Divide a vector by a scalar

        /// </summary>

        /// <param name="a">Vector operand</param>

        /// <param name="f">Scalar operand</param>

        /// <returns>Result of the division</returns>

        [Obsolete("Use static Divide() method instead.")]

        public static Vector2 Div(Vector2 a, float f)

        {

            float mult = 1.0f / f;

            a.X *= mult;

            a.Y *= mult;

            return a;

        }

        /// <summary>

        /// Divide a vector by a scalar

        /// </summary>

        /// <param name="a">Vector operand</param>

        /// <param name="f">Scalar operand</param>

        /// <param name="result">Result of the division</param>

        [Obsolete("Use static Divide() method instead.")]

        public static void Div(ref Vector2 a, float f, out Vector2 result)

        {

            float mult = 1.0f / f;

            result.X = a.X * mult;

            result.Y = a.Y * mult;

        }

        #endregion

        #endregion

        #region Add

        /// <summary>

        /// Adds two vectors.

        /// </summary>

        /// <param name="a">Left operand.</param>

        /// <param name="b">Right operand.</param>

        /// <returns>Result of operation.</returns>

        public static Vector2 Add(Vector2 a, Vector2 b)

        {

            Add(ref a, ref b, out a);

            return a;

        }

        /// <summary>

        /// Adds two vectors.

        /// </summary>

        /// <param name="a">Left operand.</param>

        /// <param name="b">Right operand.</param>

        /// <param name="result">Result of operation.</param>

        public static void Add(ref Vector2 a, ref Vector2 b, out Vector2 result)

        {

            result = new Vector2(a.X + b.X, a.Y + b.Y);

        }

        #endregion

        #region Subtract

        /// <summary>

        /// Subtract one Vector from another

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <returns>Result of subtraction</returns>

        public static Vector2 Subtract(Vector2 a, Vector2 b)

        {

            Subtract(ref a, ref b, out a);

            return a;

        }

        /// <summary>

        /// Subtract one Vector from another

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <param name="result">Result of subtraction</param>

        public static void Subtract(ref Vector2 a, ref Vector2 b, out Vector2 result)

        {

            result = new Vector2(a.X - b.X, a.Y - b.Y);

        }

        #endregion

        #region Multiply

        /// <summary>

        /// Multiplies a vector by a scalar.

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <returns>Result of the operation.</returns>

        public static Vector2 Multiply(Vector2 vector, float scale)

        {

            Multiply(ref vector, scale, out vector);

            return vector;

        }

        /// <summary>

        /// Multiplies a vector by a scalar.

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <param name="result">Result of the operation.</param>

        public static void Multiply(ref Vector2 vector, float scale, out Vector2 result)

        {

            result = new Vector2(vector.X * scale, vector.Y * scale);

        }

        /// <summary>

        /// Multiplies a vector by the components a vector (scale).

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <returns>Result of the operation.</returns>

        public static Vector2 Multiply(Vector2 vector, Vector2 scale)

        {

            Multiply(ref vector, ref scale, out vector);

            return vector;

        }

        /// <summary>

        /// Multiplies a vector by the components of a vector (scale).

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <param name="result">Result of the operation.</param>

        public static void Multiply(ref Vector2 vector, ref Vector2 scale, out Vector2 result)

        {

            result = new Vector2(vector.X * scale.X, vector.Y * scale.Y);

        }

        #endregion

        #region Divide

        /// <summary>

        /// Divides a vector by a scalar.

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <returns>Result of the operation.</returns>

        public static Vector2 Divide(Vector2 vector, float scale)

        {

            Divide(ref vector, scale, out vector);

            return vector;

        }

        /// <summary>

        /// Divides a vector by a scalar.

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <param name="result">Result of the operation.</param>

        public static void Divide(ref Vector2 vector, float scale, out Vector2 result)

        {

            Multiply(ref vector, 1 / scale, out result);

        }

        /// <summary>

        /// Divides a vector by the components of a vector (scale).

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <returns>Result of the operation.</returns>

        public static Vector2 Divide(Vector2 vector, Vector2 scale)

        {

            Divide(ref vector, ref scale, out vector);

            return vector;

        }

        /// <summary>

        /// Divide a vector by the components of a vector (scale).

        /// </summary>

        /// <param name="vector">Left operand.</param>

        /// <param name="scale">Right operand.</param>

        /// <param name="result">Result of the operation.</param>

        public static void Divide(ref Vector2 vector, ref Vector2 scale, out Vector2 result)

        {

            result = new Vector2(vector.X / scale.X, vector.Y / scale.Y);

        }

        #endregion

        #region ComponentMin

        /// <summary>

        /// Calculate the component-wise minimum of two vectors

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <returns>The component-wise minimum</returns>

        public static Vector2 ComponentMin(Vector2 a, Vector2 b)

        {

            a.X = a.X < b.X ? a.X : b.X;

            a.Y = a.Y < b.Y ? a.Y : b.Y;

            return a;

        }

        /// <summary>

        /// Calculate the component-wise minimum of two vectors

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <param name="result">The component-wise minimum</param>

        public static void ComponentMin(ref Vector2 a, ref Vector2 b, out Vector2 result)

        {

            result.X = a.X < b.X ? a.X : b.X;

            result.Y = a.Y < b.Y ? a.Y : b.Y;

        }

        #endregion

        #region ComponentMax

        /// <summary>

        /// Calculate the component-wise maximum of two vectors

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <returns>The component-wise maximum</returns>

        public static Vector2 ComponentMax(Vector2 a, Vector2 b)

        {

            a.X = a.X > b.X ? a.X : b.X;

            a.Y = a.Y > b.Y ? a.Y : b.Y;

            return a;

        }

        /// <summary>

        /// Calculate the component-wise maximum of two vectors

        /// </summary>

        /// <param name="a">First operand</param>

        /// <param name="b">Second operand</param>

        /// <param name="result">The component-wise maximum</param>

        public static void ComponentMax(ref Vector2 a, ref Vector2 b, out Vector2 result)

        {

            result.X = a.X > b.X ? a.X : b.X;

            result.Y = a.Y > b.Y ? a.Y : b.Y;

        }

        #endregion

        #region Min

        /// <summary>

        /// Returns the Vector3 with the minimum magnitude

        /// </summary>

        /// <param name="left">Left operand</param>

        /// <param name="right">Right operand</param>

        /// <returns>The minimum Vector3</returns>

        public static Vector2 Min(Vector2 left, Vector2 right)

        {

            return left.LengthSquared < right.LengthSquared ? left : right;

        }

        #endregion

        #region Max

        /// <summary>

        /// Returns the Vector3 with the minimum magnitude

        /// </summary>

        /// <param name="left">Left operand</param>

        /// <param name="right">Right operand</param>

        /// <returns>The minimum Vector3</returns>

        public static Vector2 Max(Vector2 left, Vector2 right)

        {

            return left.LengthSquared >= right.LengthSquared ? left : right;

        }

        #endregion

        #region Clamp

        /// <summary>

        /// Clamp a vector to the given minimum and maximum vectors

        /// </summary>

        /// <param name="vec">Input vector</param>

        /// <param name="min">Minimum vector</param>

        /// <param name="max">Maximum vector</param>

        /// <returns>The clamped vector</returns>

        public static Vector2 Clamp(Vector2 vec, Vector2 min, Vector2 max)

        {

            vec.X = vec.X < min.X ? min.X : vec.X > max.X ? max.X : vec.X;

            vec.Y = vec.Y < min.Y ? min.Y : vec.Y > max.Y ? max.Y : vec.Y;

            return vec;

        }

        /// <summary>

        /// Clamp a vector to the given minimum and maximum vectors

        /// </summary>

        /// <param name="vec">Input vector</param>

        /// <param name="min">Minimum vector</param>

        /// <param name="max">Maximum vector</param>

        /// <param name="result">The clamped vector</param>

        public static void Clamp(ref Vector2 vec, ref Vector2 min, ref Vector2 max, out Vector2 result)

        {

            result.X = vec.X < min.X ? min.X : vec.X > max.X ? max.X : vec.X;

            result.Y = vec.Y < min.Y ? min.Y : vec.Y > max.Y ? max.Y : vec.Y;

        }

        #endregion

        #region Normalize

        /// <summary>

        /// Scale a vector to unit length

        /// </summary>

        /// <param name="vec">The input vector</param>

        /// <returns>The normalized vector</returns>

        public static Vector2 Normalize(Vector2 vec)

        {

            float scale = 1.0f / vec.Length;

            vec.X *= scale;

            vec.Y *= scale;

            return vec;

        }

        /// <summary>

        /// Scale a vector to unit length

        /// </summary>

        /// <param name="vec">The input vector</param>

        /// <param name="result">The normalized vector</param>

        public static void Normalize(ref Vector2 vec, out Vector2 result)

        {

            float scale = 1.0f / vec.Length;

            result.X = vec.X * scale;

            result.Y = vec.Y * scale;

        }

        #endregion

        #region NormalizeFast

        /// <summary>

        /// Scale a vector to approximately unit length

        /// </summary>

        /// <param name="vec">The input vector</param>

        /// <returns>The normalized vector</returns>

        public static Vector2 NormalizeFast(Vector2 vec)

        {

            float scale = MathHelper.InverseSqrtFast(vec.X * vec.X + vec.Y * vec.Y);

            vec.X *= scale;

            vec.Y *= scale;

            return vec;

        }

        /// <summary>

        /// Scale a vector to approximately unit length

        /// </summary>

        /// <param name="vec">The input vector</param>

        /// <param name="result">The normalized vector</param>

        public static void NormalizeFast(ref Vector2 vec, out Vector2 result)

        {

            float scale = MathHelper.InverseSqrtFast(vec.X * vec.X + vec.Y * vec.Y);

            result.X = vec.X * scale;

            result.Y = vec.Y * scale;

        }

        #endregion

        #region Dot

        /// <summary>

        /// Calculate the dot (scalar) product of two vectors

        /// </summary>

        /// <param name="left">First operand</param>

        /// <param name="right">Second operand</param>

        /// <returns>The dot product of the two inputs</returns>

        public static float Dot(Vector2 left, Vector2 right)

        {

            return left.X * right.X + left.Y * right.Y;

        }

        /// <summary>

        /// Calculate the dot (scalar) product of two vectors

        /// </summary>

        /// <param name="left">First operand</param>

        /// <param name="right">Second operand</param>

        /// <param name="result">The dot product of the two inputs</param>

        public static void Dot(ref Vector2 left, ref Vector2 right, out float result)

        {

            result = left.X * right.X + left.Y * right.Y;

        }

        #endregion

        #region Lerp

        /// <summary>

        /// Returns a new Vector that is the linear blend of the 2 given Vectors

        /// </summary>

        /// <param name="a">First input vector</param>

        /// <param name="b">Second input vector</param>

        /// <param name="blend">The blend factor. a when blend=0, b when blend=1.</param>

        /// <returns>a when blend=0, b when blend=1, and a linear combination otherwise</returns>

        public static Vector2 Lerp(Vector2 a, Vector2 b, float blend)

        {

            a.X = blend * (b.X - a.X) + a.X;

            a.Y = blend * (b.Y - a.Y) + a.Y;

            return a;

        }

        /// <summary>

        /// Returns a new Vector that is the linear blend of the 2 given Vectors

        /// </summary>

        /// <param name="a">First input vector</param>

        /// <param name="b">Second input vector</param>

        /// <param name="blend">The blend factor. a when blend=0, b when blend=1.</param>

        /// <param name="result">a when blend=0, b when blend=1, and a linear combination otherwise</param>

        public static void Lerp(ref Vector2 a, ref Vector2 b, float blend, out Vector2 result)

        {

            result.X = blend * (b.X - a.X) + a.X;

            result.Y = blend * (b.Y - a.Y) + a.Y;

        }

        #endregion

        #region Barycentric

        /// <summary>

        /// Interpolate 3 Vectors using Barycentric coordinates

        /// </summary>

        /// <param name="a">First input Vector</param>

        /// <param name="b">Second input Vector</param>

        /// <param name="c">Third input Vector</param>

        /// <param name="u">First Barycentric Coordinate</param>

        /// <param name="v">Second Barycentric Coordinate</param>

        /// <returns>a when u=v=0, b when u=1,v=0, c when u=0,v=1, and a linear combination of a,b,c otherwise</returns>

        public static Vector2 BaryCentric(Vector2 a, Vector2 b, Vector2 c, float u, float v)

        {

            return a + u * (b - a) + v * (c - a);

        }