dolfin team mailing list archive

[Jake Ostien <tostien@xxxxxxxxx>]

Hi,

In Assembler::assembleCells there is and assignment/logical check that 
looks like

if (uint domain = (*domains)(*cell) < ufc.form.num_cell_integrals())
   integral = ufc.cell_integrals[domain];

Near as I can tell, if the inequality is true then domain = 1, which is 
not the intended behavior.  A simple fix might be

uint domain = (*domains)(*cell);
if (domain < ufc.form.num_cell_integrals())
   integral = ufc.cell_integrals[domain];

Now this assigns to domain the value of the MeshFunction<uint> domains 
for *cell (which may not be 1).

Jake