[IGA] Expose active CP indices and higher-order derivatives at arbitrary local coordinates for BRep Surfaces

[juancamarotti]

Motivation

Currently, Kratos allows evaluating shape functions (and first derivatives) at arbitrary local coordinates on BRepSurfaces through ShapeFunctionsValues() and ShapeFunctionsLocalGradients(). However:

• The active control-point indices at a given parameter position are not accessible (only available through QuadraturePointGeometry).
• For non-quadrature points, derivatives are limited to first order via ShapeFunctionLocalGradients(), preventing access to higher-order derivatives.

What this PR adds
Extends the BrepSurface class with a routine called ShapeFunctionsValuesAndCPIndices() which returns the following values at a certain point in the surface parameter space:

• Active control-point indices
• Shape function values
• Higher-order derivatives (up to requested order)

Comments

• A test for this new functionality was added
• This new function was exposed to the python interface

[rickyaristio]

I agree with the new function. Please check my comments.