This is the mail archive of the
fortran@gcc.gnu.org
mailing list for the GNU Fortran project.
Re: [Patch, fortran] PR24518 and PR24520 - Improvements to MOD andDOT_PRODUCT
- From: Paul Thomas <paulthomas2 at wanadoo dot fr>
- To: Tobias SchlÃter <tobias dot schlueter at physik dot uni-muenchen dot de>
- Cc: "'fortran at gcc dot gnu dot org'" <fortran at gcc dot gnu dot org>,patch <gcc-patches at gcc dot gnu dot org>
- Date: Thu, 10 Nov 2005 20:46:17 +0100
- Subject: Re: [Patch, fortran] PR24518 and PR24520 - Improvements to MOD andDOT_PRODUCT
- References: <43730120.7030808@wanadoo.fr> <437396D2.30009@physik.uni-muenchen.de>
Tobi,
I'm wondering, why are you promoting the arguments in to double precision in
your mod implementation? This will introduce double rounding problems (which
I'm not too worried about, as Fortran is only guaranteeing an approximation to
the mathematical value) and cut off precision in the REAL*10 or REAL*16 case.
To increase the range for which an answer is obtained.
I'm also wondering, why you are unconditionally using the library version of
DOT_PROD for kinds > 8.
Uh! I put that in as a backstop because I was not in a position to test
the higher precision. I meant to ask that somebody do it for me by
raising the condition. I totally forgot about it.
Can somebody do the real*16 test for me, please? I'll check out real*10.
Thanks
Paul