/*
  Copyright (C) 2006, 2007 Sony Computer Entertainment Inc.
  All rights reserved.

  Redistribution and use in source and binary forms,
  with or without modification, are permitted provided that the
  following conditions are met:
   * Redistributions of source code must retain the above copyright
     notice, this list of conditions and the following disclaimer.
   * Redistributions in binary form must reproduce the above copyright
     notice, this list of conditions and the following disclaimer in the
     documentation and/or other materials provided with the distribution.
   * Neither the name of the Sony Computer Entertainment Inc nor the names
     of its contributors may be used to endorse or promote products derived
     from this software without specific prior written permission.

  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
  AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  POSSIBILITY OF SUCH DAMAGE.
*/

#define _VECTORMATH_SOA_TEST

#include "vectormath_soa.h"
#include "test.h"

int iteration = 0;

using namespace Vectormath;
using namespace Vectormath::Soa;

void
Vector3_methods_test()
{
    Vector3 a_Vector3, b_Vector3, c_Vector3, d_Vector3, e_Vector3;
    Vector4 a_Vector4, b_Vector4, c_Vector4, d_Vector4, e_Vector4;
    Point3 a_Point3, b_Point3, c_Point3, d_Point3, e_Point3;
    Quat a_Quat, b_Quat, c_Quat, d_Quat, e_Quat;
    Aos::Vector3 aos_Vector3_0, aos_Vector3_1, aos_Vector3_2, aos_Vector3_3;
    Vector3 soa_Vector3;
    Aos::Vector4 aos_Vector4_0, aos_Vector4_1, aos_Vector4_2;
    vec_float4 rndflt1, rndflt2, rndflt3, rndflt4;
    float xyz4[12] __attribute__ ((aligned(16)));
    xyz4[0] = getfloat(randfloat());
    xyz4[1] = getfloat(randfloat());
    xyz4[2] = getfloat(randfloat());
    xyz4[3] = getfloat(randfloat());
    xyz4[4] = getfloat(randfloat());
    xyz4[5] = getfloat(randfloat());
    xyz4[6] = getfloat(randfloat());
    xyz4[7] = getfloat(randfloat());
    xyz4[8] = getfloat(randfloat());
    xyz4[9] = getfloat(randfloat());
    xyz4[10] = getfloat(randfloat());
    xyz4[11] = getfloat(randfloat());
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    print( a_Vector3, "set Vector3 with floats" );
    print( b_Vector3, "set Vector3 with floats" );
    c_Vector3 = Vector3( (vec_float4){0.0f} );
    d_Vector3 = Vector3( (vec_float4){0.0f} );
    e_Vector3 = Vector3( (vec_float4){0.0f} );
    print( c_Vector3, "set Vector3 elements to zero" );
    print( d_Vector3, "set Vector3 elements to zero" );
    print( e_Vector3, "set Vector3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Vector4, "set Vector4 with floats" );
    print( b_Vector4, "set Vector4 with floats" );
    c_Vector4 = Vector4( (vec_float4){0.0f} );
    d_Vector4 = Vector4( (vec_float4){0.0f} );
    e_Vector4 = Vector4( (vec_float4){0.0f} );
    print( c_Vector4, "set Vector4 elements to zero" );
    print( d_Vector4, "set Vector4 elements to zero" );
    print( e_Vector4, "set Vector4 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    print( a_Point3, "set Point3 with floats" );
    print( b_Point3, "set Point3 with floats" );
    c_Point3 = Point3( (vec_float4){0.0f} );
    d_Point3 = Point3( (vec_float4){0.0f} );
    e_Point3 = Point3( (vec_float4){0.0f} );
    print( c_Point3, "set Point3 elements to zero" );
    print( d_Point3, "set Point3 elements to zero" );
    print( e_Point3, "set Point3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Quat, "set Quat with floats" );
    print( b_Quat, "set Quat with floats" );
    c_Quat = Quat( (vec_float4){0.0f} );
    d_Quat = Quat( (vec_float4){0.0f} );
    e_Quat = Quat( (vec_float4){0.0f} );
    print( c_Quat, "set Quat elements to zero" );
    print( d_Quat, "set Quat elements to zero" );
    print( e_Quat, "set Quat elements to zero" );
    a_Vector3 = Vector3( a_Point3 );
    print( a_Vector3, "construct Vector3 with Point3" );
    a_Vector3 = Vector3( randfloat() );
    print( a_Vector3, "set Vector3 with float" );
    a_Vector3 = Vector3( randfloat() );
    print( a_Vector3, "set Vector3 with float" );
    aos_Vector3_0 = Aos::Vector3( 0.0f, 1.0f, 2.0f );
    aos_Vector3_1 = Aos::Vector3( 3.0f, 4.0f, 5.0f );
    aos_Vector3_2 = Aos::Vector3( 6.0f, 7.0f, 8.0f );
    aos_Vector3_3 = Aos::Vector3( 9.0f, 10.0f, 11.0f );
    soa_Vector3 = Vector3( aos_Vector3_0, aos_Vector3_1, aos_Vector3_2, aos_Vector3_3 );
    soa_Vector3.get4Aos( aos_Vector3_3, aos_Vector3_2, aos_Vector3_1, aos_Vector3_0 );
    Aos::print( aos_Vector3_0, "aos type 0" );
    Aos::print( aos_Vector3_1, "aos type 1" );
    Aos::print( aos_Vector3_2, "aos type 2" );
    Aos::print( aos_Vector3_3, "aos type 3" );
    a_Vector3 = select( a_Vector3, b_Vector3, ((vec_uint4){0,0xffffffff,0,0xffffffff}) );
    a_Vector3.get4Aos( aos_Vector3_0, aos_Vector3_1, aos_Vector3_2, aos_Vector3_3 );
    Aos::print( aos_Vector3_0, "select 0" );
    Aos::print( aos_Vector3_1, "select 1" );
    Aos::print( aos_Vector3_2, "select 2" );
    Aos::print( aos_Vector3_3, "select 3" );
    loadXYZArray( a_Vector3, (const vec_float4 *)xyz4 );
    print( a_Vector3, "load XYZ array" );
    a_Vector3 = Vector3( ( -Vector3( a_Vector3 ) ) );
    storeXYZArray( a_Vector3, (vec_float4 *)xyz4 );
    aos_Vector4_0 = Aos::Vector4( xyz4[0], xyz4[1], xyz4[2], xyz4[3] );
    aos_Vector4_1 = Aos::Vector4( xyz4[4], xyz4[5], xyz4[6], xyz4[7] );
    aos_Vector4_2 = Aos::Vector4( xyz4[8], xyz4[9], xyz4[10], xyz4[11] );
    Aos::print( aos_Vector4_0, "xyzx" );
    Aos::print( aos_Vector4_1, "yzxy" );
    Aos::print( aos_Vector4_2, "zxyz" );
    Aos::print( aos_Vector4_0, "xyzx" );
    Aos::print( aos_Vector4_1, "yzxy" );
    Aos::print( aos_Vector4_2, "zxyz" );
    printf("storeXYZ:-1.0 -2.0 -3.0 0.4\n");
    a_Vector3 = b_Vector3;
    print( a_Vector3, "assign to Vector3 from Vector3" );
    a_Vector3 = Vector3( (vec_float4){0.0f} );
    print( a_Vector3, "set Vector3 elements to zero" );
    a_Vector3 = Vector3::xAxis( );
    print( a_Vector3, "set to x axis" );
    a_Vector3 = Vector3::yAxis( );
    print( a_Vector3, "set to y axis" );
    a_Vector3 = Vector3::zAxis( );
    print( a_Vector3, "set to z axis" );
    a_Vector3.setElem( 0, randfloat() );
    print( a_Vector3, "Vector3::set( 0, float )" );
    a_Vector3[0] = randfloat();
    a_Vector3[0] = vec_mul_float( a_Vector3[0], randfloat() );
    a_Vector3[0] = divf4( a_Vector3[0], randfloat() );
    a_Vector3[0] = vec_add_float( a_Vector3[0], randfloat() );
    a_Vector3[0] = vec_sub_float( a_Vector3[0], randfloat() );
    print( a_Vector3, "Vector3::operator [](0)" );
    a_Vector3.setX( randfloat() );
    print( a_Vector3, "Vector3::setX()" );
    a_Vector3.setElem( 1, randfloat() );
    print( a_Vector3, "Vector3::set( 1, float )" );
    a_Vector3[1] = randfloat();
    a_Vector3[1] = vec_mul_float( a_Vector3[1], randfloat() );
    a_Vector3[1] = divf4( a_Vector3[1], randfloat() );
    a_Vector3[1] = vec_add_float( a_Vector3[1], randfloat() );
    a_Vector3[1] = vec_sub_float( a_Vector3[1], randfloat() );
    print( a_Vector3, "Vector3::operator [](1)" );
    a_Vector3.setY( randfloat() );
    print( a_Vector3, "Vector3::setY()" );
    a_Vector3.setElem( 2, randfloat() );
    print( a_Vector3, "Vector3::set( 2, float )" );
    a_Vector3[2] = randfloat();
    a_Vector3[2] = vec_mul_float( a_Vector3[2], randfloat() );
    a_Vector3[2] = divf4( a_Vector3[2], randfloat() );
    a_Vector3[2] = vec_add_float( a_Vector3[2], randfloat() );
    a_Vector3[2] = vec_sub_float( a_Vector3[2], randfloat() );
    print( a_Vector3, "Vector3::operator [](2)" );
    a_Vector3.setZ( randfloat() );
    print( a_Vector3, "Vector3::setZ()" );
    printf("Vector3::get( 0 ): %f\n", getfloat(a_Vector3.getElem( 0 )) );
    printf("Vector3::operator []( 0 ): %f\n", getfloat((vec_float4)a_Vector3[0]) );
    printf("Vector3::getX(): %f\n", getfloat(a_Vector3.getX( )) );
    printf("Vector3::get( 1 ): %f\n", getfloat(a_Vector3.getElem( 1 )) );
    printf("Vector3::operator []( 1 ): %f\n", getfloat((vec_float4)a_Vector3[1]) );
    printf("Vector3::getY(): %f\n", getfloat(a_Vector3.getY( )) );
    printf("Vector3::get( 2 ): %f\n", getfloat(a_Vector3.getElem( 2 )) );
    printf("Vector3::operator []( 2 ): %f\n", getfloat((vec_float4)a_Vector3[2]) );
    printf("Vector3::getZ(): %f\n", getfloat(a_Vector3.getZ( )) );
    print( ( a_Vector3 + b_Vector3 ), "Vector3 + Vector3" );
    print( ( a_Vector3 - b_Vector3 ), "Vector3 - Vector3" );
    print( ( a_Vector3 + b_Point3 ), "Vector3 + Point3" );
    print( ( a_Vector3 * randfloat() ), "Vector3 * float" );
    print( ( a_Vector3 / randfloat() ), "Vector3 / float" );
    print( ( randfloat() * a_Vector3 ), "float * Vector3" );
    print( ( -a_Vector3 ), "Vector3 negate" );
    print( mulPerElem( a_Vector3, b_Vector3 ), "mulPerElem( Vector3, Vector3 )" );
    print( divPerElem( a_Vector3, b_Vector3 ), "divPerElem( Vector3, Vector3 )" );
    print( recipPerElem( a_Vector3 ), "Vector3 recip" );
    print( sqrtPerElem( absPerElem( a_Vector3 ) ), "Vector3 sqrt" );
    print( rsqrtPerElem( absPerElem( a_Vector3 ) ), "Vector3 rsqrt" );
    print( absPerElem( a_Vector3 ), "Vector3 abs" );
    print( copySignPerElem( a_Vector3, b_Vector3 ), "Vector3 copySign" );
    print( maxPerElem( a_Vector3, b_Vector3 ), "Vector3 maximum Vector3" );
    print( minPerElem( a_Vector3, b_Vector3 ), "Vector3 minimum Vector3" );
    printf("Vector3 maximum of elements: %f\n", getfloat(maxElem( a_Vector3 )));
    printf("Vector3 minimum of elements: %f\n", getfloat(minElem( a_Vector3 )));
    printf("Vector3 sum of elements: %f\n", getfloat(sum( a_Vector3 )));
    printf("Vector3 dot Vector3: %f\n", getfloat(dot( a_Vector3, b_Vector3 )));
    printf("Vector3 lengthSqr: %f\n", getfloat(lengthSqr( a_Vector3 )));
    printf("Vector3 length: %f\n", getfloat(length( a_Vector3 )));
    print( normalize( a_Vector3 ), "Vector3 normalized" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    c_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    d_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    e_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    b_Vector3 = normalize( b_Vector3 );
    c_Vector3 = normalize( c_Vector3 );
    d_Vector3 = normalize( d_Vector3 );
    e_Vector3 = normalize( e_Vector3 );
    a_Vector3 = lerp( randfloat(), b_Vector3, c_Vector3 );
    print( a_Vector3, "Vector3 lerp" );
    a_Vector3 = slerp( randfloat(), b_Vector3, c_Vector3 );
    print( a_Vector3, "Vector3 slerp" );
}

void
Vector4_methods_test()
{
    Vector3 a_Vector3, b_Vector3, c_Vector3, d_Vector3, e_Vector3;
    Vector4 a_Vector4, b_Vector4, c_Vector4, d_Vector4, e_Vector4;
    Point3 a_Point3, b_Point3, c_Point3, d_Point3, e_Point3;
    Quat a_Quat, b_Quat, c_Quat, d_Quat, e_Quat;
    Aos::Vector4 aos_Vector4_0, aos_Vector4_1, aos_Vector4_2, aos_Vector4_3;
    Vector4 soa_Vector4;
    vec_float4 rndflt1, rndflt2, rndflt3, rndflt4;
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    print( a_Vector3, "set Vector3 with floats" );
    print( b_Vector3, "set Vector3 with floats" );
    c_Vector3 = Vector3( (vec_float4){0.0f} );
    d_Vector3 = Vector3( (vec_float4){0.0f} );
    e_Vector3 = Vector3( (vec_float4){0.0f} );
    print( c_Vector3, "set Vector3 elements to zero" );
    print( d_Vector3, "set Vector3 elements to zero" );
    print( e_Vector3, "set Vector3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Vector4, "set Vector4 with floats" );
    print( b_Vector4, "set Vector4 with floats" );
    c_Vector4 = Vector4( (vec_float4){0.0f} );
    d_Vector4 = Vector4( (vec_float4){0.0f} );
    e_Vector4 = Vector4( (vec_float4){0.0f} );
    print( c_Vector4, "set Vector4 elements to zero" );
    print( d_Vector4, "set Vector4 elements to zero" );
    print( e_Vector4, "set Vector4 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    print( a_Point3, "set Point3 with floats" );
    print( b_Point3, "set Point3 with floats" );
    c_Point3 = Point3( (vec_float4){0.0f} );
    d_Point3 = Point3( (vec_float4){0.0f} );
    e_Point3 = Point3( (vec_float4){0.0f} );
    print( c_Point3, "set Point3 elements to zero" );
    print( d_Point3, "set Point3 elements to zero" );
    print( e_Point3, "set Point3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Quat, "set Quat with floats" );
    print( b_Quat, "set Quat with floats" );
    c_Quat = Quat( (vec_float4){0.0f} );
    d_Quat = Quat( (vec_float4){0.0f} );
    e_Quat = Quat( (vec_float4){0.0f} );
    print( c_Quat, "set Quat elements to zero" );
    print( d_Quat, "set Quat elements to zero" );
    print( e_Quat, "set Quat elements to zero" );
    a_Vector4 = Vector4( a_Vector3, randfloat() );
    print( a_Vector4, "set Vector4 with Vector3, float" );
    a_Vector4 = Vector4( a_Vector3 );
    print( a_Vector4, "set Vector4 with Vector3" );
    a_Vector4 = Vector4( a_Point3 );
    print( a_Vector4, "set Vector4 with Point3" );
    a_Vector4 = Vector4( a_Quat );
    print( a_Vector4, "construct Vector4 with Quat" );
    a_Vector4 = Vector4( randfloat() );
    print( a_Vector4, "set Vector4 with float" );
    a_Vector4 = Vector4( randfloat() );
    print( a_Vector4, "set Vector4 with float" );
    aos_Vector4_0 = Aos::Vector4( 0.0f, 1.0f, 2.0f, 3.0f );
    aos_Vector4_1 = Aos::Vector4( 4.0f, 5.0f, 6.0f, 7.0f );
    aos_Vector4_2 = Aos::Vector4( 8.0f, 9.0f, 10.0f, 11.0f );
    aos_Vector4_3 = Aos::Vector4( 12.0f, 13.0f, 14.0f, 15.0f );
    soa_Vector4 = Vector4( aos_Vector4_0, aos_Vector4_1, aos_Vector4_2, aos_Vector4_3 );
    soa_Vector4.get4Aos( aos_Vector4_3, aos_Vector4_2, aos_Vector4_1, aos_Vector4_0 );
    Aos::print( aos_Vector4_0, "aos type 0" );
    Aos::print( aos_Vector4_1, "aos type 1" );
    Aos::print( aos_Vector4_2, "aos type 2" );
    Aos::print( aos_Vector4_3, "aos type 3" );
    a_Vector4 = select( a_Vector4, b_Vector4, ((vec_uint4){0,0xffffffff,0,0xffffffff}) );
    a_Vector4.get4Aos( aos_Vector4_0, aos_Vector4_1, aos_Vector4_2, aos_Vector4_3 );
    Aos::print( aos_Vector4_0, "select 0" );
    Aos::print( aos_Vector4_1, "select 1" );
    Aos::print( aos_Vector4_2, "select 2" );
    Aos::print( aos_Vector4_3, "select 3" );
    a_Vector4 = b_Vector4;
    print( a_Vector4, "assign to Vector4 from Vector4" );
    a_Vector4.setXYZ( a_Vector3 );
    print( a_Vector4, "set Vector4 xyz" );
    print( a_Vector4.getXYZ( ), "get Vector4 xyz" );
    a_Vector4 = Vector4( (vec_float4){0.0f} );
    print( a_Vector4, "set Vector4 elements to zero" );
    a_Vector4 = Vector4::xAxis( );
    print( a_Vector4, "set to x axis" );
    a_Vector4 = Vector4::yAxis( );
    print( a_Vector4, "set to y axis" );
    a_Vector4 = Vector4::zAxis( );
    print( a_Vector4, "set to z axis" );
    a_Vector4 = Vector4::wAxis( );
    print( a_Vector4, "set to w axis" );
    a_Vector4.setElem( 0, randfloat() );
    print( a_Vector4, "Vector4::set( 0, float )" );
    a_Vector4[0] = randfloat();
    a_Vector4[0] = vec_mul_float( a_Vector4[0], randfloat() );
    a_Vector4[0] = divf4( a_Vector4[0], randfloat() );
    a_Vector4[0] = vec_add_float( a_Vector4[0], randfloat() );
    a_Vector4[0] = vec_sub_float( a_Vector4[0], randfloat() );
    print( a_Vector4, "Vector4::operator [](0)" );
    a_Vector4.setX( randfloat() );
    print( a_Vector4, "Vector4::setX()" );
    a_Vector4.setElem( 1, randfloat() );
    print( a_Vector4, "Vector4::set( 1, float )" );
    a_Vector4[1] = randfloat();
    a_Vector4[1] = vec_mul_float( a_Vector4[1], randfloat() );
    a_Vector4[1] = divf4( a_Vector4[1], randfloat() );
    a_Vector4[1] = vec_add_float( a_Vector4[1], randfloat() );
    a_Vector4[1] = vec_sub_float( a_Vector4[1], randfloat() );
    print( a_Vector4, "Vector4::operator [](1)" );
    a_Vector4.setY( randfloat() );
    print( a_Vector4, "Vector4::setY()" );
    a_Vector4.setElem( 2, randfloat() );
    print( a_Vector4, "Vector4::set( 2, float )" );
    a_Vector4[2] = randfloat();
    a_Vector4[2] = vec_mul_float( a_Vector4[2], randfloat() );
    a_Vector4[2] = divf4( a_Vector4[2], randfloat() );
    a_Vector4[2] = vec_add_float( a_Vector4[2], randfloat() );
    a_Vector4[2] = vec_sub_float( a_Vector4[2], randfloat() );
    print( a_Vector4, "Vector4::operator [](2)" );
    a_Vector4.setZ( randfloat() );
    print( a_Vector4, "Vector4::setZ()" );
    a_Vector4.setElem( 3, randfloat() );
    print( a_Vector4, "Vector4::set( 3, float )" );
    a_Vector4[3] = randfloat();
    a_Vector4[3] = vec_mul_float( a_Vector4[3], randfloat() );
    a_Vector4[3] = divf4( a_Vector4[3], randfloat() );
    a_Vector4[3] = vec_add_float( a_Vector4[3], randfloat() );
    a_Vector4[3] = vec_sub_float( a_Vector4[3], randfloat() );
    print( a_Vector4, "Vector4::operator [](3)" );
    a_Vector4.setW( randfloat() );
    print( a_Vector4, "Vector4::setW()" );
    printf("Vector4::get( 0 ): %f\n", getfloat(a_Vector4.getElem( 0 )) );
    printf("Vector4::operator []( 0 ): %f\n", getfloat((vec_float4)a_Vector4[0]) );
    printf("Vector4::getX(): %f\n", getfloat(a_Vector4.getX( )) );
    printf("Vector4::get( 1 ): %f\n", getfloat(a_Vector4.getElem( 1 )) );
    printf("Vector4::operator []( 1 ): %f\n", getfloat((vec_float4)a_Vector4[1]) );
    printf("Vector4::getY(): %f\n", getfloat(a_Vector4.getY( )) );
    printf("Vector4::get( 2 ): %f\n", getfloat(a_Vector4.getElem( 2 )) );
    printf("Vector4::operator []( 2 ): %f\n", getfloat((vec_float4)a_Vector4[2]) );
    printf("Vector4::getZ(): %f\n", getfloat(a_Vector4.getZ( )) );
    printf("Vector4::get( 3 ): %f\n", getfloat(a_Vector4.getElem( 3 )) );
    printf("Vector4::operator []( 3 ): %f\n", getfloat((vec_float4)a_Vector4[3]) );
    printf("Vector4::getW(): %f\n", getfloat(a_Vector4.getW( )) );
    print( ( a_Vector4 + b_Vector4 ), "Vector4 + Vector4" );
    print( ( a_Vector4 - b_Vector4 ), "Vector4 - Vector4" );
    print( ( a_Vector4 * randfloat() ), "Vector4 * float" );
    print( ( a_Vector4 / randfloat() ), "Vector4 / float" );
    print( ( randfloat() * a_Vector4 ), "float * Vector4" );
    print( ( -a_Vector4 ), "Vector4 negate" );
    print( mulPerElem( a_Vector4, b_Vector4 ), "mulPerElem( Vector4, Vector4 )" );
    print( divPerElem( a_Vector4, b_Vector4 ), "divPerElem( Vector4, Vector4 )" );
    print( recipPerElem( a_Vector4 ), "Vector4 recip" );
    print( sqrtPerElem( absPerElem( a_Vector4 ) ), "Vector4 sqrt" );
    print( rsqrtPerElem( absPerElem( a_Vector4 ) ), "Vector4 rsqrt" );
    print( absPerElem( a_Vector4 ), "Vector4 abs" );
    print( copySignPerElem( a_Vector4, b_Vector4 ), "Vector4 copySign" );
    print( maxPerElem( a_Vector4, b_Vector4 ), "Vector4 maximum Vector4" );
    print( minPerElem( a_Vector4, b_Vector4 ), "Vector4 minimum Vector4" );
    printf("Vector4 maximum of elements: %f\n", getfloat(maxElem( a_Vector4 )));
    printf("Vector4 minimum of elements: %f\n", getfloat(minElem( a_Vector4 )));
    printf("Vector4 sum of elements: %f\n", getfloat(sum( a_Vector4 )));
    printf("Vector4 dot Vector4: %f\n", getfloat(dot( a_Vector4, b_Vector4 )));
    printf("Vector4 lengthSqr: %f\n", getfloat(lengthSqr( a_Vector4 )));
    printf("Vector4 length: %f\n", getfloat(length( a_Vector4 )));
    print( normalize( a_Vector4 ), "Vector4 normalized" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    c_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    d_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    e_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    b_Vector4 = normalize( b_Vector4 );
    c_Vector4 = normalize( c_Vector4 );
    d_Vector4 = normalize( d_Vector4 );
    e_Vector4 = normalize( e_Vector4 );
    a_Vector4 = lerp( randfloat(), b_Vector4, c_Vector4 );
    print( a_Vector4, "Vector4 lerp" );
    a_Vector4 = slerp( randfloat(), b_Vector4, c_Vector4 );
    print( a_Vector4, "Vector4 slerp" );
}

void
Point3_methods_test()
{
    Vector3 a_Vector3, b_Vector3, c_Vector3, d_Vector3, e_Vector3;
    Vector4 a_Vector4, b_Vector4, c_Vector4, d_Vector4, e_Vector4;
    Point3 a_Point3, b_Point3, c_Point3, d_Point3, e_Point3;
    Quat a_Quat, b_Quat, c_Quat, d_Quat, e_Quat;
    Aos::Point3 aos_Point3_0, aos_Point3_1, aos_Point3_2, aos_Point3_3;
    Point3 soa_Point3;
    Aos::Vector4 aos_Vector4_0, aos_Vector4_1, aos_Vector4_2;
    vec_float4 rndflt1, rndflt2, rndflt3, rndflt4;
    float xyz4[12] __attribute__ ((aligned(16)));
    xyz4[0] = getfloat(randfloat());
    xyz4[1] = getfloat(randfloat());
    xyz4[2] = getfloat(randfloat());
    xyz4[3] = getfloat(randfloat());
    xyz4[4] = getfloat(randfloat());
    xyz4[5] = getfloat(randfloat());
    xyz4[6] = getfloat(randfloat());
    xyz4[7] = getfloat(randfloat());
    xyz4[8] = getfloat(randfloat());
    xyz4[9] = getfloat(randfloat());
    xyz4[10] = getfloat(randfloat());
    xyz4[11] = getfloat(randfloat());
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    print( a_Vector3, "set Vector3 with floats" );
    print( b_Vector3, "set Vector3 with floats" );
    c_Vector3 = Vector3( (vec_float4){0.0f} );
    d_Vector3 = Vector3( (vec_float4){0.0f} );
    e_Vector3 = Vector3( (vec_float4){0.0f} );
    print( c_Vector3, "set Vector3 elements to zero" );
    print( d_Vector3, "set Vector3 elements to zero" );
    print( e_Vector3, "set Vector3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Vector4, "set Vector4 with floats" );
    print( b_Vector4, "set Vector4 with floats" );
    c_Vector4 = Vector4( (vec_float4){0.0f} );
    d_Vector4 = Vector4( (vec_float4){0.0f} );
    e_Vector4 = Vector4( (vec_float4){0.0f} );
    print( c_Vector4, "set Vector4 elements to zero" );
    print( d_Vector4, "set Vector4 elements to zero" );
    print( e_Vector4, "set Vector4 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    print( a_Point3, "set Point3 with floats" );
    print( b_Point3, "set Point3 with floats" );
    c_Point3 = Point3( (vec_float4){0.0f} );
    d_Point3 = Point3( (vec_float4){0.0f} );
    e_Point3 = Point3( (vec_float4){0.0f} );
    print( c_Point3, "set Point3 elements to zero" );
    print( d_Point3, "set Point3 elements to zero" );
    print( e_Point3, "set Point3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Quat, "set Quat with floats" );
    print( b_Quat, "set Quat with floats" );
    c_Quat = Quat( (vec_float4){0.0f} );
    d_Quat = Quat( (vec_float4){0.0f} );
    e_Quat = Quat( (vec_float4){0.0f} );
    print( c_Quat, "set Quat elements to zero" );
    print( d_Quat, "set Quat elements to zero" );
    print( e_Quat, "set Quat elements to zero" );
    a_Point3 = Point3( a_Vector3 );
    print( a_Point3, "construct Point3 with Vector3" );
    a_Point3 = Point3( randfloat() );
    print( a_Point3, "set Point3 with float" );
    a_Point3 = Point3( randfloat() );
    print( a_Point3, "set Point3 with float" );
    aos_Point3_0 = Aos::Point3( 0.0f, 1.0f, 2.0f );
    aos_Point3_1 = Aos::Point3( 3.0f, 4.0f, 5.0f );
    aos_Point3_2 = Aos::Point3( 6.0f, 7.0f, 8.0f );
    aos_Point3_3 = Aos::Point3( 9.0f, 10.0f, 11.0f );
    soa_Point3 = Point3( aos_Point3_0, aos_Point3_1, aos_Point3_2, aos_Point3_3 );
    soa_Point3.get4Aos( aos_Point3_3, aos_Point3_2, aos_Point3_1, aos_Point3_0 );
    Aos::print( aos_Point3_0, "aos type 0" );
    Aos::print( aos_Point3_1, "aos type 1" );
    Aos::print( aos_Point3_2, "aos type 2" );
    Aos::print( aos_Point3_3, "aos type 3" );
    a_Point3 = select( a_Point3, b_Point3, ((vec_uint4){0,0xffffffff,0,0xffffffff}) );
    a_Point3.get4Aos( aos_Point3_0, aos_Point3_1, aos_Point3_2, aos_Point3_3 );
    Aos::print( aos_Point3_0, "select 0" );
    Aos::print( aos_Point3_1, "select 1" );
    Aos::print( aos_Point3_2, "select 2" );
    Aos::print( aos_Point3_3, "select 3" );
    loadXYZArray( a_Point3, (const vec_float4 *)xyz4 );
    print( a_Point3, "load XYZ array" );
    a_Point3 = Point3( ( -Vector3( a_Point3 ) ) );
    storeXYZArray( a_Point3, (vec_float4 *)xyz4 );
    aos_Vector4_0 = Aos::Vector4( xyz4[0], xyz4[1], xyz4[2], xyz4[3] );
    aos_Vector4_1 = Aos::Vector4( xyz4[4], xyz4[5], xyz4[6], xyz4[7] );
    aos_Vector4_2 = Aos::Vector4( xyz4[8], xyz4[9], xyz4[10], xyz4[11] );
    Aos::print( aos_Vector4_0, "xyzx" );
    Aos::print( aos_Vector4_1, "yzxy" );
    Aos::print( aos_Vector4_2, "zxyz" );
    Aos::print( aos_Vector4_0, "xyzx" );
    Aos::print( aos_Vector4_1, "yzxy" );
    Aos::print( aos_Vector4_2, "zxyz" );
    printf("storeXYZ:-1.0 -2.0 -3.0 0.4\n");
    a_Point3 = b_Point3;
    print( a_Point3, "assign to Point3 from Point3" );
    a_Point3 = Point3( (vec_float4){0.0f} );
    print( a_Point3, "set Point3 elements to zero" );
    a_Point3.setElem( 0, randfloat() );
    print( a_Point3, "Point3::set( 0, float )" );
    a_Point3[0] = randfloat();
    a_Point3[0] = vec_mul_float( a_Point3[0], randfloat() );
    a_Point3[0] = divf4( a_Point3[0], randfloat() );
    a_Point3[0] = vec_add_float( a_Point3[0], randfloat() );
    a_Point3[0] = vec_sub_float( a_Point3[0], randfloat() );
    print( a_Point3, "Point3::operator [](0)" );
    a_Point3.setX( randfloat() );
    print( a_Point3, "Point3::setX()" );
    a_Point3.setElem( 1, randfloat() );
    print( a_Point3, "Point3::set( 1, float )" );
    a_Point3[1] = randfloat();
    a_Point3[1] = vec_mul_float( a_Point3[1], randfloat() );
    a_Point3[1] = divf4( a_Point3[1], randfloat() );
    a_Point3[1] = vec_add_float( a_Point3[1], randfloat() );
    a_Point3[1] = vec_sub_float( a_Point3[1], randfloat() );
    print( a_Point3, "Point3::operator [](1)" );
    a_Point3.setY( randfloat() );
    print( a_Point3, "Point3::setY()" );
    a_Point3.setElem( 2, randfloat() );
    print( a_Point3, "Point3::set( 2, float )" );
    a_Point3[2] = randfloat();
    a_Point3[2] = vec_mul_float( a_Point3[2], randfloat() );
    a_Point3[2] = divf4( a_Point3[2], randfloat() );
    a_Point3[2] = vec_add_float( a_Point3[2], randfloat() );
    a_Point3[2] = vec_sub_float( a_Point3[2], randfloat() );
    print( a_Point3, "Point3::operator [](2)" );
    a_Point3.setZ( randfloat() );
    print( a_Point3, "Point3::setZ()" );
    printf("Point3::get( 0 ): %f\n", getfloat(a_Point3.getElem( 0 )) );
    printf("Point3::operator []( 0 ): %f\n", getfloat((vec_float4)a_Point3[0]) );
    printf("Point3::getX(): %f\n", getfloat(a_Point3.getX( )) );
    printf("Point3::get( 1 ): %f\n", getfloat(a_Point3.getElem( 1 )) );
    printf("Point3::operator []( 1 ): %f\n", getfloat((vec_float4)a_Point3[1]) );
    printf("Point3::getY(): %f\n", getfloat(a_Point3.getY( )) );
    printf("Point3::get( 2 ): %f\n", getfloat(a_Point3.getElem( 2 )) );
    printf("Point3::operator []( 2 ): %f\n", getfloat((vec_float4)a_Point3[2]) );
    printf("Point3::getZ(): %f\n", getfloat(a_Point3.getZ( )) );
    print( ( a_Point3 - b_Point3 ), "Point3 - Point3" );
    print( ( a_Point3 + b_Vector3 ), "Point3 + Vector3" );
    print( ( a_Point3 - b_Vector3 ), "Point3 - Vector3" );
    print( mulPerElem( a_Point3, b_Point3 ), "mulPerElem( Point3, Point3 )" );
    print( divPerElem( a_Point3, b_Point3 ), "divPerElem( Point3, Point3 )" );
    print( recipPerElem( a_Point3 ), "Point3 recip" );
    print( sqrtPerElem( absPerElem( a_Point3 ) ), "Point3 sqrt" );
    print( rsqrtPerElem( absPerElem( a_Point3 ) ), "Point3 rsqrt" );
    print( absPerElem( a_Point3 ), "Point3 abs" );
    print( copySignPerElem( a_Point3, b_Point3 ), "Point3 copySign" );
    print( maxPerElem( a_Point3, b_Point3 ), "Point3 maximum Point3" );
    print( minPerElem( a_Point3, b_Point3 ), "Point3 minimum Point3" );
    printf("Point3 maximum of elements: %f\n", getfloat(maxElem( a_Point3 )));
    printf("Point3 minimum of elements: %f\n", getfloat(minElem( a_Point3 )));
    printf("Point3 sum of elements: %f\n", getfloat(sum( a_Point3 )));
    printf("Point projection: %f\n", getfloat(projection( a_Point3, b_Vector3 )));
    printf("Point distSqrFromOrigin: %f\n", getfloat(distSqrFromOrigin( a_Point3 )) );
    printf("Point distFromOrigin: %f\n", getfloat(distFromOrigin( a_Point3 )) );
    printf("Point distSqr: %f\n", getfloat(distSqr( a_Point3, b_Point3 )) );
    printf("Point dist: %f\n", getfloat(dist( a_Point3, b_Point3 )) );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    c_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    d_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    e_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    a_Point3 = lerp( randfloat(), b_Point3, c_Point3 );
    print( a_Point3, "Point3 lerp" );
}

void
Quat_methods_test()
{
    Vector3 a_Vector3, b_Vector3, c_Vector3, d_Vector3, e_Vector3;
    Vector4 a_Vector4, b_Vector4, c_Vector4, d_Vector4, e_Vector4;
    Point3 a_Point3, b_Point3, c_Point3, d_Point3, e_Point3;
    Quat a_Quat, b_Quat, c_Quat, d_Quat, e_Quat;
    Aos::Quat aos_Quat_0, aos_Quat_1, aos_Quat_2, aos_Quat_3;
    Quat soa_Quat;
    vec_float4 rndflt1, rndflt2, rndflt3, rndflt4;
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Vector3 = Vector3( rndflt1, rndflt2, rndflt3 );
    print( a_Vector3, "set Vector3 with floats" );
    print( b_Vector3, "set Vector3 with floats" );
    c_Vector3 = Vector3( (vec_float4){0.0f} );
    d_Vector3 = Vector3( (vec_float4){0.0f} );
    e_Vector3 = Vector3( (vec_float4){0.0f} );
    print( c_Vector3, "set Vector3 elements to zero" );
    print( d_Vector3, "set Vector3 elements to zero" );
    print( e_Vector3, "set Vector3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Vector4 = Vector4( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Vector4, "set Vector4 with floats" );
    print( b_Vector4, "set Vector4 with floats" );
    c_Vector4 = Vector4( (vec_float4){0.0f} );
    d_Vector4 = Vector4( (vec_float4){0.0f} );
    e_Vector4 = Vector4( (vec_float4){0.0f} );
    print( c_Vector4, "set Vector4 elements to zero" );
    print( d_Vector4, "set Vector4 elements to zero" );
    print( e_Vector4, "set Vector4 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    a_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    b_Point3 = Point3( rndflt1, rndflt2, rndflt3 );
    print( a_Point3, "set Point3 with floats" );
    print( b_Point3, "set Point3 with floats" );
    c_Point3 = Point3( (vec_float4){0.0f} );
    d_Point3 = Point3( (vec_float4){0.0f} );
    e_Point3 = Point3( (vec_float4){0.0f} );
    print( c_Point3, "set Point3 elements to zero" );
    print( d_Point3, "set Point3 elements to zero" );
    print( e_Point3, "set Point3 elements to zero" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    a_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    print( a_Quat, "set Quat with floats" );
    print( b_Quat, "set Quat with floats" );
    c_Quat = Quat( (vec_float4){0.0f} );
    d_Quat = Quat( (vec_float4){0.0f} );
    e_Quat = Quat( (vec_float4){0.0f} );
    print( c_Quat, "set Quat elements to zero" );
    print( d_Quat, "set Quat elements to zero" );
    print( e_Quat, "set Quat elements to zero" );
    a_Quat = Quat( a_Vector3, randfloat() );
    print( a_Quat, "set Quat with Vector3, float" );
    a_Quat = Quat( a_Vector4 );
    print( a_Quat, "construct Quat with Vector4" );
    a_Quat = Quat( randfloat() );
    print( a_Quat, "set Quat with float" );
    a_Quat = Quat( randfloat() );
    print( a_Quat, "set Quat with float" );
    aos_Quat_0 = Aos::Quat( 0.0f, 1.0f, 2.0f, 3.0f );
    aos_Quat_1 = Aos::Quat( 4.0f, 5.0f, 6.0f, 7.0f );
    aos_Quat_2 = Aos::Quat( 8.0f, 9.0f, 10.0f, 11.0f );
    aos_Quat_3 = Aos::Quat( 12.0f, 13.0f, 14.0f, 15.0f );
    soa_Quat = Quat( aos_Quat_0, aos_Quat_1, aos_Quat_2, aos_Quat_3 );
    soa_Quat.get4Aos( aos_Quat_3, aos_Quat_2, aos_Quat_1, aos_Quat_0 );
    Aos::print( aos_Quat_0, "aos type 0" );
    Aos::print( aos_Quat_1, "aos type 1" );
    Aos::print( aos_Quat_2, "aos type 2" );
    Aos::print( aos_Quat_3, "aos type 3" );
    a_Quat = select( a_Quat, b_Quat, ((vec_uint4){0,0xffffffff,0,0xffffffff}) );
    a_Quat.get4Aos( aos_Quat_0, aos_Quat_1, aos_Quat_2, aos_Quat_3 );
    Aos::print( aos_Quat_0, "select 0" );
    Aos::print( aos_Quat_1, "select 1" );
    Aos::print( aos_Quat_2, "select 2" );
    Aos::print( aos_Quat_3, "select 3" );
    a_Quat = b_Quat;
    print( a_Quat, "assign to Quat from Quat" );
    a_Quat.setXYZ( a_Vector3 );
    print( a_Quat, "set Quat xyz" );
    print( a_Quat.getXYZ( ), "get Quat xyz" );
    a_Quat = Quat( (vec_float4){0.0f} );
    print( a_Quat, "set Quat elements to zero" );
    a_Quat.setElem( 0, randfloat() );
    print( a_Quat, "Quat::set( 0, float )" );
    a_Quat[0] = randfloat();
    a_Quat[0] = vec_mul_float( a_Quat[0], randfloat() );
    a_Quat[0] = divf4( a_Quat[0], randfloat() );
    a_Quat[0] = vec_add_float( a_Quat[0], randfloat() );
    a_Quat[0] = vec_sub_float( a_Quat[0], randfloat() );
    print( a_Quat, "Quat::operator [](0)" );
    a_Quat.setX( randfloat() );
    print( a_Quat, "Quat::setX()" );
    a_Quat.setElem( 1, randfloat() );
    print( a_Quat, "Quat::set( 1, float )" );
    a_Quat[1] = randfloat();
    a_Quat[1] = vec_mul_float( a_Quat[1], randfloat() );
    a_Quat[1] = divf4( a_Quat[1], randfloat() );
    a_Quat[1] = vec_add_float( a_Quat[1], randfloat() );
    a_Quat[1] = vec_sub_float( a_Quat[1], randfloat() );
    print( a_Quat, "Quat::operator [](1)" );
    a_Quat.setY( randfloat() );
    print( a_Quat, "Quat::setY()" );
    a_Quat.setElem( 2, randfloat() );
    print( a_Quat, "Quat::set( 2, float )" );
    a_Quat[2] = randfloat();
    a_Quat[2] = vec_mul_float( a_Quat[2], randfloat() );
    a_Quat[2] = divf4( a_Quat[2], randfloat() );
    a_Quat[2] = vec_add_float( a_Quat[2], randfloat() );
    a_Quat[2] = vec_sub_float( a_Quat[2], randfloat() );
    print( a_Quat, "Quat::operator [](2)" );
    a_Quat.setZ( randfloat() );
    print( a_Quat, "Quat::setZ()" );
    a_Quat.setElem( 3, randfloat() );
    print( a_Quat, "Quat::set( 3, float )" );
    a_Quat[3] = randfloat();
    a_Quat[3] = vec_mul_float( a_Quat[3], randfloat() );
    a_Quat[3] = divf4( a_Quat[3], randfloat() );
    a_Quat[3] = vec_add_float( a_Quat[3], randfloat() );
    a_Quat[3] = vec_sub_float( a_Quat[3], randfloat() );
    print( a_Quat, "Quat::operator [](3)" );
    a_Quat.setW( randfloat() );
    print( a_Quat, "Quat::setW()" );
    printf("Quat::get( 0 ): %f\n", getfloat(a_Quat.getElem( 0 )) );
    printf("Quat::operator []( 0 ): %f\n", getfloat((vec_float4)a_Quat[0]) );
    printf("Quat::getX(): %f\n", getfloat(a_Quat.getX( )) );
    printf("Quat::get( 1 ): %f\n", getfloat(a_Quat.getElem( 1 )) );
    printf("Quat::operator []( 1 ): %f\n", getfloat((vec_float4)a_Quat[1]) );
    printf("Quat::getY(): %f\n", getfloat(a_Quat.getY( )) );
    printf("Quat::get( 2 ): %f\n", getfloat(a_Quat.getElem( 2 )) );
    printf("Quat::operator []( 2 ): %f\n", getfloat((vec_float4)a_Quat[2]) );
    printf("Quat::getZ(): %f\n", getfloat(a_Quat.getZ( )) );
    printf("Quat::get( 3 ): %f\n", getfloat(a_Quat.getElem( 3 )) );
    printf("Quat::operator []( 3 ): %f\n", getfloat((vec_float4)a_Quat[3]) );
    printf("Quat::getW(): %f\n", getfloat(a_Quat.getW( )) );
    print( ( a_Quat + b_Quat ), "Quat + Quat" );
    print( ( a_Quat - b_Quat ), "Quat - Quat" );
    print( ( a_Quat * b_Quat ), "Quat * Quat" );
    print( ( a_Quat * randfloat() ), "Quat * float" );
    print( ( a_Quat / randfloat() ), "Quat / float" );
    print( ( randfloat() * a_Quat ), "float * Quat" );
    print( ( -a_Quat ), "Quat negate" );
    printf("Quat dot Quat: %f\n", getfloat(dot( a_Quat, b_Quat )));
    printf("Quat lengthSqr: %f\n", getfloat(norm( a_Quat )));
    printf("Quat length: %f\n", getfloat(length( a_Quat )));
    print( normalize( a_Quat ), "Quat normalized" );
    a_Quat = Quat::identity( );
    print( a_Quat, "set to identity" );
    a_Quat = Quat::rotation( a_Vector3, b_Vector3 );
    print( a_Quat, "Quat rotation between vectors" );
    a_Quat = Quat::rotation( randfloat(), a_Vector3 );
    print( a_Quat, "Quat rotation axis angle" );
    a_Quat = Quat::rotationX( randfloat() );
    print( a_Quat, "Quat rotationX" );
    a_Quat = Quat::rotationY( randfloat() );
    print( a_Quat, "Quat rotationY" );
    a_Quat = Quat::rotationZ( randfloat() );
    print( a_Quat, "Quat rotationZ" );
    print( rotate( a_Quat, a_Vector3 ), "Quat rotate Vector3" );
    print( conj( a_Quat ), "Quat conj" );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    b_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    c_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    d_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    rndflt1 = randfloat();
    rndflt2 = randfloat();
    rndflt3 = randfloat();
    rndflt4 = randfloat();
    e_Quat = Quat( rndflt1, rndflt2, rndflt3, rndflt4 );
    b_Quat = normalize( b_Quat );
    c_Quat = normalize( c_Quat );
    d_Quat = normalize( d_Quat );
    e_Quat = normalize( e_Quat );
    a_Quat = lerp( randfloat(), b_Quat, c_Quat );
    print( a_Quat, "Quat lerp" );
    a_Quat = slerp( randfloat(), b_Quat, c_Quat );
    print( a_Quat, "Quat slerp" );
    a_Quat = squad( randfloat(), b_Quat, c_Quat, d_Quat, e_Quat );
    print( a_Quat, "Quat squad" );
}

int main()
{
    printf("\n __begin__ \n");
    for ( iteration = 0; iteration < 2; iteration++ ) {
        Vector3_methods_test();
        Vector4_methods_test();
        Point3_methods_test();
        Quat_methods_test();
    }
    printf("\n __end__ \n");
    return 0;
}