Firedrake meeting 2020 10 12

Date and time 2020-10-05 15:00UTC (16:00BST)

Action Items

1. Pick Chair and Minuter.
2. ALL: (ongoing) triage the open issues and confirm if they are indeed still open (and perhaps provide labels)
3. ??: Build master on centos to catch errors
4. RK: Report back on quadrature estimation
5. KS: Complex merge

Agenda

Present:

Apologies: Reuben Nixon-Hill (note new surname!)

JB: Numpy bug

Finally successfully installed on MacOS, we have working install script now.

Can also import Firedrake, but we need a hack:

export DYLD_INSERT_LIBRARIES=/usr/local/opt/openblas/lib/libblas.dylib:/usr/local/opt/openblas/lib/liblapack.dylib

We should discuss whether we are happy with such a hack.

KS: Subspace

Refactor tsfc/driver.py

• gem construction
• impero_c construction
• kernel construction

etc.

firedrake/tsfc_interface.py

Ex: Jacobian assembly:

all_gem_expr = []
split form according to indices (i, j) -> subform
    split subform according to test/trial subspaces (test_sub, trial_sub) -> subsubform
       gem_expr = compile_ufl(subsubform, ...)
  ->   gem_expr = transform(gem_expr, test_sub, trial_sub) (test_sub and trial_sub know right transformation rule, such as rotation)
       all_gem_expr.append(gem_expr)
    kernel = construct_kernel(all_gem_expr)

Question:

Masked is still a ufl class (now defined in Firedrake), but disappear before calling compile_ufl. This is convenient when we take action and derivative, but not sure if this is a good practice.

Date of next meeting

Discuss changing meeting to Wednesday early at Gebina Ham's request. 2020-10-19 15:00UTC (16:00BST)