[petsc-users] VecScatter betweeen Vectors with different parallel layouts

Hello PETSc users I was wondering if you've got any suggestion for the issue mentioned in the topic. I'm solving a linear system on my unstructured grid that has its own partitioning. When I set up the linear system in PETSc I use the model partitioning whose grid points might not be evenly

Re: [petsc-users] specifying vertex coordinates using DMPlexCreateFromCellListParallel

I'm interested in this topic as well Thanks Ivano 2017-04-11 16:21 GMT+02:00 Hassan Raiesi : > Hello, > > > > I’m trying to use DMPlexCreateFromCellListParallel to create a DM from an > already partitioned mesh, > > It requires an array of

Re: [petsc-users] Symmetric matrix: Setting entries below diagonal

So, as far as I understand, the only benefit of PETSc with symmetric matrices is only when Matrix values are set, by reducing the overhead of MatSetValue calls? Thanks, Ivano 2017-04-07 17:19 GMT+02:00 Barry Smith : > > > On Apr 7, 2017, at 6:40 AM, Florian Lindner

[petsc-users] Problem with IS and VecScatter

Hello users/developers I'm trying to build a vecscatter object to migrate data from a vector x to a vector x2 having same global size but different parallel layout. Prior to this, I build an Index Set using the method ISCreateStride The IS is created correctly, since the program returns ierr=0

[petsc-users] info about MatPreallocateInitialize/Finalize

Dear all How does MatPreallocateInitialize/Finalize work? I don't see any reference for examples in the online manual, soI've browsed the source code to see whetherthe routine is used and how As far as I understood it works this way: call MatPreallocateInitialize(comm, rows, cols,