4

所以这很粗糙。我有一段代码使用了 Columbia Physics System c++ 库。运行时出现段错误，很可能是因为类：CgArg 和向量 f_field_in 等未初始化。我使用了 Valgrind，发现确实各种参数都指向无效的内存。

奇怪的是，如果我在两个函数的任何位置插入 iostream 调用，段错误就会消失。当我将标志放在要调试的地方时，我发现了。如果我定义一个整数并编写一个简单的 cin >>，它也不会出现段错误。这就是为什么我认为这与 iostream 有关。

如果您知道调用 iostream 会以某种方式为参数提供代码片段的任何原因，如果您与我分享，我将非常感激。

#include <iostream>
using namespace std;

#include <config.h>
#include <util/lattice.h>
#include <util/dirac_op.h>
#include <util/gjp.h>
#include <interface.h>

#define CLOVER_MAT_SIZE 72

USING_NAMESPACE_CPS

// Same function for clover matrix and its inverse.
static
void interface(double *h_quda_clover, double *h_cps_clover)
{
  h_quda_clover[0]=h_cps_clover[0];  // c00_00_re = C0.x, A0

......还有更多这样的东西...... 

  h_quda_clover[35+36]=h_cps_clover[34+36];  // c32_31_im = C8.w, A5
}

static
void fill_h_clover_inv(Lattice &lat, int site[], double *h_quda_clover_inv_site)
{
  double h_cps_clover_inv[72];
  Vector *f_field_out, *f_field_in;
  CgArg *arg;
  CnvFrmType convert=CNV_FRM_NO;
  DiracOpClover dirac(lat,f_field_out,f_field_in,arg,convert);
  //cout << "B: " << site << endl;
  //cout << "B: site: " << site[0] << ' ' << site[1] << ' ' << site[2] << ' ' <<         site[3] << endl;
  dirac.SiteCloverMat(site,h_cps_clover_inv);
  interface(h_quda_clover_inv_site,h_cps_clover_inv);
}

void fill_h_clover_inv_all(Lattice &lat, double *h_quda_clover_inv, int parity)
{
  double *ptr=h_quda_clover_inv;
  int nsites[4];
  nsites[0]=GJP.XnodeSites();
  nsites[1]=GJP.YnodeSites();
  nsites[2]=GJP.ZnodeSites();
  nsites[3]=GJP.TnodeSites();
  int site[4];
  cout << "A: " << site << endl;
  for (site[3] = 0; site[3] < nsites[3]; ++(site[3])) {
    for (site[2] = 0; site[2] < nsites[2]; ++(site[2])) {
      for (site[1] = 0; site[1] < nsites[1]; ++(site[1])) {
        site[0] = (site[3] + site[2] + site[1] + parity)%2;
        for (; site[0] < nsites[0]; site[0] += 2) {
  //cout << "A: site: " << site[0] << ' ' << site[1] << ' ' << site[2] << ' ' <<     site[3] << endl;
          fill_h_clover_inv(lat,site,ptr);
          ptr += CLOVER_MAT_SIZE;
        }
      }
    }
  }
}
4

1 回答 1

1

问题在于一些未初始化的参数被传递给 DiracOpClover。一旦我们正确初始化了 CgArg，

CgArg cg_arg;
cg_arg.mass=1.0

一切都很棒！

感谢大家的帮助。

于 2012-12-04T20:58:09.560 回答