5 * Revision 1.1.1.1 1996/04/01 15:02:26 mclareni
12 #include "gen/imp64.inc"
14 ************************************************************************
15 * NORBAS, VERSION: 15.03.1993
16 ************************************************************************
18 * NDIGIT COMPUTES THE NUMBER OF DECIMAL DIGITS AT THE BEGINNING OF U
19 * WHICH ARE AGREE WITH CORRESPONDING DIGITS OF V
21 ************************************************************************
23 PARAMETER (SMALL = 1D-30 , TWO = 2)
24 IF(ABS(U) .LE. SMALL .AND. ABS(V) .LE. SMALL) THEN
26 ELSE IF(U*V .LT. 0) THEN
29 NDIGIT=MIN(16,ABS(INT(LOG10(MAX(SMALL,ABS((U+V))))) -
30 1 INT(LOG10(MAX(SMALL,TWO*ABS(U-V))))))