This is the mail archive of the
fortran@gcc.gnu.org
mailing list for the GNU Fortran project.
Re: [PATCH]: Integrate gfortran with MPC
- From: Thomas Koenig <tkoenig at netcologne dot de>
- To: Dominique Dhumieres <dominiq at lps dot ens dot fr>
- Cc: fortran at gcc dot gnu dot org
- Date: Mon, 01 Jun 2009 00:00:51 +0200
- Subject: Re: [PATCH]: Integrate gfortran with MPC
- References: <20090531215227.4FA533BAB9@mailhost.lps.ens.fr>
On Sun, 2009-05-31 at 23:52 +0200, Dominique Dhumieres wrote:
> Any better ideas?
Maybe not better, but...
IMHO, we should be looking at limits.
a/(c+i*d) expands to a*c/(c*c+d*d) - a*d/(c*c+d*d). For
c=0 and d=0, this is 0/0, or NaN + i * NaN.
Thomas