quadrature_rule isn't passed through to create_quadrature

[Nico Schlömer]

I'm not sure if this bug belongs in Dolfin or FFC.

When assembling a system with custom quadrature_rule in the form_compiler_parameters, e.g.,

A = assemble(a,                                                            
             form_compiler_parameters={'quadrature_rule': 'canonical',        
                                       'quadrature_degree': 1})

or

    # Solve the system.                                                         
    solve(a == L,                                                               
          sol,                                                                  
          form_compiler_parameters={'quadrature_rule': 'canonical',             
                                    'quadrature_degree': 1},                    
          bcs = bc0)

this parameter doesn't get passed through to where it is used, FFC's create_quadrature. The parameter quadrature_degree does get passed.

[Nico Schlömer]

FFC pull request in place that tentatively fixed the bug, https://bitbucket.org/fenics-project/ffc/pull-request/1/add-vertex-rule-for-quadrature.

[Jan Blechta]

This is solved, isn't it?

[Nico Schlömer]

Indeed.