Fixed BDiag bug for EBE and MFG solver types

PHASTA · May 2, 2022 · 95b6cef · 95b6cef
1 parent 779e4b5
commit 95b6cef
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/phSolver/compressible/itrdrv.f b/phSolver/compressible/itrdrv.f
@@ -177,12 +177,19 @@ subroutine itrdrv (y, ac, uold, x,
  allocate (lhsK(nflow*nflow,nnz_tot))
  endif
  endif
+ if (mod(impl(1),100)/10 .eq. 2) then
+c
+c generate the mfg data fill vectors
+c
+ allocate (BDiag(nshg,nflow,nflow))
+ endif
  if (mod(impl(1),100)/10 .eq. 3) then
 c
 c generate the ebe data fill vectors
 c
  nedof=nflow*nshape
  allocate (EGmass(numel,nedof*nedof))
+ allocate (BDiag(nshg,nflow,nflow))
  endif
 
 c..........................................