]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |