In https://github.com/SciML/NeuralOperators.jl/blob/master/src/DeepONet.jl#L122-L124, multi-dimensional input grid is handled using an Argument error. Is there a reason this was done this way?
Also, if it is a matter of contracting AbstractArray(s), we could look into OMEinsum.jl for its implemention.
Please advise.