1 /* atof_ieee.c - turn a Flonum into an IEEE floating point number
2    Copyright (C) 1987-2016 Free Software Foundation, Inc.
3 
4    This file is part of GAS, the GNU Assembler.
5 
6    GAS is free software; you can redistribute it and/or modify
7    it under the terms of the GNU General Public License as published by
8    the Free Software Foundation; either version 3, or (at your option)
9    any later version.
10 
11    GAS is distributed in the hope that it will be useful,
12    but WITHOUT ANY WARRANTY; without even the implied warranty of
13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14    GNU General Public License for more details.
15 
16    You should have received a copy of the GNU General Public License
17    along with GAS; see the file COPYING.  If not, write to the Free
18    Software Foundation, 51 Franklin Street - Fifth Floor, Boston, MA
19    02110-1301, USA.  */
20 
21 #include "as.h"
22 
23 /* Flonums returned here.  */
24 extern FLONUM_TYPE generic_floating_point_number;
25 
26 /* Precision in LittleNums.  */
27 /* Don't count the gap in the m68k extended precision format.  */
28 #define MAX_PRECISION  5
29 #define F_PRECISION    2
30 #define D_PRECISION    4
31 #define X_PRECISION    5
32 #define P_PRECISION    5
33 
34 /* Length in LittleNums of guard bits.  */
35 #define GUARD          2
36 
37 #ifndef TC_LARGEST_EXPONENT_IS_NORMAL
38 #define TC_LARGEST_EXPONENT_IS_NORMAL(PRECISION) 0
39 #endif
40 
41 static const unsigned long mask[] =
42 {
43   0x00000000,
44   0x00000001,
45   0x00000003,
46   0x00000007,
47   0x0000000f,
48   0x0000001f,
49   0x0000003f,
50   0x0000007f,
51   0x000000ff,
52   0x000001ff,
53   0x000003ff,
54   0x000007ff,
55   0x00000fff,
56   0x00001fff,
57   0x00003fff,
58   0x00007fff,
59   0x0000ffff,
60   0x0001ffff,
61   0x0003ffff,
62   0x0007ffff,
63   0x000fffff,
64   0x001fffff,
65   0x003fffff,
66   0x007fffff,
67   0x00ffffff,
68   0x01ffffff,
69   0x03ffffff,
70   0x07ffffff,
71   0x0fffffff,
72   0x1fffffff,
73   0x3fffffff,
74   0x7fffffff,
75   0xffffffff,
76 };
77 
78 static int bits_left_in_littlenum;
79 static int littlenums_left;
80 static LITTLENUM_TYPE *littlenum_pointer;
81 
82 static int
83 next_bits (int number_of_bits)
84 {
85   int return_value;
86 
87   if (!littlenums_left)
88     return 0;
89 
90   if (number_of_bits >= bits_left_in_littlenum)
91     {
92       return_value = mask[bits_left_in_littlenum] & *littlenum_pointer;
93       number_of_bits -= bits_left_in_littlenum;
94       return_value <<= number_of_bits;
95 
96       if (--littlenums_left)
97 	{
98 	  bits_left_in_littlenum = LITTLENUM_NUMBER_OF_BITS - number_of_bits;
99 	  --littlenum_pointer;
100 	  return_value |=
101 	    (*littlenum_pointer >> bits_left_in_littlenum)
102 	    & mask[number_of_bits];
103 	}
104     }
105   else
106     {
107       bits_left_in_littlenum -= number_of_bits;
108       return_value =
109 	mask[number_of_bits] & (*littlenum_pointer >> bits_left_in_littlenum);
110     }
111   return return_value;
112 }
113 
114 /* Num had better be less than LITTLENUM_NUMBER_OF_BITS.  */
115 
116 static void
117 unget_bits (int num)
118 {
119   if (!littlenums_left)
120     {
121       ++littlenum_pointer;
122       ++littlenums_left;
123       bits_left_in_littlenum = num;
124     }
125   else if (bits_left_in_littlenum + num > LITTLENUM_NUMBER_OF_BITS)
126     {
127       bits_left_in_littlenum =
128 	num - (LITTLENUM_NUMBER_OF_BITS - bits_left_in_littlenum);
129       ++littlenum_pointer;
130       ++littlenums_left;
131     }
132   else
133     bits_left_in_littlenum += num;
134 }
135 
136 static void
137 make_invalid_floating_point_number (LITTLENUM_TYPE *words)
138 {
139   as_bad (_("cannot create floating-point number"));
140   /* Zero the leftmost bit.  */
141   words[0] = (LITTLENUM_TYPE) ((unsigned) -1) >> 1;
142   words[1] = (LITTLENUM_TYPE) -1;
143   words[2] = (LITTLENUM_TYPE) -1;
144   words[3] = (LITTLENUM_TYPE) -1;
145   words[4] = (LITTLENUM_TYPE) -1;
146   words[5] = (LITTLENUM_TYPE) -1;
147 }
148 
149 /* Warning: This returns 16-bit LITTLENUMs.  It is up to the caller to
150    figure out any alignment problems and to conspire for the
151    bytes/word to be emitted in the right order.  Bigendians beware!  */
152 
153 /* Note that atof-ieee always has X and P precisions enabled.  it is up
154    to md_atof to filter them out if the target machine does not support
155    them.  */
156 
157 /* Returns pointer past text consumed.  */
158 
159 char *
160 atof_ieee (char *str,			/* Text to convert to binary.  */
161 	   int what_kind,		/* 'd', 'f', 'x', 'p'.  */
162 	   LITTLENUM_TYPE *words)	/* Build the binary here.  */
163 {
164   /* Extra bits for zeroed low-order bits.
165      The 1st MAX_PRECISION are zeroed, the last contain flonum bits.  */
166   static LITTLENUM_TYPE bits[MAX_PRECISION + MAX_PRECISION + GUARD];
167   char *return_value;
168   /* Number of 16-bit words in the format.  */
169   int precision;
170   long exponent_bits;
171   FLONUM_TYPE save_gen_flonum;
172 
173   /* We have to save the generic_floating_point_number because it
174      contains storage allocation about the array of LITTLENUMs where
175      the value is actually stored.  We will allocate our own array of
176      littlenums below, but have to restore the global one on exit.  */
177   save_gen_flonum = generic_floating_point_number;
178 
179   return_value = str;
180   generic_floating_point_number.low = bits + MAX_PRECISION;
181   generic_floating_point_number.high = NULL;
182   generic_floating_point_number.leader = NULL;
183   generic_floating_point_number.exponent = 0;
184   generic_floating_point_number.sign = '\0';
185 
186   /* Use more LittleNums than seems necessary: the highest flonum may
187      have 15 leading 0 bits, so could be useless.  */
188 
189   memset (bits, '\0', sizeof (LITTLENUM_TYPE) * MAX_PRECISION);
190 
191   switch (what_kind)
192     {
193     case 'f':
194     case 'F':
195     case 's':
196     case 'S':
197       precision = F_PRECISION;
198       exponent_bits = 8;
199       break;
200 
201     case 'd':
202     case 'D':
203     case 'r':
204     case 'R':
205       precision = D_PRECISION;
206       exponent_bits = 11;
207       break;
208 
209     case 'x':
210     case 'X':
211     case 'e':
212     case 'E':
213       precision = X_PRECISION;
214       exponent_bits = 15;
215       break;
216 
217     case 'p':
218     case 'P':
219       precision = P_PRECISION;
220       exponent_bits = -1;
221       break;
222 
223     default:
224       make_invalid_floating_point_number (words);
225       return (NULL);
226     }
227 
228   generic_floating_point_number.high
229     = generic_floating_point_number.low + precision - 1 + GUARD;
230 
231   if (atof_generic (&return_value, ".", EXP_CHARS,
232 		    &generic_floating_point_number))
233     {
234       make_invalid_floating_point_number (words);
235       return NULL;
236     }
237   gen_to_words (words, precision, exponent_bits);
238 
239   /* Restore the generic_floating_point_number's storage alloc (and
240      everything else).  */
241   generic_floating_point_number = save_gen_flonum;
242 
243   return return_value;
244 }
245 
246 /* Turn generic_floating_point_number into a real float/double/extended.  */
247 
248 int
249 gen_to_words (LITTLENUM_TYPE *words, int precision, long exponent_bits)
250 {
251   int return_value = 0;
252 
253   long exponent_1;
254   long exponent_2;
255   long exponent_3;
256   long exponent_4;
257   int exponent_skippage;
258   LITTLENUM_TYPE word1;
259   LITTLENUM_TYPE *lp;
260   LITTLENUM_TYPE *words_end;
261 
262   words_end = words + precision;
263 #ifdef TC_M68K
264   if (precision == X_PRECISION)
265     /* On the m68k the extended precision format has a gap of 16 bits
266        between the exponent and the mantissa.  */
267     words_end++;
268 #endif
269 
270   if (generic_floating_point_number.low > generic_floating_point_number.leader)
271     {
272       /* 0.0e0 seen.  */
273       if (generic_floating_point_number.sign == '+')
274 	words[0] = 0x0000;
275       else
276 	words[0] = 0x8000;
277       memset (&words[1], '\0',
278 	      (words_end - words - 1) * sizeof (LITTLENUM_TYPE));
279       return return_value;
280     }
281 
282   /* NaN:  Do the right thing.  */
283   if (generic_floating_point_number.sign == 0)
284     {
285       if (TC_LARGEST_EXPONENT_IS_NORMAL (precision))
286 	as_warn (_("NaNs are not supported by this target\n"));
287       if (precision == F_PRECISION)
288 	{
289 	  words[0] = 0x7fff;
290 	  words[1] = 0xffff;
291 	}
292       else if (precision == X_PRECISION)
293 	{
294 #ifdef TC_M68K
295 	  words[0] = 0x7fff;
296 	  words[1] = 0;
297 	  words[2] = 0xffff;
298 	  words[3] = 0xffff;
299 	  words[4] = 0xffff;
300 	  words[5] = 0xffff;
301 #else /* ! TC_M68K  */
302 #ifdef TC_I386
303 	  words[0] = 0xffff;
304 	  words[1] = 0xc000;
305 	  words[2] = 0;
306 	  words[3] = 0;
307 	  words[4] = 0;
308 #else /* ! TC_I386  */
309 	  abort ();
310 #endif /* ! TC_I386  */
311 #endif /* ! TC_M68K  */
312 	}
313       else
314 	{
315 	  words[0] = 0x7fff;
316 	  words[1] = 0xffff;
317 	  words[2] = 0xffff;
318 	  words[3] = 0xffff;
319 	}
320       return return_value;
321     }
322   else if (generic_floating_point_number.sign == 'P')
323     {
324       if (TC_LARGEST_EXPONENT_IS_NORMAL (precision))
325 	as_warn (_("Infinities are not supported by this target\n"));
326 
327       /* +INF:  Do the right thing.  */
328       if (precision == F_PRECISION)
329 	{
330 	  words[0] = 0x7f80;
331 	  words[1] = 0;
332 	}
333       else if (precision == X_PRECISION)
334 	{
335 #ifdef TC_M68K
336 	  words[0] = 0x7fff;
337 	  words[1] = 0;
338 	  words[2] = 0;
339 	  words[3] = 0;
340 	  words[4] = 0;
341 	  words[5] = 0;
342 #else /* ! TC_M68K  */
343 #ifdef TC_I386
344 	  words[0] = 0x7fff;
345 	  words[1] = 0x8000;
346 	  words[2] = 0;
347 	  words[3] = 0;
348 	  words[4] = 0;
349 #else /* ! TC_I386  */
350 	  abort ();
351 #endif /* ! TC_I386  */
352 #endif /* ! TC_M68K  */
353 	}
354       else
355 	{
356 	  words[0] = 0x7ff0;
357 	  words[1] = 0;
358 	  words[2] = 0;
359 	  words[3] = 0;
360 	}
361       return return_value;
362     }
363   else if (generic_floating_point_number.sign == 'N')
364     {
365       if (TC_LARGEST_EXPONENT_IS_NORMAL (precision))
366 	as_warn (_("Infinities are not supported by this target\n"));
367 
368       /* Negative INF.  */
369       if (precision == F_PRECISION)
370 	{
371 	  words[0] = 0xff80;
372 	  words[1] = 0x0;
373 	}
374       else if (precision == X_PRECISION)
375 	{
376 #ifdef TC_M68K
377 	  words[0] = 0xffff;
378 	  words[1] = 0;
379 	  words[2] = 0;
380 	  words[3] = 0;
381 	  words[4] = 0;
382 	  words[5] = 0;
383 #else /* ! TC_M68K  */
384 #ifdef TC_I386
385 	  words[0] = 0xffff;
386 	  words[1] = 0x8000;
387 	  words[2] = 0;
388 	  words[3] = 0;
389 	  words[4] = 0;
390 #else /* ! TC_I386  */
391 	  abort ();
392 #endif /* ! TC_I386  */
393 #endif /* ! TC_M68K  */
394 	}
395       else
396 	{
397 	  words[0] = 0xfff0;
398 	  words[1] = 0x0;
399 	  words[2] = 0x0;
400 	  words[3] = 0x0;
401 	}
402       return return_value;
403     }
404 
405   /* The floating point formats we support have:
406      Bit 15 is sign bit.
407      Bits 14:n are excess-whatever exponent.
408      Bits n-1:0 (if any) are most significant bits of fraction.
409      Bits 15:0 of the next word(s) are the next most significant bits.
410 
411      So we need: number of bits of exponent, number of bits of
412      mantissa.  */
413   bits_left_in_littlenum = LITTLENUM_NUMBER_OF_BITS;
414   littlenum_pointer = generic_floating_point_number.leader;
415   littlenums_left = (1
416 		     + generic_floating_point_number.leader
417 		     - generic_floating_point_number.low);
418 
419   /* Seek (and forget) 1st significant bit.  */
420   for (exponent_skippage = 0; !next_bits (1); ++exponent_skippage);
421   exponent_1 = (generic_floating_point_number.exponent
422 		+ generic_floating_point_number.leader
423 		+ 1
424 		- generic_floating_point_number.low);
425 
426   /* Radix LITTLENUM_RADIX, point just higher than
427      generic_floating_point_number.leader.  */
428   exponent_2 = exponent_1 * LITTLENUM_NUMBER_OF_BITS;
429 
430   /* Radix 2.  */
431   exponent_3 = exponent_2 - exponent_skippage;
432 
433   /* Forget leading zeros, forget 1st bit.  */
434   exponent_4 = exponent_3 + ((1 << (exponent_bits - 1)) - 2);
435 
436   /* Offset exponent.  */
437   lp = words;
438 
439   /* Word 1.  Sign, exponent and perhaps high bits.  */
440   word1 = ((generic_floating_point_number.sign == '+')
441 	   ? 0
442 	   : (1 << (LITTLENUM_NUMBER_OF_BITS - 1)));
443 
444   /* Assume 2's complement integers.  */
445   if (exponent_4 <= 0)
446     {
447       int prec_bits;
448       int num_bits;
449 
450       unget_bits (1);
451       num_bits = -exponent_4;
452       prec_bits =
453 	LITTLENUM_NUMBER_OF_BITS * precision - (exponent_bits + 1 + num_bits);
454 #ifdef TC_I386
455       if (precision == X_PRECISION && exponent_bits == 15)
456 	{
457 	  /* On the i386 a denormalized extended precision float is
458 	     shifted down by one, effectively decreasing the exponent
459 	     bias by one.  */
460 	  prec_bits -= 1;
461 	  num_bits += 1;
462 	}
463 #endif
464 
465       if (num_bits >= LITTLENUM_NUMBER_OF_BITS - exponent_bits)
466 	{
467 	  /* Bigger than one littlenum.  */
468 	  num_bits -= (LITTLENUM_NUMBER_OF_BITS - 1) - exponent_bits;
469 	  *lp++ = word1;
470 	  if (num_bits + exponent_bits + 1
471 	      > precision * LITTLENUM_NUMBER_OF_BITS)
472 	    {
473 	      /* Exponent overflow.  */
474 	      make_invalid_floating_point_number (words);
475 	      return return_value;
476 	    }
477 #ifdef TC_M68K
478 	  if (precision == X_PRECISION && exponent_bits == 15)
479 	    *lp++ = 0;
480 #endif
481 	  while (num_bits >= LITTLENUM_NUMBER_OF_BITS)
482 	    {
483 	      num_bits -= LITTLENUM_NUMBER_OF_BITS;
484 	      *lp++ = 0;
485 	    }
486 	  if (num_bits)
487 	    *lp++ = next_bits (LITTLENUM_NUMBER_OF_BITS - (num_bits));
488 	}
489       else
490 	{
491 	  if (precision == X_PRECISION && exponent_bits == 15)
492 	    {
493 	      *lp++ = word1;
494 #ifdef TC_M68K
495 	      *lp++ = 0;
496 #endif
497 	      *lp++ = next_bits (LITTLENUM_NUMBER_OF_BITS - num_bits);
498 	    }
499 	  else
500 	    {
501 	      word1 |= next_bits ((LITTLENUM_NUMBER_OF_BITS - 1)
502 				  - (exponent_bits + num_bits));
503 	      *lp++ = word1;
504 	    }
505 	}
506       while (lp < words_end)
507 	*lp++ = next_bits (LITTLENUM_NUMBER_OF_BITS);
508 
509       /* Round the mantissa up, but don't change the number.  */
510       if (next_bits (1))
511 	{
512 	  --lp;
513 	  if (prec_bits >= LITTLENUM_NUMBER_OF_BITS)
514 	    {
515 	      int n = 0;
516 	      int tmp_bits;
517 
518 	      n = 0;
519 	      tmp_bits = prec_bits;
520 	      while (tmp_bits > LITTLENUM_NUMBER_OF_BITS)
521 		{
522 		  if (lp[n] != (LITTLENUM_TYPE) - 1)
523 		    break;
524 		  --n;
525 		  tmp_bits -= LITTLENUM_NUMBER_OF_BITS;
526 		}
527 	      if (tmp_bits > LITTLENUM_NUMBER_OF_BITS
528 		  || (lp[n] & mask[tmp_bits]) != mask[tmp_bits]
529 		  || (prec_bits != (precision * LITTLENUM_NUMBER_OF_BITS
530 				    - exponent_bits - 1)
531 #ifdef TC_I386
532 		      /* An extended precision float with only the integer
533 			 bit set would be invalid.  That must be converted
534 			 to the smallest normalized number.  */
535 		      && !(precision == X_PRECISION
536 			   && prec_bits == (precision * LITTLENUM_NUMBER_OF_BITS
537 					    - exponent_bits - 2))
538 #endif
539 		      ))
540 		{
541 		  unsigned long carry;
542 
543 		  for (carry = 1; carry && (lp >= words); lp--)
544 		    {
545 		      carry = *lp + carry;
546 		      *lp = carry;
547 		      carry >>= LITTLENUM_NUMBER_OF_BITS;
548 		    }
549 		}
550 	      else
551 		{
552 		  /* This is an overflow of the denormal numbers.  We
553                      need to forget what we have produced, and instead
554                      generate the smallest normalized number.  */
555 		  lp = words;
556 		  word1 = ((generic_floating_point_number.sign == '+')
557 			   ? 0
558 			   : (1 << (LITTLENUM_NUMBER_OF_BITS - 1)));
559 		  word1 |= (1
560 			    << ((LITTLENUM_NUMBER_OF_BITS - 1)
561 				- exponent_bits));
562 		  *lp++ = word1;
563 #ifdef TC_I386
564 		  /* Set the integer bit in the extended precision format.
565 		     This cannot happen on the m68k where the mantissa
566 		     just overflows into the integer bit above.  */
567 		  if (precision == X_PRECISION)
568 		    *lp++ = 1 << (LITTLENUM_NUMBER_OF_BITS - 1);
569 #endif
570 		  while (lp < words_end)
571 		    *lp++ = 0;
572 		}
573 	    }
574 	  else
575 	    *lp += 1;
576 	}
577 
578       return return_value;
579     }
580   else if ((unsigned long) exponent_4 > mask[exponent_bits]
581 	   || (! TC_LARGEST_EXPONENT_IS_NORMAL (precision)
582 	       && (unsigned long) exponent_4 == mask[exponent_bits]))
583     {
584       /* Exponent overflow.  Lose immediately.  */
585 
586       /* We leave return_value alone: admit we read the
587 	 number, but return a floating exception
588 	 because we can't encode the number.  */
589       make_invalid_floating_point_number (words);
590       return return_value;
591     }
592   else
593     {
594       word1 |= (exponent_4 << ((LITTLENUM_NUMBER_OF_BITS - 1) - exponent_bits))
595 	| next_bits ((LITTLENUM_NUMBER_OF_BITS - 1) - exponent_bits);
596     }
597 
598   *lp++ = word1;
599 
600   /* X_PRECISION is special: on the 68k, it has 16 bits of zero in the
601      middle.  Either way, it is then followed by a 1 bit.  */
602   if (exponent_bits == 15 && precision == X_PRECISION)
603     {
604 #ifdef TC_M68K
605       *lp++ = 0;
606 #endif
607       *lp++ = (1 << (LITTLENUM_NUMBER_OF_BITS - 1)
608 	       | next_bits (LITTLENUM_NUMBER_OF_BITS - 1));
609     }
610 
611   /* The rest of the words are just mantissa bits.  */
612   while (lp < words_end)
613     *lp++ = next_bits (LITTLENUM_NUMBER_OF_BITS);
614 
615   if (next_bits (1))
616     {
617       unsigned long carry;
618       /* Since the NEXT bit is a 1, round UP the mantissa.
619 	 The cunning design of these hidden-1 floats permits
620 	 us to let the mantissa overflow into the exponent, and
621 	 it 'does the right thing'. However, we lose if the
622 	 highest-order bit of the lowest-order word flips.
623 	 Is that clear?  */
624 
625       /* #if (sizeof(carry)) < ((sizeof(bits[0]) * BITS_PER_CHAR) + 2)
626 	 Please allow at least 1 more bit in carry than is in a LITTLENUM.
627 	 We need that extra bit to hold a carry during a LITTLENUM carry
628 	 propagation. Another extra bit (kept 0) will assure us that we
629 	 don't get a sticky sign bit after shifting right, and that
630 	 permits us to propagate the carry without any masking of bits.
631 	 #endif */
632       for (carry = 1, lp--; carry; lp--)
633 	{
634 	  carry = *lp + carry;
635 	  *lp = carry;
636 	  carry >>= LITTLENUM_NUMBER_OF_BITS;
637 	  if (lp == words)
638 	    break;
639 	}
640       if (precision == X_PRECISION && exponent_bits == 15)
641 	{
642 	  /* Extended precision numbers have an explicit integer bit
643 	     that we may have to restore.  */
644 	  if (lp == words)
645 	    {
646 #ifdef TC_M68K
647 	      /* On the m68k there is a gap of 16 bits.  We must
648 		 explicitly propagate the carry into the exponent.  */
649 	      words[0] += words[1];
650 	      words[1] = 0;
651 	      lp++;
652 #endif
653 	      /* Put back the integer bit.  */
654 	      lp[1] |= 1 << (LITTLENUM_NUMBER_OF_BITS - 1);
655 	    }
656 	}
657       if ((word1 ^ *words) & (1 << (LITTLENUM_NUMBER_OF_BITS - 1)))
658 	{
659 	  /* We leave return_value alone: admit we read the number,
660 	     but return a floating exception because we can't encode
661 	     the number.  */
662 	  *words &= ~(1 << (LITTLENUM_NUMBER_OF_BITS - 1));
663 	}
664     }
665   return return_value;
666 }
667 
668 #ifdef TEST
669 char *
670 print_gen (gen)
671      FLONUM_TYPE *gen;
672 {
673   FLONUM_TYPE f;
674   LITTLENUM_TYPE arr[10];
675   double dv;
676   float fv;
677   static char sbuf[40];
678 
679   if (gen)
680     {
681       f = generic_floating_point_number;
682       generic_floating_point_number = *gen;
683     }
684   gen_to_words (&arr[0], 4, 11);
685   memcpy (&dv, &arr[0], sizeof (double));
686   sprintf (sbuf, "%x %x %x %x %.14G   ", arr[0], arr[1], arr[2], arr[3], dv);
687   gen_to_words (&arr[0], 2, 8);
688   memcpy (&fv, &arr[0], sizeof (float));
689   sprintf (sbuf + strlen (sbuf), "%x %x %.12g\n", arr[0], arr[1], fv);
690 
691   if (gen)
692     generic_floating_point_number = f;
693 
694   return (sbuf);
695 }
696 #endif
697 
698 #define MAX_LITTLENUMS 6
699 
700 /* This is a utility function called from various tc-*.c files.  It
701    is here in order to reduce code duplication.
702 
703    Turn a string at input_line_pointer into a floating point constant
704    of type TYPE (a character found in the FLT_CHARS macro), and store
705    it as LITTLENUMS in the bytes buffer LITP.  The number of chars
706    emitted is stored in *SIZEP.  BIG_WORDIAN is TRUE if the littlenums
707    should be emitted most significant littlenum first.
708 
709    An error message is returned, or a NULL pointer if everything went OK.  */
710 
711 const char *
712 ieee_md_atof (int type,
713 	      char *litP,
714 	      int *sizeP,
715 	      bfd_boolean big_wordian)
716 {
717   LITTLENUM_TYPE words[MAX_LITTLENUMS];
718   LITTLENUM_TYPE *wordP;
719   char *t;
720   int prec = 0;
721 
722   if (strchr (FLT_CHARS, type) != NULL)
723     {
724       switch (type)
725 	{
726 	case 'f':
727 	case 'F':
728 	case 's':
729 	case 'S':
730 	  prec = F_PRECISION;
731 	  break;
732 
733 	case 'd':
734 	case 'D':
735 	case 'r':
736 	case 'R':
737 	  prec = D_PRECISION;
738 	  break;
739 
740 	case 't':
741 	case 'T':
742 	  prec = X_PRECISION;
743 	  type = 'x';		/* This is what atof_ieee() understands.  */
744 	  break;
745 
746 	case 'x':
747 	case 'X':
748 	case 'p':
749 	case 'P':
750 #ifdef TC_M68K
751 	  /* Note: on the m68k there is a gap of 16 bits (one littlenum)
752 	     between the exponent and mantissa.  Hence the precision is
753 	     6 and not 5.  */
754 	  prec = P_PRECISION + 1;
755 #else
756 	  prec = P_PRECISION;
757 #endif
758 	  break;
759 
760 	default:
761 	  break;
762 	}
763     }
764   /* The 'f' and 'd' types are always recognised, even if the target has
765      not put them into the FLT_CHARS macro.  This is because the 'f' type
766      can come from the .dc.s, .dcb.s, .float or .single pseudo-ops and the
767      'd' type from the .dc.d, .dbc.d or .double pseudo-ops.
768 
769      The 'x' type is not implicitly recongised however, even though it can
770      be generated by the .dc.x and .dbc.x pseudo-ops because not all targets
771      can support floating point values that big.  ie the target has to
772      explicitly allow them by putting them into FLT_CHARS.  */
773   else if (type == 'f')
774     prec = F_PRECISION;
775   else if (type == 'd')
776     prec = D_PRECISION;
777 
778   if (prec == 0)
779     {
780       *sizeP = 0;
781       return _("Unrecognized or unsupported floating point constant");
782     }
783 
784   gas_assert (prec <= MAX_LITTLENUMS);
785 
786   t = atof_ieee (input_line_pointer, type, words);
787   if (t)
788     input_line_pointer = t;
789 
790   *sizeP = prec * sizeof (LITTLENUM_TYPE);
791 
792   if (big_wordian)
793     {
794       for (wordP = words; prec --;)
795 	{
796 	  md_number_to_chars (litP, (valueT) (* wordP ++), sizeof (LITTLENUM_TYPE));
797 	  litP += sizeof (LITTLENUM_TYPE);
798 	}
799     }
800   else
801     {
802       for (wordP = words + prec; prec --;)
803 	{
804 	  md_number_to_chars (litP, (valueT) (* -- wordP), sizeof (LITTLENUM_TYPE));
805 	  litP += sizeof (LITTLENUM_TYPE);
806 	}
807     }
808 
809   return NULL;
810 }
811