Does pdepe accept 'Vectorized', 'on' as an option? If so, what does it do?

I am writing a code that uses pdepe. It would be easy to write the pdefun in vectorized form. That is, if there are N variables and M position points, if the position x is input as a vector, then c, f, and s would be output as NxM arrays. It seems like this would be more efficient than only allowing x to be a scalar. However, odeset does have a "Vectorized" option. Does this option apply to pdepe? If so, what does it do? Does it allow the input variable x to be a vector rather than a scalar?

Thanks for your help!

Answer 1

[Torsten]

If you look into the documentation for "pdepe", you will find that the "Vectorized" option cannot be chosen:

Cited from the documentation of "pdepe":

options — Option structure

structure array

Option structure, specified as a structure array. Use the odeset function to create or modify the option structure. pdepe supports these options:

• 'RelTol'
• 'AbsTol'
• 'NormControl'
• 'InitialStep'
• 'MaxStep'
• 'Events'

In most cases, default values for these options provide satisfactory solutions.

If I remember correctly, this variant of "pdepe" has the choice to vectorize the calls to the function where you provide c, s and f:

https://uk.mathworks.com/matlabcentral/fileexchange/97437-pde1dm

[Bill Greene]

Yes, pde1dm provides a "Vectorized=on" option that can substantially improve performance for refined meshes. I don't know why pdepe does not also support this. (There are several other ode15s options that could also be supported by pdepe but, unfortunately, are not.)