I did not see VecRestoreArrayRead_Nest. On VecNest, VecRestoreArrayRead defaults to VecRestoreArray(x,a) , which copies a[] to x and is expensive. Is it an oversight?
--Junchao Zhang
Zhang, Junchao via petsc-dev Mon, 11 Feb 2019 11:43:37 -0800
I did not see VecRestoreArrayRead_Nest. On VecNest, VecRestoreArrayRead defaults to VecRestoreArray(x,a) , which copies a[] to x and is expensive. Is it an oversight?
--Junchao Zhang