qemu/roms/ipxe/src/core/isqrt.c

   1 /*
   2  * Copyright (C) 2014 Michael Brown <mbrown@fensystems.co.uk>.
   3  *
   4  * This program is free software; you can redistribute it and/or
   5  * modify it under the terms of the GNU General Public License as
   6  * published by the Free Software Foundation; either version 2 of the
   7  * License, or any later version.
   8  *
   9  * This program is distributed in the hope that it will be useful, but
  10  * WITHOUT ANY WARRANTY; without even the implied warranty of
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12  * General Public License for more details.
  13  *
  14  * You should have received a copy of the GNU General Public License
  15  * along with this program; if not, write to the Free Software
  16  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
  17  * 02110-1301, USA.
  18  */
  19
  20 FILE_LICENCE ( GPL2_OR_LATER );
  21
  22 /** @file
  23  *
  24  * Integer square root
  25  *
  26  */
  27
  28 #include <ipxe/isqrt.h>
  29
  30 /**
  31  * Find integer square root
  32  *
  33  * @v value             Value
  34  * @v isqrt             Integer square root of value
  35  */
  36 unsigned long isqrt ( unsigned long value ) {
  37         unsigned long result = 0;
  38         unsigned long bit = ( 1UL << ( ( 8 * sizeof ( bit ) ) - 2 ) );
  39
  40         while ( bit > value )
  41                 bit >>= 2;
  42         while ( bit ) {
  43                 if ( value >= ( result + bit ) ) {
  44                         value -= ( result + bit );
  45                         result = ( ( result >> 1 ) + bit );
  46                 } else {
  47                         result >>= 1;
  48                 }
  49                 bit >>= 2;
  50         }
  51         return result;
  52 }