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