On 10/07/14 22:40, Zach Davis wrote:
> I was curious whether writing L1 or L2 norm information for each
> equation in the system being solved along with any user-defined
> force/moment integration information could be output and/or collected
> for purposes of plotting rates of convergence once the simulation has
> completed running.  I imagine this might incur a small performance hit,
> but is this something easily implemented and/or on the horizon?  I’m
> working on a short write-up about running PyFR on Rescale, and it would
> be nice to have that sort of information available.

Hi Zach,

This is currently not possible with the current version of PyFR however
it is on the roadmap.  The functionality will build on top of the
"completed step handler" mechanism that already exists within PyFR.  An
example of such a handler is the NaN checker:

    def nansweep(intg):
        if intg.nsteps % args.nansweep == 0:
            if any(np.isnan(np.sum(s)) for s in intg.soln):
                raise RuntimeError('NaNs detected at t = {}'
                                   .format(intg.tcurr))
    solver.completed_step_handlers.append(nansweep)

with the idea being that after each time step a set of user-defined
functions can be invoked which can interrogate the solution without
having to write it out to disk.

Regards, Freddie.

Attachment: signature.asc
Description: OpenPGP digital signature

Reply via email to