很小数相除的浮点型精度问题，求指教

[iam1210900506]

编程时算法遇到问题，分析原因可能是由于出现类似 很小数/很大数 类型的算式，导致机器截断误差，精度不够导致。本人是摄影测量专业，精度要求就很高，求解决办法。发现double型不够用，有什么精度更高的浮点型？奇怪的类型请给出使用方法，谢谢

[derekrose]

用字符串？

[mujiok2003]

//因为浮点数在计算机中都误差，如果对精度要求很高，就不要用了。

#include <stdio.h>

typedef struct tagFloat
{
  int exponent; //指数
  int mantisa;  //尾数
}Float;


int main()
{
   
   double a = 0.0001, b = 200000, c = a / b;
   printf("%lf\n", c);
   Float x = {-4, 1}, y = {5, 2}, z = {-5-5, 10/2}; // z = x / y
   printf("%dE%d", z.mantisa, z.exponent);
   return 0;
}
/*
0.000000
5E-10
*/

[bravery36]

如果真的需要高精度的话，可以考虑下实现大数的加减乘除（假如说这样的运算就足够了的话），实现方法一般是用数组存储，具体的可以百度。不过这样做有很大局限性，最明显的问题就是只能做简单的运算，用不了其他库函数。

[赵4老师]

用10进制小数不能精确表示某些三进制小数0.1(3)=0.33333333333……(10) 同理，用二进制小数也不能精确表示某些10进制小数。 include\float.h

...
#define DBL_DIG         15                      /* # of decimal digits of precision */
#define DBL_EPSILON     2.2204460492503131e-016 /* smallest such that 1.0+DBL_EPSILON != 1.0 */
#define DBL_MANT_DIG    53                      /* # of bits in mantissa */
#define DBL_MAX         1.7976931348623158e+308 /* max value */
#define DBL_MAX_10_EXP  308                     /* max decimal exponent */
#define DBL_MAX_EXP     1024                    /* max binary exponent */
#define DBL_MIN         2.2250738585072014e-308 /* min positive value */
#define DBL_MIN_10_EXP  (-307)                  /* min decimal exponent */
#define DBL_MIN_EXP     (-1021)                 /* min binary exponent */
#define _DBL_RADIX      2                       /* exponent radix */
#define _DBL_ROUNDS     1                       /* addition rounding: near */
...

[dingqiang107]

C++标准有long double型，但不一定比double精度高，你可以先在你的系统中试一下。 x86架构上多数编译器支持80位的long double，包括Intel和GCC。但是VC++中long double和double是一样的。 GCC上有非标准的__float128。 如果不行，可以试试　GNU MP，任意精度整数、分数、浮点数　（http://gmplib.org/）。

[oN5GrzoN]

我在做地球椭球投影的时候(也是测量学的)，公司也是需要高精度，网上有很多成熟的大数类。 可惜我忘了那个有官网的大数类。只现成找了这个 http://bbs.csdn.net/topics/190173310

[caoh110]

一般double类型都够用了，我是搞科学计算的，我想你摄影专业的应该也够了