拉帕克例程DSTEV中的分割故障
Segmentation Fault in Lapack routine dstev
我目前只是编写一个小程序来使用C 中的有限差方法求解微分方程。问题是我正在使用拉帕克(Lapack)来求解我的tridiagonal矩阵,但是在调用例程之后,我会遇到分段故障。例程仍然带有信息= 0退出,所以我不知道问题是什么。这是程序。
#include <iostream>
#include <cmath>
#include <vector>
#define pi 3.14159265358979323846
using namespace std;
extern "C" void dstev_(char* job, int* N, double* D, double* OFFD, double* EV, int* VDIM, double* WORK, int* INFO);
int main(){
int i = 0;
// systems parameters
double a = 2e10-6;
double k = pi/1.55e-6;
double n1_sq = pow(3.0, 2), n2_sq = pow(3.1, 2);
double step = 1e-6;
double factor = 1/(k*k*step*step);
// parameters for lapack
char job = 'V';
int N = 80, vdim = 100, info;
// create and initialize vector arrays for lapack routine
vector<double> d, offd, work;
vector< vector<double> > ev(vdim, vector<double>(N));
// set up diagonal elements
d.push_back(n1_sq - factor);
for(i=1; i<N-1; i++){
if(i<N/4 || i>=3*N/4){
d.push_back(n1_sq - 2*factor);
}
else{
d.push_back(n2_sq - 2*factor);
}
}
d.push_back(n1_sq - factor);
// set up off diagonal elements
for(i=0; i<N-1; i++){
offd.push_back(factor);
}
// initialize other arrays
for(i=0; i<vdim; i++){
work.push_back(0.0);
}
cout << "Before routine" << endl;
dstev_(&job, &N, &*d.begin(), &*offd.begin(), &*(ev.begin()->begin()), &vdim, &*work.begin(), &info);
cout << "After routine"<< endl;
return 0;
}
这是我第一次使用Lapack,所以我不知道该代码有什么问题/正确,但是对任何帮助都表示赞赏。另外,要将矢量传递到例程中,为什么我们需要传递begin()迭代器?这就是我学习的方式,但我认为我完全不了解它。
我不知道这是否是此代码的唯一问题,但是std::vector< std::vector< double> >
不分配连续的内存...因此,Lapack无法访问Matrix的所有元素ev
。您应该给出一个大小 n*n
的向量,该矩阵的条目(i,j)
为 i+j*n
。
相关文章:
- 分割故障C++矩阵
- C ,二进制树的高度,而不是检查我的子树是否为空,而是在检查我的子树是否是叶子节点.抛出分割故障
- 我正在为此代码分割故障
- 在执行strcpy时获取分割故障错误
- 空隙指针分割故障
- Xcode会产生分割故障,而不是例外
- OPENCV C 中的分割故障误差
- 分割故障:核心倾倒C 矢量对字符串:
- C 2D数组导致分割故障(核心倾倒)
- 分割故障/错误读数变量-C 列表
- DLSYM导致分割故障
- 使用C 中的操作员删除的分割故障
- 使用Linux X64下的Libmozjs-52(SpidermonKey)进行分割故障
- 使用长变量时分割故障(11)错误
- C - 返回字符串时的分割故障
- 创建VertexArray时的分割故障
- 向量分割故障的向量
- 从std :: fileSystem ::路径对象的类中退出时的分割故障
- C 矢量分割故障
- 基于向量的数据格式的分割故障误差