.
On Wed, Feb 27, 2013 at 2:52 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Feb 26, 2013, at 4:44 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
These calls compute default values for normtype and pcside based on
partial KSP options,
which can cause an unwarranted error to be thrown
On Tue, Feb 26, 2013 at 8:59 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Feb 26, 2013, at 9:54 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Here's what I get with this simple example (source and makefile also
attached):
#include petscksp.h
#undef __FUNCT__
#define
These calls compute default values for normtype and pcside based on partial
KSP options,
which can cause an unwarranted error to be thrown.
For example, in code there may be
KSPSetPCSide(ksp, PC_RIGHT); /* hardwired default, e.g., in Moose. */
Now, when I'm using command-line arguments
-ksp_type
On Thu, Feb 21, 2013 at 12:54 PM, Richard Tran Mills rtm at eecs.utk.eduwrote:
Hi Jed,
The reason is that in the PFLOTRAN code we currently have some code like
this:
select case(discretization%itype)
case(STRUCTURED_GRID,STRUCTURED_GRID_MIMETIC)
call
Done:
https://bitbucket.org/petsc/petsc-dev/commits/aee59bf3b744565306435f5b76084937
Dmitry.
On Tue, Feb 12, 2013 at 10:06 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
I'll remove it.
Dmitry
On Feb 12, 2013 9:50 PM, Barry Smith bsmith at mcs.anl.gov wrote:
How hard will it be to get
I'll remove it.
Dmitry
On Feb 12, 2013 9:50 PM, Barry Smith bsmith at mcs.anl.gov wrote:
How hard will it be to get rid of the sys/shell directory? All that
code that needs to be constantly updated with Jed's crazy schemes is
driving me nuts.
Barry
-- next part
On Sat, Feb 9, 2013 at 12:39 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 8, 2013 at 7:48 PM, Satish Balay balay at mcs.anl.gov wrote:
- *never* sync with petsc-dev until this work is done.
- merge with petsc-dev when its done. [either merge - or rebase it
over]
The
On Tue, Jan 15, 2013 at 2:31 PM, Barry Smith bsmith at mcs.anl.gov wrote:
So who wrote the PetscOpFList() code and who will refactor it? Is this
a high priority or just something to be listed on that website and then
ignored for years?
I wrote it when I was planning to implement extra
On Tue, Jan 15, 2013 at 2:14 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Jan 15, 2013 at 1:50 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Do we really want mallocs for EVERY single multiple-dispatch function
call in PETSc? I suggest not having MatQueryOp() but calling directly
On Tue, Jan 15, 2013 at 2:51 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Jan 15, 2013 at 2:47 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Yes. Jed is too in love with Git to ever do this so you have to.
What does this have to do with choice of VCS? They all support commit
hooks.
On Tue, Jan 15, 2013 at 2:55 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Jan 15, 2013 at 2:51 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
I think passing PetscObject is natural here, in which case using
obj-type_id instead of obj-type_name can be done later. But on that
front
On Tue, Jan 15, 2013 at 2:59 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Jan 15, 2013 at 2:55 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
Some projects (e.g., moose) reject a noncompliant push. The committer
can then rewrite the patch locally before it is published and resubmit
On Tue, Jan 15, 2013 at 3:20 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Jan 15, 2013 at 3:14 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
There is, by the way, the 'type' field in the header, which the docs say
should no longer be used, and which is typically set to 0 or -1
I personally find git (and its branches) rather cumbersome and wish libmesh
used mercurial instead :-)
And if hgsubversion actually worked there would no need for git :-)
Dmitry
On Jan 9, 2013 7:17 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
Libmesh just moved to github as well.
I think if
On Jan 9, 2013 9:51 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jan 9, 2013, at 7:17 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
Libmesh just moved to github as well.
I think if you carefully consider the branching model, it has a clear
advantage over everything else. Dusty Phillips
On Jan 9, 2013 10:10 PM, Sean Farley sean.michael.farley at gmail.com
wrote:
On Wed, Jan 9, 2013 at 9:13 PM, Dmitry Karpeev karpeev at mcs.anl.gov
wrote:
I personally find git (and its branches) rather cumbersome and wish
libmesh
used mercurial instead :-)
Not to mention git's atrocious
branches.
On Jan 9, 2013 10:12 PM, Dmitry Karpeev karpeev at gmail.com wrote:
On Jan 9, 2013 10:10 PM, Sean Farley sean.michael.farley at gmail.com
wrote:
On Wed, Jan 9, 2013 at 9:13 PM, Dmitry Karpeev karpeev at mcs.anl.gov
wrote:
I personally find git (and its branches) rather
On Jan 9, 2013 10:13 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jan 9, 2013, at 10:03 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
On Jan 9, 2013 9:51 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Then I would consider switching over to git and proceeding to make
fun of hg
This appears to be similar to GASM with the whole index range in the IS on
a single rank: GASM will use MatGetSubMatices to gather the matrix, do the
solve and scatter back.
On Jan 7, 2013 8:55 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Mon, Jan 7, 2013 at 8:51 PM, Barry Smith bsmith at
until the beginning of
January, but I can help you with getting an account and getting started.
Btw: I just got a call regarding Altera hardware, we might have chances to
get our fingers on their OpenCL-enabled hardware.
Best regards,
Karli
--
Dmitry Karpeev, Ph.D.
Assistant
sys/shell is mine. It's a badly named attempt at simple language
interoperability. Since it's not really being used, as far as I know, it
can go, although the corresponding parts of petsc4py would have to be
removed, too. I can clean it out this week along with some other long dead
experimental
in the code about its purpose and current state.
Barry
On Dec 11, 2012, at 7:43 AM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
sys/shell is mine. It's a badly named attempt at simple language
interoperability. Since it's not really being used, as far as I know, it
can go, although
/ctetgen-dev
Any thoughts?
thanks,
Satish
--
Dmitry Karpeev, Ph.D.
Assistant Computational Mathematician
Mathematics and Computer Science
Argonne National Laboratory
Argonne, Illinois, USA
and
Fellow
Computation Institute
On Sat, Sep 29, 2012 at 12:16 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Sep 29, 2012, at 11:03 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Sep 28, 2012 at 9:24 PM, Barry Smith bsmith at mcs.anl.gov wrote:
I don't like having two typedefs
I tried
typedef const
On Thu, Sep 13, 2012 at 9:12 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
Constructor tab completion is better for DMCreateDA().
But will clash with all the other methods of the form DMDAGetXXX() etc.
I would vote for converting MatCreateSeqAIJ to MatSeqAIJCreate().
Note also that Mat*Create
the misbehaving thing.
Okay, I'll do that shortly. I've just tested that this will work with
recursive splits (at least the ones I obtain using libMesh/Moose).
Dmitry.
On Wed, Aug 22, 2012 at 6:40 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Wed, Aug 22, 2012 at 1:16 AM, Dmitry Karpeev karpeev
On Tue, Aug 21, 2012 at 8:46 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Aug 21, 2012 at 2:53 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Mon, Aug 20, 2012 at 8:36 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Mon, Aug 20, 2012 at 4:02 AM, Dmitry Karpeev karpeev
On Mon, Aug 20, 2012 at 8:36 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Mon, Aug 20, 2012 at 4:02 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Sun, Aug 19, 2012 at 1:17 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sun, Aug 19, 2012 at 12:48 PM, Dmitry Karpeev karpeev
AM, Jed Brown jedbrown at mcs.anl.govwrote:
On Mon, Aug 20, 2012 at 4:02 AM, Dmitry Karpeev karpeev at
mcs.anl.govwrote:
On Sun, Aug 19, 2012 at 1:17 PM, Jed Brown jedbrown at
mcs.anl.govwrote:
On Sun, Aug 19, 2012 at 12:48 PM, Dmitry Karpeev
karpeev at mcs.anl.gov wrote:
The main
On Sun, Aug 19, 2012 at 1:17 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sun, Aug 19, 2012 at 12:48 PM, Dmitry Karpeev karpeev at
mcs.anl.govwrote:
The main reason for this being a patch in 3.3 is that recursive
FieldSplit is broken there (which is why I couldn't enable it in
libMesh
I just pushed a couple of fixes to PCFIELDSPLIT in petsc-3.3 (
http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/943dde820f7f,
http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/5a9ebf885615) In
particular, the second changeset
modifies the way the splits' nullspaces are configured. In my
The main reason for this being a patch in 3.3 is that recursive FieldSplit
is broken there (which is why I couldn't enable it in libMesh/Moose
correctly). More comments below.
On Sun, Aug 19, 2012 at 9:32 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sun, Aug 19, 2012 at 2:31 AM, Dmitry
On Wed, Aug 1, 2012 at 5:01 PM, Todd Munson tmunson at mcs.anl.gov wrote:
On Aug 1, 2012, at 4:13 PM, Jed Brown wrote:
On Wed, Aug 1, 2012 at 11:10 AM, Todd Munson tmunson at mcs.anl.gov
wrote:
Yes, I am all in favor of actually specifying a VI externally and then
internally
);CHKERRQ(ierr);
That works for me. Matt, I thought that didn't work for SIMPLE?
Dmitry.
On Mon, Jul 9, 2012 at 1:56 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 9, 2012, at 1:17 AM, Dmitry Karpeev wrote:
On Sun, Jul 8, 2012 at 7:46 PM, Barry Smith bsmith at mcs.anl.gov wrote
I have a couple of questions about the usage DMGetSubDM() to define the
splits:
1) Should this be DMGetFieldSubDM(), since I can imagine other types of
subDMs (e.g., for a DD subdomain)?
In fact, I would advocate something like DMGetFieldSplitSubDM(), because
term field is getting abused and means
On Sun, Jul 8, 2012 at 7:46 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 8, 2012, at 3:36 PM, Dmitry Karpeev wrote:
As currently designed, there are always two KSPs for A00 - one is
created as part of MatSchurComplement and the other one sits in the FS
split.
Really? Isn't
On Sat, Jul 7, 2012 at 4:01 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 7, 2012, at 3:23 AM, Dmitry Karpeev wrote:
On Fri, Jul 6, 2012 at 9:56 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 6, 2012, at 9:35 PM, Matthew Knepley wrote:
On Fri, Jul 6, 2012 at 8:18 PM
I wanted to revive this (month-old) thread to try to finalize the API
change related to pivoting in PCFIELDSPLIT.
On Thu, Jun 7, 2012 at 6:42 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Jun 7, 2012 at 6:34 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
Well, the above defines
in
mind?
Dmitry.
On Sunday, July 8, 2012, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 8, 2012, at 2:03 AM, Dmitry Karpeev wrote:
On Sat, Jul 7, 2012 at 4:01 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jul 7, 2012, at 3:23 AM, Dmitry Karpeev wrote:
On Fri, Jul 6, 2012 at 9:56 PM
Barry
On Jul 6, 2012, at 8:39 AM, Matthew Knepley wrote:
On Fri, Jul 6, 2012 at 7:28 AM, Dmitry Karpeev karpeev at mcs.anl.gov
wrote:
On Fri, Jul 6, 2012 at 8:17 AM, Matthew Knepley knepley at gmail.com
wrote:
On Fri, Jul 6, 2012 at 5:06 AM, Dmitry Karpeev karpeev
I have the following problem with the prefix choice for the
MatSchurComplement KSP introduced in this changeset (
http://petsc.cs.iit.edu/petsc/petsc-dev/rev/0d4ccb990bb8).
I'm talking about the inner KSP for A00, effecting inv(A00) in the
definition S = A11 - A10 inv(A00) A01.
We also have the
Here's the line in question (also see the immediately preceding code):
http://petsc.cs.iit.edu/petsc/petsc-dev/rev/0d4ccb990bb8#l1.127
Dmitry.
-- Forwarded message --
From: Dmitry Karpeev karp...@mcs.anl.gov
Date: Fri, Jul 6, 2012 at 6:04 AM
Subject: Re: [petsc-dev] Problematic
On Fri, Jul 6, 2012 at 8:17 AM, Matthew Knepley knepley at gmail.com wrote:
On Fri, Jul 6, 2012 at 5:06 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
Here's the line in question (also see the immediately preceding code):
http://petsc.cs.iit.edu/petsc/petsc-dev/rev/0d4ccb990bb8#l1.127
This code should be disabled for now, in part, because DMKKT is incomplete.
The reason for the dependence on PC, though, is that DMKKT is built out of
pieces of GAMG, which is currently a PC, but ought to be a DM serving up
prolongations etc. to PCMG. I am planning to sort out the longer-term
I guess that will work: I can print a separate header before viewing the
Seq matrix.
Dmitry.
On Wed, Jun 6, 2012 at 7:30 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Wed, Jun 6, 2012 at 6:59 PM, Barry Smith bsmith at mcs.anl.gov wrote:
A better solution is to fix the parallel views for
If you built a nested fieldsplit and wanted to look at the ksp monitor
output from all the levels of nesting, before this fix
the output would be extremely confusing. This is because setting the tab
level on a PCFIELDSPLIT's inner KSPs doesn't
actually set the tab level on its PC. If one of
In the short run, however, we need to agree on an API for defining splits
that have an input and output
space. In the code the splits already carry row and col ISs, but there is
no way to set them separately,
at least not for IS-defined splits. Should we add something like
On Thu, Jun 7, 2012 at 6:21 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jun 7, 2012, at 6:17 PM, Dmitry Karpeev wrote:
In the short run, however, we need to agree on an API for defining
splits that have an input and output
space. In the code the splits already carry row and col ISs
When can we expect that (I mean 4.0)?
Moose now rejects commits with trailing white space.
Dmitry.
On Thu, Jun 7, 2012 at 6:24 PM, Barry Smith bsmith at mcs.anl.gov wrote:
For PETSc 4.0 we adopt an appropriate pre-commit source code pretty
print system; that prevents incorrectly formatted
On Thu, Jun 7, 2012 at 6:24 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Jun 7, 2012 at 6:17 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
DMCreateFieldDecomposition(DM dm, PetscInt *numsplits, char **splitnames,
IS *rowislist, IS *colislist, DM *dmlist)
How can this be right
On Wed, Jun 6, 2012 at 4:58 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Mar 23, 2012 at 10:15 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Mar 23, 2012, at 6:22 PM, Jed Brown wrote:
But what does using B as 00 block actually mean? It seems to me that
you pivot, producing P' P
When viewing a MATMPIAIJ on a communicator of size1, one can see something
like this (in an ASCII viewer):
Matrix Object: 1 MPI processes
type: mpiaij
row 0: (0, -8.32667e-17 + 0.4 i) (1, -1) (2, -1)
row 1: (0, -1) (1, -8.32667e-17 + 0.4 i) (3, -1) (8, -1)
row 2: (0, -1) (2, -8.32667e-17 +
On Mon, May 28, 2012 at 5:22 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Mon, May 28, 2012 at 5:10 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Mon, May 28, 2012 at 5:02 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
This is an enormous amount of new code and changes to put in during
, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I think I have a fix and want to reenable DM to serve up decompositions
to PCASM.
Was there a particular test that was failing? I don't see any
decomposition-related errors in alltest.
Dmitry.
On Sat, May 26, 2012 at 10:38 PM, Jed Brown jedbrown
. That fix will have to wait for the next release, I guess.
Dmitry.
On Mon, May 28, 2012 at 4:32 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
I refactored the code to get rid of field/domain decomposition mixing.
The use of DM in ASM is reenabled here:
http://petsc.cs.iit.edu/petsc/petsc-dev
On Saturday, May 26, 2012, Jed Brown jedbrown at mcs.anl.gov wrote:
Dmitry, I have problems with this commit.
http://petsc.cs.iit.edu/petsc/petsc-dev/rev/440ec9350d51
If I try running -pc_type asm with a vector problem defined on a DM,
PCASM calls DMCreateDecomposition which splits fields on
On Sat, May 26, 2012 at 9:07 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, May 26, 2012 at 8:04 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Saturday, May 26, 2012, Jed Brown jedbrown at mcs.anl.gov wrote:
Dmitry, I have problems with this commit.
http://petsc.cs.iit.edu
The problem is most likely in the creation routine or, rather, in the
arguments passed into it:
is it possible that the Vecs fed into VecScatterCreate have different block
sizes on different procs?
Dmitry.
On Fri, May 18, 2012 at 7:20 AM, Mark F. Adams mark.adams at
columbia.eduwrote:
I've
On Fri, May 11, 2012 at 11:40 AM, Barry Smith bsmith at mcs.anl.gov wrote:
Ok, in order to get the release out the door, please do not push
development work to petsc-dev. Only push fixes and removal of dead code
(DMMG for example). Also please run extensive tests and check the nightly
More accurately, the PetscViewerASCII subcomm viewer mechanism appears
broken.
I'm ready to push a fix, but since this is a heavily used and sensitive
part of the code
I wanted to make sure I'm not missing anything this close to a release.
The problem seems to start here:
Mark,
PetscTypeCompare got renamed to PetscObjectTypeCompare in line with the
PETSc naming guidelines:
http://petsc.cs.iit.edu/petsc/petsc-dev/rev/3fc8f665421d
The gamg source was part of this update.
Any chance you introduced a new instance of PetscTypeCompare() after that
changeset?
If so,
On Wed, Apr 25, 2012 at 4:09 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Apr 25, 2012, at 3:48 PM, Barry Smith wrote:
On Apr 25, 2012, at 2:44 PM, Mark F. Adams wrote:
On Apr 25, 2012, at 12:22 PM, Hong Zhang wrote:
Mark,
I attach it with this code:
/* attach
On Wed, Apr 25, 2012 at 4:41 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Apr 25, 2012, at 4:33 PM, Dmitry Karpeev wrote:
On Wed, Apr 25, 2012 at 4:09 PM, Mark F. Adams mark.adams at columbia.edu
wrote:
On Apr 25, 2012, at 3:48 PM, Barry Smith wrote:
or is it a more vague
On Wed, Apr 25, 2012 at 4:58 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Wed, Apr 25, 2012 at 16:41, Barry Smith bsmith at mcs.anl.gov wrote:
I admit MatSetVirtualColumnBlockSize() is clunky and would be happy to
hear alternatives.
What about having a separate row and column block
On Wed, Apr 25, 2012 at 5:58 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Apr 25, 2012, at 6:10 PM, Barry Smith wrote:
On Apr 25, 2012, at 4:58 PM, Jed Brown wrote:
On Wed, Apr 25, 2012 at 16:41, Barry Smith bsmith at mcs.anl.gov wrote:
I admit
Should this be renamed PetscObjectTypeCompare, since it operates on
PetscObjects?
Dmitry.
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120424/427ba991/attachment.html
Are you thinking of using PCSetModifyMatrices()? Actually, this brings up
another issue (now on petsc-dev).
Should this callback be attached to a Mat, instead of a PC? Then
MatGetSubMatrices() can use this callback on each of the extracted
submatrices, if necessary. In the case of a KKT system,
, 2012, at 9:35 PM, Dmitry Karpeev wrote:
Now it should work the way you intended.
Let me know if it doesn't, though.
Dmitry.
On Thu, Mar 29, 2012 at 7:49 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Mar 29, 2012, at 11:38 AM, Dmitry Karpeev wrote:
On Thu, Mar 29, 2012 at 10
Barry,
This looks like petsc-dev from about a week ago -- I fixed this bug today.
It runs on my laptop right now.
Dmitry.
On Fri, Mar 30, 2012 at 9:00 PM, Barry Smith bsmith at mcs.anl.gov wrote:
[0]PETSC ERROR:
I don't think so.
On Fri, Mar 30, 2012 at 9:39 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Do I have to rerun configure ?
On Mar 30, 2012, at 9:35 PM, Dmitry Karpeev wrote:
Maybe that arch didn't get rebuilt? The HG time stamp in your stack
trace says 'March 23'.
Hong alerted me
On Thu, Mar 29, 2012 at 9:04 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Mar 29, 2012 at 07:58, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I am not sure how to solve the problem of the relationship between
KSPReset() and PCReset() when the pc is being shared.
How useful
On Thu, Mar 29, 2012 at 10:13 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Mar 29, 2012 at 08:28, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Back to the original problem: how shall we prevent a shared pc from being
gutted by an unprotected PCReset() cascading from a KSPDestroy
(), but in this case building a whole new XXX from scratch
might be the right
thing to do: we are no longer reusing the configuration.
Dmitry.
On Thu, Mar 29, 2012 at 11:45 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
Pushed: http://petsc.cs.iit.edu/petsc/petsc-dev/rev/76f5543c7e89
On Thu, Mar 29
Okay, I see: refcnt is incremented in KSPGetPC().
On Thu, Mar 29, 2012 at 1:12 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
And refct is decremented from 2 to 1.
On Mar 29, 2012 12:04 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Are you sure this is correct?
KSPDestroy(eksp) will still
Now it should work the way you intended.
Let me know if it doesn't, though.
Dmitry.
On Thu, Mar 29, 2012 at 7:49 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Mar 29, 2012, at 11:38 AM, Dmitry Karpeev wrote:
On Thu, Mar 29, 2012 at 10:13 AM, Jed Brown jedbrown at mcs.anl.gov
Mark,
You should pass your subdomain ISs through the first IS list argument, and
leave the second as PETSC_NULL.
In that case GASM assumes that the user passed in a nonoverlapping
decomposition through the first list.
This is the behavior that GASM inherited from ASM and it isn't very well
(or at
On Mon, Mar 12, 2012 at 1:57 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Mar 12, 2012, at 1:50 PM, Sean Farley wrote:
What you prefer doesn't matter.
Clearly. I should have known Pedantic McGee over here would want to type
more and drop a simple name. Why not just put all the
On Tue, Feb 28, 2012 at 2:06 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Feb 28, 2012 at 13:52, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I don't understand the need for multiple function evaluation routines and
their explicit naming.
I think the situation here is completely
On Tue, Feb 28, 2012 at 2:21 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Feb 28, 2012 at 14:16, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Okay. I would still refrain from calling these things by the class of
the object consuming them. DMFormFunction() and DMFormIFunction() seem
On Mon, Feb 27, 2012 at 10:48 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
I'm getting lost with too many ways to do the same thing. The DM has only
one slot for a residual function, I'm not wild about having TS overwrite
that. So TS uses SNESSetFunction(), but now I need it on levels and I
On Sat, Feb 25, 2012 at 12:18 PM, Matthew Knepley knepley at gmail.com wrote:
On Sat, Feb 25, 2012 at 12:13 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
[migrating from petsc-maint]
On Sat, Feb 25, 2012 at 11:45, Barry Smith bsmith at mcs.anl.gov wrote:
I agree, but it would be quite
wrote:
The current API guarantees the FormFunction is called before FormJacobian.
Do you have a use case where that is not satisfied?
On Feb 17, 2012 11:13 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
The need for computation reuse (between J and J_pre, for example) might
require a change
make my interaction with libMesh a bit easier, but
isn't strictly necessary. Still, it would be good to know any thoughts on
this.
Dmitry.
On Sun, Feb 12, 2012 at 10:32 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Sat, Feb 11, 2012 at 8:36 AM, Jed Brown jedbrown at mcs.anl.gov wrote
On Sat, Feb 11, 2012 at 8:36 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Feb 11, 2012 at 08:32, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
This sounds fine to me, but what's the API for doing this? You want to do
KSPSetDM() and PCSetDM(),
so that the user has to explicitly pull
On Sat, Feb 11, 2012 at 8:24 AM, Matthew Knepley knepley at gmail.com wrote:
Barry is clearly right that each DM should produce a single Mat. However,
the proposed 2 DM (or
multiple DM) interface is just as clearly wrong. The 2 Mat interface
exists so that a user can specify
the system and
On Thu, Feb 9, 2012 at 11:56 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Feb 9, 2012 at 23:51, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Good point. It is still strange, however, that SNES insists on using the
same Mat obtained from DMGetMatrix for J and P.
Oh, you're
On Fri, Feb 10, 2012 at 12:01 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Thu, Feb 9, 2012 at 23:55, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
In a somewhat related matter, it appears that I cannot duplicate a
preallocated MATXXXAIJ until it has been assembled:
if my DM implementation
On Fri, Feb 10, 2012 at 12:11 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 00:02, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Yes. And, by the way, preallocation for J and P could be different (e.g.,
mass matrix vs lumped mass matrix for projection problems
On Fri, Feb 10, 2012 at 12:13 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 00:05, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I don't think so. Preallocating doesn't set any values and seems to
leave the matrix marked !assembled.
MatDuplicate for such a matrix
On Fri, Feb 10, 2012 at 12:13 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 00:05, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I don't think so. Preallocating doesn't set any values and seems to
leave the matrix marked !assembled.
MatDuplicate for such a matrix
On Fri, Feb 10, 2012 at 3:38 PM, Matthew Knepley knepley at gmail.com wrote:
On Fri, Feb 10, 2012 at 3:24 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Fri, Feb 10, 2012 at 12:13 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 00:05, Dmitry Karpeev karpeev
On Fri, Feb 10, 2012 at 3:40 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:24, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
So, as I understand, this behavior is reasonable for the use with SNES,
since the matrix is assembled with a call to SNES(DM)ComputeJacobian
On Fri, Feb 10, 2012 at 3:53 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:50, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Why not simply allow DM to return two matrices as SNES and KSP would want
anyway?
How do you specify the reduced stencil to use
On Fri, Feb 10, 2012 at 3:58 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 3:53 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:50, Dmitry Karpeev karpeev at mcs.anl.govwrote:
Why not simply allow DM to return two matrices as SNES and KSP
On Fri, Feb 10, 2012 at 3:59 PM, Matthew Knepley knepley at gmail.com wrote:
On Fri, Feb 10, 2012 at 3:58 PM, Dmitry Karpeev karpeev at mcs.anl.govwrote:
On Fri, Feb 10, 2012 at 3:53 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:50, Dmitry Karpeev karpeev
On Fri, Feb 10, 2012 at 4:00 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:58, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
It seems to me that the sparsity pattern (and even the type of the
matrices) is up to the particular DM type,
so it is, in principle, opaque
On Fri, Feb 10, 2012 at 4:02 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:59, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
As a separate matter, however, should preallocated, but not assembled
matrices be able to be duplicated or copied?
It doesn't even make
On Fri, Feb 10, 2012 at 4:05 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 15:59, Matthew Knepley knepley at gmail.com wrote:
This makes no sense to me unless the two matrices being returned have
different sparsity patterns, and we have no
good way for specifying that,
On Fri, Feb 10, 2012 at 4:28 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 16:13, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
Well, two is the number that's used in KSP and SNES.
DM can be used outside of KSP and SNES.
It may well make sense to have J and Jpre
On Fri, Feb 10, 2012 at 4:26 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Feb 10, 2012 at 16:06, Dmitry Karpeev karpeev at mcs.anl.gov wrote:
I'm not sure I follow. If a matrix has been created, its type has been
set and it has been preallocated (e.g., MATMPIAIJ in
ksp/ksp/examples
1 - 100 of 222 matches
Mail list logo