/*
 * Copyright (c) 2013, Linz Center of Mechatronics GmbH (LCM), web: www.lcm.at
 * All rights reserved.
 */
/*
 * This file is licensed according to the BSD 3-clause license as follows:
 * 
 * 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 "Linz Center of Mechatronics GmbH" and "LCM" 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 "Linz Center of Mechatronics GmbH" 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.
 */
/*
 * This file is part of X2C. web: x2c.lcm.at
 * $LastChangedRevision: 4111 $
 * $LastChangedDate:: 2026-03-09 14:42:15 +0100#$
 */
/* USERCODE-BEGIN:Description                                                                                         */
/*      Description: With lookup table approximated arc tangent function      */
/*                      Calculation:                                          */
/*                          theta = atan2(y,x)                                */
/*                      Range:                                                */
/*                          y:      -1..1                                     */
/*                          x:      -1..1                                     */
/*                          theta:  -pi..pi                                   */
/*                                      --> with -pi..pi ~ -1..1              */
/*                                                                            */
/* USERCODE-END:Description                                                                                           */
#include <stddef.h>
#include "Atan2_FiP32.h"

/* all used update functions to ram for c2000                                                                         */
#if defined(__ALL_UPDATE_FUNC_2_RAM_C2000__)
    #pragma CODE_SECTION(Atan2_FiP32_Update, "ramfuncs")
#endif

/* USERCODE-BEGIN:PreProcessor                                                                                        */
#include "Atan_Data.h"

/* Inputs */
#define IN1                 (*block->y)
#define IN2                 (*block->x)

/* Outputs */
#define THETA               (block->Out)

/* Quadrant Offset Macros */
#define ATAN2_OFFSET_Q2     0x40000000
#define ATAN2_OFFSET_Q34    0x80000000u
#define ATAN2_OFFSET_COT    0x3FFFFFFF

#define QUARTER_PI          0x20000000      /* 45� */

/* USERCODE-END:PreProcessor                                                                                          */

/**********************************************************************************************************************/
/** Update                                                                                                           **/
/**********************************************************************************************************************/
void Atan2_FiP32_Update(ATAN2_FIP32 *block)
{
/* USERCODE-BEGIN:UpdateFnc                                                                                           */
    int32   iFact, iY, iX;
    int64   arg  = (int64)0;
    int8    iTan =  (int8)1;        /* 1 -> tangent for Q1 & Q3, 0 -> cotangent for Q2 & Q4 */
    int8    Q13  =  (int8)1;        /* 1 -> Q1 or Q3 */
    int32   out1 = (int32)0;

    iY = IN1;
    iX = IN2;

    /* check input range (negating lowest negative number could lead to division by zero) */
    if (iX == (int32)0x80000000u)
    {
        iX = (int32)0x80000001u;
    }
    if (iY == (int32)0x80000000u)
    {
        iY = (int32)0x80000001u;
    }

    if ((iY == 0) && (iX == 0))
    {
        out1 = 0;
    }
    else
    {
        /* Check quadrants */
        if (iY < (int32)0)
        {
            iY   = -iY;
            out1 = (int32)ATAN2_OFFSET_Q34;
            if (iX > (int32)0)  /* 4th quadrant */
            {
                out1 += ATAN2_OFFSET_Q2;
                Q13   = (int8)0;
            }
            else
            {           /* 3rd quadrant */
                iX   = -iX;
            }
        }
        else
        {
            if (iX < (int32)0)  /* 2nd quadrant */
            {
                iX   = -iX;
                out1 = ATAN2_OFFSET_Q2;
                Q13  = (int8)0;
            }
        }

        /* Check region in quadrant */
        if (iY > iX)    /* angle > 45� */
        {
            iFact = iX;         /* swap input arguments  (iFact is used as temporary variable) */
            iX    = iY;         /* swap input arguments */
            iY    = iFact;      /* swap input arguments */
            iTan  = (int8)0;    /* use cotangent */
        }

        if (iY != iX)
        {
            /* Calculation of y/x */
            /* RULECHECKER_comment(1:0, 1:0, check_inappropriate_int, "see X2C MISRA Deviations, DEV0004", true_no_defect) */
            arg  = (int64)iY << 32; /* Q63,64 */
            arg /= iX;              /* Q32,64 */

            /* Interpolation Previous: LOOKUP32(Atan_Table32, (uint32)arg, iFact); */
            iFact = lookupTableInt32(Atan_Table32, (int32)((uint32)arg));
        }
        else
        {
            iFact = QUARTER_PI;     /* angle = 45� */
        }



        /* Check for need of cotangent correction */
        if (Q13 != iTan)
        {
            out1 += ((int32)ATAN2_OFFSET_COT - iFact);
        }
        else
        {
            out1 += iFact;
        }
    }

    THETA = out1;           /* assign output */

/* USERCODE-END:UpdateFnc                                                                                             */

}

/**********************************************************************************************************************/
/** Initialization                                                                                                   **/
/**********************************************************************************************************************/
void Atan2_FiP32_Init(ATAN2_FIP32 *block)
{
    block->identifier = ATAN2_FIP32_ID;
    block->Out = 0;
/* USERCODE-BEGIN:InitFnc                                                                                             */
/* USERCODE-END:InitFnc                                                                                               */
}

/**********************************************************************************************************************/
/** No Load function                                                                                                 **/
/**********************************************************************************************************************/

/**********************************************************************************************************************/
/** No Save function                                                                                                 **/
/**********************************************************************************************************************/

/**********************************************************************************************************************/
/** Get block element address                                                                                        **/
/**********************************************************************************************************************/
#if !defined(ATAN2_FIP32_ISLINKED)
void* Atan2_FiP32_GetAddress(ATAN2_FIP32* block, uint16 elementId)
{
    void* addr;
    switch (elementId)
    {
        case 1:
            addr = block->y;
            break;
        case 2:
            addr = block->x;
            break;
        case 3:
            addr = &block->Out;
            break;
        default:
            addr = NULL;
            break;
    }
    return (addr);
}
#endif
