[petsc-users] DM structures of PetscInt

[Matthew Knepley]

On Sun, Nov 8, 2015 at 7:37 AM, Wee Beng Tay <zonexo at gmail.com> wrote:
>
> Sent using CloudMagic Email
> <https://cloudmagic.com/k/d/mailapp?ct=pa&cv=7.4.10&pv=5.0.2&source=email_footer_2>
> 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. Based on my description, can I use Is? Is there an
> example I can follow?
>

Yes, use an IS. If you need unstructured communication, you will need to
use the PetscSF structure.

  Matt


> 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
>
>


-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151108/cc248f5c/attachment.html>