not actually L-p norms, since it doesn't take into account the dx of the system. A user can take that into account or not.For p != 0, returns pth root of sum of pth powers over all points in all fabs and all components in the interval:
( sum [ |A[i][pt,var]|^p : FArrayBox A[i], point pt in A[i].box(), var in interval ] )^(1/p)
To turn into an L-p norm, one needs to multiply this by dx^(SpaceDim/p).
For p == 0, returns global max over all points in all fabs and all components in the interval:
max [ |A[i][pt,var]| : FArrayBox A[i], point pt in A[i].box(), var in interval ]
Some people don't like that this norm is not normalized based on number of points in A. Normalization is your problem.
This software is copyright (C) by the Lawrence Berkeley National Laboratory. Permission is granted to reproduce this software for non-commercial purposes provided that this notice is left intact.
It is acknowledged that the U.S. Government has rights to this software under Contract DE-AC03-765F00098 between the U.S. Department of Energy and the University of California.
This software is provided as a professional and academic contribution for joint exchange. Thus it is experimental, is provided ``as is'', with no warranties of any kind whatsoever, no support, no promise of updates, or printed documentation. By using this software, you acknowledge that the Lawrence Berkeley National Laboratory and Regents of the University of California shall have no liability with respect to the infringement of other copyrights by any part of this software.