Roc Wang <[email protected]> writes: > The coordinates of nodes are needed to plot the 3-D distributions of > the solution. The matrix and the vector are managed with DMDA, so the > array for the coordinates can be obtained by > DMDAGetCoordinateDA(da,&cda) and DMDAVecGetArray(cda,gc,&coors). > Here, coors is defined as DMDACoor3d ***coors. Since the function > VecVeiw can only output one vector, Vec sol, to a binary file , the > array of coordinates must be output to another file and thus must be > read separately.
Not true, just read them in the same order you wrote them. > The problem is that the array of coordinates is local array and the > number of processes in the post-processing program must be same as > that in the solver. DMGetCoordinates() returns a global Vec. Store that, not the local Vec. >From the names, you must have an older version of PETSc. Please upgrade to petsc-3.4.
