我定义了这个结构:
struct s_molecule
{
std::string res_name;
std::vector<t_particle> my_particles;
std::vector<t_bond> my_bonds;
std::vector<t_angle> my_angles;
std::vector<t_dihedral> my_dihedrals;
s_molecule& operator=(const s_molecule &to_assign)
{
res_name = to_assign.res_name;
my_particles = to_assign.my_particles;
my_bonds = to_assign.my_bonds;
my_angles = to_assign.my_angles;
my_dihedrals = to_assign.my_dihedrals;
return *this;
}
};
和这些结构:
typedef struct s_particle
{
t_coordinates position;
double charge;
double mass;
std::string name;
std::vector<t_lj_param>::iterator my_particle_kind_iter;
s_particle& operator=(const s_particle &to_assign)
{
position = to_assign.position;
charge = to_assign.charge;
mass = to_assign.mass;
name = to_assign.name;
my_particle_kind_iter = to_assign.my_particle_kind_iter;
return *this;
}
} t_particle;
struct s_bond
{
t_particle * particle_1;
t_particle * particle_2;
std::vector<t_bond_param>::iterator my_bond_kind_iter;
s_bond& operator=(const s_bond &to_assign)
{
particle_1 = to_assign.particle_1;
particle_2 = to_assign.particle_2;
my_bond_kind_iter = to_assign.my_bond_kind_iter;
return *this;
}
};
然后在我的代码中我返回一个指向s_molecule的指针(typedef'd to t_molecule,但仍然)。
使用此指针我可以使用此代码:
for (unsigned int i = 0;
i < current_molecule->my_particles.size();
i++)
{
std::cout << "Particle "
<< current_molecule->my_particles[i].name << std::endl
<< "Charge: "
<< current_molecule->my_particles[i].charge << std::endl
<< "Mass: "
<< current_molecule->my_particles[i].mass << std::endl
<< "Particle Kind Name: "
<< (*current_molecule->my_particles[i].my_particle_kind_iter).atom_kind_name
<< std::endl
<< "x: " << current_molecule->my_particles[i].position.x
<< " y: " << current_molecule->my_particles[i].position.y
#ifdef USE_3D_GEOM
<< "z: " << current_molecule->my_particles[i].position.z
#endif
<< std::endl;
}
如果我将其替换为:
for (std::vector<t_particle>::iterator it = current_molecule->my_particles.begin();
it !=current_molecule->my_particles.end();
it++)
{
std::cout << "Particle "
<< (*it).name << std::endl
<< "Charge: "
<< (*it).charge << std::endl
<< "Mass: "
<< (*it).mass << std::endl
<< "Particle Kind Name: "
<< (*(*it).my_particle_kind_iter).atom_kind_name
<< std::endl
<< "x: " << (*it).position.x
<< " y: " << (*it).position.y
#ifdef USE_3D_GEOM
<< "z: " << (*it).position.z
#endif
<< std::endl;
}
我现在得到了令人讨厌的段错误......
不要在这里放太多,但是当我尝试这样做时,我也会遇到段错误:
std::cout << "Bond ATOMS : "
<< (*current_molecule).my_bonds[0].particle_1->name
<< std::endl
同样,current_molecule是一个指向s_molecule结构的指针,该结构包含结构数组,而结构数组又直接包含vars或指针。我不能让这些多层间接工作。关于修复这些段错误的建议。
仅供参考我正在使用g ++和使用自定义makefile系统在Linux Centos 5.4上进行编译。
答案 0 :(得分:0)
@sbi感谢您的好建议!我相信你是对的 - 分配重载运算符是不必要的,应该废弃。
我已经按照评论出来的方法而且非常困惑。基本上在将指向我的特定分子的指针传递给主要功能进行打印的功能中,我可以完美地看到该分子中的所有数据(键,粒子,名称等),并与cout打印。
一旦我将它作为ptr传递给main,如果我将该ptr与迭代器一起使用,我会得到一个段错误。换一种说法。也是出于某种原因,债券数据(我可以在我的函数中自由打印,返回到指针)也会出现段错误,如果我尝试打印它,即使我使用[]索引键的矢量(适用于粒子)向量)。
这是我现在可以提供的最佳信息。
答案 1 :(得分:0)
一个疯狂的猜测:您使用的是共享库吗?我记得在共享库边界之间来回传递STL容器时遇到了困难。
答案 2 :(得分:0)
杰森回答说:您是否返回指向局部变量的指针?
没有它是一个类变量的ptr。该类非常存在(它包含返回分子的函数)。
除非你在谈论一个真正的类变量(限定为static
),否则类存在的事实与它没有太大关系。存在一个类的 Instances ,即使你只是在它们上面调用了一个函数,它们也可能已经不存在了。
因此,问题是:
current_molecule
的类的实例是否仍然存在?current_molecule
是否为static
,即是真正的类变量?如果两个问题的答案都是“否”,那么你就在Undefined County。
此时,您发布可供我们使用的源代码以实际重现问题变得非常重要;它可能位于源代码中,但您没有向我们展示。
答案 3 :(得分:0)
同样,这个问题在这里得到了回答: Weird Pointer issue in C++ 由DeadMG提供。对不起双重帖子。