Landau ex1 should work. I will test.

On Wed, Jun 23, 2021 at 10:47 AM Matthew Knepley <knepley@gmail.com> wrote:
On Wed, Jun 23, 2021 at 10:44 AM Junchao Zhang <junchao.zhang@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 <mfadams@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