1 files changed, 49 insertions, 0 deletions
diff --git a/reference/C/CONTRIB/SNIP/scaldate.c b/reference/C/CONTRIB/SNIP/scaldate.c
new file mode 100755
index 0000000..ca0f00d
--- /dev/null
+++ b/reference/C/CONTRIB/SNIP/scaldate.c
@@ -0,0 +1,49 @@
+/*
+** scalar date routines    --    public domain by Ray Gardner
+** These will work over the range 1/01/01 thru 14699/12/31
+*/
+
+#include "scaldate.h"
+
+int isleap (unsigned yr)
+{
+   return yr % 400 == 0 || (yr % 4 == 0 && yr % 100 != 0);
+}
+
+static unsigned months_to_days (unsigned month)
+{
+   return (month * 3057 - 3007) / 100;
+}
+
+static long years_to_days (unsigned yr)
+{
+   return yr * 365L + yr / 4 - yr / 100 + yr / 400;
+}
+
+long ymd_to_scalar (unsigned yr, unsigned mo, unsigned day)
+{
+   long scalar;
+   scalar = day + months_to_days(mo);
+   if ( mo > 2 )                         /* adjust if past February */
+      scalar -= isleap(yr) ? 1 : 2;
+   yr--;
+   scalar += years_to_days(yr);
+   return scalar;
+}
+
+void scalar_to_ymd (long scalar, unsigned *yr, unsigned *mo, unsigned *day)
+{
+   unsigned n;                /* compute inverse of years_to_days() */
+
+   for ( n = (unsigned)((scalar * 400L) / 146097); years_to_days(n) < scalar;)
+      n++;                          /* 146097 == years_to_days(400) */
+   *yr = n;
+   n = (unsigned)(scalar - years_to_days(n-1));
+   if ( n > 59 ) {                       /* adjust if past February */
+      n += 2;
+      if ( isleap(*yr) )
+         n -= n > 62 ? 1 : 2;
+   }
+   *mo = (n * 100 + 3007) / 3057;    /* inverse of months_to_days() */
+   *day = n - months_to_days(*mo);
+}