Merge branch 'master' of https://github.com/lehner/gpt

lehner · Oct 26, 2023 · 6800f15 · 6800f15
2 parents bd22300 + d005514
commit 6800f15
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/benchmarks/stencil_tensor.py b/benchmarks/stencil_tensor.py
@@ -2,9 +2,9 @@
 import gpt as g
 
 #grid = g.grid([64,64,64,64], g.double)
-#grid = g.grid([32,32,32,32], g.double)
+grid = g.grid([32,32,32,32], g.double)
 #grid = g.grid([32,16,16,16], g.double)
-grid = g.grid([16,16,16,32], g.double)
+#grid = g.grid([16,16,16,32], g.double)
 #grid = g.grid([2*4,4*3,3*4,3*3*4], g.double)
 m1 = g.mcolor(grid)
 m2 = g.mcolor(grid)

diff --git a/lib/cgpt/lib/eval.cc b/lib/cgpt/lib/eval.cc
@@ -299,6 +299,8 @@ void eval_general(std::vector<cgpt_Lattice_base*>& dst, std::vector<_eval_term_>
 
 static inline void simplify(_eval_term_& term) {
   auto& factors = term.factors;
+  if (!factors.size())
+    return;
   for (size_t i=0;i<factors.size()-1;) {
     if ((factors[i  ].type == _eval_factor_::GAMMA) &&
 	(factors[i+1].type == _eval_factor_::GAMMA)) {
@@ -345,7 +347,7 @@ EXPORT(eval,{
     // if (expr_class_prod()) 
     //   eval_prod()
     // else
-    
+
     // General code path:
     eval_general(dst,terms,unary,ac);