ugly init but working

ref idaholab#15204
tophmatthews · May 27, 2020 · 2bac9b6 · 2bac9b6
1 parent 7961ede
commit 2bac9b6
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/framework/src/problems/FEProblemBase.C b/framework/src/problems/FEProblemBase.C
@@ -1658,6 +1658,8 @@ FEProblemBase::reinitDirac(const Elem * elem, THREAD_ID tid)
         _scalar_zero[tid].resize(FORTYTHIRD, 0);
         _zero[tid].resize(max_qpts, 0);
         _grad_zero[tid].resize(max_qpts, RealGradient(0.));
+        _grad_phi_zero[tid].resize(max_qpts, {0., 0., 0.});
+        _grad_test_zero[tid].resize(max_qpts, {0., 0., 0.});
         _second_zero[tid].resize(max_qpts, RealTensor(0.));
         _vector_zero[tid].resize(max_qpts, RealGradient(0.));
         _vector_curl_zero[tid].resize(max_qpts, RealGradient(0.));
@@ -4642,6 +4644,8 @@ FEProblemBase::createQRules(QuadratureType type, Order order, Order volume_order
     _zero[tid].resize(max_qpts, 0);
     _ad_zero[tid].resize(max_qpts, 0);
     _grad_zero[tid].resize(max_qpts, RealGradient(0.));
+    _grad_phi_zero[tid].resize(max_qpts, {0., 0., 0.});
+    _grad_test_zero[tid].resize(max_qpts, {0., 0., 0.});
     _ad_grad_zero[tid].resize(max_qpts, ADRealGradient(0));
     _second_zero[tid].resize(max_qpts, RealTensor(0.));
     _ad_second_zero[tid].resize(max_qpts, ADRealTensorValue(0));