Subversion Repositories neo430

// #################################################################################################

// #  < neo430_muldiv.c - Multiplier/Divider function >                                            #

// # ********************************************************************************************* #

// # BSD 3-Clause License                                                                          #

// #                                                                                               #

// # Copyright (c) 2020, Stephan Nolting. All rights reserved.                                     #

// #                                                                                               #

// # Redistribution and use in source and binary forms, with or without modification, are          #

// # permitted provided that the following conditions are met:                                     #

// #                                                                                               #

// # 1. Redistributions of source code must retain the above copyright notice, this list of        #

// #    conditions and the following disclaimer.                                                   #

// #                                                                                               #

// # 2. 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.                                                            #

// #                                                                                               #

// # 3. Neither the name of the copyright holder 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 HOLDER 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.                                                            #

// # ********************************************************************************************* #

// # The NEO430 Processor - https://github.com/stnolting/neo430                                    #

// #################################################################################################

#include "neo430.h"

#include "neo430_muldiv.h"

// Macros

#define muldiv_processing_delay {asm volatile("nop"); asm volatile("nop"); asm volatile("nop");}

/* ------------------------------------------------------------

 * INFO Unsigned 16x16-bit multiplication

 * PARAM 16-bit factor a

 * PARAM 16-bit factor b

 * RETURN 32-bit product

 * ------------------------------------------------------------ */

uint32_t neo430_umul32(uint16_t a, uint16_t b) {

  MULDIV_OPA_RESX = a;

  MULDIV_OPB_UMUL_RESY = b;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  return MULDIV_R32bit;

}

/* ------------------------------------------------------------

 * INFO Unsigned 32x32-bit multiplication (large!)

 * PARAM 32-bit factor a

 * PARAM 32-bit factor b

 * RETURN 32-bit product

 * ------------------------------------------------------------ */

uint32_t neo430_umul32_32(uint32_t a, uint32_t b) {

  union uint32_u tmp32;

  // get input words (16-bit)

  tmp32.uint32 = a;

  register uint16_t a_lo = tmp32.uint16[0];

  register uint16_t a_hi = tmp32.uint16[1];

  tmp32.uint32 = b;

  register uint16_t b_lo = tmp32.uint16[0];

  register uint16_t b_hi = tmp32.uint16[1];

  // compute partial results and arrange in 64-bit word

  MULDIV_OPA_RESX      = a_lo;

  MULDIV_OPB_UMUL_RESY = b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp32.uint16[0] = MULDIV_OPA_RESX;

  tmp32.uint16[1] = MULDIV_OPB_UMUL_RESY;

  uint32_t res32 = tmp32.uint32;

//MULDIV_OPA_RESX      = a_lo;

  MULDIV_OPB_UMUL_RESY = b_hi;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp32.uint16[0] = 0;

  tmp32.uint16[1] = MULDIV_OPA_RESX;

  res32 += tmp32.uint32;

  MULDIV_OPA_RESX      = a_hi;

  MULDIV_OPB_UMUL_RESY = b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp32.uint16[0] = 0;

  tmp32.uint16[1] = MULDIV_OPA_RESX;

  res32 += tmp32.uint32;

  return res32;

}

/* ------------------------------------------------------------

 * INFO Unsigned 32x32-bit multiplication (large!)

 * PARAM 32-bit factor a

 * PARAM 32-bit factor b

 * RETURN 64-bit product

 * ------------------------------------------------------------ */

uint64_t neo430_umul64(uint32_t a, uint32_t b) {

  union uint32_u tmp32;

  union uint64_u tmp64;

  // get input words (16-bit)

  tmp32.uint32 = a;

  register uint16_t a_lo = tmp32.uint16[0];

  register uint16_t a_hi = tmp32.uint16[1];

  tmp32.uint32 = b;

  register uint16_t b_lo = tmp32.uint16[0];

  register uint16_t b_hi = tmp32.uint16[1];

  // compute partial results and arrange in 64-bit word

  MULDIV_OPA_RESX      = a_lo;

  MULDIV_OPB_UMUL_RESY = b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.uint16[0] = MULDIV_OPA_RESX;

  tmp64.uint16[1] = MULDIV_OPB_UMUL_RESY;

  tmp64.uint16[2] = 0;

  tmp64.uint16[3] = 0;

  uint64_t res64 = tmp64.uint64;

//MULDIV_OPA_RESX      = a_lo;

  MULDIV_OPB_UMUL_RESY = b_hi;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.uint16[0] = 0;

  tmp64.uint16[1] = MULDIV_OPA_RESX;

  tmp64.uint16[2] = MULDIV_OPB_UMUL_RESY;

  tmp64.uint16[3] = 0;

  res64 += tmp64.uint64;

  MULDIV_OPA_RESX      = a_hi;

  MULDIV_OPB_UMUL_RESY = b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.uint16[0] = 0;

  tmp64.uint16[1] = MULDIV_OPA_RESX;

  tmp64.uint16[2] = MULDIV_OPB_UMUL_RESY;

  tmp64.uint16[3] = 0;

  res64 += tmp64.uint64;

//MULDIV_OPA_RESX      = a_hi;

  MULDIV_OPB_UMUL_RESY = b_hi;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.uint16[0] = 0;

  tmp64.uint16[1] = 0;

  tmp64.uint16[2] = MULDIV_OPA_RESX;

  tmp64.uint16[3] = MULDIV_OPB_UMUL_RESY;

  res64 += tmp64.uint64;

  return res64;

}

/* ------------------------------------------------------------

 * INFO Signed 16x16-bit multiplication

 * PARAM 16-bit factor a

 * PARAM 16-bit factor b

 * RETURN 32-bit product

 * ------------------------------------------------------------ */

int32_t neo430_mul32(int16_t a, int16_t b) {

  MULDIV_OPA_RESX = (uint16_t)a;

  MULDIV_OPB_SMUL = (uint16_t)b;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  int32_t r = (int32_t)MULDIV_R32bit;

  return r;

}

/* ------------------------------------------------------------

 * INFO Signed 32x32-bit multiplication (large!)

 * PARAM 32-bit factor a

 * PARAM 32-bit factor b

 * RETURN 64-bit product

 * ------------------------------------------------------------ */

int64_t neo430_mul64(int32_t a, int32_t b) {

  union int32_u tmp32;

  union int64_u tmp64;

  // get input words (16-bit)

  tmp32.int32 = a;

  register int16_t a_lo = tmp32.int16[0];

  register int16_t a_hi = tmp32.int16[1];

  tmp32.int32 = b;

  register int16_t b_lo = tmp32.int16[0];

  register int16_t b_hi = tmp32.int16[1];

  // compute partial results and arrange in 64-bit word

  MULDIV_OPA_RESX = (uint16_t)a_lo;

  MULDIV_OPB_SMUL = (uint16_t)b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.int16[0] = MULDIV_OPA_RESX;

  tmp64.int16[1] = MULDIV_OPB_UMUL_RESY;

  tmp64.int16[2] = 0;

  tmp64.int16[3] = 0;

  int64_t res64 = tmp64.int64;

//MULDIV_OPA_RESX = (uint16_t)a_lo;

  MULDIV_OPB_SMUL = (uint16_t)b_hi;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.int16[0] = 0;

  tmp64.int16[1] = MULDIV_OPA_RESX;

  tmp64.int16[2] = MULDIV_OPB_UMUL_RESY;

  if (tmp64.int16[2] < 0) // sign extension

    tmp64.int16[3] = 0xffff;

  else

    tmp64.int16[3] = 0;

  res64 += tmp64.int64;

  MULDIV_OPA_RESX = (uint16_t)a_hi;

  MULDIV_OPB_SMUL = (uint16_t)b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.int16[0] = 0;

  tmp64.int16[1] = MULDIV_OPA_RESX;

  tmp64.int16[2] = MULDIV_OPB_UMUL_RESY;

  if (tmp64.int16[2] < 0) // sign extension

    tmp64.int16[3] = 0xffff;

  else

    tmp64.int16[3] = 0;

  res64 += tmp64.int64;

//MULDIV_OPA_RESX = (uint16_t)a_hi;

  MULDIV_OPB_SMUL = (uint16_t)b_lo;

#if NEO430_HWMUL_DSP

  // no processind delay required for DSP (MUL) blocks

#else

  muldiv_processing_delay // HW processing delay

#endif

  tmp64.int16[0] = 0;

  tmp64.int16[1] = 0;

  tmp64.int16[2] = MULDIV_OPA_RESX;

  tmp64.int16[3] = MULDIV_OPB_UMUL_RESY;

  res64 += tmp64.int64;

  return res64;

}

/* ------------------------------------------------------------

 * INFO Unsigned 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * RETURN 16-bit quotient

 * ------------------------------------------------------------ */

uint16_t neo430_udiv16(uint16_t dividend, uint16_t divisor) {

  MULDIV_OPA_RESX = dividend;

  MULDIV_OPB_UDIV = divisor;

  muldiv_processing_delay // HW processing delay

  return MULDIV_OPA_RESX;

}

/* ------------------------------------------------------------

 * INFO Signed 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * RETURN 16-bit quotient

 * ------------------------------------------------------------ */

int16_t neo430_div16(int16_t dividend, int16_t divisor) {

  // make positive / unsigned

  if (dividend < 0)

    dividend = 0 - dividend;

  if (divisor < 0)

    divisor = 0 - divisor;

  MULDIV_OPA_RESX = (uint16_t)dividend;

  MULDIV_OPB_UDIV = (uint16_t)divisor;

  muldiv_processing_delay // HW processing delay

  int16_t r = (int16_t)MULDIV_OPA_RESX;

  if (dividend < 0)

    return 0 - r;

  else

    return r;

}

/* ------------------------------------------------------------

 * INFO Unsigned 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * RETURN 16-bit remainder

 * ------------------------------------------------------------ */

uint16_t neo430_umod16(uint16_t dividend, uint16_t divisor) {

  MULDIV_OPA_RESX = dividend;

  MULDIV_OPB_UDIV = divisor;

  muldiv_processing_delay // HW processing delay

  return MULDIV_OPB_UMUL_RESY;

}

/* ------------------------------------------------------------

 * INFO Signed 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * RETURN 16-bit remainder

 * ------------------------------------------------------------ */

int16_t neo430_mod16(int16_t dividend, int16_t divisor) {

  int16_t dividend_int = dividend;

  // make positive / unsigned

  if (dividend_int < 0)

    dividend_int = 0 - dividend_int;

  if (divisor < 0)

    divisor = 0 - divisor;

  MULDIV_OPA_RESX = (uint16_t)dividend_int;

  MULDIV_OPB_UDIV = (uint16_t)divisor;

  muldiv_processing_delay // HW processing delay

  int16_t r = (int16_t)MULDIV_OPB_UMUL_RESY;

  if (dividend < 0)

    return 0 - r;

  else

    return r;

}

/* ------------------------------------------------------------

 * INFO Unsigned 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * PARAM Pointer to store 16-bit remainder

 * RETURN 16-bit quotient

 * ------------------------------------------------------------ */

uint16_t neo430_umoddiv16(uint16_t *remainder, uint16_t dividend, uint16_t divisor) {

  MULDIV_OPA_RESX = dividend;

  MULDIV_OPB_UDIV = divisor;

  muldiv_processing_delay // HW processing delay

  *remainder = MULDIV_OPB_UMUL_RESY;

  return MULDIV_OPA_RESX;

}

/* ------------------------------------------------------------

 * INFO Signed 16 by 16-bit division

 * PARAM 16-bit dividend

 * PARAM 16-bit divisor

 * PARAM Pointer to store 16-bit remainder

 * RETURN 16-bit quotient

 * ------------------------------------------------------------ */

int16_t neo430_moddiv16(int16_t *remainder, int16_t dividend, int16_t divisor) {

  int16_t sign = dividend ^ divisor;

  int16_t dividend_int = dividend;

  // make positive / unsigned

  if (dividend_int < 0)

    dividend_int = 0 - dividend_int;

  if (divisor < 0)

    divisor = 0 - divisor;

  MULDIV_OPA_RESX = (uint16_t)dividend_int;

  MULDIV_OPB_UDIV = (uint16_t)divisor;

  muldiv_processing_delay // HW processing delay

  int16_t q = (int16_t)MULDIV_OPA_RESX;

  int16_t r = (int16_t)MULDIV_OPB_UMUL_RESY;

  if (dividend < 0)

    *remainder = 0 - r;

  else

    *remainder = r;

  if (sign < 0)

    return 0 - q;

  else

    return q;

}

// *****************************************************************************

// * DANGER ZONE!!!                                                            *

// * Override primitives for multiplication to use the MULDIV hardware unit    *

// *****************************************************************************

#if NEO430_HWMUL_ABI_OVERRIDE

int16_t __mulhi2(int16_t x, int16_t y) {

  return (int16_t)neo430_umul32((uint16_t)x,(uint16_t)y);

}

int32_t __mulhisi2(int16_t x, int16_t y) {

  return neo430_mul32(x,y);

}

uint32_t __umulhisi2(uint16_t x, uint16_t y) {

  return neo430_umul32(x,y);

}

int32_t __mulsi2(int32_t x, int32_t y) {

  return (int32_t)neo430_umul32_32((uint32_t)x,(uint32_t)y);

}

int32_t __mulsidi2(int32_t x, int32_t y) {

  return neo430_mul64(x,y);

}

uint64_t __umulsidi2(uint32_t x, uint32_t y) {

  return neo430_umul64(x,y);

}

#endif

// *****************************************************************************