Question

我为拉格朗日插值编写了一个编译程序。不幸的是，当我在代码块中运行程序时，它只显示警告＆＃34;程序停止工作＆＃34; 。

在Dev-C ++中，它返回的值超过32 000 000。

你能检查代码并告诉我原因吗？你有代码：

#include <iostream>
#include<stdlib.h>
using namespace std;

struct wielo
{
    double a, b, c, d, e, f;
};

double *mnozenie(double A[], double B[], int m, int n)
{
    double *prod = new double[m + n - 1];

    for (int i = 0; i < m + n - 1; i++)
        prod[i] = 0;

    for (int i = 0; i < m; i++)
    {
        for (int j = 0; j < n; j++)
            prod[i + j] += A[i] * B[j];
    }

    return prod;
    delete prod;
}

int main()
{
    int wierz = 3; // max 6
    int sizeW = wierz - 1;
    double wynik = 0;
    double szukana = 5;
    double tabX[wierz];
    double tabY[wierz];
    double *tabTT;
    double wynikowy[wierz];

    tabX[0] = 1;
    tabX[1] = 2;
    tabX[2] = 4;

    tabY[0] = 3;
    tabY[1] = 4;
    tabY[2] = 6;

    wielo w[] =
    {
    { 0, 0, 0, 0, 0, 0 },
    { 0, 0, 0, 0, 0, 0 },
    { 0, 0, 0, 0, 0, 0 } };

    wynikowy[0] = 0.0;
    wynikowy[1] = 0.0;
    wynikowy[2] = 0.0;
    int k = 0;
    for (int i = 0; i < wierz; i++)
    {
        double dzialanie = 0;
        dzialanie = tabY[i];

        for (int j = 0; j < wierz; j++)
        {
            if (tabX[i] != tabX[j])
            {
                dzialanie *= (szukana - tabX[j]);
            }
        }

        for (int j = 0; j < wierz; j++)
        {
            if (tabX[i] != tabX[j])
            {
                dzialanie /= (tabX[i] - tabX[j]);
            }
            if (tabX[i] != tabX[j] && i < j)
            {
                double *tabT1;
                double *tabT2;
                tabT1[0] = 1;
                tabT1[1] = tabX[j] * (-1);
                tabT2[0] = 1;
                tabT2[1] = tabX[i] * (-1);
                tabTT = mnozenie(tabT1, tabT2, sizeW, sizeW);
                w[k].a = tabTT[0];
                w[k].b = tabTT[1];
                w[k].c = tabTT[2];
                k++;
            }
        }
        wynik += dzialanie;
        cout << dzialanie << endl;
    }
    double x, y, z;
    x = tabY[0] / ((tabX[0] - tabX[1]) * (tabX[0] - tabX[2]));
    y = tabY[1] / ((tabX[1] - tabX[0]) * (tabX[1] - tabX[2]));
    z = tabY[2] / ((tabX[2] - tabX[0]) * (tabX[2] - tabX[1]));
    for (int i = 0; i < wierz; i++)
    {
        if (i == 0)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * x);
            wynikowy[1] = wynikowy[1] + (w[i].b * x);
            wynikowy[2] = wynikowy[2] + (w[i].c * x);
        }
        if (i == 1)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * y);
            wynikowy[1] = wynikowy[1] + (w[i].b * y);
            wynikowy[2] = wynikowy[2] + (w[i].c * y);
        }
        if (i == 2)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * z);
            wynikowy[1] = wynikowy[1] + (w[i].b * z);
            wynikowy[2] = wynikowy[2] + (w[i].c * z);
        }

    }
    cout << "Wzor to: " << wynikowy[0] << "X^2 " << wynikowy[1] << "X " << wynikowy[2] << endl;
    cout << "Wynik dla x=" << szukana << " " << "y=" << wynik << endl;

    return 0;
}

Answer 1

这里有你的主要问题。

double *tabT1;
double *tabT2;
tabT1[0]=1;
tabT1[1]=tabX[j]*(-1);
tabT2[0]=1;
tabT2[1]=tabX[i]*(-1);

你还没有分配内存，相反，你刚刚宣布了double ptrs tabT1和tabT2并通过假装你已经分配来访问它们。

double *tabT1 = new double[2];
double *tabT2 = new double[2];

会解决此问题，但是，我强烈建议您使用 smart pointers ，这在您的情况下会更安全，因为看起来您有很多指针数组。 / p>

编辑：上面的修复将导致内存泄漏，因为在使用它们之后，每次在for循环结束时都不会想要删除这些指针数组。一个可能的智能指针修复可能是：

std::unique_ptr<double[]> tabTT = nullptr;  // change
.....
.....

if(tabX[i] != tabX[j] && i<j)
{
   std::unique_ptr<double[]> tabT1 = std::unique_ptr<double[]>(new double[2]);  // change
   std::unique_ptr<double[]> tabT2 = std::unique_ptr<double[]>(new double[2]);  // change
   .....
   .....
   tabTT = mnozenie(std::move(tabT1), std::move(tabT2),sizeW,sizeW);
   .....
}

在功能中：

auto mnozenie(std::unique_ptr<double[]> A, std::unique_ptr<double[]> B, int m, int n)
{
   std::unique_ptr<double[]> prod = std::unique_ptr<double[]>(new double[m+n-1]);
   // calculations....    
   return std::move(prod);
}

在ISO C ++标准中禁止使用：double tabX[wierz];类型的声明（可变长度数组）。替代方案是std::vector<>，它还将为您提供动态分配的连续数组。

试图运行时程序停止工作

1 个答案: