From: <cli...@li...> - 2005-02-26 04:19:17
|
Send clisp-cvs mailing list submissions to cli...@li... To subscribe or unsubscribe via the World Wide Web, visit https://lists.sourceforge.net/lists/listinfo/clisp-cvs or, via email, send a message with subject or body 'help' to cli...@li... You can reach the person managing the list at cli...@li... When replying, please edit your Subject line so it is more specific than "Re: Contents of clisp-cvs digest..." CLISP CVS commits for today Today's Topics: 1. clisp/src ChangeLog,1.4289,1.4290 comptran.d,1.32,1.33 lisparit.d,1.75,1.76 (Bruno Haible) 2. clisp/src realtran.d,1.18,1.19 comptran.d,1.33,1.34 lisparit.d,1.76,1.77 ChangeLog,1.4290,1.4291 (Bruno Haible) 3. clisp/src realtran.d,1.19,1.20 comptran.d,1.34,1.35 ChangeLog,1.4291,1.4292 (Bruno Haible) --__--__-- Message: 1 From: Bruno Haible <ha...@us...> To: cli...@li... Subject: clisp/src ChangeLog,1.4289,1.4290 comptran.d,1.32,1.33 lisparit.d,1.75,1.76 Date: Fri, 25 Feb 2005 13:14:08 +0000 Reply-To: cli...@li... Update of /cvsroot/clisp/clisp/src In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv2972/src Modified Files: ChangeLog comptran.d lisparit.d Log Message: Remove start_p argument - always true. Index: lisparit.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/lisparit.d,v retrieving revision 1.75 retrieving revision 1.76 diff -u -d -r1.75 -r1.76 --- lisparit.d 29 Jan 2005 14:55:06 -0000 1.75 +++ lisparit.d 25 Feb 2005 13:14:05 -0000 1.76 @@ -930,7 +930,8 @@ if (complexp(STACK_0)) pushSTACK(R_R_contagion_R(TheComplex(STACK_0)->c_real, TheComplex(STACK_0)->c_imag)); - else pushSTACK(STACK_0); + else + pushSTACK(STACK_0); VALUES1(N_exp_N(STACK_1,true,&STACK_0)); skipSTACK(2); } @@ -953,7 +954,7 @@ arg = STACK_1; } else STACK_0 = arg; - VALUES1(N_log_N(arg,true,&STACK_0)); + VALUES1(N_log_N(arg,&STACK_0)); } else { /* LOG with two arguments */ base = check_number(base); VALUES1(N_N_log_N(STACK_1,base)); Index: comptran.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/comptran.d,v retrieving revision 1.32 retrieving revision 1.33 diff -u -d -r1.32 -r1.33 --- comptran.d 21 Feb 2005 12:36:56 -0000 1.32 +++ comptran.d 25 Feb 2005 13:14:05 -0000 1.33 @@ -81,33 +81,32 @@ } } -/* N_log_N(x) liefert (log x), wo x eine Zahl ist. +/* N_log_N(x,&end_precision) liefert (log x), wo x eine Zahl ist. can trigger GC Methode: (complex (log (abs x)) (phase x)) */ -local maygc object N_log_N (object x, bool start_p, gcv_object_t *end_p) +local maygc object N_log_N (object x, gcv_object_t *end_p) { pushSTACK(x); /* save x */ pushSTACK(N_abs_R(x)); /* (abs x) */ if (R_zerop(STACK_0)) /* (abs x) = 0 -> Error */ divide_0(); - STACK_0 = R_ln_R(STACK_0,start_p,end_p); /* (log (abs x)) */ - if (start_p) { /* increase precision */ - if (floatp(STACK_1)) - STACK_1 = F_extend_F(STACK_1); - else if (complexp(STACK_1) - && (floatp(TheComplex(STACK_1)->c_real) - || floatp(TheComplex(STACK_1)->c_imag))) { - var object realpart = TheComplex(STACK_1)->c_real; - if (floatp(realpart)) - realpart = F_extend_F(realpart); - pushSTACK(realpart); - var object imagpart = TheComplex(STACK_(1+1))->c_imag; - if (floatp(imagpart)) - imagpart = F_extend_F(imagpart); - realpart = popSTACK(); - STACK_1 = R_R_complex_C(realpart,imagpart); - } + STACK_0 = R_ln_R(STACK_0,true,end_p); /* (log (abs x)) */ + /* Increase precision: */ + if (floatp(STACK_1)) + STACK_1 = F_extend_F(STACK_1); + else if (complexp(STACK_1) + && (floatp(TheComplex(STACK_1)->c_real) + || floatp(TheComplex(STACK_1)->c_imag))) { + var object realpart = TheComplex(STACK_1)->c_real; + if (floatp(realpart)) + realpart = F_extend_F(realpart); + pushSTACK(realpart); + var object imagpart = TheComplex(STACK_(1+1))->c_imag; + if (floatp(imagpart)) + imagpart = F_extend_F(imagpart); + realpart = popSTACK(); + STACK_1 = R_R_complex_C(realpart,imagpart); } STACK_1 = N_phase_R(STACK_1,true); /* (phase x) */ if (end_p != NULL && floatp(STACK_1)) @@ -206,8 +205,8 @@ } } else { /* normal complex case */ pushSTACK(a); pushSTACK(b); - STACK_1 = N_log_N(STACK_1,true,&STACK_1); /* (log a) */ - STACK_0 = N_log_N(STACK_0,true,&STACK_0); /* (log b) */ + STACK_1 = N_log_N(STACK_1,&STACK_1); /* (log a) */ + STACK_0 = N_log_N(STACK_0,&STACK_0); /* (log b) */ a = N_N_durch_N(STACK_1,STACK_0); /* divide */ skipSTACK(2); return a; } @@ -463,7 +462,7 @@ var uintL x_prec = R_float_digits(x); if (x_prec < F_float_digits(STACK_0)) STACK_2 = N_N_float_N(STACK_2,STACK_0); /* extend precision of x */ - STACK_2 = N_log_N(STACK_2,true,NULL); /* (log x) */ + STACK_2 = N_log_N(STACK_2,NULL); /* (log x) */ STACK_2 = N_N_float_N(STACK_2,STACK_0); /* rounded (log x) */ STACK_4 = N_N_float_N(STACK_4,STACK_0); /* rounded y */ var object temp = N_N_mal_N(STACK_2,STACK_4); /* (* (log x) y) */ Index: ChangeLog =================================================================== RCS file: /cvsroot/clisp/clisp/src/ChangeLog,v retrieving revision 1.4289 retrieving revision 1.4290 diff -u -d -r1.4289 -r1.4290 --- ChangeLog 25 Feb 2005 03:29:03 -0000 1.4289 +++ ChangeLog 25 Feb 2005 13:13:57 -0000 1.4290 @@ -1,3 +1,9 @@ +2005-02-20 Bruno Haible <br...@cl...> + + * comptran.d (N_log_N): Remove start_p argument. + (N_N_log_N, N_N_expt_N): Update. + * lisparit.d (LOG): Update. + 2005-02-24 Sam Steingold <sd...@gn...> * lispbibl.d (HAVE_DISASSEMBLER): define always --__--__-- Message: 2 From: Bruno Haible <ha...@us...> To: cli...@li... Subject: clisp/src realtran.d,1.18,1.19 comptran.d,1.33,1.34 lisparit.d,1.76,1.77 ChangeLog,1.4290,1.4291 Date: Fri, 25 Feb 2005 13:15:48 +0000 Reply-To: cli...@li... Update of /cvsroot/clisp/clisp/src In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv3474/src Modified Files: realtran.d comptran.d lisparit.d ChangeLog Log Message: Remove start_p argument - always true. Index: lisparit.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/lisparit.d,v retrieving revision 1.76 retrieving revision 1.77 diff -u -d -r1.76 -r1.77 --- lisparit.d 25 Feb 2005 13:14:05 -0000 1.76 +++ lisparit.d 25 Feb 2005 13:15:42 -0000 1.77 @@ -1807,7 +1807,7 @@ > 1 , so that it is not recalculated to often: */ oldlen += floor(oldlen,2); /* oldlen * 3/2 */ var uintC newlen = (d < oldlen ? oldlen : d); - ln_x = *objptr = LF_shorten_LF(R_ln_R(I_to_LF(x,newlen,true),true,NULL), + ln_x = *objptr = LF_shorten_LF(R_ln_R(I_to_LF(x,newlen,true),NULL), newlen); /* (ln x) - LF of len newlen */ return (d < newlen ? LF_shorten_LF(ln_x,d) : ln_x); } else if (d > FF_mant_len+1) /* a double-float is sufficient */ Index: realtran.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/realtran.d,v retrieving revision 1.18 retrieving revision 1.19 diff -u -d -r1.18 -r1.19 --- realtran.d 21 Feb 2005 12:36:56 -0000 1.18 +++ realtran.d 25 Feb 2005 13:15:42 -0000 1.19 @@ -778,7 +778,7 @@ ); } -# R_ln_R(x) liefert zu einer reellen Zahl x>0 die Zahl ln(x). +# R_ln_R(x,&end_precision) liefert zu einer reellen Zahl x>0 die Zahl ln(x). # can trigger GC # Methode: # x rational -> bei x=1 0 als Ergebnis, sonst x in Float umwandeln. @@ -788,7 +788,7 @@ # (m,e) := (decode-float x), so dass 1/2 <= m < 1. # m<2/3 -> m:=2m, e:=e-1, so dass 2/3 <= m <= 4/3. # ln(m) errechnen, ln(x)=ln(m)+e*ln(2) als Ergebnis. -local maygc object R_ln_R (object x, bool start_p, gcv_object_t* end_p) +local maygc object R_ln_R (object x, gcv_object_t* end_p) { if (R_rationalp(x)) { if (eq(x,Fixnum_1)) { return Fixnum_0; } /* x=1 -> return 0 */ @@ -796,8 +796,7 @@ } /* x -- float */ pushSTACK(x); /* save x */ - if (start_p) /* increase computational precision */ - x = F_extend2_F(x); + x = F_extend2_F(x); /* increase computational precision */ F_decode_float_F_I_F(x); /* compute m,e,s */ /* Stack layout: x, m, e, s. */ if (F_F_comp(STACK_2, @@ -1006,8 +1005,8 @@ } } # Nun a,b beide Floats. - pushSTACK(R_ln_R(STACK_1,true,NULL)); /* (ln a) */ - pushSTACK(R_ln_R(STACK_1,true,NULL)); /* (ln b) */ + pushSTACK(R_ln_R(STACK_1,NULL)); /* (ln a) */ + pushSTACK(R_ln_R(STACK_1,NULL)); /* (ln b) */ STACK_0 = F_F_durch_F(STACK_1,STACK_0); /* (/ (ln a) (ln b)) */ STACK_1 = R_R_contagion_R(STACK_2,STACK_3); var object ret = F_R_float_F(STACK_0,STACK_1); Index: comptran.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/comptran.d,v retrieving revision 1.33 retrieving revision 1.34 diff -u -d -r1.33 -r1.34 --- comptran.d 25 Feb 2005 13:14:05 -0000 1.33 +++ comptran.d 25 Feb 2005 13:15:42 -0000 1.34 @@ -91,7 +91,7 @@ pushSTACK(N_abs_R(x)); /* (abs x) */ if (R_zerop(STACK_0)) /* (abs x) = 0 -> Error */ divide_0(); - STACK_0 = R_ln_R(STACK_0,true,end_p); /* (log (abs x)) */ + STACK_0 = R_ln_R(STACK_0,end_p); /* (log (abs x)) */ /* Increase precision: */ if (floatp(STACK_1)) STACK_1 = F_extend_F(STACK_1); @@ -160,7 +160,7 @@ b = STACK_1; if (R_rationalp(b)) b = RA_F_float_F(b,angle,true); - b = F_ln_F(b,true,&STACK_1); STACK_0 = F_F_durch_F(STACK_0,b); + b = F_ln_F(b,&STACK_1); STACK_0 = F_F_durch_F(STACK_0,b); } /* Stackaufbau: a, b, Imaginärteil. Realteil (/ (log (abs a)) (log b)) errechnen: */ @@ -190,12 +190,12 @@ } } /* Keine Chance für rationalen Realteil */ - pushSTACK(F_ln_F(N_abs_R(a),true,&STACK_3)); /* (log (abs a)), a float */ + pushSTACK(F_ln_F(N_abs_R(a),&STACK_3)); /* (log (abs a)), a float */ /* durch (log b) dividieren, liefert den Realteil: */ b = STACK_2; if (R_rationalp(b)) b = RA_F_float_F(b,STACK_0,true); - b = F_ln_F(b,true,&STACK_2); STACK_0 = F_F_durch_F(STACK_0,b); + b = F_ln_F(b,&STACK_2); STACK_0 = F_F_durch_F(STACK_0,b); real_ok: /* stack layout: a, b, imagpart, realpart. */ { @@ -877,7 +877,7 @@ STACK_0 = temp; } /* stack layout: |(1+x)/(1-x)| (>0), Im. */ - STACK_1 = F_I_scale_float_F(R_ln_R(STACK_1,true,&STACK_1),Fixnum_minus1); /* ln / 2 */ + STACK_1 = F_I_scale_float_F(R_ln_R(STACK_1,&STACK_1),Fixnum_minus1); /* ln / 2 */ return; } pushSTACK(x); pushSTACK(y); @@ -916,7 +916,7 @@ temp = F_F_durch_F(STACK_0,temp); /* ((1+x)^2+y^2)/((1-x)^2+y^2), a float >=0 */ if (R_zerop(temp)) /* should be >0 */ divide_0(); - temp = R_ln_R(temp,true,NULL); /* ln(temp), a float */ + temp = R_ln_R(temp,NULL); /* ln(temp), a float */ STACK_6 = F_I_scale_float_F(temp,sfixnum(-2)); /* .../4 =: u */ } } @@ -1075,7 +1075,7 @@ else temp = F_F_plus_F(temp,y); /* temp = sqrt(y^2-1)+|y|, ein Float >1 */ - STACK_1 = R_ln_R(temp,true,&STACK_0); /* ln(|y|+sqrt(y^2-1)), Float >0 */ + STACK_1 = R_ln_R(temp,&STACK_0); /* ln(|y|+sqrt(y^2-1)), Float >0 */ temp = F_I_scale_float_F(pi(STACK_1),Fixnum_minus1); /* (scale-float pi -1) = pi/2 */ if (!R_minusp(STACK_0)) { /* Vorzeichen von y */ /* y>1 -> v = pi/2 */ @@ -1102,9 +1102,9 @@ STACK_1 = F_atanhx_F(F_F_durch_F(x,temp)); /* u = atanh(x/sqrt(1+x^2)) */ } else { /* |x| >= 1/2 */ if (!R_minusp(x)) /* x >= 1/2 */ - STACK_1 = R_ln_R(F_F_plus_F(temp,x),true,&STACK_1); /* u = ln(x+sqrt(1+x^2)) */ + STACK_1 = R_ln_R(F_F_plus_F(temp,x),&STACK_1); /* u = ln(x+sqrt(1+x^2)) */ else /* x <= -1/2 */ - STACK_1 = F_minus_F(R_ln_R(F_F_minus_F(temp,x),true,&STACK_1)); /* u = -ln(-x+sqrt(1+x^2)) */ + STACK_1 = F_minus_F(R_ln_R(F_F_minus_F(temp,x),&STACK_1)); /* u = -ln(-x+sqrt(1+x^2)) */ } return; } @@ -1203,7 +1203,7 @@ temp = R_R_minus_R(F_square_F(temp),Fixnum_1); /* z^2-1, ein Float >=0 */ temp = F_sqrt_F(temp); /* sqrt(z^2-1), ein Float >=0 */ temp = F_F_plus_F(STACK_0,temp); /* z+sqrt(z^2-1), float >1 */ - temp = R_ln_R(temp,true,&STACK_0); /* ln(z+sqrt(z^2-1)), float >=0 */ + temp = R_ln_R(temp,&STACK_0); /* ln(z+sqrt(z^2-1)), float >=0 */ skipSTACK(1); return R_R_complex_C(Fixnum_0,temp); } @@ -1272,7 +1272,7 @@ STACK_0 = z = RA_float_F(z); /* z Float <= -1 */ z = F_sqrt_F(R_R_minus_R(F_square_F(z),Fixnum_1)); /* sqrt(z^2-1), ein Float >=0 */ - STACK_0 = R_ln_R(F_F_minus_F(z,STACK_0),true,&STACK_0); /* log(sqrt(z^2-1)-z), ein Float >=0 */ + STACK_0 = R_ln_R(F_F_minus_F(z,STACK_0),&STACK_0); /* log(sqrt(z^2-1)-z), ein Float >=0 */ z = pi(STACK_0); /* and imaginary part == pi */ return R_R_complex_C(popSTACK(),z); } Index: ChangeLog =================================================================== RCS file: /cvsroot/clisp/clisp/src/ChangeLog,v retrieving revision 1.4290 retrieving revision 1.4291 diff -u -d -r1.4290 -r1.4291 --- ChangeLog 25 Feb 2005 13:13:57 -0000 1.4290 +++ ChangeLog 25 Feb 2005 13:15:42 -0000 1.4291 @@ -1,5 +1,13 @@ 2005-02-20 Bruno Haible <br...@cl...> + * realtran.d (R_ln_R): Remove start_p argument. + (R_R_log_R): Update. + * comptran.d (N_log_N, N_N_log_N, R_R_atanh_R_R, R_R_asinh_R_R, + N_acos_N, N_acosh_N): Update. + * lisparit.d (log_digits): Update. + +2005-02-20 Bruno Haible <br...@cl...> + * comptran.d (N_log_N): Remove start_p argument. (N_N_log_N, N_N_expt_N): Update. * lisparit.d (LOG): Update. --__--__-- Message: 3 From: Bruno Haible <ha...@us...> To: cli...@li... Subject: clisp/src realtran.d,1.19,1.20 comptran.d,1.34,1.35 ChangeLog,1.4291,1.4292 Date: Fri, 25 Feb 2005 13:16:53 +0000 Reply-To: cli...@li... Update of /cvsroot/clisp/clisp/src In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv3702/src Modified Files: realtran.d comptran.d ChangeLog Log Message: Remove start_p argument - always true. Index: comptran.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/comptran.d,v retrieving revision 1.34 retrieving revision 1.35 diff -u -d -r1.34 -r1.35 --- comptran.d 25 Feb 2005 13:15:42 -0000 1.34 +++ comptran.d 25 Feb 2005 13:16:50 -0000 1.35 @@ -495,7 +495,7 @@ if (R_rationalp(STACK_0)) /* b */ STACK_0 = RA_float_F(STACK_0); var bool same_precision = (F_float_digits(STACK_1) == F_float_digits(STACK_0)); - R_cosh_sinh_R_R(STACK_0,true,NULL); /* cosh(b) sinh(b) */ + R_cosh_sinh_R_R(STACK_0,NULL); /* cosh(b) sinh(b) */ /* stack layout: a, b, cosh(b), sinh(b). b != Fixnum_0 ==> sinh(b) != Fixnum_0. */ R_cos_sin_R_R(STACK_3,true,NULL); /* cos(a)!=0, sin(a) */ @@ -540,7 +540,7 @@ if (R_rationalp(STACK_0)) /* b */ STACK_0 = RA_float_F(STACK_0); var bool same_precision = (F_float_digits(STACK_1) == F_float_digits(STACK_0)); - R_cosh_sinh_R_R(STACK_0,true,NULL); /* cosh(b), sinh(b) */ + R_cosh_sinh_R_R(STACK_0,NULL); /* cosh(b), sinh(b) */ /* stack layout: a, b, cosh(b), sinh(b). */ R_cos_sin_R_R(STACK_3,true,NULL); /* cos(a), sin(a) */ /* stack layout: a, b, cosh(b), sinh(b), cos(a), sin(a). */ @@ -586,7 +586,7 @@ if (R_rationalp(STACK_0)) /* b */ STACK_0 = RA_float_F(STACK_0); var bool same_precision = (F_float_digits(STACK_1) == F_float_digits(STACK_0)); - R_cosh_sinh_R_R(STACK_0,true,NULL); /* cosh(b), sinh(b) */ + R_cosh_sinh_R_R(STACK_0,NULL); /* cosh(b), sinh(b) */ /* stack layout: a, b, cosh(b), sinh(b). */ R_cos_sin_R_R(STACK_3,true,NULL); /* cos(a), sin(a) */ /* stack layout: a, b, cosh(b), sinh(b), cos(a), sin(a). */ @@ -686,7 +686,7 @@ R_cos_sin_R_R(STACK_0,true,NULL); /* cos(b), sin(b) */ /* stack layout: a, b, cos(b), sin(b). b != Fixnum_0 ==> sin(b) != Fixnum_0. */ - R_cosh_sinh_R_R(STACK_3,true,NULL); /* cosh(a), sinh(a); cosh(a) != Fixnum 0 */ + R_cosh_sinh_R_R(STACK_3,NULL); /* cosh(a), sinh(a); cosh(a) != Fixnum 0 */ /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a). */ pushSTACK(R_R_contagion_R(STACK_4,STACK_5)); /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a), resfloat. */ @@ -730,7 +730,7 @@ var bool same_precision = (F_float_digits(STACK_1) == F_float_digits(STACK_0)); R_cos_sin_R_R(STACK_0,true,NULL); /* cos(b), sin(b) */ /* stack layout: a, b, cos(b), sin(b). */ - R_cosh_sinh_R_R(STACK_3,true,NULL); /* cosh(a), sinh(a) */ + R_cosh_sinh_R_R(STACK_3,NULL); /* cosh(a), sinh(a) */ /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a). */ pushSTACK(R_R_contagion_R(STACK_4,STACK_5)); /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a), resfloat. */ @@ -776,7 +776,7 @@ var bool same_precision = (F_float_digits(STACK_1) == F_float_digits(STACK_0)); R_cos_sin_R_R(STACK_0,true,NULL); /* cos(b), sin(b) */ /* stack layout: a, b, cos(b), sin(b). */ - R_cosh_sinh_R_R(STACK_3,true,NULL); /* cosh(a), sinh(a) */ + R_cosh_sinh_R_R(STACK_3,NULL); /* cosh(a), sinh(a) */ /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a). */ pushSTACK(R_R_contagion_R(STACK_4,STACK_5)); /* stack layout: a, b, cos(b), sin(b), cosh(a), sinh(a), resfloat. */ Index: realtran.d =================================================================== RCS file: /cvsroot/clisp/clisp/src/realtran.d,v retrieving revision 1.19 retrieving revision 1.20 diff -u -d -r1.19 -r1.20 --- realtran.d 25 Feb 2005 13:15:42 -0000 1.19 +++ realtran.d 25 Feb 2005 13:16:50 -0000 1.20 @@ -1229,7 +1229,7 @@ # falls e>0: y:=exp(x) errechnen, # (scale-float (+ y (/ y)) -1) und (scale-float (- y (/ y)) -1) bilden. # Genauigkeit wieder verringern. -local maygc void R_cosh_sinh_R_R (object x, bool start_p, gcv_object_t* end_p) +local maygc void R_cosh_sinh_R_R (object x, gcv_object_t* end_p) { if (R_rationalp(x)) { /* x rational */ if (eq(x,Fixnum_0)) /* x=0 -> return (1,0) */ @@ -1241,7 +1241,7 @@ if (e > 0) { /* e>0 -> use exp(x) */ var object temp; pushSTACK(x); - pushSTACK(temp = R_exp_R(x,start_p,NULL)); /* y:=exp(x) */ + pushSTACK(temp = R_exp_R(x,true,NULL)); /* y:=exp(x) */ pushSTACK(temp = F_durch_F(temp)); /* (/ y) */ /* stack layout: x, exp(x), exp(-x). */ temp = F_F_plus_F(STACK_1,temp); /* + y */ @@ -1255,8 +1255,7 @@ } else { /* e<=0 */ if (R_zerop(x) /* e <= (1-d)/2 <==> e <= -ceiling((d-1)/2) ? */ || (e <= (sintL)(1-F_float_digits(x))>>1)) { - if (start_p) /* increase computational precision */ - x = F_extend_F(x); + x = F_extend_F(x); /* increase computational precision */ pushSTACK(x); pushSTACK(x); if (end_p != NULL) { STACK_1 = RA_R_float_F(Fixnum_1,*end_p); /* cosh=1 */ @@ -1266,7 +1265,7 @@ return; } pushSTACK(x); - { var object temp = (start_p ? F_extend_F(x) : x); + { var object temp = F_extend_F(x); pushSTACK(temp); pushSTACK(F_square_F(temp)); /* x*x */ pushSTACK(temp = F_sinhx_F(STACK_1)); /* y:=(sinh(x)/x)^2 */ @@ -1290,7 +1289,7 @@ local maygc object R_tanh_R (object x) { pushSTACK(x); - R_cosh_sinh_R_R(x,true,NULL); + R_cosh_sinh_R_R(x,NULL); /* stack layout: x, cosh(x), sinh(x). */ var object result = R_R_durch_R(STACK_0,STACK_1); if (floatp(STACK_0) || floatp(STACK_1)) Index: ChangeLog =================================================================== RCS file: /cvsroot/clisp/clisp/src/ChangeLog,v retrieving revision 1.4291 retrieving revision 1.4292 diff -u -d -r1.4291 -r1.4292 --- ChangeLog 25 Feb 2005 13:15:42 -0000 1.4291 +++ ChangeLog 25 Feb 2005 13:16:50 -0000 1.4292 @@ -1,5 +1,12 @@ 2005-02-20 Bruno Haible <br...@cl...> + * realtran.d (R_cosh_sinh_R_R): Remove start_p argument. + (R_tanh_R): Update. + * comptran.d (N_sin_N, N_cos_N, N_tan_N, N_sinh_N, N_cosh_N, N_tanh_N): + Update. + +2005-02-20 Bruno Haible <br...@cl...> + * realtran.d (R_ln_R): Remove start_p argument. (R_R_log_R): Update. * comptran.d (N_log_N, N_N_log_N, R_R_atanh_R_R, R_R_asinh_R_R, --__--__-- _______________________________________________ clisp-cvs mailing list cli...@li... https://lists.sourceforge.net/lists/listinfo/clisp-cvs End of clisp-cvs Digest |