/* Standard C source for D1MACH -- remove the * in column 1 */
#include <stdio.h>
#include <float.h>
#include <math.h>

#if    LINKTYP == 1

#define d1mach d1mach_

#endif

#if  LINKTYP == 2

#define d1mach D1MACH

#endif

double d1mach(int *i)
{
	switch(*i){
	  case 1: return DBL_MIN;
	  case 2: return DBL_MAX;
	  case 3: return DBL_EPSILON/FLT_RADIX;
	  case 4: return DBL_EPSILON;
	  case 5: return log10((double)FLT_RADIX);
	  }
	fprintf(stderr, "invalid argument: d1mach(%d)\n", *i);
	exit(1); 
	return 0; /* some compilers demand return values */
}