- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
由于这是我的第一篇文章,我希望我不会犯太多错误......
我面临以下问题。当我运行这段代码时(很抱歉,相当长),我在 LAPACKE 调用期间的“解决”函数中遇到了段错误。
问题是当我不尝试在主函数中打开文件“dataf”时,此代码运行时没有任何(可见)问题。
我想我弄乱了我的鲍勃类型和指针,但我想听听你的观点。
致以诚挚的问候
蒂博
#include <stdio.h>
#include <stdlib.h>
#include <lapacke/lapacke.h>
typedef struct{
double *data;
double *tmp;
double **A;
int size;
} Bob;
void allocate(Bob *b);
void init_data(Bob *b);
void init_A(Bob* b);
void allocate(Bob *b){
b->data = malloc(b->size*sizeof(double));
b->tmp = malloc(b->size*sizeof(double));
b->A = malloc(4*sizeof(double*));
b->A[0] = malloc(b->size*sizeof(double));
b->A[1] = malloc(b->size*sizeof(double));
b->A[2] = malloc(b->size*sizeof(double));
b->A[3] = malloc(b->size*sizeof(double));
}
void init_data(Bob *b){
int i;
for (i=0;i<b->size;i++){
b->data[i] = 0.0;
}
}
void init_A(Bob* b){
int i;
double *diag = b->A[0];
double *down = b->A[1];
double *up = b->A[2];
double *up2 = b->A[3];
for (i=0;i<b->size;i++){
diag[i] = 2.0;
down[i] = -1.0;
up[i] = -1.0;
up2[i] = 0.0;
}
}
void solve(Bob *b){
lapack_int info;
lapack_int *ipiv = malloc(b->size*(sizeof(lapack_int)));
lapack_int lsize = b->size;
lapack_int lone = 1;
int i;
for (i=0;i<b->size;i++){
b->tmp[i] = b->data[i];
}
info = LAPACKE_dgttrs(LAPACK_COL_MAJOR,'N',
lsize,lone,\
b->A[1],b->A[0],b->A[2],b->A[3],\
ipiv, b->tmp, lsize);
printf("%d\n",info);
free(ipiv);
}
int main(void){
Bob b;
b.size = 10000;
allocate(&b);
init_A(&b);
init_data(&b);
// No Segfault if dataf is not opened
FILE *dataf;
dataf = fopen("data","w");
fclose(dataf);
solve(&b);
free(b.data);
free(b.tmp);
free(b.A[0]);
free(b.A[1]);
free(b.A[2]);
free(b.A[3]);
free(b.A);
return 0;
}
使用 Makefile:
P=mini
CC=gcc
CFLAGS=-g -Wall
LFLAGS=-lm -llapacke -llapack -lblas
FILES=min.c
all: $(FILES)
$(CC) $(CFLAGS) $(LFLAGS) $(FILES) -o $(P)
编辑:我在原来的帖子中犯了一个错误。 A[0] -> A[3] 的分配是 double ,而不是指向 double 的指针。
解决后更改 fclose 的位置似乎可以解决问题。但在我的真实代码中它不会改变任何东西......
编辑2:该问题似乎与 LAPACK 函数有关。我更改了此函数(更改为更适合我的问题的函数),现在内存问题消失了。
这不是我原来问题的解决方案,但我以某种方式解决了它。感谢所有非常有用的回复和评论:)
编辑3:我找到了问题的根源。与LAPACK有关。在函数solve中,我为ipiv分配空间。但是 LAPACKE_dgttrs 将 init_A 中的函数 LAPACKE_dgttrf 给出的 ipiv 作为输入。
因此,我在结构中为 ipiv 分配了一个指针,并在程序的其余部分中使用它。
最佳答案
错误是由于fclose函数造成的,在fclose之前检查dataf是否为NULL或noy
FILE *dataf;
dataf = fopen("data","w"); //If file is not opened for writing, the dataf is NULL
fclose(dataf);
solve(&b);
更改为
FILE *dataf;
dataf = fopen("data","w");
if(dataf != NULL)
fclose(dataf);
solve(&b);
关于c - 在 C 中的结构内部使用指针时出现段错误(可能与 LAPACKE 有关),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/25035532/
这是我的测试用例。 http://tobeythorn.com/isi/dummy2.svg http://tobeythorn.com/isi/isitest.html 如果我自己打开 svg,内部
这是我的测试用例。 http://tobeythorn.com/isi/dummy2.svg http://tobeythorn.com/isi/isitest.html 如果我自己打开 svg,内部
我正在尝试做类似的事情: SELECT SUM( CASE WHEN ( AND EXISTS(SELECT 1
我想问如何在外部 ng-repeat 内部正确使用内部 ng-repeat: 这意味着你想使用这样的东西: {{milestone.id}} {{
我希望在 wordpress 的仪表板内编辑 css 样式并且如果可能的话不必编辑 php 文件。 我知道至少可以编辑一些属性,所以我希望我可以直接在仪表板中编辑所有属性。 更具体地说如何更改自定义类
我在安装在 windows10 上的 vmware 中的 Ubuntu 上安装了伪分布式独立 hadoop 版本。 我从网上下载了一个文件,复制到ubuntu本地目录/lab/data 我在 ubun
我有一个如下所示的 WHERE 语句: WHERE ((@Value1 IS NULL AND [value1_id] IS NULL) OR [value1_id] = ISNULL(@Va
我有一个如下所示的 WHERE 语句: WHERE ((@Value1 IS NULL AND [value1_id] IS NULL) OR [value1_id] = ISNULL(@Va
在我的一些测试帮助程序代码中,我有一个名为 FakeDbSet(Of T) 的 IDbSet(Of T) 实现,它模拟了许多 EF 行为,但没有实际的数据库。我将类声明为 Friend ,因为我想强制
我正在寻找 Cassandra/CQL 的常见 SQL 习语 INSERT INTO ... SELECT ... FROM ... 的表亲。并且一直无法找到任何以编程方式或在 CQL 中执行此类操作
如何防止内部 while 循环无限运行?问题是,如果没有外部 while 循环,内部循环将毫无问题地运行。我知道它必须对外循环执行某些操作,但我无法弄清楚是什么导致了问题。 import java.u
我正在努力学习更多有关 C++ 的知识,但在国际象棋程序中遇到了一些代码,需要帮助才能理解。我有一个 union ,例如: union b_union { Bitboard b; st
这是我项目网页中的代码片段。这里我想显示用户选择的类别,然后想显示属于该类别的主题。在那里,用户可以拥有多个类别,这没有问题。我可以在第一个 while 循环中打印所有这些类别。问题是当我尝试打印主题
我想知道如何在 swing 中显示内部框架。这意味着,当需要 JFrame 时,通常我所做的是, new MyJFrame().setVisible(true); 假设之前的表单也应该显示。当显示这个
我最近发现了一些有趣的行为,这让我想知道对象如何知道存在哪些全局变量。例如,假设我有一个文件“test.py”: globalVar = 1 toDelete = 2 class Test(objec
我知道它已经在这里得到回答: google maps drag and drop objects into google maps from outside the Map ,但这并不完全是我所需要的
我目前正在学习Javascript DOM和innerHTML,发现在理解innerHTML方面存在一些问题。 这是我的代码:http://jsfiddle.net/hphchan/bfjx1w70/
我构建了一个布局如下的库: lib/ private_class_impl.cc private_class_decl.h public_class_impl.cc include/
我有一个使用 bootstrap 3 的组合 wordpress 网站。它基本上是一个图像网格。当屏幕展开时,它会从三列变为四列。移动时它是一列。 我想出了如何调整图像的顶部和底部边距,但我希望图像的
我正在试用 MSP-EXP430G2 的教程程序,使用 Code Composer Studio 使 LED 闪烁。最初,它有一个闪烁的无限循环: for(;;) // This emp
我是一名优秀的程序员,十分优秀!