MINICERN/mathlib/gen/e/ndigit.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:26  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       FUNCTION NDIGIT(U,V)
  11
  12 #include "gen/imp64.inc"
  13
  14 ************************************************************************
  15 *   NORBAS, VERSION: 15.03.1993
  16 ************************************************************************
  17 *
  18 *   NDIGIT  COMPUTES THE NUMBER OF DECIMAL DIGITS AT THE BEGINNING OF  U
  19 *   WHICH ARE AGREE WITH CORRESPONDING DIGITS OF  V
  20 *
  21 ************************************************************************
  22
  23       PARAMETER (SMALL = 1D-30 , TWO = 2)
  24       IF(ABS(U) .LE. SMALL .AND. ABS(V) .LE. SMALL) THEN
  25        NDIGIT=16
  26       ELSE IF(U*V .LT. 0) THEN
  27        NDIGIT=0
  28       ELSE
  29        NDIGIT=MIN(16,ABS(INT(LOG10(MAX(SMALL,ABS((U+V))))) -
  30      1                   INT(LOG10(MAX(SMALL,TWO*ABS(U-V))))))
  31       ENDIF
  32       RETURN
  33       END
  34