From 2935c17f781ef43c20f34fc6e2cac4f30d98bcd5 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 10 Oct 2024 05:19:32 -0400
Subject: [PATCH] avoid segfault in msm/dielectric

---
 src/DIELECTRIC/msm_dielectric.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/DIELECTRIC/msm_dielectric.cpp b/src/DIELECTRIC/msm_dielectric.cpp
index fd09e40acc..46bb9c0415 100644
--- a/src/DIELECTRIC/msm_dielectric.cpp
+++ b/src/DIELECTRIC/msm_dielectric.cpp
@@ -113,9 +113,11 @@ void MSMDielectric::compute(int eflag, int vflag)
   if (atom->nmax > nmax) {
     memory->destroy(part2grid);
     memory->destroy(efield);
+    memory->destroy(phi);
     nmax = atom->nmax;
     memory->create(part2grid,nmax,3,"msm:part2grid");
     memory->create(efield,nmax,3,"msm:efield");
+    memory->create(phi,nmax,"msm:phi");
   }
 
   // find grid points for all my particles