MercuryMath.cpp

Go to the documentation of this file.

#include "MercuryMath.h"

#if !defined( _MSC_VER ) || defined( SLOW_CODE )

//Generic Math functions. Compile these if you can not use optimized functions.

void Mul4f(const float* first, const float* second, float* out)
{
    out[0] = first[0] * second[0];
    out[1] = first[1] * second[1];
    out[2] = first[2] * second[2];
    out[3] = first[3] * second[3];
}

void Div4f(const float* first, const float* second, float* out)
{
    out[0] = first[0] / second[0];
    out[1] = first[1] / second[1];
    out[2] = first[2] / second[2];
    out[3] = first[3] / second[3];
}

void Add4f(const float* first, const float* second, float* out)
{
    out[0] = first[0] + second[0];
    out[1] = first[1] + second[1];
    out[2] = first[2] + second[2];
    out[3] = first[3] + second[3];
}

void Sub4f(const float* first, const float* second, float* out)
{
    out[0] = first[0] - second[0];
    out[1] = first[1] - second[1];
    out[2] = first[2] - second[2];
    out[3] = first[3] - second[3];
}

void Copy4f( void * dest, const void * source )
{
    ((float*)dest)[0] = ((float*)source)[0];
    ((float*)dest)[1] = ((float*)source)[1];
    ((float*)dest)[2] = ((float*)source)[2];
    ((float*)dest)[3] = ((float*)source)[3];
}

void Copy8f( void * dest, const void * source )
{
    ((float*)dest)[0] = ((float*)source)[0];
    ((float*)dest)[1] = ((float*)source)[1];
    ((float*)dest)[2] = ((float*)source)[2];
    ((float*)dest)[3] = ((float*)source)[3];

    ((float*)dest)[4] = ((float*)source)[4];
    ((float*)dest)[5] = ((float*)source)[5];
    ((float*)dest)[6] = ((float*)source)[6];
    ((float*)dest)[7] = ((float*)source)[7];
}

void Copy16f( void * dest, const void * source )
{
    ((float*)dest)[0] = ((float*)source)[0];
    ((float*)dest)[1] = ((float*)source)[1];
    ((float*)dest)[2] = ((float*)source)[2];
    ((float*)dest)[3] = ((float*)source)[3];

    ((float*)dest)[4] = ((float*)source)[4];
    ((float*)dest)[5] = ((float*)source)[5];
    ((float*)dest)[6] = ((float*)source)[6];
    ((float*)dest)[7] = ((float*)source)[7];

    ((float*)dest)[8] = ((float*)source)[8];
    ((float*)dest)[9] = ((float*)source)[9];
    ((float*)dest)[10] = ((float*)source)[10];
    ((float*)dest)[11] = ((float*)source)[11];

    ((float*)dest)[12] = ((float*)source)[12];
    ((float*)dest)[13] = ((float*)source)[13];
    ((float*)dest)[14] = ((float*)source)[14];
    ((float*)dest)[15] = ((float*)source)[15];
}

void R_ConcatTransforms4 ( const float* in1, const float* in2, float* out)
{
    out[0] = in1[0] * in2[0] + in1[1] * in2[4] +
                in1[2] * in2[8] + in1[3] * in2[12];
    out[1] = in1[0] * in2[1] + in1[1] * in2[5] +
                in1[2] * in2[9] + in1[3] * in2[13];
    out[2] = in1[0] * in2[2] + in1[1] * in2[6] +
                in1[2] * in2[10] + in1[3] * in2[14];
    out[3] = in1[0] * in2[3] + in1[1] * in2[7] +
                in1[2] * in2[11] + in1[3] * in2[15];

    out[4] = in1[4] * in2[0] + in1[5] * in2[4] +
                in1[6] * in2[8] + in1[7] * in2[12];
    out[5] = in1[4] * in2[1] + in1[5] * in2[5] +
                in1[6] * in2[9] + in1[7] * in2[13];
    out[6] = in1[4] * in2[2] + in1[5] * in2[6] +
                in1[6] * in2[10] + in1[7] * in2[14];
    out[7] = in1[4] * in2[3] + in1[5] * in2[7] +
                in1[6] * in2[11] + in1[7] * in2[15];

    out[8] = in1[8] * in2[0] + in1[9] * in2[4] +
                in1[10] * in2[8] + in1[11] * in2[12];
    out[9] = in1[8] * in2[1] + in1[9] * in2[5] +
                in1[10] * in2[9] + in1[11] * in2[13];
    out[10] = in1[8] * in2[2] + in1[9] * in2[6] +
                in1[10] * in2[10] + in1[11] * in2[14];
    out[11] = in1[8] * in2[3] + in1[9] * in2[7] +
                in1[10] * in2[11] + in1[11] * in2[15];

    out[12] = in1[12] * in2[0] + in1[13] * in2[4] +
                in1[14] * in2[8] + in1[15] * in2[12];
    out[13] = in1[12] * in2[1] + in1[13] * in2[5] +
                in1[14] * in2[9] + in1[15] * in2[13];
    out[14] = in1[12] * in2[2] + in1[13] * in2[6] +
                in1[14] * in2[10] + in1[15] * in2[14];
    out[15] = in1[12] * in2[3] + in1[13] * in2[7] +
                in1[14] * in2[11] + in1[15] * in2[15];
}

#endif

void VectorMultiply4f( const float * m, float *p, float *out )
{
    out[0] = p[0] * m[0] + p[1] * m[1] + p[2] * m[2] + p[3] * m[3];
    out[1] = p[0] * m[4] + p[1] * m[5] + p[2] * m[6] + p[3] * m[7];
    out[2] = p[0] * m[8] + p[1] * m[9] + p[2] * m[10] + p[3] * m[11];
    out[3] = p[0] * m[12] + p[1] * m[13] + p[2] * m[14] + p[3] * m[15];
}

/*
 * (c) 2006 Joshua Allen, Charles Lohr
 * All rights reserved.
 * 
 * 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, and/or sell copies of the Software, and to permit persons to
 * whom the Software is furnished to do so, provided that the above
 * copyright notice(s) and this permission notice appear in all copies of
 * the Software and that both the above copyright notice(s) and this
 * permission notice appear in supporting documentation.
 * 
 * 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 OF
 * THIRD PARTY RIGHTS. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR HOLDERS
 * INCLUDED IN THIS NOTICE BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL INDIRECT
 * OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS
 * OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR
 * OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
 * PERFORMANCE OF THIS SOFTWARE.
 */

---

Hosted by