How to implement multi-dimensional independent variables?

[burakatakan]

Dear all,
I just wanted to play around with chaospy. The problems, we are interested in are typically functions of 2-3 independent variables/ coordinates.
In the examples I could only find 1D examples (time).
Could you perhaps provide an example for fitting samples from a function like:

y[:] = sin (alpha * x[0, :] * exp(beta * x[1, :])

to obtain a surrogate model?

Thanks in advance!
Burak

[jonathf]

I agree that examples of multivariate coordinates are lacking.

Following the quick tutorial, how about redefining the model solver. E.g. if you want a 2D mesh, you can do something like this:

coordinates = numpy.meshgrid(numpy.linspace(0, 10, 300), numpy.linspace(30, 40, 300))

def model_solver(parameters, coordinates=coordinates):
    """Simple ordinary differential equation solver."""
    alpha, beta = parameters
    return numpy.sin(alpha * coordinates[0, :]) * numpy.exp(beta * coordinates[1, :])

[burakatakan]

Thanks a lot!