Thank you very much Andrew. That fixed my problems!

Anna.

P.S. Happy holidays and Merry Christmas to all of you :)


On 22 Dec 2014, at 14:22, Eike Mueller <e.mueller@bath.ac.uk<mailto:e.mueller@bath.ac.uk>> wrote:

Oh, I should have done what Andrew suggests before sending my email...

op2.init()
parameters["coffee"]["O2"] = False

seems to fix at least part of my problem :-)

Sorry for the noise,

Eike

--

Dr Eike Hermann Mueller
Research Associate (PostDoc)

Department of Mathematical Sciences
University of Bath
Bath BA2 7AY, United Kingdom

+44 1225 38 5803
e.mueller@bath.ac.uk<mailto:e.mueller@bath.ac.uk>
http://people.bath.ac.uk/em459/

On 22 Dec 2014, at 13:58, Andrew McRae <a.mcrae12@imperial.ac.uk<mailto:a.mcrae12@imperial.ac.uk>> wrote:

Ugh, one of the recent changes we made was to turn on COFFEE's optimisations by default.  Unfortunately, this has led to some teething problems.

I assume that all components -- COFFEE, Firedrake, PyOP2, FFC, UFL, and FIAT -- are up to date, and that you're running with a clean cache (from the Firedrake base directory, run "./scripts/firedrake-clean").

If the problem persists, as a short-term workaround, can you add

op2.init()
parameters["coffee"]["O2"] = False

to the start of your file, just after the imports?  This will result in sub-optimal performance during the assembly stage, so isn't ideal in the long-term.

Fabio, our resident barista, will have some fixes after the new year, I'm sure.

Andrew

On 22 December 2014 at 13:47, Anna Kalogirou <A.Kalogirou@leeds.ac.uk<mailto:A.Kalogirou@leeds.ac.uk>> wrote:
Dear all,

I very recently updated firedrake, pyop2 and installed coffee, but since then none of my codes include nonlinear solvers don’t work (they used to work fine before the update).

I even tried to run the burgers.py code from the demo and even that diverges.

I get something like:

Traceback (most recent call last):
  File "benney-luke.py", line 57, in <module>
    phi_solver = solver_phi(phi1, phi1, phi0, eta0, etaR, phi, gamma, dt, mu, epsilon, is_linear);
  File "/Users/matak/Documents/Firedrake/modules/solvers.py", line 26, in solver_phi
    phi_solver = NonlinearVariationalSolver(phi_problem,solver_parameters={'snes_monitor': True,'ksp_monitor': True,'snes_linesearch_monitor': True})
  File "/Users/matak/firedrake/firedrake/variational_solver.py", line 97, in __init__
    form_compiler_parameters=self._problem.form_compiler_parameters)
  File "/Users/matak/firedrake/firedrake/assemble.py", line 58, in assemble
    form_compiler_parameters=form_compiler_parameters)
  File "/Users/matak/firedrake/firedrake/assemble.py", line 79, in _assemble
    kernels = ffc_interface.compile_form(f, "form", parameters=form_compiler_parameters)
  File "/Users/matak/firedrake/firedrake/ffc_interface.py", line 268, in compile_form
    ffc_kernel = FFCKernel(f, name + str(i) + str(j), parameters)
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/caching.py", line 202, in __new__
    obj = make_obj()
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/caching.py", line 192, in make_obj
    obj.__init__(*args, **kwargs)
  File "/Users/matak/firedrake/firedrake/ffc_interface.py", line 194, in __init__
    (name, it.integral_type(), it.subdomain_id()), opts, inc),
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/backends.py", line 118, in __call__
    return t(*args, **kwargs)
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/caching.py", line 202, in __new__
    obj = make_obj()
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/caching.py", line 192, in make_obj
    obj.__init__(*args, **kwargs)
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/base.py", line 3632, in __init__
    self._code = self._ast_to_c(code, opts)
  File "/usr/local/Cellar/python/2.7.8_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/PyOP2-0.11.0_148_g8a1955c-py2.7-macosx-10.9-x86_64.egg/pyop2/host.py", line 63, in _ast_to_c
    ast_handler.plan_cpu(opts)
  File "/usr/local/lib/python2.7/site-packages/coffee/plan.py", line 405, in plan_cpu
    loop_opts = _generate_cpu_code(self, **params)
  File "/usr/local/lib/python2.7/site-packages/coffee/plan.py", line 241, in _generate_cpu_code
    loop_opt.rewrite(rewrite)
  File "/usr/local/lib/python2.7/site-packages/coffee/optimizer.py", line 110, in rewrite
    ew.expand()
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 168, in expand
    ee.expand(self.stmt.children[1], self.stmt, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 565, in expand
    return self.expand(node.children[0], node, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 567, in expand
    l_node, l_type = self.expand(node.children[0], node, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 565, in expand
    return self.expand(node.children[0], node, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 599, in expand
    l_node, l_type = self.expand(node.children[0], node, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 565, in expand
    return self.expand(node.children[0], node, itvars_occs, itvar_exp)
  File "/usr/local/lib/python2.7/site-packages/coffee/rewriter.py", line 608, in expand
    raise RuntimeError("Expansion error: unknown node: %s" % str(node))
RuntimeError: Expansion error: unknown node: (LI_IPK_1_0[k] * LI_IPJ_1_1[j]) / LI_C_1_0


Thank you.

Anna.



_______________________________________________
firedrake mailing list

firedrake@imperial.ac.uk<mailto:firedrake@imperial.ac.uk>
https://mailman.ic.ac.uk/mailman/listinfo/firedrake

_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk<mailto:firedrake@imperial.ac.uk>
https://mailman.ic.ac.uk/mailman/listinfo/firedrake

_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk<mailto:firedrake@imperial.ac.uk>

https://mailman.ic.ac.uk/mailman/listinfo/firedrake


_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake