initial field to dataout

[e-dinesh]

@rloekvh @exapde suggestion to print the initial field data. The following commands start printing from t=0+dt onwards, not t=0. << pde.dt = dt*ones(1,nt); % time step sizes pde.soltime = 1:length(pde.dt);

[exapde]

Exasim does not store the initial solution. I think we may want to have an option flag for users to decide whether they want to get the initial solution.

[e-dinesh]

@exapde any possible workaround (at least temporarily) ? I feel this is important for the users, who start using Exasim, to check whether their initialization is correct.

[exapde]

That's a good suggestion. We'll update the code to store the initial solution.

[exapde]

Temporarily, you can get the mesh structure and call the initu function in the pdeapp script to compute the initial solution. Alternatively, you can visualize the sequence of time-dependent solutions using paraview to see if they looks okay.