[petsc-users] DM structures of PetscInt

[TAY wee-beng]

On 9/11/2015 2:08 AM, Barry Smith wrote:
>> 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.
>
Hi,

Ya, thought that this is simpler. It's working now. Thanks.
>> 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