- firedrake - imperial.ac.uk

Reduce number of output files for parallel runs
by Buesing, Henrik 24 Mar '17

24 Mar '17

Dear all, Is it possible to reduce the number of output files in Firedrake when running in parallel? I am seeing an increase in the number of output files when increasing the number of processes used. Is this normal behavior? Thank you! Henrik -- Dipl.-Math. Henrik Büsing Institute for Applied Geophysics and Geothermal Energy E.ON Energy Research Center RWTH Aachen University ------------------------------------------------------ Mathieustr. 10 | Tel +49 (0)241 80 49907 52074 Aachen, Germany | Fax +49 (0)241 80 49889 ------------------------------------------------------ http://www.eonerc.rwth-aachen.de/GGE hbuesing(a)eonerc.rwth-aachen.de ------------------------------------------------------

1 0

Firedrake PETSc options_left
by Buesing, Henrik 23 Mar '17

23 Mar '17

Dear all, PETSc options "-options_left" and "-log_view" seem to have no effect when setting as solver_parameters. I can prescribe them from the command line, but then I do not get the Firedrake/PETSc "options_left". How can I print the Firedrake/PETSc options_left. Thank you! Henrik -- Dipl.-Math. Henrik Büsing Institute for Applied Geophysics and Geothermal Energy E.ON Energy Research Center RWTH Aachen University ------------------------------------------------------ Mathieustr. 10 | Tel +49 (0)241 80 49907 52074 Aachen, Germany | Fax +49 (0)241 80 49889 ------------------------------------------------------ http://www.eonerc.rwth-aachen.de/GGE hbuesing(a)eonerc.rwth-aachen.de ------------------------------------------------------

2 3

Deactivate UFL Warning
by Buesing, Henrik 23 Mar '17

23 Mar '17

Dear all, I am getting a lot of UFL:WARNING Discontinuous Lagrange element requested on quadrilateral * interval, creating DQ element Should I care? How can I deactivate these warnings? Thank you! Henrik -- Dipl.-Math. Henrik Büsing Institute for Applied Geophysics and Geothermal Energy E.ON Energy Research Center RWTH Aachen University ------------------------------------------------------ Mathieustr. 10 | Tel +49 (0)241 80 49907 52074 Aachen, Germany | Fax +49 (0)241 80 49889 ------------------------------------------------------ http://www.eonerc.rwth-aachen.de/GGE hbuesing(a)eonerc.rwth-aachen.de ------------------------------------------------------

2 6

NonLinearVariationalProblem: How to say nest=False?
by Buesing, Henrik 22 Mar '17

22 Mar '17

Dear all, In November 2015 I did problem = NonlinearVariationalProblem(F,u,nest=False,bcs=[bc0, bc1, bc2, bc3, bc4, bc5]) This is not the correct syntax anymore. How do I say nest=False today? Thank you! Henrik -- Dipl.-Math. Henrik Büsing Institute for Applied Geophysics and Geothermal Energy E.ON Energy Research Center RWTH Aachen University ------------------------------------------------------ Mathieustr. 10 | Tel +49 (0)241 80 49907 52074 Aachen, Germany | Fax +49 (0)241 80 49889 ------------------------------------------------------ http://www.eonerc.rwth-aachen.de/GGE hbuesing(a)eonerc.rwth-aachen.de ------------------------------------------------------

2 2

Re: [firedrake] NonLinearVariationalProblem: How to say nest=False?
by Andrew McRae 22 Mar '17

22 Mar '17

Per http://www.firedrakeproject.org/solving-interface.html, set "mat_type": "aij" in the solver_parameters dictionary. On 22 March 2017 at 16:05, Buesing, Henrik <HBuesing(a)eonerc.rwth-aachen.de> wrote: > Dear all, > > > > In November 2015 I did > > > > *problem = NonlinearVariationalProblem(F,u,nest=False,bcs=[bc0, bc1, bc2, > bc3, bc4, bc5])* > > > > This is not the correct syntax anymore. How do I say nest=False today? > > Thank you! > > Henrik > > > > > > > > -- > > Dipl.-Math. Henrik Büsing > > Institute for Applied Geophysics and Geothermal Energy > > E.ON Energy Research Center > > RWTH Aachen University > > ------------------------------------------------------ > > Mathieustr. 10 | Tel +49 (0)241 80 49907 > <+49%20241%208049907> > > 52074 Aachen, Germany | Fax +49 (0)241 80 49889 > <+49%20241%208049889> > > ------------------------------------------------------ > > http://www.eonerc.rwth-aachen.de/GGE > > hbuesing(a)eonerc.rwth-aachen.de > > ------------------------------------------------------ > > >

1 0

Guidelines for using/running jobs on larger HPC machines
by Justin Chang 22 Mar '17

22 Mar '17

Hi all, Generally speaking, what should the PYOP2/FIREDRAKE_TSFC_KERNEL_CACHE_DIR environment variables be set to? Is it sufficient to have something like: export PYOP2/FIREDRAKE_TSFC_KERNEL_CACHE_DIR=$HOME/fd-cache This works on my smaller university cluster but I wonder if for a system like Edison at NERSC ( http://www.nersc.gov/users/computational-systems/edison/) if there is a better directory for this. Also, from what I read, Edison's SLURM scheduler loads the executable to the allocated compute nodes from the current working directory, which supposedly can be really slow - they recommend using something like: srun --bcast=/tmp/$SLURM_JOB_ID --compress=lz4 ... if 2000 or more nodes are needed. But even on jobs that require no more than a single compute node (24 cores), the firedrake/python modules seem to load very slowly. This is the description of the file storage system on Edison ( http://www.nersc.gov/users/computational-systems/edison/file-storage-and-i-…) btw Any help or thoughts appreciated. Justin

2 1

Re: [firedrake] firedrake-install on macOS Sierra
by David Ham 21 Mar '17

21 Mar '17

I think we've seen this occasionally on new macs. I think it's caused by the virtualenv getting installed in a different python installation from the one which is being used to install Firedrake. Try the following: Open a new terminal (this is to try to avoid any bash path caching effects) Type "which python". The answer you want is /usr/local/bin/python run python and check that you can import virtualenv If all the above works, try running firedrake-install in the new shell. If anything above fails, tell us. Regards, David On Tue, 21 Mar 2017 at 10:28 Anna Kalogirou <A.Kalogirou(a)leeds.ac.uk> wrote: > Dear all, > > > I am installing firedrake on my new laptop from scratch (running macOS > Sierra). > > > The installation script runs and installs virtualenv, up to the point > where I get this: > > Virtual env installed. Please run firedrake-install again. > > > But when I run > > python firedrake-install > > again, it just does the same thing (similar to the previous step) and does > not install firedrake (and petsc etc..). What shall I do? > > > Best, Anna. > -- Dr David Ham Department of Mathematics Imperial College London

2 2

firedrake-install on macOS Sierra
by Anna Kalogirou 21 Mar '17

21 Mar '17

Dear all, I am installing firedrake on my new laptop from scratch (running macOS Sierra). The installation script runs and installs virtualenv, up to the point where I get this: Virtual env installed. Please run firedrake-install again. But when I run python firedrake-install again, it just does the same thing (similar to the previous step) and does not install firedrake (and petsc etc..). What shall I do? Best, Anna.

1 0

Re: [firedrake] circles in firedrake?
by David Ham 21 Mar '17

21 Mar '17

Hi Francis, There appears to be something called vext.vtk which claims to make system vtk accessible from inside a virtualenv. There is also this email: http://vtk.1045678.n5.nabble.com/Re-VTK-in-virtualenv-td5716725.html I have to confess that I have not really tried either. Regards, David On Mon, 20 Mar 2017 at 01:40 Francis Poulin <fpoulin(a)uwaterloo.ca> wrote: > > Thanks for suggesting PyGmsh! I bumped into it a while ago but never > really figured it out. Now I'm keen since it should help building and > modifying meshes more easily. > > I installed pygmsh easily using the instructions on the website. It > imports correctly. There is even an example that builds a circle, I think. > > When I try running this, as well as other examples, I get an error because > no module named vtk is found. See below. I get this comes from meshio? > > Traceback (most recent call last): > > File "circle.py", line 24, in <module> > meshio.write('circle.vtu', *out) > File > "/home/fpoulin/software/firedrake/local/lib/python2.7/site-packages/meshio/__init__.py", > line 187, in write > field_data=field_data > File > "/home/fpoulin/software/firedrake/local/lib/python2.7/site-packages/meshio/vtk_io.py", > line 204, in write > import vtk > ImportError: No module named vtk > > When I do a pip install it tells me that it is already installed and up to > date. > > > (firedrake) fpoulin@fpoulin-Gazelle:~/software/pygmsh/test/examples$ pip > install -U meshio > Requirement already up-to-date: meshio in > /home/fpoulin/software/firedrake/lib/python2.7/site-packages > Requirement already up-to-date: pipdated in > /home/fpoulin/software/firedrake/lib/python2.7/site-packages (from meshio) > Requirement already up-to-date: numpy in > /home/fpoulin/software/firedrake/lib/python2.7/site-packages (from meshio) > Requirement already up-to-date: requests in > /home/fpoulin/software/firedrake/lib/python2.7/site-packages (from > pipdated->meshio) > Requirement already up-to-date: appdirs in > /home/fpoulin/software/firedrake/lib/python2.7/site-packages (from > pipdated->meshio) > > > I realize this isn't necessarily a firedrake question but if anyone had an > idea on how to get this working I would be keen to try it. > > > Cheers, Francis > ------------------ > Francis Poulin > Associate Professor > Department of Applied Mathematics > University of Waterloo > > email: fpoulin(a)uwaterloo.ca > Web: https://uwaterloo.ca/poulin-research-group/ > Telephone: +1 519 888 4567 x32637 <(519)%20888-4567> > > ------------------------------ > *From:* firedrake-bounces(a)imperial.ac.uk [firedrake-bounces(a)imperial.ac.uk] > on behalf of G. D. McBain [gdmcbain(a)protonmail.com] > *Sent:* Sunday, March 19, 2017 6:16 PM > > *To:* firedrake(a)imperial.ac.uk > *Subject:* Re: [firedrake] circles in firedrake? > > I would expect a user to run gmsh separately (at the command line or > through the GUI), rather than from within the Python script. I don't think > it's possible to interact with gmsh from Python nowadays without jumping > through considerable hoops. > > > Have you tried the (free third-party) PyGmsh? It generates .geo code from > Python. From there one calls Gmsh using Python subprocess and reads in the > resulting .msh. I've found it handy. > > https://github.com/nschloe/pygmsh > > -- Dr David Ham Department of Mathematics Imperial College London

3 2

Mismatching function spaces error
by Anna Kalogirou 20 Mar '17

20 Mar '17

Dear all, I am trying to solve a set of PDEs and ODEs together as a mixed system. The ODE variables only depend on time, so I define them as functions from the real function space R. First question: Is this a correct approach? Second question: I am running into a “Mismatching function spaces” error; can someone please take a look at my code (found here<https://bitbucket.org/annakalog/buoy2d/src/957e88739223e11d84826c76923a9b07…>) and let me know what I’m doing wrong? The error arises when the value of one of the mixed-system variables is assigned to another variable which belongs to the same function space. [I remember having the same issue some time ago, but I thought that it was fixed.] Thank you. Best, Anna.

3 6