# HG changeset patch
# User jwe
# Date 1109790854 0
# Node ID 10172b51969552ddd96d02e7b2eea2662bedddff
# Parent  7a25bb94ef65a65dc9a3fb578264a774904ac26b
[project @ 2005-03-02 19:14:13 by jwe]

diff -r 7a25bb94ef65 -r 10172b519695 src/ChangeLog
--- a/src/ChangeLog	Wed Mar 02 17:38:32 2005 +0000
+++ b/src/ChangeLog	Wed Mar 02 19:14:14 2005 +0000
@@ -1,5 +1,9 @@
 2005-03-02  John W. Eaton  <jwe@octave.org>
 
+	* oct-stream.cc (octave_scan (std::istream&, const
+	scanf_format_elt&, double*)): New specialization to handle Inf,
+	NaN, and NA.
+
 	* parse.y (looks_like_copyright):
 	Rename from looks_like_octave_copyright.  Change all uses.
 	Simply match "Copyright".
diff -r 7a25bb94ef65 -r 10172b519695 src/oct-stream.cc
--- a/src/oct-stream.cc	Wed Mar 02 17:38:32 2005 +0000
+++ b/src/oct-stream.cc	Wed Mar 02 19:14:14 2005 +0000
@@ -1142,6 +1142,181 @@
 octave_scan (std::istream&, const scanf_format_elt&, float*);
 #endif
 
+template <>
+std::istream&
+octave_scan (std::istream& is, const scanf_format_elt& fmt, double* valptr)
+{
+  double& ref = *valptr;
+
+  switch (fmt.type)
+    {
+    case 'e':
+    case 'f':
+    case 'g':
+      {
+	int c1;
+
+	while (is && (c1 = is.get ()) != EOF && isspace (c1))
+	  /* skip whitespace */;
+
+	if (c1 != EOF)
+	  {
+	    if (c1 == 'N')
+	      {
+		int c2 = is.get ();
+
+		if (c2 != EOF)
+		  {
+		    if (c2 == 'A')
+		      {
+			int c3 = is.get ();
+
+			if (c3 != EOF)
+			  {
+			    is.putback (c3);
+
+			    if (isspace (c3) || ispunct (c3))
+			      ref = octave_NA;
+			    else
+			      {
+				is.putback (c2);
+				is.putback (c1);
+
+				is >> ref;
+			      }
+			  }
+			else
+			  {
+			    is.clear ();
+
+			    ref = octave_NA;
+			  }
+		      }
+		    else if (c2 == 'a')
+		      {
+			int c3 = is.get ();
+
+			if (c3 != EOF)
+			  {
+			    if (c3 == 'N')
+			      {
+				int c4 = is.get ();
+
+				if (c4 != EOF)
+				  {
+				    is.putback (c4);
+
+				    if (isspace (c4) || ispunct (c4))
+				      ref = octave_NaN;
+				    else
+				      {
+					is.putback (c3);
+					is.putback (c2);
+					is.putback (c1);
+
+					is >> ref;
+				      }
+				  }
+				else
+				  {
+				    is.clear ();
+
+				    ref = octave_NaN;
+				  }
+			      }
+			    else
+			      {
+				is.putback (c3);
+				is.putback (c2);
+				is.putback (c1);
+
+				is >> ref;
+			      }
+			  }
+		      }
+		    else
+		      {
+			is.putback (c2);
+			is.putback (c1);
+
+			is >> ref;
+		      }
+		  }
+	      }
+	    else if (c1 == 'I')
+	      {
+		int c2 = is.get ();
+
+		if (c2 != EOF)
+		  {
+		    if (c2 == 'n')
+		      {
+			int c3 = is.get ();
+
+			if (c3 != EOF)
+
+			  if (c3 == 'f')
+			    {
+			      int c4 = is.get ();
+
+			      if (c4 != EOF)
+				{
+				  is.putback (c4);
+
+				  if (isspace (c4) || ispunct (c4))
+				    ref = octave_Inf;
+				  else
+				    {
+				      is.putback (c3);
+				      is.putback (c2);
+				      is.putback (c1);
+
+				      is >> ref;
+				    }
+				}
+			      else
+				{
+				  is.clear ();
+
+				  ref = octave_Inf;
+				}
+			    }
+			  else
+			    {
+			      is.putback (c3);
+			      is.putback (c2);
+			      is.putback (c1);
+
+			      is >> ref;
+			    }
+		      }
+		    else
+		      {
+			is.putback (c2);
+			is.putback (c1);
+
+			is >> ref;
+		      }
+		  }
+	      }
+	    else
+	      {
+		is.putback (c1);
+
+		is >> ref;
+	      }
+	  }
+      }
+      break;
+
+    default:
+      panic_impossible ();
+      break;
+    }
+
+  return is;
+}
+
 template std::istream&
 octave_scan (std::istream&, const scanf_format_elt&, double*);