[petsc-users] DM structures of PetscInt

[Barry Smith]

> On Nov 8, 2015, at 7:37 AM, Wee Beng Tay <zonexo at gmail.com> wrote:
> 
> 
> 
> Sent using CloudMagic Email
> On Sun, Nov 08, 2015 at 7:41 PM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Sat, Nov 7, 2015 at 11:27 PM, TAY wee-beng <zonexo at gmail.com> 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