Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-06 Thread David Knezevic
On Mon, Jun 6, 2016 at 4:16 PM, Barry Smith wrote: > > David, > >Actually all the routines for accessing the MUMPS information already > exist. I have added this to the manual page for MATSOLVERMUMPS > > Notes: When a MUMPS factorization fails inside a KSP solve,

Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-06 Thread Barry Smith
David, Actually all the routines for accessing the MUMPS information already exist. I have added this to the manual page for MATSOLVERMUMPS Notes: When a MUMPS factorization fails inside a KSP solve, for example with a KSP_DIVERGED_PCSETUP_FAILED, one can find the MUMPS information

Re: [petsc-users] How to create symmetric matrix?

2016-06-06 Thread Barry Smith
> On Jun 6, 2016, at 1:02 PM, Faraz Hussain wrote: > > I am learning petsc and using the following code to insert values into my > matrix: > > ierr = MatSetValue(A, row, col, value,INSERT_VALUES); > ierr = MatSetValue(A, col, row, value,INSERT_VALUES); > > Since my

Re: [petsc-users] Functionality in Petsc

2016-06-06 Thread Patrick Sanan
Assuming that you don't want to store this dense rank-1 matrix explicitly, one way to go is to define your own MATSHELL which holds your two vectors and can apply the outer product to a vector. > Am 06.06.2016 um 19:49 schrieb Sophie Léger : > > Hello, > > I’m

[petsc-users] Functionality in Petsc

2016-06-06 Thread Sophie Léger
Hello, I'm trying to multiply a column vector with a line vector (which would result in a matrix) and am wondering if this is possible with Petsc? I have looked, but have not found anything yet. Would it be possible to let me know if such a functionality exists and if so guide me in the right

Re: [petsc-users] Parallel assembly FEM matrix

2016-06-06 Thread Lawrence Mitchell
> On 6 Jun 2016, at 01:43, Nestor Cerpa wrote: > > Hello, > > I am learning to use Petsc in order to modify my sequential FEM code to work > in parallel. First, I am trying to understand how a FEM assembly loop would > work in parallel using Petsc but I run into some