[petsc-users] How Can I scatter a global vector to a local vector without using DMGlobaltoLocalBegin/DMGlobaltoLocalEnd?

2019-01-15 Thread leejearl via petsc-users
Hi, Matt > > Hi, all Petscer: > > I am very sorry for disturbing you for my questions. > > > I like your questions. > > > > I have a > > problems as follow. > > 1. A DMPlex object dm is created. > > 2. A global vector and a local vector are created using the > following > > routines DMCr

Re: [petsc-users] PETSc turtorial example generating unexpected result

2019-01-15 Thread Maahi Talukder via petsc-users
Hi, Thank you for your reply. Yeah, as you suspected, I printed these numbers after that statement. Now it gives the expected values. Regards, Maahi Talukder On Mon, Jan 14, 2019 at 9:33 PM Mark Adams wrote: > This code puts the error in x: > > call VecAXPY(x,neg_one,u,ierr) > > I suspect you

[petsc-users] How Can I scatter a global vector to a local vector without using DMGlobaltoLocalBegin/DMGlobaltoLocalEnd?

2019-01-15 Thread leejearl via petsc-users
Hi Matt, Thanks for your helps. I have another question, but I think I should ask the question in a new subject. leejearl > > Hi, Matt > > > > > > Hi, all Petscer: > > > > I am very sorry for disturbing you for my questions. > > > > > > > > > I like your questions. > > > > > > > > > > I have

[petsc-users] Is there easy to update the ghost value in the vector created using DMCreatLocalVector?

2019-01-15 Thread leejearl via petsc-users
Hi all Petscer, I have a question about how to update the ghost value in the vector. The prolblem is as follow. 1. A dmplex dm object is created using DMPlexCreateFromFile 2. The dm is distributed using DMPlexDistribute, and the overlap is set as 1. 3. The ghosted dm is constructed by DMConstruct