> On Nov 8, 2015, at 7:37 AM, Wee Beng Tay <[email protected]> wrote: > > > > Sent using CloudMagic Email > On Sun, Nov 08, 2015 at 7:41 PM, Matthew Knepley <[email protected]> wrote: > > On Sat, Nov 7, 2015 at 11:27 PM, TAY wee-beng <[email protected]> wrote: > Hi, > > I need to use DM structures of type PETScInt. > > I tried to use: > > DM da_cu_types > > Vec cu_types_local,cu_types_global > > PetscInt,pointer :: cu_types_array(:,:,:) > > Is this allowed? The cu_types_array is to be of integer type. Because I got > into problems compiling. > > We do not support integer Vecs. You should use IS. What are you trying to do? > > Matt > > Hi, > > I've a cfd code and I need to separate cells inside or outside a body. I > denote the inside /boundary /outside cells as type =2/1/0. I've a subroutine > which determines the type in parallel and hence I need a Vec with integer > type.
The simplest thing is just to use a real value of 2. 1. or 0. Memoywise it is not a big deal since you will have dozens of vectors plus memory heavy matrix so you won't even notice this extra memory. > Based on my description, can I use Is? Is there an example I can follow? > > Thanks > > > -- > Thank you. > > Yours sincerely, > > TAY wee-beng > > > > > -- > What most experimenters take for granted before they begin their experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener
