Home

Resume

Blog

Teikitu


/* =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-==-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= */
/*  »Project«   Teikitu Gaming System (TgS) (∂)
    »File«      TgS Common - Math API [Vector].c
    »Author«    Andrew Aye (EMail: mailto:andrew.aye@gmail.com, Web: http://www.andrewaye.com)
    »Version«   4.51 / »GUID« A9981407-3EC9-42AF-8B6F-8BE6DD919615                                                                                                        */
/*   -------------------------------------------------------------------------------------------------------------------------------------------------------------------- */
/*  Copyright: © 2002-2017, Andrew Aye.  All Rights Reserved.
    This software is free for non-commercial use.  Redistribution and use in source and binary forms, with or without modification, are permitted provided that the
      following conditions are met:
        Redistribution of source code must retain this copyright notice, this list of conditions and the following disclaimers.
        Redistribution in binary form must reproduce this copyright notice, this list of conditions and the following disclaimers in the documentation and other materials
          provided with the distribution.
    The name of the author may not be used to endorse or promote products derived from this software without specific prior written permission.
    The intellectual property rights of the algorithms used reside with Andrew Aye.
    You may not use this software, in whole or in part, in support of any commercial product without the express written consent of the author.
    There is no warranty or other guarantee of fitness of this software for any purpose. It is provided solely "as is".                                                   */
/* =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-==-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= */
/* == Common ============================================================================================================================================================ */

/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */
/*  File Local Constants - Random Number Generator - Mersenne Twister */
/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */

enum { KTgMERSENNE_TWISTER_PERIOD_N         = 624 };




/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */
/*  File Local Functions And Data - Random Number Generator - Mersenne Twister */
/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */

TgEXTN TgTLS TgUINT32                       tls_uiMT_Index;
TgEXTN TgTLS union
{
    TgUINT32                                    m_uiData[KTgMERSENNE_TWISTER_PERIOD_N];
    TgVEC_U32_04                                m_vData[KTgMERSENNE_TWISTER_PERIOD_N >> 2];
}                                           tls_uiMT_State;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
TgEXTN STg1_MP_CS                           g_csMT_Lock;
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

TgVOID                                      tgCM_Random_MT_Reload( TgVOID );




/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */
/*  Public Functions - Random Number Generator - Mersenne Twister */
/* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. */

/* ---- F_RAND_F32_04 --------------------------------------------------------------------------------------------------------------------------------------------------- */
/* ---------------------------------------------------------------------------------------------------------------------------------------------------------------------- */
TgVEC_F32_04 F_RAND_F32_04( TgVOID )
{
    TgUINT32                                ui0_0, ui1_0, ui2_0, ui3_0;
    TgVEC_F32_04                            tvRand_Val;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Enter_Block( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    if (tls_uiMT_Index + 4 > KTgMERSENNE_TWISTER_PERIOD_N)
    {
        tgCM_Random_MT_Reload();
    }

    ui0_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 0];
    ui1_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 1];
    ui2_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 2];
    ui3_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 3];

    tls_uiMT_Index += 4;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Exit( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    {
        C_TgUINT32                          ui0_1 = ui0_0 ^ (ui0_0 >> 11);
        C_TgUINT32                          ui1_1 = ui1_0 ^ (ui1_0 >> 11);
        C_TgUINT32                          ui2_1 = ui2_0 ^ (ui2_0 >> 11);
        C_TgUINT32                          ui3_1 = ui3_0 ^ (ui3_0 >> 11);

        C_TgUINT32                          ui0_2 = ui0_1 ^ ((ui0_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui1_2 = ui1_1 ^ ((ui1_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui2_2 = ui2_1 ^ ((ui2_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui3_2 = ui3_1 ^ ((ui3_1 << 7) & 0x9d2c5680UL);

        C_TgUINT32                          ui0_3 = ui0_2 ^ ((ui0_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui1_3 = ui1_2 ^ ((ui1_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui2_3 = ui2_2 ^ ((ui2_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui3_3 = ui3_2 ^ ((ui3_2 << 15) & 0xefc60000UL);

        C_TgUINT32                          ui0_4 = ui0_3 ^ (ui0_3 >> 18);
        C_TgUINT32                          ui1_4 = ui1_3 ^ (ui1_3 >> 18);
        C_TgUINT32                          ui2_4 = ui2_3 ^ (ui2_3 >> 18);
        C_TgUINT32                          ui3_4 = ui3_3 ^ (ui3_3 >> 18);

        tvRand_Val.m.x = ui0_4 / 4294967296.0F;
        tvRand_Val.m.y = ui1_4 / 4294967296.0F;
        tvRand_Val.m.z = ui2_4 / 4294967296.0F;
        tvRand_Val.m.w = ui3_4 / 4294967296.0F;
    }

    return (tvRand_Val);
}


/* ---- F_RAND_F64_04 --------------------------------------------------------------------------------------------------------------------------------------------------- */
/* ---------------------------------------------------------------------------------------------------------------------------------------------------------------------- */
TgVEC_F64_04 F_RAND_F64_04( TgVOID )
{
    TgUINT32                                ui0_0, ui1_0, ui2_0, ui3_0;
    TgVEC_F64_04                            tvRand_Val;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Enter_Block( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    if (tls_uiMT_Index + 4 > KTgMERSENNE_TWISTER_PERIOD_N)
    {
        tgCM_Random_MT_Reload();
    }

    ui0_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 0];
    ui1_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 1];
    ui2_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 2];
    ui3_0 = tls_uiMT_State.m_uiData[tls_uiMT_Index + 3];

    tls_uiMT_Index += 4;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Exit( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    {
        C_TgUINT32                          ui0_1 = ui0_0 ^ (ui0_0 >> 11);
        C_TgUINT32                          ui1_1 = ui1_0 ^ (ui1_0 >> 11);
        C_TgUINT32                          ui2_1 = ui2_0 ^ (ui2_0 >> 11);
        C_TgUINT32                          ui3_1 = ui3_0 ^ (ui3_0 >> 11);

        C_TgUINT32                          ui0_2 = ui0_1 ^ ((ui0_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui1_2 = ui1_1 ^ ((ui1_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui2_2 = ui2_1 ^ ((ui2_1 << 7) & 0x9d2c5680UL);
        C_TgUINT32                          ui3_2 = ui3_1 ^ ((ui3_1 << 7) & 0x9d2c5680UL);

        C_TgUINT32                          ui0_3 = ui0_2 ^ ((ui0_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui1_3 = ui1_2 ^ ((ui1_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui2_3 = ui2_2 ^ ((ui2_2 << 15) & 0xefc60000UL);
        C_TgUINT32                          ui3_3 = ui3_2 ^ ((ui3_2 << 15) & 0xefc60000UL);

        C_TgUINT32                          ui0_4 = ui0_3 ^ (ui0_3 >> 18);
        C_TgUINT32                          ui1_4 = ui1_3 ^ (ui1_3 >> 18);
        C_TgUINT32                          ui2_4 = ui2_3 ^ (ui2_3 >> 18);
        C_TgUINT32                          ui3_4 = ui3_3 ^ (ui3_3 >> 18);

        tvRand_Val.m.x = ui0_4 / 4294967296.0;
        tvRand_Val.m.y = ui1_4 / 4294967296.0;
        tvRand_Val.m.z = ui2_4 / 4294967296.0;
        tvRand_Val.m.w = ui3_4 / 4294967296.0;
    }

    return (tvRand_Val);
}


/* ---- M_RAND_F32_04 --------------------------------------------------------------------------------------------------------------------------------------------------- */
/* ---------------------------------------------------------------------------------------------------------------------------------------------------------------------- */
TgVEC_M_F32_04 M_RAND_F32_04( TgVOID )
{
    TgUINT32                                uiInc = 4 + ((4 - (tls_uiMT_Index & 0x3)) & 0x3);
    TgVECTOR_ALIAS_UNION                    tvRand_Val;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Enter_Block( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    if (tls_uiMT_Index + uiInc > KTgMERSENNE_TWISTER_PERIOD_N)
    {
        tgCM_Random_MT_Reload();
    }

    TgERROR(0 == ((tls_uiMT_Index + uiInc - 4) & 0x3));

    tvRand_Val.m_u32_v04 = tls_uiMT_State.m_vData[(tls_uiMT_Index + uiInc - 4) >> 2];

    tls_uiMT_Index += uiInc;

#if !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD
    tgCM_MP_CS_Exit( &g_csMT_Lock );
/*# !TgCOMPILE_THREAD_LOCAL && TgCOMPILE_THREAD */
#endif

    return (tvRand_Val.m_f32_v04.m_mData);
}