Dear petsc-users,
I have a 3x3 block system built as a single MatNest (with 9 Mats in it). I'd
like to treat this as a 2x2 system:
[ A B
C D ]
where A is 2x2 and precondition the A block with a schur complement. Since I
have a mat nest, the ISes for the three fields are just
On Mon, Dec 2, 2013 at 2:26 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
Dear petsc-users,
I have a 3x3 block system built as a single MatNest (with 9 Mats in it).
I'd like to treat this as a 2x2 system:
[ A B
C D ]
where A is 2x2 and precondition the A block with a
On 2 Dec 2013, at 12:23, Matthew Knepley knep...@gmail.com wrote:
On Mon, Dec 2, 2013 at 2:26 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
Dear petsc-users,
I have a 3x3 block system built as a single MatNest (with 9 Mats in it). I'd
like to treat this as a 2x2 system:
if (i == 0) {
for (f = 0; f numFields; ++f) {
ierr = PCFieldSplitSetIS(pc, fieldNames[f],
fields[f]);CHKERRQ(ierr);
ierr = PetscFree(fieldNames[f]);CHKERRQ(ierr);
ierr = ISDestroy(fields[f]);CHKERRQ(ierr);
}
numFields should be 2.
Mark Adams mfad...@lbl.gov writes:
if (i == 0) {
for (f = 0; f numFields; ++f) {
ierr = PCFieldSplitSetIS(pc, fieldNames[f],
fields[f]);CHKERRQ(ierr);
ierr = PetscFree(fieldNames[f]);CHKERRQ(ierr);
ierr = ISDestroy(fields[f]);CHKERRQ(ierr);
So this is a different KSP; not the one you want to use fieldsplit on?
Yes,
Please always use different prefixes for different solvers.
We use prefixes. I left them out of the example. So things like
'-s2_pc_type fieldsplit' cause an error. The 's2_' solver is a one field
(Laplacian)
Mark Adams mfad...@lbl.gov writes:
We use prefixes. I left them out of the example. So things like
'-s2_pc_type fieldsplit' cause an error. The 's2_' solver is a one field
(Laplacian) solve.
So when you say to use fieldsplit for a solve with only one field, you
get an error saying that
So when you say to use fieldsplit for a solve with only one field, you
get an error saying that fieldsplit requires at least two fields?
Yes
1. Don't do that.
My point was that this user got this message and so I suggesting a possible
problem: giving fieldsplit parameters to a single
On Mon, Dec 2, 2013 at 6:31 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
On 2 Dec 2013, at 12:23, Matthew Knepley knep...@gmail.com wrote:
On Mon, Dec 2, 2013 at 2:26 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
Dear petsc-users,
I have a 3x3 block
On 2 Dec 2013, at 15:20, Matthew Knepley knep...@gmail.com wrote:
On Mon, Dec 2, 2013 at 6:31 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
On 2 Dec 2013, at 12:23, Matthew Knepley knep...@gmail.com wrote:
On Mon, Dec 2, 2013 at 2:26 AM, Lawrence Mitchell
Mark Adams mfad...@lbl.gov writes:
My point was that this user got this message and so I suggesting a possible
problem: giving fieldsplit parameters to a single field matrix/solver.
Okay, I interpreted your original statement as being that this error
message has something to do with using
On Mon, Dec 2, 2013 at 9:47 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
On 2 Dec 2013, at 15:20, Matthew Knepley knep...@gmail.com wrote:
On Mon, Dec 2, 2013 at 6:31 AM, Lawrence Mitchell
lawrence.mitch...@imperial.ac.uk wrote:
On 2 Dec 2013, at 12:23, Matthew
Dear Developers of SLEPc,
I have found that SLEPc has now two eigenvalues solver that are based on
contour integration, rather than on Krylov space method.
They are EPSFEAST and EPSCISS.
I'm eager to try these solvers in my code both for Hermitian and
non-Hermitian problems.
My questions are
13 matches
Mail list logo