c++ - 程序在尝试运行时停止工作
问题描述
我为拉格朗日插值编写了一个可以编译的程序。不幸的是,当我在代码块中运行程序时,它只显示警告“程序停止工作”。
在 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;
}
解决方案
你的主要问题在这里。
double *tabT1;
double *tabT2;
tabT1[0]=1;
tabT1[1]=tabX[j]*(-1);
tabT2[0]=1;
tabT2[1]=tabX[i]*(-1);
您还没有分配内存,相反,您刚刚声明了double
ptrstabT1
并tabT2
通过假装您已分配来访问它们。
double *tabT1 = new double[2];
double *tabT2 = new double[2];
将解决此问题,但是,我强烈建议您改用智能指针,这在您的情况下更安全,因为看起来您有很多指针数组。
编辑:上述修复将导致内存泄漏,因为您不必担心在每次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);
}
附带说明:double tabX[wierz];
ISO C++ 标准中禁止声明类型(可变长度数组)。另一种方法是std::vector<>
,它还将为您动态分配一个连续的数组。
推荐阅读
- typescript - 是否有禁止前缀类型断言语法的 TSLint 规则?
- filemaker - FileMaker 错误:无法在此磁盘上创建 PDF
- java - 错误无法打开属性文件:来自 Ant 脚本的 AppManage.tra
- python - 这段代码有什么问题?它没有显示任何错误,但也没有显示任何输出
- report - 将DataBand放在某个空间
- javascript - else 语句 deadend 后无法从头开始启动路由
- java - 如何将所有 ASCII 值放入数组中?
- oracle-apex - APEX 行选择器第 2 部分
- sql - sql查询删除重复行
- php - TinyMce 加长