Interpolating vs Evaluating an Expression on the Mesh Coordinates
it will output an array taking the form of [..., x, x, nan, 0, ..., 0, x, x, x, ...] but when I evaluate the expression, eg:
numpy.array([j(x) for x in mesh.coordinates()])
it'll output [..., x, x, x, 0, ..., 0, nan, x, x, ...]. It just so happens that these NaNs are showing up on each of the internal boundaries depending on how I evaluate the expression. It seems a bit odd that to get the output I expect that I'd have to first interpolate, then replace the nth value with the evaluated value. I'd appreciate any insight.