Hi Piotr, On 27/04/2016 05:06, Piotr Prusiński wrote: > I have just found out a curious thing that when changing the nstep value > to 1 (unity) in [soln-plugin-residual] closure, the code does not > produce residuals at all in the output file, a *.csv file is created but > stays empty.
This is due to the design of the residual plugin which estimates the residual through: residual = (prev_soln - curr_soln) / delta_t with the plugin itself having two phases: the first phase is run one time step before we are due to save out the residual and its purpose is to assign: self._prev = curr and then in the second phase (run on the next time step) we actually compute an approximation to the residual itself. The problem when nsteps = 1 is that in a single time step we have to perform both phases. This confuses the code as it is currently written. I can certainly update the plugin to support this. However, generally speaking nsteps = 1 is extremely frequent in terms of residual output. Usually, we would not run below nsteps = 50 given that with an explicit code such as PyFR the time steps themselves are usually extremely small. Regards, Freddie. -- You received this message because you are subscribed to the Google Groups "PyFR Mailing List" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send an email to [email protected]. Visit this group at https://groups.google.com/group/pyfrmailinglist. For more options, visit https://groups.google.com/d/optout.
signature.asc
Description: OpenPGP digital signature
