主页 / user-217660

AVK's questions

Martin Hope
AVK
Asked: 2022-03-17 18:42:18 +0000 UTC
  • 0

我需要具有任意精度的计算，因此我决定使用 Boost Multiprecision。我下载了 Boost 1.77 并尝试使用 mingw 编译这个程序：

#include <iostream>
#include <boost/multiprecision/cpp_bin_float.hpp>

using namespace std;

int main()
{
    cout << "Hello world!" << endl;
    return 0;
}

这个程序没有编译。编译器报错：

C:\include\boost\math\special_functions\detail\bernoulli_details.hpp|559|错误：命名空间'std'中的'mutex'没有命名类型| C:\include\boost\math\special_functions\detail\bernoulli_details.hpp|21|注意：'std::mutex' 定义在头文件''中；你忘了'#include mutex'吗？

等等。看着 bernoulli_details.hpp 我看到了这个：

#ifdef BOOST_HAS_THREADS
std::mutex m_mutex;
atomic_counter_type m_counter, m_current_precision;
#else
int m_counter;
int m_current_precision;
#endif // BOOST_HAS_THREADS

有点惊讶的是，理论上应该计算伯努利数的文件说这个，我决定我需要找到配置文件并删除#define BOOST_HAS_THREADS. 我不需要多线程。顺便说一句，我也不需要伯努利数。正如您可能猜到的那样，我失败了……BOOST_HAS_THREADS定义了 18 个文件。其中一些是平台配置文件，其余的是来自各种 Boost 数学库的文件，其中该宏被无条件地重新定义。修改所有这些文件可能不值得，而且没有任何意义。我再次感到惊讶，我想问：是否有可能以某种方式在 mingw 中编译我的程序，还是放弃这个想法并使用 Visual Studio 更好？理想情况下，我想以某种方式禁用多线程，因为对此有一个宏并且它在 mingw 中被破坏了。

c++
Martin Hope
AVK
Asked: 2021-11-29 03:53:26 +0000 UTC
  • 1

亲爱的参与者，您对这段代码有什么想说的：

#include <iostream>

class A {
public:
    int x = 2020;

    class B& GetInterface()
    {
            // ...
            // Тут на самом деле делается нечто.
            // Но это не важно для вопроса
        return (B&)(*this);
    }
};

class B : public A {
public:
    void f()
    {
        std::cout << x << std::endl;
    }
};

int main()
{
    A a;
    a.GetInterface().f();
    return 0;
}

这里有未定义的行为吗？B如果类没有被覆盖的虚拟类函数A、数据——除了非虚拟函数之外什么都没有，这样的代码是否有效？

c++
Martin Hope
AVK
Asked: 2020-07-06 02:01:51 +0000 UTC
  • 2

我有一个用于积分微分方程系统的代码。它是一个带有参数的模板函数——系统的维度，类型为int。它看起来像这样：

template<int n> std::array<double, n> operator+(const std::array<double, n>& a, 
    const std::array<double, n>& b)
{       // Вспомогательная функция
    std::array<double, n> ret;
    for (int i = 0; i < n; i++) ret[i] = a[i] + b[i];
    return ret; // NRVO
}

template<int n> std::array<double, n> operator*(double h,const std::array<double, n>& b)
{        // Вспомогательная функция
    std::array<double, n> ret;
    for (int i = 0; i < n; i++) ret[i] = h * b[i];
    return ret; // NRVO
}

template<int n> std::vector<std::array<double, n>>
    rk4(std::function<std::array<double, n>(double, const std::array<double, n>&)> f,
        double t0, const std::array<double, n>& x0, double tend, double h)
{        // Функция для решения системы
    if (tend < t0) return {};
    size_t points = static_cast<size_t>(1 + floor((tend - t0) / h));
    std::vector<std::array<double, n>> ret{ points };
    ret[0] = x0;
    double t = 0;
    for (size_t i = 1; i < points; i++)
    {       // Метод Рунге-Кутты 4 порядка с пост. шагом
        auto k1 = f(t, ret[i-1]);
        auto k2 = f(t + h / 2, ret[i - 1] + (h / 2 ) * k1);
        auto k3 = f(t + h / 2, ret[i - 1] + (h / 2 ) * k2);
        auto k4 = f(t + h, ret[i - 1] + h * k3);
        ret[i] = ret[i - 1] + (h / 6) * (k1 + 2*k2 + 2*k3 + k4);
    }
    return ret;
}

这是我的使用方法：

std::array<double, 2> test1(double t, const std::array<double, 2>& x)
{    // функция, считающая правую часть системы размерности 2
     return { x[1], -x[0] };  
}

代码中的某处：

auto result = rk4<2>(test1, 0, { 1,0 }, 6.28*10, 0.05);

您不能删除尖括号中的两个，编译器将无法理解什么n等于两个。尽管这里显然是多余的，因为只有 2n才能相同。我试图为rk4演绎指南编写，以便从第三个参数的大小推断维度：

template <typename F,typename T0, typename X0, typename Tend, typename H> 
    rk4(F, T0, X0, Tend, H)->rk4<std::tuple_size<X0>::value>;

但它不起作用。我的问题是：是否可以在这里写一个演绎指南，如果可以，怎么写？

c++
Martin Hope
AVK
Asked: 2020-04-09 03:31:22 +0000 UTC
  • 3

是否有一些技巧可以让您相对快速地计算或至少估计两个给定数字之间包含的浮点数的数量？在 C++ 中，可以这样完成：

#include <cmath>
#include <iostream>

using namespace std;

int main()
{
    double a= 1, b= 1.00000001; // a - первое число, b - второе число
    long long count = 0;
    for (double x = a; x < b; x = nextafter(x, numeric_limits<double>::max()))
    {
        count++;
    }
    cout << count << endl;
    return 0;
}

但是太慢了。不需要计算确切的金额，-50 + 200% 的误差是完全可以接受的

c++