为什么 boost 数值库使用调试或发布模式会给出不同的结果？

Question

我遇到过这样的现象，当我使用调试模式或发布模式时，我的代码给了我不同的结果。我已将问题简化为下面的代码。我正在使用 Microsoft Visual Studio Professional 2013 和 libeary boost 1.62

#include "stdafx.h"
#include <iostream>
#include <math.h>

#include <boost/numeric/interval.hpp>
#include <boost/numeric/interval/rounded_arith.hpp>

using namespace std;
using namespace boost::numeric::interval_lib;
using namespace boost::numeric;

typedef interval<double, policies<save_state<rounded_transc_std<double> >,
checking_base<double> > > Interval;

int _tmain(int argc, _TCHAR* argv[])
{   
    Interval result = (Interval(3.15, 4.6) - Interval(-0.6, 2.1))*sqrt(Interval(2, 2) + Interval(-2, -2)*Interval(10.022631612535406, 10.031726559552226));
    cout << "result: " << result.lower() << " " << result.upper();

    return 0;
}

在调试模式下的结果是 1.#QNAN 1.#QNAN

在释放模式下的结果是 0 0

我想知道导致此问题的原因以及如何解决此问题。因为如果我不能依赖结果，这会给我的项目带来严重的问题。

Answer 1

sqrt 的负数是一个艰难的命题。问题是Interval(-2, -2)。产生 0、0 仍然是 VisualStudio 的魔力。:)。 nan 是对sqrt(-x) 最合适的答案。你可以sqrt 或std::complex<T>。