eikehmueller

from firedrake import *

n = 32

mesh = UnitSquareMesh(n,n)
V_u = FunctionSpace(mesh,'RT',1)
V_phi = FunctionSpace(mesh,'CG',1)
expression = Expression('(1.-cos(2.*pi*x[0]))*(1.-cos(2.*pi*x[1]))')

phi = Function(V_phi)

eikehmueller

Traceback (most recent call last):
  File "divide.py", line 23, in <module>
    par_loop('for(int i=0; i<f3.dofs; i++) { f3[i][0] = f1[i][0]/f2[i][0]; }',dS,{'f3':(f3,WRITE), 'f1':(f1,READ), 'f2':(f2,READ)})
  File "/Users/eikemueller/PostDocBath/EllipticSolvers/firedrake/firedrake/parloops.py", line 137, in par_loop
    op2args.append(_map['itspace'](mesh, measure))
  File "/Users/eikemueller/PostDocBath/EllipticSolvers/firedrake/firedrake/parloops.py", line 18, in <lambda>
    'itspace': lambda mesh, measure: mesh.interior_facets.measure_set(measure)
TypeError: measure_set() takes exactly 3 arguments (2 given)