Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ft07 does not seem to work on Linux #69

Open
scopatz opened this issue Feb 6, 2020 · 3 comments
Open

ft07 does not seem to work on Linux #69

scopatz opened this issue Feb 6, 2020 · 3 comments

Comments

@scopatz
Copy link

scopatz commented Feb 6, 2020

I with conda install -c conda-forge fenics jupyterlab matplotlib, I am seeing the following:

$ python ft07_navier_stokes_channel.py                  
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Traceback (most recent call last):
  File "ft07_navier_stokes_channel.py", line 118, in <module>
    error = np.abs(u_e.vector().array() - u_.vector().array()).max()
AttributeError: 'dolfin.cpp.la.PETScVector' object has no attribute 'array'
@gdmcbain
Copy link

gdmcbain commented Feb 7, 2020

    error = np.abs(u_e.compute_vertex_values() -
                   u_.compute_vertex_values()).max()

There's a minimally modifying version of the whole script that works in FEniCS 2019.1.0 at:

https://github.com/gdmcbain/fenics-tuto-in-skfem/blob/119258655559ab7ac6b563bf75be6923d494d2be/07_navier_stokes_channel/ft07_navier_stokes_channel1.py#L100-L101

@scopatz
Copy link
Author

scopatz commented Feb 12, 2020

Thanks. That should probably be in master

@gdmcbain
Copy link

You're welcome.

I'm still waiting for the merger of #49 (it'll be a year on Friday fortnight), so I'm not sure of the protocol for getting it into master.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants