decide upon Lagrange monomial evaluation

currently the Lagrange polynomial class allows evaluating all monomials (= all individual "Lagrange polynomials" of the basis) separately.

This is required for the current version of the optimiser in order to perform a "per monomial" gradient descent (e.g. with additional uncertainty). However this use case might not be general enough to keep this functionality so closely connected to the core classes/code.

Would it be useful to implement a "monomial evaluation" function for the polynomial classes of all the different bases, e.g. evaluate Lagrange/Newton/classical monomials depending on the basis in use?