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