[chriss85]

I'm wondering if it is possible to have pyFoam plot residuals (or other values) during each iteration step. So far, I have only found out how to plot the values once each time step (depending on the accumulation property in customRegexp).

Alternatively, I can modify the solver so it prints time steps instead of iteration numbers, but this is not really a proper solution.

I know there are various other methods for plotting residuals, I'd just prefer to have one solution that works for all cases.

[gschaider]

Originally Posted by chriss85

I'm wondering if it is possible to have pyFoam plot residuals (or other values) during each iteration step. So far, I have only found out how to plot the values once each time step (depending on the accumulation property in customRegexp).

Alternatively, I can modify the solver so it prints time steps instead of iteration numbers, but this is not really a proper solution.

I know there are various other methods for plotting residuals, I'd just prefer to have one solution that works for all cases.

You're a bit vague about "each iteration step". Iteration of the linear solver? SIMPLE/PISO-iteration?

If you mean "I have multiple instances of solving psi per timestep and want to visualize them all": online plotting won't work. But you can write out files (from the help of pyFoamPlotRunner.py)

Data files
----------
How data files are written

--single-data-files-only
                        Don't create consecutive data files 'value',
                        'value_2', 'value_3' etc but put all the data into a
                        single file
--write-files           Writes the parsed data to files

One of the reasons why there is no plotting of multiple residuals is: how do you assign t to them on the t-axis?
- all the same time? Doesn't look good and gives no indication of progress (that is what single-data-files will give you)
- distribute evenly between this and next time? Problem is that from "outside" (and that is from where PyFoam looks) there is no way to tell how many iterations there are going to be

[chriss85]

In my case I'm talking about PIMPLE/PISO loop.

You're right about the axis of course, but I think in such a case a simple counter would suffice. Using gnuplot, one could probably go more advaned, draw vertical lines and such, but I don't think this is really needed.

It would also be useful if it was possible to specify the regular expression of the loop iteration log file output manually. This way, plotting iterations of custom loops would also be possible.

[gschaider]

Originally Posted by chriss85

In my case I'm talking about PIMPLE/PISO loop.

You're right about the axis of course, but I think in such a case a simple counter would suffice. Using gnuplot, one could probably go more advaned, draw vertical lines and such, but I don't think this is really needed.

As I said: the data can be written to file. And then you can do what you want with it.

I won't do it in online plotting because the algorithm that collects data for that and "coarsens" it (to make sure that runs with thousands of timesteps don't use all the CPU just for plotting) assumes that there is only one data item for each timestep

Originally Posted by chriss85

It would also be useful if it was possible to specify the regular expression of the loop iteration log file output manually. This way, plotting iterations of custom loops would also be possible.

I'm not quite sure what you mean with "specify manually". Of course you can target the output of the linear solver with your own customRegexp too. Or do you mean something else? Two things that MIGHT help you there are "phases" (see the ReleaseNotes) or modifying in the configuration SolverOutput/timeregexp but I'll stay as vague with the answer as you are in your question

[chriss85]

Originally Posted by gschaider

I'm not quite sure what you mean with "specify manually". Of course you can target the output of the linear solver with your own customRegexp too. Or do you mean something else? Two things that MIGHT help you there are "phases" (see the ReleaseNotes) or modifying in the configuration SolverOutput/timeregexp but I'll stay as vague with the answer as you are in your question

Indeed, the timeregexp might be what I'm looking for, as it would allow me to plot during each iteration step if I add a counting output during each iteration. I basically meant this setting, but specified for a single plot. This way should be fine though because it prevents having to modify the solver time output.

I didn't know about the configuration file before, thanks for pointing me to it!

[gschaider]

Originally Posted by chriss85

Indeed, the timeregexp might be what I'm looking for, as it would allow me to plot during each iteration step if I add a counting output during each iteration. I basically meant this setting, but specified for a single plot. This way should be fine though because it prevents having to modify the solver time output.

I didn't know about the configuration file before, thanks for pointing me to it!

OK. But don't override it globally but in LocalConfigPyFoam in the case