function [da,err] = PetscDMDACreate2d(wrap,stentype,M,N,mm,nn,dof,s) da = PetscDM(); err = da.SetType('da'); da.SetPeriodicity(wrap); da.SetSizes([M N 1]); da.SetDof(dof); da.SetStencilType(stentype); da.SetStencilWidth(s); da.SetDim(2); da.SetFromOptions(); err = da.SetUp();