0

所以我有一个包含三列的文件,例如:

1 1 750

这些是 x,y,强度值。

然后我尝试将这个长文件读入数组。这是到目前为止的代码:

using std::cout;
using std::endl;
using std::string;

string fOutFileName("gaintest.root"); 


int main()
{

  std::ifstream file1("data_p30.dat");

  double intensity;
  int i;
  int j ;
  double c[1000][1000];
if (file1.is_open()) {
    file1.seekg(0);
    while (!file1.eof()) {
      file1 >> i >> j >> intensity;
      c[i][j]=intensity;
      cout<<c[i][j]<<'/n'; 
    }
    file1.close();
  } else cout << "Error, cannot open file 1";
}

所以最终我希望能够拥有一个与强度相关的二维数组。关于我为什么失败的任何想法?它编译得很好,但是当它运行时,它看起来像这样:

root [0] 
 *** Break *** segmentation violation



===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0  0x0000003701098c05 in waitpid () from /lib64/libc.so.6
#1  0x000000370103c481 in do_system () from /lib64/libc.so.6
#2  0x00002b036f5ebc6a in TUnixSystem::StackTrace() ()
   from /batchsoft/root/root528-x86_64-slc5-43/lib/libCore.so
#3  0x00002b036f5eb63c in TUnixSystem::DispatchSignals(ESignals) ()
   from /batchsoft/root/root528-x86_64-slc5-43/lib/libCore.so
#4  <signal handler called>
#5  0x00002b0370acd515 in TRint::Run(bool) ()
   from /batchsoft/root/root528-x86_64-slc5-43/lib/libRint.so
#6  0x000000000040106d in main ()
===========================================================


The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5  0x00002b0370acd515 in TRint::Run(bool) ()
   from /batchsoft/root/root528-x86_64-slc5-43/lib/libRint.so
#6  0x000000000040106d in main ()
===========================================================
4

2 回答 2

2

在 C++ 中

double c[1000, 1000];

前 1000 个被丢弃,编译器生成以下代码:

double c[1000];

你写的 [1000, 1000] 是 Pascal 中的多维数组,但不是 C++。在 C++ 中,您将使用:

double c[1000][1000];
于 2013-03-06T13:37:31.630 回答
0

您显示的错误是由 ROOT 生成的。由于这不是 ROOT 论坛,ROOT 是由 CERN 的科学家开发的 C++ 数据分析框架。对于好奇的人,他们的网站在这里。但是,您的问题与 ROOT 无关。

二维数组由

double a[2][2];

这对你有用。

也许更安全的是使用模板类

std::vector<double> v;

对于 2D 应用程序,这看起来像

std::vector<std::vector<double> > v2;

这样做的好处是可以根据需要调整大小,

v.push_back(d);

将在向量的末尾添加一个元素,并在必要时将其拉长。也可以使用数组语法访问元素,例如

v[1] 

或者

v2[1][2]
于 2013-04-09T00:18:56.377 回答