Building operator to extract solution at points, discrepancy?

[abhijit-c]

My apologies, I belatedly found this question which I believe resolves my issue.

How to construct a sparse matrix used for calculating function values at multiple points

Hi everyone. Recently, I try to convert my program based on FEniCS 2019 to the new version FEniCSx. I meet a problem for transforming the following code to the current version due to lack of some functions, e.g., dolfin_element.evaluate_basis_all import fenics as fe def construct_measurement_matrix(xs, V): # This function generate measurement matrix # xs: measurement points # V: function space generated by FEniCS # Example: # Let V be a function space generated by FEniCS…