[petsc-users] Fwd: what is the equivalent DMDAVecRestoreArray() function in petsc4py?
HeeHo Park
heeho.park at gmail.com
Wed Feb 14 17:05:37 CST 2018
I just found a user group on PETSc website. Can someone please answer the
question below?
Thanks!
---------- Forwarded message ----------
From: HeeHo Park <heeho.park at gmail.com>
Date: Wed, Feb 14, 2018 at 5:04 PM
Subject: what is the equivalent DMDAVecRestoreArray() function in petsc4py?
To: dalcinl at gmail.com
Hi Lisandro,
I cannot find DMDAVecRestoreArray() equivalent in petsc4py.
I'm trying to set a 1D initial condition like this.
def initial_conditions(ts, U, appctx):
da = ts.getDM()
mstart,xm = da.getCorners()
mstart = mstart[0]
xm = xm[0]
M = da.getSizes()[0]
h = 1.0/M
mend = mstart + xm
u = da.getVecArray(U)
for i in range(mstart, mend):
u[i] = np.sin(np.pi*i*6.*h) + 3.*np.sin(np.pi*i*2.*h)
da.getVecRestoreArray(u)
Also, is there a better way to ask questions about petsc4py? a forum? or
google-group?
Thanks,
--
HeeHo Daniel Park
--
HeeHo Daniel Park
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180214/0fcd7762/attachment.html>
More information about the petsc-users
mailing list