Changeset 70e9af9 in rtems


Ignore:
Timestamp:
Nov 9, 1999, 10:45:35 PM (20 years ago)
Author:
Joel Sherrill <joel.sherrill@…>
Branches:
4.10, 4.11, 4.8, 4.9, master
Children:
3d67661
Parents:
e18ba7f
Message:

Changed B, H, and W variables to BVar, HVar, and WVar to avoid conflicts
with ITRON 3.0 data types regardless of how dumb it is to name
types with a simple letter.

Files:
2 edited

Legend:

Unmodified
Added
Removed
  • c/src/tests/samples/paranoia/paranoia.c

    re18ba7f r70e9af9  
    274274FLOAT   Third;
    275275FLOAT   F6, F9;
    276 FLOAT   H, HInvrse;
     276FLOAT   HVar, HInvrse;
    277277int     I;
    278278FLOAT   StickyBit, J;
     
    284284FLOAT   OneUlp, UfThold, U1, U2;
    285285FLOAT   V, V0, V9;
    286 FLOAT   W;
     286FLOAT   WVar;
    287287FLOAT   X, X1, X2, X8, Random1;
    288288FLOAT   Y, Y1, Y2, Random2;
     
    446446    }
    447447    printf ("Searching for Radix and Precision.\n");
    448     W = One;
     448    WVar = One;
    449449    do {
    450         W = W + W;
    451         Y = W + One;
    452         Z = Y - W;
     450        WVar = WVar + WVar;
     451        Y = WVar + One;
     452        Z = Y - WVar;
    453453        Y = Z - One;
    454454    }
    455455    while (MinusOne + FABS (Y) < Zero);
    456     /*.. now W is just big enough that |((W+1)-W)-1| >= 1 ...*/
     456    /*.. now WVar is just big enough that |((WVar+1)-WVar)-1| >= 1 ...*/
    457457    Precision = Zero;
    458458    Y = One;
    459459    do {
    460         Radix = W + Y;
     460        Radix = WVar + Y;
    461461        Y = Y + Y;
    462         Radix = Radix - W;
     462        Radix = Radix - WVar;
    463463    }
    464464    while (Radix == Zero);
     
    467467    printf ("Radix = %f .\n", Radix);
    468468    if (Radix != 1) {
    469         W = One;
     469        WVar = One;
    470470        do {
    471471            Precision = Precision + One;
    472             W = W * Radix;
    473             Y = W + One;
    474         }
    475         while ((Y - W) == One);
    476     }
    477     /*... now W == Radix^Precision is barely too big to satisfy (W+1)-W == 1
     472            WVar = WVar * Radix;
     473            Y = WVar + One;
     474        }
     475        while ((Y - WVar) == One);
     476    }
     477    /*... now WVar == Radix^Precision is barely too big to satisfy (WVar+1)-WVar == 1
    478478                                                      ...*/
    479     U1 = One / W;
     479    U1 = One / WVar;
    480480    U2 = Radix * U1;
    481481    printf ("Closest relative separation found is U1 = %.7e .\n\n", U1);
     
    529529    else
    530530        printf ("gets better closest relative separation U1 = %.7e .\n", U1);
    531     W = One / U1;
     531    WVar = One / U1;
    532532    F9 = (Half - U1) + Half;
    533533    Radix = FLOOR (0.01 + U2 / U1);
     
    653653    /*=============================================*/
    654654    if (Radix >= Two) {
    655         X = W / (Radix * Radix);
     655        X = WVar / (Radix * Radix);
    656656        Y = X + One;
    657657        Z = Y - X;
     
    11951195                    }
    11961196                    while (!(U2 * D >= F9));
    1197                     if (D * Radix - D != W - D)
     1197                    if (D * Radix - D != WVar - D)
    11981198                        Anomaly = True;
    11991199                    else {
     
    12081208                        SR3750 ();
    12091209                        NewD ();
    1210                         if (D - Z2 != W - Z2)
     1210                        if (D - Z2 != WVar - Z2)
    12111211                            Anomaly = True;
    12121212                        else {
     
    12261226        if ((I == 0) || Anomaly) {
    12271227            BadCond (Failure, "Anomalous arithmetic with Integer < ");
    1228             printf ("Radix^Precision = %.7e\n", W);
     1228            printf ("Radix^Precision = %.7e\n", WVar);
    12291229            printf (" fails test whether sqrt rounds or chops.\n");
    12301230            SqRWrng = True;
     
    12831283    N = 0;
    12841284    Z = A1;
    1285     M = (int) FLOOR (Two * LOG (W) / LOG (A1));
     1285    M = (int) FLOOR (Two * LOG (WVar) / LOG (A1));
    12861286    Break = False;
    12871287    do {
     
    13591359    else
    13601360        HInvrse = Radix;
    1361     H = One / HInvrse;
    1362     /* ... 1/HInvrse == H == Min(1/Radix, 1/2) */
     1361    HVar = One / HInvrse;
     1362    /* ... 1/HInvrse == HVar == Min(1/Radix, 1/2) */
    13631363    CInvrse = One / C;
    13641364    E0 = C;
    1365     Z = E0 * H;
     1365    Z = E0 * HVar;
    13661366    /* ...1/Radix^(BIG Integer) << 1 << CInvrse == 1/C */
    13671367    do {
    13681368        Y = E0;
    13691369        E0 = Z;
    1370         Z = E0 * H;
     1370        Z = E0 * HVar;
    13711371    }
    13721372    while ((E0 > Z) && (Z + Z > Z));
     
    13901390    } else {
    13911391        Underflow = D;
    1392         PseudoZero = Underflow * H;
     1392        PseudoZero = Underflow * HVar;
    13931393        UfThold = Zero;
    13941394        do {
     
    14021402                    UfThold = Y1;
    14031403            }
    1404             PseudoZero = PseudoZero * H;
     1404            PseudoZero = PseudoZero * HVar;
    14051405        }
    14061406        while ((Underflow > PseudoZero)
     
    14331433    /*=============================================*/
    14341434    if (CInvrse * Y > CInvrse * Y1) {
    1435         S = H * S;
     1435        S = HVar * S;
    14361436        E0 = Underflow;
    14371437    }
     
    15031503        if (setjmp (ovfl_buf)) {
    15041504            printf ("Underflow / UfThold failed!\n");
    1505             R = H + H;
     1505            R = HVar + HVar;
    15061506        } else
    15071507            R = SQRT (Underflow / UfThold);
    15081508        sigsave = 0;
    1509         if (R <= H) {
     1509        if (R <= HVar) {
    15101510            Z = R * UfThold;
    1511             X = Z * (One + R * H * (One + H));
     1511            X = Z * (One + R * HVar * (One + HVar));
    15121512        } else {
    15131513            Z = UfThold;
    1514             X = Z * (One + H * H * (One + H));
     1514            X = Z * (One + HVar * HVar * (One + HVar));
    15151515        }
    15161516        if (!((X == Z) || (X - Z != Zero))) {
     
    17801780        if (((V9 < (One - Two * Radix * E9) * Z) || (V9 > Z))) {
    17811781            Y = V9;
    1782             if (X < W)
     1782            if (X < WVar)
    17831783                BadCond (Serious, "");
    17841784            else
     
    20672067SR3750 ()
    20682068{
    2069     if (!((X - Radix < Z2 - Radix) || (X - Z2 > W - Z2))) {
     2069    if (!((X - Radix < Z2 - Radix) || (X - Z2 > WVar - Z2))) {
    20702070        I = I + 1;
    20712071        X2 = SQRT (X * D);
     
    21122112        X = Z * X;
    21132113    }
    2114     while (X < W);
     2114    while (X < WVar);
    21152115}
    21162116
  • testsuites/samples/paranoia/paranoia.c

    re18ba7f r70e9af9  
    274274FLOAT   Third;
    275275FLOAT   F6, F9;
    276 FLOAT   H, HInvrse;
     276FLOAT   HVar, HInvrse;
    277277int     I;
    278278FLOAT   StickyBit, J;
     
    284284FLOAT   OneUlp, UfThold, U1, U2;
    285285FLOAT   V, V0, V9;
    286 FLOAT   W;
     286FLOAT   WVar;
    287287FLOAT   X, X1, X2, X8, Random1;
    288288FLOAT   Y, Y1, Y2, Random2;
     
    446446    }
    447447    printf ("Searching for Radix and Precision.\n");
    448     W = One;
     448    WVar = One;
    449449    do {
    450         W = W + W;
    451         Y = W + One;
    452         Z = Y - W;
     450        WVar = WVar + WVar;
     451        Y = WVar + One;
     452        Z = Y - WVar;
    453453        Y = Z - One;
    454454    }
    455455    while (MinusOne + FABS (Y) < Zero);
    456     /*.. now W is just big enough that |((W+1)-W)-1| >= 1 ...*/
     456    /*.. now WVar is just big enough that |((WVar+1)-WVar)-1| >= 1 ...*/
    457457    Precision = Zero;
    458458    Y = One;
    459459    do {
    460         Radix = W + Y;
     460        Radix = WVar + Y;
    461461        Y = Y + Y;
    462         Radix = Radix - W;
     462        Radix = Radix - WVar;
    463463    }
    464464    while (Radix == Zero);
     
    467467    printf ("Radix = %f .\n", Radix);
    468468    if (Radix != 1) {
    469         W = One;
     469        WVar = One;
    470470        do {
    471471            Precision = Precision + One;
    472             W = W * Radix;
    473             Y = W + One;
    474         }
    475         while ((Y - W) == One);
    476     }
    477     /*... now W == Radix^Precision is barely too big to satisfy (W+1)-W == 1
     472            WVar = WVar * Radix;
     473            Y = WVar + One;
     474        }
     475        while ((Y - WVar) == One);
     476    }
     477    /*... now WVar == Radix^Precision is barely too big to satisfy (WVar+1)-WVar == 1
    478478                                                      ...*/
    479     U1 = One / W;
     479    U1 = One / WVar;
    480480    U2 = Radix * U1;
    481481    printf ("Closest relative separation found is U1 = %.7e .\n\n", U1);
     
    529529    else
    530530        printf ("gets better closest relative separation U1 = %.7e .\n", U1);
    531     W = One / U1;
     531    WVar = One / U1;
    532532    F9 = (Half - U1) + Half;
    533533    Radix = FLOOR (0.01 + U2 / U1);
     
    653653    /*=============================================*/
    654654    if (Radix >= Two) {
    655         X = W / (Radix * Radix);
     655        X = WVar / (Radix * Radix);
    656656        Y = X + One;
    657657        Z = Y - X;
     
    11951195                    }
    11961196                    while (!(U2 * D >= F9));
    1197                     if (D * Radix - D != W - D)
     1197                    if (D * Radix - D != WVar - D)
    11981198                        Anomaly = True;
    11991199                    else {
     
    12081208                        SR3750 ();
    12091209                        NewD ();
    1210                         if (D - Z2 != W - Z2)
     1210                        if (D - Z2 != WVar - Z2)
    12111211                            Anomaly = True;
    12121212                        else {
     
    12261226        if ((I == 0) || Anomaly) {
    12271227            BadCond (Failure, "Anomalous arithmetic with Integer < ");
    1228             printf ("Radix^Precision = %.7e\n", W);
     1228            printf ("Radix^Precision = %.7e\n", WVar);
    12291229            printf (" fails test whether sqrt rounds or chops.\n");
    12301230            SqRWrng = True;
     
    12831283    N = 0;
    12841284    Z = A1;
    1285     M = (int) FLOOR (Two * LOG (W) / LOG (A1));
     1285    M = (int) FLOOR (Two * LOG (WVar) / LOG (A1));
    12861286    Break = False;
    12871287    do {
     
    13591359    else
    13601360        HInvrse = Radix;
    1361     H = One / HInvrse;
    1362     /* ... 1/HInvrse == H == Min(1/Radix, 1/2) */
     1361    HVar = One / HInvrse;
     1362    /* ... 1/HInvrse == HVar == Min(1/Radix, 1/2) */
    13631363    CInvrse = One / C;
    13641364    E0 = C;
    1365     Z = E0 * H;
     1365    Z = E0 * HVar;
    13661366    /* ...1/Radix^(BIG Integer) << 1 << CInvrse == 1/C */
    13671367    do {
    13681368        Y = E0;
    13691369        E0 = Z;
    1370         Z = E0 * H;
     1370        Z = E0 * HVar;
    13711371    }
    13721372    while ((E0 > Z) && (Z + Z > Z));
     
    13901390    } else {
    13911391        Underflow = D;
    1392         PseudoZero = Underflow * H;
     1392        PseudoZero = Underflow * HVar;
    13931393        UfThold = Zero;
    13941394        do {
     
    14021402                    UfThold = Y1;
    14031403            }
    1404             PseudoZero = PseudoZero * H;
     1404            PseudoZero = PseudoZero * HVar;
    14051405        }
    14061406        while ((Underflow > PseudoZero)
     
    14331433    /*=============================================*/
    14341434    if (CInvrse * Y > CInvrse * Y1) {
    1435         S = H * S;
     1435        S = HVar * S;
    14361436        E0 = Underflow;
    14371437    }
     
    15031503        if (setjmp (ovfl_buf)) {
    15041504            printf ("Underflow / UfThold failed!\n");
    1505             R = H + H;
     1505            R = HVar + HVar;
    15061506        } else
    15071507            R = SQRT (Underflow / UfThold);
    15081508        sigsave = 0;
    1509         if (R <= H) {
     1509        if (R <= HVar) {
    15101510            Z = R * UfThold;
    1511             X = Z * (One + R * H * (One + H));
     1511            X = Z * (One + R * HVar * (One + HVar));
    15121512        } else {
    15131513            Z = UfThold;
    1514             X = Z * (One + H * H * (One + H));
     1514            X = Z * (One + HVar * HVar * (One + HVar));
    15151515        }
    15161516        if (!((X == Z) || (X - Z != Zero))) {
     
    17801780        if (((V9 < (One - Two * Radix * E9) * Z) || (V9 > Z))) {
    17811781            Y = V9;
    1782             if (X < W)
     1782            if (X < WVar)
    17831783                BadCond (Serious, "");
    17841784            else
     
    20672067SR3750 ()
    20682068{
    2069     if (!((X - Radix < Z2 - Radix) || (X - Z2 > W - Z2))) {
     2069    if (!((X - Radix < Z2 - Radix) || (X - Z2 > WVar - Z2))) {
    20702070        I = I + 1;
    20712071        X2 = SQRT (X * D);
     
    21122112        X = Z * X;
    21132113    }
    2114     while (X < W);
     2114    while (X < WVar);
    21152115}
    21162116
Note: See TracChangeset for help on using the changeset viewer.