In line 139 you have an Expression, but H0 is undefined. As just demonstrated, making mistakes in Expression C strings produces unhelpful error messages. Therefore users are strongly encouraged to use UFL expressions instead. For example, Expression("0.5*(1+copysign(1.0,x[0]-xb))*slope*(x[0]-xb)",xb=xb, slope=sb) can be rewritten as 0.5*(1+sign(x[0]-xb))*sb*(x[0]-xb) Cheers, Miklos On 05/05/17 12:18, Floriane Gidel [RPG] wrote:
Dear Miklos,
You can find my code attached. The code to run is 'main.py'.
Firedrake is up-to-date (from this morning) but it doesn't fix the error.
I solve equations for h_n1 and hat_psi_star simultaneously (in the h_solver) with the nonlinear solver, where hat_psi_star is a vector function. Then I solve the equations for psi_1_n1 with the linear solver psi_solver.
I also solve an linear equation for hat_psi_n1 (alone), and that one works (hat_psi_solver).
Thanks,
Floriane
------------------------------------------------------------------------ *De :* firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> de la part de Miklós Homolya <m.homolya14@imperial.ac.uk> *Envoyé :* vendredi 5 mai 2017 12:08:12 *À :* firedrake@imperial.ac.uk *Objet :* Re: [firedrake] Compilation error
Dear Floriane,
You need to send the code that triggers this problem since we cannot immediately see from this error message what is going wrong.
Also make sure that your Firedrake installation is up-to-date, and that you can reproduce the problem with the latest Firedrake.
Thanks,
Miklos
On 05/05/17 12:01, Floriane Gidel [RPG] wrote:
Dear all,
I try to solve a weak form with the NonlinearVariationalSolver but get the error below and couldn't figure out what causes this error. I get the same error when solving another weak form with the LinearVariationalSolver, but my third weak form (also solved with LinearVariationalSolver) works. I have updated firedrake but it didn't solve the issue.
Does it ring a bell to someone ? If it is not an "obvious" error, I can also send my code.
Thanks a lot,
Floriane
Traceback (most recent call last):
File "main.py", line 366, in <module>
h_solver.solve()
File "/Users/mmfg/firedrake/src/firedrake/firedrake/variational_solver.py", line 220, in solve
self.snes.solve(None, v)
File "PETSc/SNES.pyx", line 537, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:172544)
File "PETSc/petscsnes.pxi", line 261, in petsc4py.PETSc.SNES_Function (src/petsc4py.PETSc.c:33927)
File "/Users/mmfg/firedrake/src/firedrake/firedrake/solving_utils.py", line 433, in form_function
with ctx._F.dat.vec_ro as v:
File "/usr/local/Cellar/python/2.7.13/Frameworks/Python.framework/Versions/2.7/lib/python2.7/contextlib.py", line 17, in __enter__
return self.gen.next()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/petsc_base.py", line 395, in vecscatter
with acc(d) as v:
File "/usr/local/Cellar/python/2.7.13/Frameworks/Python.framework/Versions/2.7/lib/python2.7/contextlib.py", line 17, in __enter__
return self.gen.next()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/petsc_base.py", line 328, in vec_context
self._force_evaluation(read=True, write=not readonly)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/base.py", line 1681, in _force_evaluation
_trace.evaluate(reads, writes)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/base.py", line 193, in evaluate
comp._run()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/base.py", line 4085, in _run
return self.compute()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/base.py", line 4127, in compute
fun = self._jitmodule
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/utils.py", line 64, in __get__
obj.__dict__[self.__name__] = result = self.fget(obj)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/sequential.py", line 805, in _jitmodule
pass_layer_arg=self._pass_layer_arg)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/caching.py", line 200, in __new__
obj = make_obj()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/caching.py", line 190, in make_obj
obj.__init__(*args, **kwargs)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/sequential.py", line 635, in __init__
self.compile()
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/sequential.py", line 716, in compile
comm=self.comm)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/compilation.py", line 385, in load
dll = compiler.get_so(src, extension)
File "/Users/mmfg/firedrake/src/PyOP2/pyop2/compilation.py", line 227, in get_so
Compile errors in %s""" % (e.cmd, e.returncode, logfile, errfile))
pyop2.exceptions.CompilationError: Command "['mpicc', '-std=c99', '-fPIC', '-Wall', '-framework', 'Accelerate', '-march=native', '-O3', '-I/Users/mmfg/firedrake/lib/python2.7/site-packages/petsc/include', '-I/Users/mmfg/firedrake/src/PyOP2/pyop2', '-msse', '-o', '/tmp/pyop2-cache-uid1010350243/c993a8a607879888f3be0da8f3d4009d_p95228.so.tmp', '/tmp/pyop2-cache-uid1010350243/c993a8a607879888f3be0da8f3d4009d_p95228.c', '-dynamiclib', '-L/Users/mmfg/firedrake/lib/python2.7/site-packages/petsc/lib', '-Wl,-rpath,/Users/mmfg/firedrake/lib/python2.7/site-packages/petsc/lib', '-lpetsc', '-lm']" return error status 1.
Unable to compile code
Compile log in /tmp/pyop2-cache-uid1010350243/c993a8a607879888f3be0da8f3d4009d_p95228.log
Compile errors in /tmp/pyop2-cache-uid1010350243/c993a8a607879888f3be0da8f3d4009d_p95228.err
_______________________________________________ firedrake mailing list 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