Landau ex1 should work. I will test. On Wed, Jun 23, 2021 at 10:47 AM Matthew Knepley <knep...@gmail.com> wrote:
> On Wed, Jun 23, 2021 at 10:44 AM Junchao Zhang <junchao.zh...@gmail.com> > wrote: > >> Use VecGetArrayRead/Write() to get up-to-date host pointers to the vector >> array. >> > > I think Mark is saying that those are not working. We do call > VecGetArrayRead() in the HDF5 code. > > Mark, it seem like a small broken code is necessary. > > Thanks, > > Matt > > >> --Junchao Zhang >> >> >> On Wed, Jun 23, 2021 at 9:15 AM Mark Adams <mfad...@lbl.gov> wrote: >> >>> First, there seem to be two pages for VecGetArrayAndMemType (one has a >>> pointer to the other). >>> >>> So I need to get a CPU array for HDF5 viewing. Totally broken for >>> devices. >>> >>> I don't find a VecGetArrayCpu[HOST] that does the right thing. >>> >>> Perhaps have VecGetArrayAndMemType return a valid CPU pointer when >>> "mtype==NULL"? >>> >>> Mark >>> >> > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > > https://www.cse.buffalo.edu/~knepley/ > <http://www.cse.buffalo.edu/~knepley/> >