lib_umul32x64.c 3.3 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586
  1. /****************************************************************************
  2. * libs/libc/misc/lib_umul32x64.c
  3. *
  4. * Copyright (C) 2016 Gregory Nutt. All rights reserved.
  5. * Author: Gregory Nutt <gnutt@nuttx.org>
  6. *
  7. * Redistribution and use in source and binary forms, with or without
  8. * modification, are permitted provided that the following conditions
  9. * are met:
  10. *
  11. * 1. Redistributions of source code must retain the above copyright
  12. * notice, this list of conditions and the following disclaimer.
  13. * 2. Redistributions in binary form must reproduce the above copyright
  14. * notice, this list of conditions and the following disclaimer in
  15. * the documentation and/or other materials provided with the
  16. * distribution.
  17. * 3. Neither the name NuttX nor the names of its contributors may be
  18. * used to endorse or promote products derived from this software
  19. * without specific prior written permission.
  20. *
  21. * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
  22. * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
  23. * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
  24. * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
  25. * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
  26. * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
  27. * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS
  28. * OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED
  29. * AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  30. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
  31. * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  32. * POSSIBILITY OF SUCH DAMAGE.
  33. *
  34. ****************************************************************************/
  35. /****************************************************************************
  36. * Included Files
  37. ****************************************************************************/
  38. #include <nuttx/lib/math32.h>
  39. /****************************************************************************
  40. * Public Functions
  41. ****************************************************************************/
  42. /****************************************************************************
  43. * Name: umul32x64
  44. *
  45. * Description:
  46. * Multiply one 32-bit and one 64-bit values, factor1 and factor2,
  47. * respectively, and return the truncated 64-bit product.
  48. *
  49. * Input Parameters:
  50. * factor1 and factor2 - The values to be multiplied
  51. * product - The location to return the product of the two values.
  52. *
  53. ****************************************************************************/
  54. void umul32x64(uint32_t factor1, FAR const struct uint64_s *factor2,
  55. FAR struct uint64_s *product)
  56. {
  57. struct uint64_s part1;
  58. struct uint64_s part2;
  59. /* factor2 = factor2->ms << 32 + factor2->ls
  60. *
  61. * Full 128-bit product:
  62. * factor1 * factor2 = factor1 * (factor2->ms << 32) +
  63. * factor1 * factor2->ls
  64. */
  65. /* Get part1 = factor1 * factor2->ms, shifting left by 32-bits
  66. * (truncating to 64-bits)
  67. */
  68. part1.ms = factor1 * factor2->ms;
  69. part1.ls = 0;
  70. /* Get the full 64-bit part2 = factor1 * factor2->ls */
  71. umul32(factor1, factor2->ls, &part2);
  72. /* The product is then the sum */
  73. uadd64(&part1, &part2, product);
  74. }