Fenicsx: use an array of data for boundary conditions

[dokken]

@belmont

def f(x):
   xT = x.T
   n=x.shape[1]
   values = np.zeros(n)
   for i in range(n):
    values[i]=random.gauss(1,0.01)
    print(xT[i][0],xT[i][1])
   return values

u_bc.interpolate(f, bc_dofs)

You should remove this code, as the whole point of

@belmont

for facet in boundary_facets:
    # Find local index for each facet
    cells = f_to_c.links(facet)
    assert len(cells) == 1
    facets = c_to_f.links(cells[0])
    local_index = np.flatnonzero(facets == facet)
    # Find local index for dofs in cell
    closure_dofs = dof_layout.entity_closure_dofs(
        mesh.topology.dim-1,  local_index)
    cell_dofs = V.dofmap.cell_dofs(cells[0])
    for dof in closure_dofs:
        local_dof = cell_dofs[dof]
        dof_coordinate = coords[local_dof]
        print(local_dof, dof_coordinate)
        for b in range(V.dofmap.bs):
            num_dofs += 1
            u_bc.x.array[local_dof * V.dofmap.bs + b] = dof_coordinate[b]
        bc_dofs.append(local_dof) 

was to replace the interpolate command