On Dec 1, 2010, at 8:00 AM, Dave May wrote: > Hi Barry, > I have updated the DMGetInterpolation_DA_XXX_Q1() family of > functions to enable them to > work with DMDA's which are nested, but have general coordinates > defined by the user. > The original used code which would compare two floating point numbers > with "==" and in general it > would seg fault with 1) non-regular coordinates and 2) even with > uniform coordinates with non constant dx,dy.dz. > > I have confirmed that the modified GetInterpolation() routines produce > second order interpolants, but before > pushing anything (using Jed as my proxy) I wanted to test the routines > with some galerkin multigrid examples / tests. > Grepping the source code I could not find any tests which use the flag > -pc_mg_galerkin. > > Can someone advise which tests actually test the functionality of the > interpolants constructed > from DMGetInterpolation_DA_XXX_Q1()?
Unlikely to be such tests, just push the stuff and let the chips fall where they May. Barry > > Cheers, > Dave
