[petsc-users] Fwd: Interpolation in staggered grid

2018-10-15 Thread Manuel Valera
Thanks Matthew, i have made some progress but i am still unsure on how to proceed to make the DMDAs work as intended, i will try to lay out what i am trying now: I was able to port the interpolation into the DMDA model but it works only in serial runs, and it becomes unstable for more than one pro

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Matthew Knepley
On Mon, Oct 15, 2018 at 12:03 PM Yingjie Wu wrote: > Thanks,Matt > > I'm sorry I didn't describe the problem clearly. > I'm still wondering how, although VecLoad supports parallel vector > reading, how can I guarantee that the local vector on the processor > corresponds to the DM's IS? > For exam

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Yingjie Wu
Thanks,Matt I'm sorry I didn't describe the problem clearly. I'm still wondering how, although VecLoad supports parallel vector reading, how can I guarantee that the local vector on the processor corresponds to the DM's IS? For example, when in FormInitialGuess, I want to make sure that the array

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Dave May
On Mon, 15 Oct 2018 at 16:54, Matthew Knepley wrote: > On Mon, Oct 15, 2018 at 10:42 AM Yingjie Wu wrote: > >> Dear Petsc developer: >> Hi, >> Thank you very much for your previous reply. >> I recently wanted to modify my program to parallel version, and >> encountered some problems in modifying

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Matthew Knepley
On Mon, Oct 15, 2018 at 11:13 AM Yingjie Wu wrote: > Thank you very much, Matt. > > My problem is actually a two-dimensional problem, and the matrix in my > binary file is actually a two-dimensional array, each point corresponding > to a value on the two-dimensional grid. > These are still vecto

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Yingjie Wu
Thank you very much, Matt. My problem is actually a two-dimensional problem, and the matrix in my binary file is actually a two-dimensional array, each point corresponding to a value on the two-dimensional grid. Your suggestion is very useful to me. I can try to convert the matrix binary file into

Re: [petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Matthew Knepley
On Mon, Oct 15, 2018 at 10:42 AM Yingjie Wu wrote: > Dear Petsc developer: > Hi, > Thank you very much for your previous reply. > I recently wanted to modify my program to parallel version, and > encountered some problems in modifying it. > 1. There are functions (read matrix) in the program that

[petsc-users] Some Problems in Modifying Parallel Programs

2018-10-15 Thread Yingjie Wu
Dear Petsc developer: Hi, Thank you very much for your previous reply. I recently wanted to modify my program to parallel version, and encountered some problems in modifying it. 1. There are functions (read matrix) in the program that reads files,will they affect my parallelism? The codes are as f