[petsc-users] DM structures of PetscInt

Wee Beng Tay zonexo at gmail.com
Sun Nov 8 07:37:26 CST 2015


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 
[knepley at gmail.com] > wrote:
On Sat, Nov 7, 2015 at 11:27 PM, TAY wee-beng < zonexo at gmail.com 
[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?
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151108/1904bd7a/attachment.html>


More information about the petsc-users mailing list