c - binning/jack-knife c 程序中的段错误-6ren

c - binning/jack-knife c 程序中的段错误

转载作者：太空宇宙更新时间：2023-11-03 23:24:52

为什么会出现段错误(核心已转储)？我不知道问题出在哪里。

该程序应该读取包含热化数据的数据文件，对数据进行装箱，然后使用折刀法重新采样。

无论使用多少数据，我都会收到段错误。

/* Ben Marshall 12/03/15 
Bin data & Jack-Knife  */

#include <stdio.h>
#include <stdlib.h>


int main(void)
{
  int binSize, i=0, j=0, k=0, s=0, dataSize=0;
  double *data, *binData, *jackData, mean,  test; //variance

  /* open files, therm.dat is the data to be jack-knifed, 
     jack-knife.dat will store the new data  */
  FILE *thermalised = fopen("therm.dat","r");
  FILE *jack = fopen("jack-knife.dat","w");
  if(thermalised==NULL || jack==NULL) {
    return(-1);
  }

  // check amount of data to be used
  while(fscanf(thermalised, "%lf", &test) != EOF) {
    i++;
  }

  // close and reopen file to reset to start of file
  fclose(thermalised);
  FILE *therm = fopen("therm.dat","r");
  if(therm==NULL) {
    return(-1);
  }

  dataSize=i;
  data=malloc(dataSize*sizeof(double));
  binData=malloc(dataSize*sizeof(double));
  jackData=malloc(dataSize*sizeof(double));

  i=0;
  while(fscanf(therm, "%lf", &data[i]) != EOF) {
    i++;
  }

  fclose(therm);

  /* increasing bin size in multiples of 2  */
  for(binSize=1; (dataSize/binSize)>20; (binSize=binSize*2)) {

    /* binning the data */
    for(j=0; j<dataSize; j=(j+binSize)) {
      for(k=0; k<binSize; k++) {
    binData[s] = binData[s] + data[j+k];
      }
      binData[s] = binData[s]/binSize;
      s++;
    }

    /* jack-knifing the binned data */
    for(i=0; i<(dataSize/binSize); i++) {
      for(s=0; s<(dataSize/binSize); s++) {
    jackData[i] = jackData[i] + binData[s] - binData[i];
      }
      jackData[i] = jackData[i]/((dataSize/binSize)-1); 
    }

    /* calculate mean and variance(FINISH!!) of jackData */
    for(i=0; i<(dataSize/binSize); i++) {
      mean = mean + jackData[i];
    }
    mean = mean/(dataSize/binSize);
    fprintf(jack, "%d\t%lf\n",binSize,mean);
  }

  free(data);
  free(binData);
  free(jackData);
  fclose(jack);
  return(0);
}

包含 512 个数据点的 Therm.dat 文件

最佳答案

您的代码存在一些问题:

while(fscanf(thermalised, "%lf", &test) != EOF)

如果文件包含无法解析为 float 的内容，这将导致无限循环。改用这个:

while(fscanf(thermalised, "%lf", &test) == 1)

为什么然后关闭并重新打开文件而不是调用 rewind() ？

在第二个读取循环中，你应该更加小心，以防另一个进程正在写入 "them.dat"异步归档。您应该尝试解析尽可能多的数字:

for (i = 0; i < dataSize; i++) {
    if (fscanf(therm, "%lf", &data[i]) != 1) {
        printf("error: them.dat changed\n");
        exit(1);
    }
}

你不应该返回 -1来自 main() , 大多数系统只支持 0 之间的退出状态和 255 .

在装箱阶段，您有一个潜在的缓冲区溢出:

for(binSize=1; (dataSize/binSize)>20; (binSize=binSize*2)) {

    /* binning the data */
    for(j=0; j<dataSize; j=(j+binSize)) {
        for(k=0; k<binSize; k++) {
            binData[s] = binData[s] + data[j+k];
        }
        binData[s] = binData[s]/binSize;
        s++;
    }

如果dataSize不是 2 的幂，最后一个 bin 可能会超出 data 的末尾数组，调用未定义的行为，可能是您遇到的崩溃。修改边界测试以避免这种情况。

此外，binData和 jackData在初始化之前使用。您应该使用 calloc 分配这些数组所以他们从0.0开始(在具有 IEEE float 的系统上 ;-)，或使用循环将它们初始化为适当的值。

当您计算 mean 时:

    /* calculate mean and variance(FINISH!!) of jackData */
    for(i=0; i<(dataSize/binSize); i++) {
        mean = mean + jackData[i];
    }

mean没有初始化!

打印时:

fprintf(jack, "%d\t%lf\n",binSize,mean);

"%lf"未定义，你应该使用 "%f"打印 double或 float .

关于c - binning/jack-knife c 程序中的段错误，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/29721976/

文章推荐： c# - 将列表或数组从 Arden MLM 传递到 C#

文章推荐： node.js - 如何使用 ES6 在 Koa.js 中创建中间件

文章推荐： c - 为什么我会遇到段错误(二进制数据解析)c

git - .gitignore 语法 : bin vs bin/vs. bin/* 与 bin/**
bin、bin/、bin/*和bin/**有什么区别我的 .gitignore 文件？我一直在使用 bin/，但正在查看 other .gitignore files (在 eclipse file
optimization - Bin Packing : Set amount on bins, 想要最小化最大 bin 重量
给定 n 个无限容量的箱子，我想将 m 件元素装入其中(每个元素都有特定的重量)，同时最小化最重箱子的重量。这不是传统的垃圾箱包装/背包问题，其中垃圾箱的容量有限，而您试图尽量减少使用的垃圾箱数量；
ssh - bin/sh，bin/bash，sbin/nologin，bin/tcsh等有什么区别？
Closed. This question is off-topic。它当前不接受答案。想改善这个问题吗？ Update the question，所以它是用于堆栈溢出的on-topic。 7年前关
java - bin packing 的具体变化(n bins 优化以最大化 bins 的最小值)
我一直在寻找一个问题的解决方案，这个问题显然比我最初想象的更不寻常。 Optaplanner 看起来很有前途，但由于我对 Java 相对缺乏经验，所以我想在深入研究之前调查一下这是否完全不可能。我正
python - 在所有数据框列上应用具有不同 bin 大小的 binning
我有一个琐碎的问题。我有一个非常大的 df 有很多列。我正在尝试找到最有效的方法来对具有不同 bin 大小的所有列进行 bin 并创建一个新的 df。这是一个仅对单个列进行分箱的示例: import
sql - 在每个 bin 中创建具有唯一值的 bin
我想以这样的方式对数字列 (var) 进行分类，使每个分类中的行数大致相同。我的附加要求是该列中的一个(唯一)值不能分配给多个 bin。例如，如果 var 列中的值 1 分配给 bin 1，则不允许将
linux -/bin/sh :/bin: Is a directory
我在 centos (rhel 7)_ 上使用 make 命令和 gcc 7.1 通过代码进行编译，但遇到了以下错误: /bin/sh:/bin: 是一个目录 which sh 返回/usr/bin/
c++ -/usr/bin/ld : cannot find -ldlib/usr/bin/ld: cannot find -lcblas/usr/bin/ld: cannot find -llapack
/usr/bin/ld: cannot find -ldlib /usr/bin/ld: cannot find -lcblas /usr/bin/ld: cannot find -llapack 在
r - 如何将 span bin 信息转换为单独的 bins？
我进行了一项眼动追踪实验，试图检测两个区域(感兴趣区域，名为“代理”和“患者”的 AOI)的注视分布。我将整个时间划分为时间段，例如得到以下列表: Stimulus Participant A
linux - ~/.local/bin 中的文件优先于/usr/bin
我正在尝试替换 whoami带有 ~/.local/bin/ 中的脚本的命令.有没有办法让我的 whoami 获得优先权，这样当我运行 whoami 时，我的脚本就会运行？最佳答案这就是我的 ~/
import - 如何在另一个 bin 中重用主 bin 中的代码？
我的项目结构是这样的: . ├── Cargo.lock ├── Cargo.toml └── src ├── bin │ └── other.rs ├── main.rs
Java - 从 BIN 文件读取数组或将数组写入 BIN 文件
我正在开发一个小型图书馆应用程序，它以这种格式存储技术手册: 目前，我正在尝试将库的内容保存并根据用户的需要加载到 bin 文件中。但是，当我尝试加载库文件时，唯一的变化是“-1”被打印到控制台。此
GitIgnore - 忽略 bin/但包括 bin/*.refresh
我试图忽略 Web 项目的 bin 文件夹，但包含该 bin 文件夹中的 .refresh 文件。这是我的 .gitignore 中的内容: [Bb]in/ #Allow .refresh file
verilog - Functional Coverage - 收集所有未在其他 bin 中收集的值的 bin
ipv6_hdr__f_next_header_cp: coverpoint this.ipv6_hdr.ipv6_f_next_header iff (this.has_ipv6_header){
bash - #!/bin/sh 和 #:/bin/sh 之间的区别
今天有人给我发了一个以 #: 开头的脚本，谷歌搜索后我没有找到任何答案。即使脚本有效，我想知道那是什么意思。最佳答案哇!这让我想起了很多记忆! 回到 1980 年代和 90 年代初期，有两种基本
linux - 将#!/bin/ksh 日期转换函数翻译成#!/bin/sh
我使用这个 ksh 函数将“1-Jan-2011”格式转换为“1.1.2011”。 #!/bin/ksh ##---- function to convert 3 char month into nu
c# - Bin/Debug 和 Bin/Release 的区别
我在 c# winform 项目中使用 sql-server compact 数据库。如果我查看构成解决方案的文件，我可以在以下两个文件夹中看到 exe 和 sdf 文件的副本: /bin/Debu
python - Seaborn 直方图 bin 宽度未扩展到 bin 标签
这个问题与我上一个问题不同。我正在通过以下代码使用 facetgrid 打印直方图。 import numpy as np import pandas as pd import seaborn as
python - 给定一系列 bin 概率，如何生成 bin 计数的随机样本？
我有一个整数需要根据概率分布分成 bin。例如，如果我有 N=100 对象进入 [0.02, 0.08, 0.16, 0.29, 0.45] 那么你可能会得到 [1, 10, 20 , 25, 44]
linux - #!/bin/bash --login vs #!/bin/bash
我在寻找使用 Jenkins 运行 Appium 的脚本时阅读了以下脚本 #!/bin/bash --login killall -9 "iPhone Simulator" &> /dev/null

太空宇宙

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

c - binning/jack-knife c 程序中的段错误