saving parallel vectors

2006-11-24 Thread Satish Balay
On Fri, 24 Nov 2006, Matteo Semplice wrote: > Thansk for all the replies. > > I always thought that binary format would be best, but I couldn't find any > description of the binary format. Could you point me to a relevant piece of > the docs? This is documented in the man page for VecLoad() and

saving parallel vectors

2006-11-24 Thread Matteo Semplice
Thansk for all the replies. I always thought that binary format would be best, but I couldn't find any description of the binary format. Could you point me to a relevant piece of the docs? Can I load the PETSC binary files into MatLab or DataExplorer for further processing and visualization? M

saving parallel vectors

2006-11-23 Thread Satish Balay
To complete your query - you'll also need a Vec with the same communicator - i.e PETSC_COMM_SELF - for your approach to work. One way to do this is: VecGetArray(vec,&v) VecGetLocalSize(vec,size) VecCreateSeqWithArray(MPI_COMM_SELF,size,v,&seqvec) VecView(seqvec,seqviewer) VecDestroy(seqvec) Vec

saving parallel vectors

2006-11-23 Thread Matteo Semplice
This is really a newbie question but I am struggling to solve out of memory troubles on my first processor. I have a parallel global vector, say v, that I obtain with a call to DAGetGlobalVector. I save it to disk opening an ASCII standard viewer on PETSC_COMM_WORLD and calling VecView. As I un

saving parallel vectors

2006-11-23 Thread Barry Smith
NEVER use ascii for large data sets. Use the binary viewer to save them. See PetscViewerBinaryOpen(). Barry On Thu, 23 Nov 2006, Matteo Semplice wrote: > This is really a newbie question but I am struggling to solve out of memory > troubles on my first processor. > > I have a parallel g