#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <stdint.h>
#include <limits.h>
#include <float.h>

#define _CRT_SECURE_NO_WARNINGS
#define LN10    2.302585092994045684018
#define RLN10  0.43429448190325182765

#define  MAX_LEN  128
......................

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>

#define MAX_Q  1000000000 
// 返回最大下标
int calc_cfrac(double x, int a[], int p[], int q[], int max_idx)
{
   double ipart;
   int i;
   
......................

#include <stdio.h>
#include <stdlib.h>
#include <float.h>
#include <stdint.h>
#include <math.h>
#include <mpfr.h>
#include <mpreal.h>
#include <time.h>
#include <inttypes.h>

#define LG2  0.30102999566398119521
#define RLN2 1.44269504088896340735
......................

#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>

typedef union
{
   double f;
   uint32_t i[2];
}Packet;

#define M  ((1<<20)-1)
void test( double f)
......................