【发布时间】：2011-08-26 14:24:06
【问题描述】：
在以下代码中：
void ResourceFitter::copyToLarvalResources( const gsl_vector* input ){
assert( input->size == invLarvalResources.size() );
if( invLarvalResources.size() != 365 ){
cout<<"error: iLR.size(): "<<invLarvalResources.size()<<endl;
exit(21);
}
// inverting larval resources may help fitting algorithm, so we do that here:
for( size_t i=0; i<invLarvalResources.size(); ++i ){
if( i >= 365 ){
cout<<"error: i="<<i<<endl;
exit(22);
}
double val = gsl_vector_get( input, i );
invLarvalResources[i] = 1.0 / val;
}
}
这是来自 ResourceFitter.cpp 的几行代码。第 380 行是最后一行代码（在 invLarvalResources 中赋值）。 invLarvalResources 是vector<double>。
valgrind 抱怨：
==30152== Invalid write of size 8
==30152== at 0x8ED3E5: OM::Transmission::Vector::ResourceFitter::copyToLarvalResources(gsl_vector const*) (ResourceFitter.cpp:380)
==30152== by 0x8ECE99: OM::Transmission::Vector::ResourceFitter::sampler(gsl_vector const*) (ResourceFitter.cpp:334)
==30152== by 0x8EB64D: OM::Transmission::Vector::ResourceFitter_minimise_sampler(gsl_vector const*, void*) (ResourceFitter.cpp:88)
==30152== by 0x4F3249A: ??? (in /usr/lib/libgsl.so.0.16.0)
==30152== by 0x8ED5C4: OM::util::MultidimMinimiser::MultidimMinimiser(gsl_multimin_fminimizer_type const*, unsigned long, double (*)(gsl_vector const*, void*), void*, gsl_vector*, gsl_vector*) (MultidimSolver.h:71)
==30152== by 0x8EC45E: OM::Transmission::Vector::ResourceFitter::fit(unsigned long, OM::Transmission::Vector::ResourceFitter::FitMethod, unsigned long) (ResourceFitter.cpp:217)
==30152== by 0x8EC308: OM::Transmission::Vector::ResourceFitter::fit() (ResourceFitter.cpp:183)
==30152== by 0x8E474F: OM::Transmission::Vector::SpeciesModel::init2(unsigned long, std::list<OM::Host::Human, std::allocator<OM::Host::Human> > const&, int, double) (SpeciesModel.cpp:393)
==30152== by 0x8DA29E: OM::Transmission::VectorModel::init2(std::list<OM::Host::Human, std::allocator<OM::Host::Human> > const&, int) (VectorModel.cpp:163)
==30152== by 0x81C08F: OM::Population::createInitialHumans() (Population.cpp:165)
==30152== by 0x8172E8: OM::Simulation::start() (Simulation.cpp:120)
==30152== by 0x816615: main (openMalaria.cpp:53)
==30152== Address 0x8ca3d10 is 0 bytes inside a block of size 2,920 free'd
==30152== at 0x4C2658C: operator delete(void*) (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==30152== by 0x822515: __gnu_cxx::new_allocator<double>::deallocate(double*, unsigned long) (new_allocator.h:98)
==30152== by 0x81FDBB: std::_Vector_base<double, std::allocator<double> >::_M_deallocate(double*, unsigned long) (stl_vector.h:156)
==30152== by 0x81FC98: std::_Vector_base<double, std::allocator<double> >::~_Vector_base() (stl_vector.h:142)
==30152== by 0x81E82B: std::vector<double, std::allocator<double> >::~vector() (stl_vector.h:351)
==30152== by 0x8DB5D2: OM::Transmission::Vector::MosqLifeCycleParams::~MosqLifeCycleParams() (in /home/dhardy/code/openmalaria/build-debug/openMalaria)
==30152== by 0x8DBA24: OM::Transmission::Vector::MosquitoTransmission::~MosquitoTransmission() (in /home/dhardy/code/openmalaria/build-debug/openMalaria)
==30152== by 0x8E4719: OM::Transmission::Vector::SpeciesModel::init2(unsigned long, std::list<OM::Host::Human, std::allocator<OM::Host::Human> > const&, int, double) (SpeciesModel.cpp:391)
==30152== by 0x8DA29E: OM::Transmission::VectorModel::init2(std::list<OM::Host::Human, std::allocator<OM::Host::Human> > const&, int) (VectorModel.cpp:163)
==30152== by 0x81C08F: OM::Population::createInitialHumans() (Population.cpp:165)
==30152== by 0x8172E8: OM::Simulation::start() (Simulation.cpp:120)
==30152== by 0x816615: main (openMalaria.cpp:53)
调试行 (cout
那么 valgrind 在这方面是完全错误的还是发生了其他事情？在没有 valgrind 的情况下运行我没有遇到这个问题，但是我确实得到了一个奇怪的数值结果，这在独立的单元测试中是找不到的。在解决另一个错误之前，我还遇到了一个 SIGSEGV（我有把握地确定不应该导致 SIGSEGV）。
Build-flags 是 -g -rdynamic，编译器是 gcc 4.6.1 via ccache。如果你想要完整的资源，我很乐意为你指出：svn co http://openmalaria.googlecode.com/svn/branches/vec-lifecycle
任何帮助将不胜感激！
【问题讨论】：
-
我个人从未见过 Valgrind 错了
-
可以
gsl_vector_get修改i的值吗？这段代码是多线程的吗？ -
也许没有错误，但它可能会产生误导，因为它不一定会检测到对分配内存的错误写入。
-
@Matthieu：没有（i 是按值传递的），没有。
-
你能显示 SpeciesModel.cpp:391 周围的线条吗？ valgrind 声称
OM::Transmission::Vector::MosquitoTransmission::~MosquitoTransmission是从那里调用的（表明 MosquitoTransmission 类及其内容的破坏）。