First of all thanks for the wonderful framework. I am able to run the material constants (such as elastic modulus, thermal conductivity, etc) as a function of space and time only in the functional form (sympy expression). Can you please provide any suggestion on whether can we give this data in a discrete form or some tabular format? Thanks.

For discrete data I think the easiest route is to define another constraint from numpy dictionaries that would contain the discrete inputs and outputs. Most of the constraints you’re using are usually designed for a continuous representation via SymPy (which is what you have).

This is a pretty useful utility to build constraints from a fixed dataset. Some of the wave examples use this function. You will need to create an input/target dataset of sampled boundary points but this should be achievable manually or using the geometry module.