[petsc-users] Provide Matrix Factorization to EPS for Generalized Eigenvalue Problem

2016-10-31 Thread Peetz, Darin T
Hello, I'm wondering how I could go about providing a matrix factorization calculated in Petsc to the eigenvalue routines in Slepc. I'm trying to solve the eigenvalue problem for stability, where the solution to KU=F is needed to construct the A-matrix (K_sigma) in the eigenvalue problem.

[petsc-users] --download-metis and build of stand-alone tools

2016-10-31 Thread Klaij, Christiaan
Satish, I've noticed that SuperLU depends on metis and parmetis and that PETSc downloads the versions 5.1.0-p3 and 4.0.3-p3. These are different from the Karypis latest stable versions (without the -p3). Do I really need these -p3 versions? If so, after configure, compilation and installation by

Re: [petsc-users] --download-metis and build of stand-alone tools

2016-10-31 Thread Jed Brown
"Klaij, Christiaan" writes: > Satish, > > I've noticed that SuperLU depends on metis and parmetis and that > PETSc downloads the versions 5.1.0-p3 and 4.0.3-p3. These are > different from the Karypis latest stable versions (without the > -p3). Do I really need these -p3

Re: [petsc-users] --download-metis and build of stand-alone tools

2016-10-31 Thread Satish Balay
On Mon, 31 Oct 2016, Klaij, Christiaan wrote: > Satish, > > I've noticed that SuperLU depends on metis and parmetis and that > PETSc downloads the versions 5.1.0-p3 and 4.0.3-p3. These are > different from the Karypis latest stable versions (without the > -p3). Do I really need these -p3

Re: [petsc-users] GAMG

2016-10-31 Thread Matthew Knepley
On Mon, Oct 31, 2016 at 10:29 AM, Kong, Fande wrote: > On Mon, Oct 31, 2016 at 8:44 AM, Jed Brown wrote: > >> Jeremy Theler writes: >> >> > Hi again >> > >> > I have been wokring on these issues. Long story short: it is about the >> >

Re: [petsc-users] GAMG

2016-10-31 Thread Kong, Fande
On Mon, Oct 31, 2016 at 8:44 AM, Jed Brown wrote: > Jeremy Theler writes: > > > Hi again > > > > I have been wokring on these issues. Long story short: it is about the > > ordering of the unknown fields in the vector. > > > > Long story: > > The physics

Re: [petsc-users] GAMG

2016-10-31 Thread Jed Brown
Jeremy Theler writes: > Hi again > > I have been wokring on these issues. Long story short: it is about the > ordering of the unknown fields in the vector. > > Long story: > The physics is linear elastic problem, you can see it does work with LU > over a simple cube (warp

Re: [petsc-users] GAMG

2016-10-31 Thread Jed Brown
"Kong, Fande" writes: > If the boundary values are not zero, no way to maintain symmetry unless we > reduce the extra part of the matrix. Not updating the columns is better in > this situation. The inhomogeneity of the boundary condition has nothing to do with operator

Re: [petsc-users] --download-metis and build of stand-alone tools

2016-10-31 Thread Klaij, Christiaan
Jed, Thanks, that line in the cmake file is exactly what I needed to know. A petsc configure option would be nice to have, but it's too difficult for me to do right now, I'll just hack the file instead. Chris dr. ir. Christiaan Klaij | CFD Researcher | Research & Development MARIN | T +31 317

Re: [petsc-users] GAMG

2016-10-31 Thread Jeremy Theler
On Mon, 2016-10-31 at 08:44 -0600, Jed Brown wrote: > > After understanding Matt's point about the near nullspace (and reading > > some interesting comments from Jed on scicomp stackexchange) I did built > > my own vectors (I had to take a look at MatNullSpaceCreateRigidBody() > > because I found