Is there any way to figure out what NDSolve default spatial discretization is for a given set of initial conditions? The notes state that it is "based on the initial conditions", but I cannot figure out a set of options that will let me see spatial locations at which the PDE is actually being solved.
Motivation: The default discretization is yielding results that I cannot replicate using "Coordinates" to fix my discretization.
Any assistance is appreciated.
First to address your question: try to have a look at Menu/Help/Advanced Numerical Differential Equation Solving in Mathematica/Partial Differential Equations/Controlling the Spatial Grid Selection You may fund some useful hints there.
I would like, however, to ask a further question.
At the Mathematica conference in London this Summer Mr. Stephen Wolfram informed us that PDE solving by NDSolve will be included in Mathematica 9. Mr. Wolfram highlighted this point in his talk, when he explained how much M9 differs from M8. I understood this message the way that some important advancements are achieved in this area.
Just to remind, Mathematica 8 was already able to solve certain types of PDEs. It used the so-called, Method of Lines. The latter may only be applied, provided the PDE may be formulated as initial+boundary problem. A natural example of such a formulation is the time-dependent diffusion equation. M8 was unable to solve pure boundary problems, such as, for example, the Laplace equation.
My question is: if something has changed in M9 with the PDE solution and what the changes are? Does it apply something other method(s) in addition to the Method of Lines?
I checked the help. On one hand I find no indication that only the Method of Lines may be applied. On the other hand, all PDE examples in the Help/NDSolve are those I have already seen in M8.
Besides, one can directly check that Laplace equation is still not understood by NDSolve: