Sorry, I sent the wrong file. Please find attached the file that gave me the error: Traceback (most recent call last): File "<stdin>", line 21, in <module> File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 110, in __lshift__ self._file << data File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 524, in __lshift__ self._update_PVD(data) File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 542, in _update_PVD new_vtk << function File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 180, in __lshift__ mesh = function.function_space().mesh() AttributeError: 'Sum' object has no attribute 'function_space' Angwenyi

On 19 Feb 2016, at 15:04, Andrew McRae <A.T.T.McRae@bath.ac.uk> wrote:

Whoops, ignore, email sync failure.

On 19 Feb 2016 14:03, "Andrew McRae" <A.T.T.McRae@bath.ac.uk <mailto:A.T.T.McRae@bath.ac.uk>> wrote: What is in the .err file and the corresponding .c file?

On 19 Feb 2016 13:15, "Angwenyi David" <kipkoej@gmail.com <mailto:kipkoej@gmail.com>> wrote: Hi all,

I still have a problem with the time loop. I get the following error.

Traceback (most recent call last): File "<stdin>", line 3, in <module> File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/constant.py", line 116, in assign self.dat.data = value File "<decorator-gen-21>", line 2, in data File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/versioning.py", line 115, in modifies retval = method(self, *args, **kwargs) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/base.py", line 2761, in data _trace.evaluate(set(), set([self])) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/base.py", line 170, in evaluate comp._run() File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/base.py", line 4070, in _run return self.compute() File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/base.py", line 4111, in compute fun = self._jitmodule File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/utils.py", line 64, in __get__ obj.__dict__[self.__name__] = result = self.fget(obj) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/sequential.py", line 151, in _jitmodule direct=self.is_direct, iterate=self.iteration_region) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 203, in __new__ obj = make_obj() File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/caching.py", line 193, in make_obj obj.__init__(*args, **kwargs) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/host.py", line 700, in __init__ self.compile() File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/host.py", line 801, in compile compiler=compiler.get('name')) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/compilation.py", line 292, in load dll = compiler.get_so(src, extension) File "/Users/master/firedrake/lib/python2.7/site-packages/pyop2/compilation.py", line 146, 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/master/firedrake/lib/python2.7/site-packages/petsc/include', '-I/Users/master/firedrake/lib/python2.7/site-packages/pyop2', '-msse', '-o', '/var/folders/8q/tv2c50fj0n3_27j_3yjmyv7c0000gp/T/pyop2-cache-uid502/59ddbd69d6d404a58876bfb1da21738a_p675.so.tmp', '/var/folders/8q/tv2c50fj0n3_27j_3yjmyv7c0000gp/T/pyop2-cache-uid502/59ddbd69d6d404a58876bfb1da21738a_p675.c', '-dynamiclib', '-L/Users/master/firedrake/lib/python2.7/site-packages/petsc/lib', '-Wl,-rpath,/Users/master/firedrake/lib/python2.7/site-packages/petsc/lib', '-lpetsc', '-lm']" return error status 1. Unable to compile code Compile log in /var/folders/8q/tv2c50fj0n3_27j_3yjmyv7c0000gp/T/pyop2-cache-uid502/59ddbd69d6d404a58876bfb1da21738a_p675.log Compile errors in /var/folders/8q/tv2c50fj0n3_27j_3yjmyv7c0000gp/T/pyop2-cache-uid502/59ddbd69d6d404a58876bfb1da21738a_p675.err

How can I fix it?

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

Hi Lawrence, Phin is the value of phi at time step n. Perhaps I should do phi = phin within the time loop and output the result as phi. D.

On 19 Feb 2016, at 16:20, Lawrence Mitchell <lawrence.mitchell@imperial.ac.uk> wrote:

Hi David,

...

On 19 Feb 2016, at 15:05, Angwenyi David <kipkoej@gmail.com> wrote:

<WaveEquation.py> Sorry, I sent the wrong file. Please find attached the file that gave me the error:

Traceback (most recent call last): File "<stdin>", line 21, in <module> File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 110, in __lshift__ self._file << data File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 524, in __lshift__ self._update_PVD(data) File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 542, in _update_PVD new_vtk << function File "/Users/master/firedrake/lib/python2.7/site-packages/firedrake/io.py", line 180, in __lshift__ mesh = function.function_space().mesh() AttributeError: 'Sum' object has no attribute 'function_space'

You're trying to output "phin", but inside the time loop you reassign the name phin to an expression:

phin = phin_1 - (dt * theta * pn + dt * (1 - theta) * pn_1)

This is now no longer a Function.

I'm confused about what you're trying to do here. Do you want phin to take the point wise value of the right hand side?

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

Thanks a lot. It is now working. Cheers. David

On 19 Feb 2016, at 16:36, Lawrence Mitchell <lawrence.mitchell@imperial.ac.uk> wrote:

phin.assign(phin_1 - (dt * theta * pn + dt * (1 - theta) * pn_1))