[Bug fortran/33197] Fortran 2008: math functions

2014-06-09 Thread fxcoudert at gcc dot gnu.org
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=33197 Bug 33197 depends on bug 36096, which changed state. Bug 36096 Summary: F2008 Bessel: Documentation/diagnostic errors https://gcc.gnu.org/bugzilla/show_bug.cgi?id=36096 What|Removed |Added --

[Bug fortran/33197] Fortran 2008: math functions

2010-08-27 Thread burnus at gcc dot gnu dot org
--- Comment #41 from burnus at gcc dot gnu dot org 2010-08-27 19:20 --- FIXED on the 4.6 trunk. Thanks to all involved in implementing the math intrinsics. -- burnus at gcc dot gnu dot org changed: What|Removed |Added -

[Bug fortran/33197] Fortran 2008: math functions

2010-08-27 Thread burnus at gcc dot gnu dot org
--- Comment #40 from burnus at gcc dot gnu dot org 2010-08-27 19:18 --- Subject: Bug 33197 Author: burnus Date: Fri Aug 27 19:17:45 2010 New Revision: 163595 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=163595 Log: gcc/fortran/ 2010-08-27 Tobias Burnus PR fortran/3

[Bug fortran/33197] Fortran 2008: math functions

2010-08-22 Thread burnus at gcc dot gnu dot org
--- Comment #39 from burnus at gcc dot gnu dot org 2010-08-22 21:02 --- (In reply to comment #37) > PS: NORM2 is described as "careful calculation of Euclidean norm" in the BCS > slides and in the what's new in F2008 article. Currently, I use the trivial > brute-force method. Maybe somet

[Bug fortran/33197] Fortran 2008: math functions

2010-08-22 Thread burnus at gcc dot gnu dot org
--- Comment #38 from burnus at gcc dot gnu dot org 2010-08-22 19:34 --- (In reply to comment #37) > - PARITY: ... and the same issues as for NORM2. And PARITY documentation is missing. Reminder: For NORM2's inline version, think also about hypot. -- http://gcc.gnu.org/bugzilla/show

[Bug fortran/33197] Fortran 2008: math functions

2010-08-22 Thread burnus at gcc dot gnu dot org
--- Comment #37 from burnus at gcc dot gnu dot org 2010-08-22 19:32 --- Created an attachment (id=21545) --> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=21545&action=view) Draft patch for NORM2 -- and draft run-time test, early PARITY patch Draft patch for NORM2; as PARITY is essent

[Bug fortran/33197] Fortran 2008: math functions

2010-08-21 Thread burnus at gcc dot gnu dot org
--- Comment #36 from burnus at gcc dot gnu dot org 2010-08-21 10:35 --- TODO (carried on): - NORM2 (L2 norm). For size-2, one should simply call HYPOT. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33197

[Bug fortran/33197] Fortran 2008: math functions

2010-08-21 Thread burnus at gcc dot gnu dot org
--- Comment #35 from burnus at gcc dot gnu dot org 2010-08-21 10:13 --- Subject: Bug 33197 Author: burnus Date: Sat Aug 21 10:12:53 2010 New Revision: 163440 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=163440 Log: 2010-08-21 Tobias Burnus PR fortran/36158

[Bug fortran/33197] Fortran 2008: math functions

2010-08-19 Thread burnus at gcc dot gnu dot org
--- Comment #34 from burnus at gcc dot gnu dot org 2010-08-19 07:28 --- Subject: Bug 33197 Author: burnus Date: Thu Aug 19 07:28:17 2010 New Revision: 163364 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=163364 Log: 2010-08-19 Tobias Burnus PR fortran/36158

[Bug fortran/33197] Fortran 2008: math functions

2010-05-02 Thread dfranke at gcc dot gnu dot org
--- Comment #33 from dfranke at gcc dot gnu dot org 2010-05-02 08:47 --- TODO (carried on): - BESSEL_JN and BESSEL_YN: Transitional form is missing (c.f. PR36096, PR36158) - NORM2 is missing -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33197

[Bug fortran/33197] Fortran 2008: math functions

2009-10-01 Thread ghazi at gcc dot gnu dot org
--- Comment #32 from ghazi at gcc dot gnu dot org 2009-10-02 03:52 --- Subject: Bug 33197 Author: ghazi Date: Fri Oct 2 03:52:05 2009 New Revision: 152394 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=152394 Log: PR fortran/33197 * gfortran.h (HAVE_mpc_arc): De

[Bug fortran/33197] Fortran 2008: math functions

2009-07-26 Thread burnus at gcc dot gnu dot org
--- Comment #31 from burnus at gcc dot gnu dot org 2009-07-26 17:26 --- Subject: Bug 33197 Author: burnus Date: Sun Jul 26 17:25:56 2009 New Revision: 150100 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=150100 Log: 2009-07-26 Tobias Burnus PR fortran/33197

[Bug fortran/33197] Fortran 2008: math functions

2009-07-25 Thread burnus at gcc dot gnu dot org
--- Comment #30 from burnus at gcc dot gnu dot org 2009-07-25 19:44 --- (Patch was approved by Steve on IRC.) TODO (carried on): - Support inverse functions as initialization expressions, when MPC is ready - TAN(X,Y) (= TAN2(X,Y)) - and fix wrong TAN2 description in intrinsic.texi - BE

[Bug fortran/33197] Fortran 2008: math functions

2009-07-25 Thread burnus at gcc dot gnu dot org
--- Comment #29 from burnus at gcc dot gnu dot org 2009-07-25 19:39 --- Subject: Bug 33197 Author: burnus Date: Sat Jul 25 19:39:07 2009 New Revision: 150087 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=150087 Log: 2009-07-25 Tobias Burnus Francois-Xavier Couder

[Bug fortran/33197] Fortran 2008: math functions

2009-07-14 Thread burnus at gcc dot gnu dot org
--- Comment #28 from burnus at gcc dot gnu dot org 2009-07-15 06:08 --- > In a quick scan of the patch, I note that the mpfr versions > of the simplifications weren't in the patch. MPFR or MPC? The MPFR version should be there since years; the MPC version of tan/sinh/cosh/tanh should be

[Bug fortran/33197] Fortran 2008: math functions

2009-07-14 Thread kargl at gcc dot gnu dot org
--- Comment #27 from kargl at gcc dot gnu dot org 2009-07-15 02:36 --- (In reply to comment #26) > Patch: http://gcc.gnu.org/ml/gcc-patches/2009-07/msg00739.html > - MPC compile-time evaluation of complex tan and sinh/cosh/tanh > - Run-time complex a(sin,cos,tan)(h) with C99 fallback (fo

[Bug fortran/33197] Fortran 2008: math functions

2009-07-14 Thread burnus at gcc dot gnu dot org
--- Comment #26 from burnus at gcc dot gnu dot org 2009-07-14 07:43 --- Patch: http://gcc.gnu.org/ml/gcc-patches/2009-07/msg00739.html - MPC compile-time evaluation of complex tan and sinh/cosh/tanh - Run-time complex a(sin,cos,tan)(h) with C99 fallback (for finite values only) TODO: -

[Bug fortran/33197] Fortran 2008: math functions

2009-07-10 Thread burnus at gcc dot gnu dot org
--- Comment #25 from burnus at gcc dot gnu dot org 2009-07-11 00:03 --- Subject: Bug 33197 Author: burnus Date: Sat Jul 11 00:03:07 2009 New Revision: 149503 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=149503 Log: 2009-07-09 Tobias Burnus PR fortran/33197

[Bug fortran/33197] Fortran 2008: math functions

2009-07-10 Thread burnus at gcc dot gnu dot org
--- Comment #24 from burnus at gcc dot gnu dot org 2009-07-10 16:17 --- Regarding expressions with Inf with x being finite according to Mathematica 7: ArcSin (+/-Inf,x) = (0, -/+Inf) ArcSin (x,+/-Inf) = (0, +/-Inf) ArcCos (+/-Inf,x) = (0, +/-Inf) ArcCos (x,+/-Inf) = (0, -/+Inf) ArcTan (+

[Bug fortran/33197] Fortran 2008: math functions

2009-07-10 Thread burnus at gcc dot gnu dot org
--- Comment #23 from burnus at gcc dot gnu dot org 2009-07-10 15:45 --- Created an attachment (id=18173) --> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=18173&action=view) Test implementation in C for complex arc{tan,sin,cos}{,h} > Fallback implementations via complex logarithm Att

[Bug fortran/33197] Fortran 2008: math functions

2009-07-10 Thread burnus at gcc dot gnu dot org
--- Comment #22 from burnus at gcc dot gnu dot org 2009-07-10 11:46 --- For z a complex number: Patch for tan(z), sinh(z), cosh(z), tanh(z) - see at http://gcc.gnu.org/ml/fortran/2009-07/msg00071.html Remark: atan(z), asin(z) and acos(z) are also missing besides a{sin,cos,tan}h. Additi

[Bug fortran/33197] Fortran 2008: math functions

2009-05-16 Thread fxcoudert at gcc dot gnu dot org
--- Comment #21 from fxcoudert at gcc dot gnu dot org 2009-05-16 18:09 --- ERFC_SCALED compile-time simplification was committed. TODO (carried on from above): - TAN(X,Y) (= TAN2(X,Y)) - BESSEL_JN and BESSEL_YN: Transitional form is missing - NORM2 - Finish complex (A)TAN(H), (A)COSH,

[Bug fortran/33197] Fortran 2008: math functions

2009-05-16 Thread fxcoudert at gcc dot gnu dot org
--- Comment #20 from fxcoudert at gcc dot gnu dot org 2009-05-16 17:33 --- Subject: Bug 33197 Author: fxcoudert Date: Sat May 16 17:33:23 2009 New Revision: 147621 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=147621 Log: PR fortran/33197 * intrinsic.c (add_fu

[Bug fortran/33197] Fortran 2008: math functions

2009-05-02 Thread fxcoudert at gcc dot gnu dot org
--- Comment #19 from fxcoudert at gcc dot gnu dot org 2009-05-02 13:46 --- Patch submitted for review for the ERFC_SCALED compile-time version here: http://gcc.gnu.org/ml/fortran/2009-05/msg00012.html -- fxcoudert at gcc dot gnu dot org changed: What|Removed

[Bug fortran/33197] Fortran 2008: math functions

2009-02-17 Thread burnus at gcc dot gnu dot org
--- Comment #18 from burnus at gcc dot gnu dot org 2009-02-17 20:58 --- Created an attachment (id=17315) --> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=17315&action=view) Complex (A)TAN(H), (A)COSH, (A)SINH, without no-C99 fallback & w/o tests "Implementation" of the run-time versi

[Bug fortran/33197] Fortran 2008: math functions

2008-06-15 Thread fxcoudert at gcc dot gnu dot org
-- fxcoudert at gcc dot gnu dot org changed: What|Removed |Added AssignedTo|fxcoudert at gcc dot gnu dot|unassigned at gcc dot gnu |org

[Bug fortran/33197] Fortran 2008: math functions

2008-03-11 Thread fxcoudert at gcc dot gnu dot org
--- Comment #17 from fxcoudert at gcc dot gnu dot org 2008-03-11 16:53 --- Adding a note: fallback formulas for systems that lack C99 complex functions can be found in Cephes (http://www.netlib.org/cephes/) -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33197

[Bug fortran/33197] Fortran 2008: math functions

2008-03-06 Thread burnus at gcc dot gnu dot org
--- Comment #16 from burnus at gcc dot gnu dot org 2008-03-06 13:50 --- (In reply to comment #15) > > - Implementation of NORM2 > There are quite a few other new F2008 intrinsics, including other > transformational intrinsics, so I think it's better to focus this PR on math > intrinsics.

[Bug fortran/33197] Fortran 2008: math functions

2008-03-06 Thread fxcoudert at gcc dot gnu dot org
--- Comment #15 from fxcoudert at gcc dot gnu dot org 2008-03-06 13:15 --- (In reply to comment #14) > I think also missing are: > - Compile-time evaluation of erfc_scaled Right. A though one. > - Implementation of NORM2 There are quite a few other new F2008 intrinsics, including oth

[Bug fortran/33197] Fortran 2008: math functions

2008-03-06 Thread burnus at gcc dot gnu dot org
--- Comment #14 from burnus at gcc dot gnu dot org 2008-03-06 13:10 --- > Remaining tasks in this PR: [...] I think also missing are: - Compile-time evaluation of erfc_scaled - Implementation of NORM2 -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33197

[Bug fortran/33197] Fortran 2008: math functions

2008-03-06 Thread fxcoudert at gcc dot gnu dot org
--- Comment #13 from fxcoudert at gcc dot gnu dot org 2008-03-06 12:45 --- Remaining tasks in this PR: taking care of math functions that now accept complex args (in particular, A{SIN,COS,TAN}{,H}), taking care of the two arguments form of ATAN, and deal with the 3 args form of BESSEL_J

[Bug fortran/33197] Fortran 2008: math functions and other small changes

2008-03-06 Thread fxcoudert at gcc dot gnu dot org
--- Comment #12 from fxcoudert at gcc dot gnu dot org 2008-03-06 12:41 --- Subject: Bug 33197 Author: fxcoudert Date: Thu Mar 6 12:40:28 2008 New Revision: 132970 URL: http://gcc.gnu.org/viewcvs?root=gcc&view=rev&rev=132970 Log: PR fortran/33197 * intrinsic.c (add_fun

[Bug fortran/33197] Fortran 2008: math functions and other small changes

2008-03-04 Thread fxcoudert at gcc dot gnu dot org
--- Comment #11 from fxcoudert at gcc dot gnu dot org 2008-03-04 14:26 --- On ppc longdouble, convergence requires argument >= 8.6 and, again, less than 100 terms. Why does convergence radius depend on the format of the floating-point type, I wonder... For the record, here's the small C

[Bug fortran/33197] Fortran 2008: math functions and other small changes

2008-03-04 Thread fxcoudert at gcc dot gnu dot org
--- Comment #10 from fxcoudert at gcc dot gnu dot org 2008-03-04 14:21 --- A bit of experimentation reveals that we can use the asymptotic expression at http://en.wikipedia.org/wiki/Error_function#Asymptotic_expansion for arguments larger than 7, and it does converge with less than 100

[Bug fortran/33197] Fortran 2008: math functions and other small changes

2008-03-03 Thread fxcoudert at gcc dot gnu dot org
--- Comment #9 from fxcoudert at gcc dot gnu dot org 2008-03-03 23:51 --- Tobias noted the following also need updating: 3999: /* Fortran 2008 draft allows BIND(C) for internal procedures. */ 4000- if (gfc_current_state () == COMP_CONTAINS 4001- && sym->ns->proc_name->at