为什么 Valgrind 在 PETSc 中分配 char* 时报告大小为 8 的无效写入？答案

【问题标题】：Why Valgrind reports an invalid write of size 8 on the assignment of a char* in PETSc?为什么 Valgrind 在 PETSc 中分配 char* 时报告大小为 8 的无效写入？
【发布时间】：2016-02-23 20:00:24
【问题描述】：

我在 PETSc 中的申请被终止了。我检查了 Valgrind 的情况，但我不明白它的报告：

==97331== Invalid write of size 8
==97331==    at 0x10007FED5: PetscHeaderCreate_Private (inherit.c:40)
==97331==    by 0x1013EFE23: TSResilCreate (tsresil.c:525)
==97331==    by 0x10139A877: TSGetResil (ts.c:5252)
==97331==    by 0x10138D17E: TSSetFromOptions (ts.c:421)
==97331==    by 0x100011A47: SolveODE (in ./ex31)
==97331==    by 0x100012E98: main (in ./ex31)
==97331==  Address 0x102dde080 is 1,776 bytes inside a block of size 1,780 alloc'd
==97331==    at 0x100025EA1: malloc (vg_replace_malloc.c:303)
==97331==    by 0x10010E8C1: PetscMallocAlign (mal.c:34)
==97331==    by 0x10011063D: PetscTrMallocDefault (mtr.c:188)
==97331==    by 0x1013EFDA0: TSResilCreate (tsresil.c:525)
==97331==    by 0x10139A877: TSGetResil (ts.c:5252)
==97331==    by 0x10138D17E: TSSetFromOptions (ts.c:421)
==97331==    by 0x100011A47: SolveODE (in ./ex31)
==97331==    by 0x100012E98: main (in ./ex31)

在 tsresil.c:525（我在 PETSc 中的实现）：

ierr = PetscHeaderCreate(resil,TSRESIL_CLASSID,"TSResil","Time stepping resilience","TS",comm,TSResilDestroy,TSResilView);CHKERRQ(ierr);

在inherit.c中：

 22: /*
 23:    PetscHeaderCreate_Private - Creates a base PETSc object header and fills
 24:    in the default values.  Called by the macro PetscHeaderCreate().
 25: */
 26: PetscErrorCode  PetscHeaderCreate_Private(PetscObject h,PetscClassId classid,const char class_name[],const char descr[],const char mansec[],
 27:                                           MPI_Comm comm,PetscObjectDestroyFunction destroy,PetscObjectViewFunction view)
 28: {
 29:   static PetscInt idcnt = 1;
 30:   PetscErrorCode  ierr;
 31: #if defined(PETSC_USE_LOG)
 32:   PetscObject     *newPetscObjects;
 33:   PetscInt         newPetscObjectsMaxCounts,i;
 34: #endif

 37:   h->classid               = classid;
 38:   h->type                  = 0;
 39:   h->class_name            = (char*)class_name;
 40:   h->description           = (char*)descr;
 41:   h->mansec                = (char*)mansec;

 ...

一个 PetscObject 是一个指向：

typedef struct _p_PetscObject {
  ...
  char                 *class_name;    /*  for example, "Vec" */
  char                 *description;
  ...

因此，问题似乎与char* descr的分配有关，但我不明白为什么。大小 8 会通知指针错误不是吗？

【问题讨论】：

很难从上下文中分辨出来，但看起来 Valgrind 正在抱怨，因为您只为 PetscObject 分配了 1780 个字节，而这 8 个字节的写入（从字节 1776 开始）足以写入进入你没有分配的内存。请更新问题以显示 PetscObject 的分配方式。
有道理！ PetscObject 的分配肯定是太小了。等一下我去查
你把我引向了正确的方向。我不知道，但“所有 PETSc 对象都以 PETSCHEADER 中定义的字段开头。”。将 PETSCHEADER 添加到我的 PetscObject 中解决了我的问题。
听起来不错。然后我会发表我的评论作为答案。

标签： c valgrind petsc

【解决方案1】：

Valgrind 似乎在抱怨，因为您只为 PetscObject 分配了 1780 个字节，而这 8 个字节的写入（从字节 1776 开始）足以写入您未分配的内存。

您需要增加分配的大小，或者增加结构的大小，具体取决于分配的完整定义和方法。

【讨论】：