[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