<div dir="ltr">Hello everyone,<div><br></div><div>I have a questions about vectors in an DMDA with DOF>1. For example, in 1d with number of grid N and DOF=2 (two fields u and v), the length of the global vector is 2*N. </div>
<div><br></div><div>What is the best way to restrict this vector (length 2*N) to a vector (length N) corresponding to the field u only?  This will help me obtain the properties of field u by using VecSum and other vec functions.</div>
<div><br></div><div>With DMDAVecGetArray and looping over only u field can do the job. However, I am just wondering whether any petsc function can provide either the restriction matrix or the vector restricted to a single field.</div>
<div><br></div><div>Thank you.</div><div><br></div><div>Best,</div><div>Xiangdong </div></div>