I did do the updates and all went through well it seems, but both with the codes of Fryderyk and Will matters crash at the solve stage; I also attache d F's code he just send me. I run on an macbook and they run under linux. The code runs for Fryderyk, he said. mtob$ python ns_plane_poiseuille.py COFFEE finished in 0.00192809 seconds (flops: 0 -> 0) COFFEE finished in 0.00323701 seconds (flops: 0 -> 0) COFFEE finished in 0.00282502 seconds (flops: 300 -> 300) pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done L2norm = 0.182574185835 compute_form_data finished in 0.203755 seconds. compile_integral finished in 0.122444 seconds. compile_integral finished in 0.047055 seconds. compile_integral finished in 0.0981019 seconds. compile_integral finished in 0.0492978 seconds. TSFC finished in 0.520766 seconds. COFFEE finished in 0.126961 seconds (flops: 27327 -> 24555) COFFEE finished in 0.0907421 seconds (flops: 3882 -> 1722) COFFEE finished in 0.0888281 seconds (flops: 3882 -> 1722) COFFEE finished in 0.087405 seconds (flops: 3882 -> 1722) compute_form_data finished in 0.177042 seconds. compile_integral finished in 0.077399 seconds. compile_integral finished in 0.028666 seconds. compile_integral finished in 0.0420549 seconds. compile_integral finished in 0.028265 seconds. TSFC finished in 0.353534 seconds. COFFEE finished in 0.046375 seconds (flops: 1196 -> 1196) compute_form_data finished in 0.196338 seconds. compile_integral finished in 0.0776689 seconds. compile_integral finished in 0.0208788 seconds. compile_integral finished in 0.034956 seconds. compile_integral finished in 0.0276241 seconds. TSFC finished in 0.357583 seconds. COFFEE finished in 0.0538301 seconds (flops: 1322 -> 1250) compute_form_data finished in 0.187414 seconds. compile_integral finished in 0.0602942 seconds. compile_integral finished in 0.0215549 seconds. compile_integral finished in 0.030211 seconds. compile_integral finished in 0.0244799 seconds. TSFC finished in 0.324083 seconds. COFFEE finished in 0.00157189 seconds (flops: 0 -> 0) COFFEE finished in 0.00145602 seconds (flops: 0 -> 0) pyop2:INFO Solving nonlinear variational problem... COFFEE finished in 0.00165892 seconds (flops: 0 -> 0) pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done COFFEE finished in 0.00231886 seconds (flops: 0 -> 0) COFFEE finished in 0.001647 seconds (flops: 0 -> 0) pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done compute_form_data finished in 0.166336 seconds. compile_integral finished in 0.0847039 seconds. compile_integral finished in 0.055573 seconds. compile_integral finished in 0.056812 seconds. compile_integral finished in 0.065223 seconds. TSFC finished in 0.428837 seconds. COFFEE finished in 0.122875 seconds (flops: 4752 -> 1616) COFFEE finished in 0.0609031 seconds (flops: 594 -> 270) COFFEE finished in 0.0570309 seconds (flops: 594 -> 270) COFFEE finished in 0.058157 seconds (flops: 594 -> 270) compute_form_data finished in 0.139532 seconds. compile_integral finished in 0.0625341 seconds. compile_integral finished in 0.022583 seconds. compile_integral finished in 0.0320051 seconds. compile_integral finished in 0.024895 seconds. TSFC finished in 0.281678 seconds. COFFEE finished in 0.0402429 seconds (flops: 842 -> 410) COFFEE finished in 0.00145102 seconds (flops: 0 -> 0) COFFEE finished in 0.00148582 seconds (flops: 0 -> 0) COFFEE finished in 0.00240588 seconds (flops: 0 -> 0) pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done compute_form_data finished in 0.167164 seconds. compile_integral finished in 0.101004 seconds. compile_integral finished in 0.053689 seconds. compile_integral finished in 0.106863 seconds. compile_integral finished in 0.0500929 seconds. TSFC finished in 0.478959 seconds. COFFEE finished in 0.11501 seconds (flops: 26991 -> 24219) COFFEE finished in 0.0758209 seconds (flops: 3786 -> 1626) COFFEE finished in 0.0740719 seconds (flops: 3786 -> 1626) COFFEE finished in 0.074147 seconds (flops: 3786 -> 1626) compute_form_data finished in 0.147847 seconds. compile_integral finished in 0.0788112 seconds. compile_integral finished in 0.028301 seconds. compile_integral finished in 0.0441599 seconds. compile_integral finished in 0.028101 seconds. TSFC finished in 0.327342 seconds. compute_form_data finished in 0.157457 seconds. compile_integral finished in 0.072794 seconds. compile_integral finished in 0.0211651 seconds. compile_integral finished in 0.0256679 seconds. compile_integral finished in 0.0287292 seconds. TSFC finished in 0.305935 seconds. COFFEE finished in 0.0498259 seconds (flops: 1178 -> 1178) compute_form_data finished in 0.14918 seconds. compile_integral finished in 0.0586369 seconds. compile_integral finished in 0.0208559 seconds. compile_integral finished in 0.0232871 seconds. compile_integral finished in 0.0236769 seconds. TSFC finished in 0.275756 seconds. pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done pyop2:INFO Compiling wrapper... pyop2:INFO Compiling wrapper...done Traceback (most recent call last): File "ns_plane_poiseuille.py", line 104, in <module> solve(F_u == 0, u, bcs = [bc0,bc1,bc2],solver_parameters={'pc_type': 'lu', 'pc_factor_mat_solver_package': 'mumps'} ) # File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving.py", line 120, in solve _solve_varproblem(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving.py", line 164, in _solve_varproblem solver.solve() File "<decorator-gen-295>", line 2, in solve File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/profiling.py", line 203, in wrapper return f(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/variational_solver.py", line 190, in solve solving_utils.check_snes_convergence(self.snes) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving_utils.py", line 62, in check_snes_convergence %s""" % (snes.getIterationNumber(), msg)) RuntimeError: Nonlinear solve failed to converge after 0 nonlinear iterations. Reason: Inner linear solve failed to converge after 0 iterations with reason: unknown reason (petsc4py enum incomplete?) Similar problem for Will's code: pyop2:INFO Compiling wrapper...done compute_form_data finished in 0.171724 seconds. compile_integral finished in 0.0246358 seconds. compile_integral finished in 0.105376 seconds. compile_integral finished in 0.145405 seconds. TSFC finished in 0.447226 seconds. Traceback (most recent call last): File "incompressible_acoustic.py", line 92, in <module> solve(DIVUdagger == DIVU, udagger, solver_parameters={'ksp_rtol': 1e-20, 'pc_type' : 'none', 'ksp_type': 'lsqr'} File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving.py", line 120, in solve _solve_varproblem(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving.py", line 147, in _solve_varproblem solver.solve() File "<decorator-gen-295>", line 2, in solve File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/profiling.py", line 203, in wrapper return f(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/variational_solver.py", line 188, in solve self.snes.solve(None, v) File "PETSc/SNES.pyx", line 520, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:172525) File "PETSc/petscsnes.pxi", line 253, in petsc4py.PETSc.SNES_Function (src/petsc4py.PETSc.c:33588) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/solving_utils.py", line 187, in form_function nest=problem._nest) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/assemble.py", line 66, in assemble inverse=inverse, nest=nest) File "<decorator-gen-296>", line 2, in _assemble File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/utils.py", line 62, in wrapper return f(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/assemble.py", line 100, in _assemble inverse=inverse) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/tsfc_interface.py", line 243, in compile_form number_map).kernels File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 203, in __new__ obj = make_obj() File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 193, in make_obj obj.__init__(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/firedrake/tsfc_interface.py", line 174, in __init__ kernels.append(KernelInfo(kernel=Kernel(ast, ast.name, opts=opts), File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/backends.py", line 118, in __call__ return t(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 203, in __new__ obj = make_obj() File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 193, in make_obj obj.__init__(*args, **kwargs) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/base.py", line 3847, in __init__ self._code = self._ast_to_c(self._ast, self._opts) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/pyop2/host.py", line 60, in _ast_to_c ast_handler.plan_cpu(self._opts) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/plan.py", line 184, in plan_cpu _generate_cpu_code(self, kernel, **params) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/plan.py", line 110, in _generate_cpu_code loop_opt.eliminate_zeros() File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/optimizer.py", line 163, in eliminate_zeros self.nz_syms = zls.reschedule(self.header) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/scheduler.py", line 844, in reschedule self.exprs.update(elf.fission(stmt, expr_info)) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/scheduler.py", line 463, in fission split, splittable = self.cutter.cut(*splittable) File "/Users/amtob/werk/firedrake/lib/python2.7/site-packages/coffee/scheduler.py", line 300, in cut index = expr_info.parent.children.index(node) ValueError: <coffee.base.Incr object at 0x10e4186d0> is not in list