[petsc-dev] easyinstall for petsc4py

2010-08-05 Thread Matthew Knepley
On Thu, Aug 5, 2010 at 9:56 AM, Lisandro Dalcin dalcinl at gmail.com wrote: On 5 August 2010 00:39, Lisandro Dalcin dalcinl at gmail.com wrote: On 4 August 2010 12:44, Barry Smith bsmith at mcs.anl.gov wrote: Lisandro, We would like to help you make it possible for petsc4py to be

[petsc-dev] PetscDLSym called on an empty handle

2010-08-05 Thread Matthew Knepley
On Thu, Aug 5, 2010 at 11:40 AM, Dmitry Karpeev karpeev at mcs.anl.gov wrote: What is the intended behavior of PetscDLSym when called on an empty handle: PetscDLSym(handle=PETSC_NULL, symbol, value) My understanding was that this should look for symbol in the symbol table of the main

[petsc-dev] PetscDLSym called on an empty handle

2010-08-05 Thread Matthew Knepley
On Thu, Aug 5, 2010 at 11:47 AM, Dmitry Karpeev karpeev at mcs.anl.gov wrote: On Thu, Aug 5, 2010 at 11:42 AM, Matthew Knepley knepley at gmail.comwrote: On Thu, Aug 5, 2010 at 11:40 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote: What is the intended behavior of PetscDLSym when called

[petsc-dev] easyinstall for petsc4py

2010-08-05 Thread Matthew Knepley
On Thu, Aug 5, 2010 at 3:48 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On 5 August 2010 11:59, Matthew Knepley knepley at gmail.com wrote: What I was hoping for is something along the lines of Atlas in numpy. However, looking for Atlas availability is way simpler that asking numpy

[petsc-dev] something to consider adding to PETSc

2010-08-12 Thread Matthew Knepley
Prob will not help that much over Metis for MatMult, but fill-reducing orderings might be a big win for LU. Matt On Thu, Aug 12, 2010 at 1:38 PM, Barry Smith bsmith at mcs.anl.gov wrote: From: Bas Fagginger Auer B.O.FaggingerAuer at uu.nl Date: Wed, 4 Aug 2010 15:02:44 +0200 Subject:

[petsc-dev] something to consider adding to PETSc

2010-08-12 Thread Matthew Knepley
(the original paper by Harald), but what is meant here? Matt is cache-oblivious, this would be a great feature to have it worked. Barry On Aug 12, 2010, at 3:07 PM, Matthew Knepley wrote: Prob will not help that much over Metis for MatMult, but fill-reducing orderings might be a big

[petsc-dev] Triangular solve on the GPU

2010-08-16 Thread Matthew Knepley
I am no longer as pessimistic as I once was. To start, this is a good report on prefix-sum in parallel: http://www.cs.cmu.edu/~blelloch/papers/Ble93.pdf It can solve first-order recurrences optimally. Higher order recurrences can be reduced to first-order, however it is not work optimal (too

[petsc-dev] Stop breaking my build

2010-08-19 Thread Matthew Knepley
libfast in: /PETSc3/petsc/petsc-dev/src/snes/impls/lsvi lsvi.c: In function 'PetscErrorCode SNESLSVIComputeMeritFunction(_p_Vec*, PetscScalar*)': lsvi.c:137: error: no matching function for call to 'VecNormBegin(_p_Vec*, NormType, PetscScalar*)' /PETSc3/petsc/petsc-dev/include/petscvec.h:468:

[petsc-dev] Stop breaking my build

2010-08-19 Thread Matthew Knepley
On Thu, Aug 19, 2010 at 5:37 PM, Satish Balay balay at mcs.anl.gov wrote: On Thu, 19 Aug 2010, Jed Brown wrote: On Thu, 19 Aug 2010 12:33:05 -0500 (CDT), Satish Balay balay at mcs.anl.gov wrote: Matt, Nobody knows how you build stuff. Atleast mention the configure command used

[petsc-dev] using consistant username for commits

2010-08-19 Thread Matthew Knepley
I am going to vote for consistent formatting the way Jed and I have it now. Is that alright? Matt On Thu, Aug 19, 2010 at 6:56 PM, Satish Balay balay at mcs.anl.gov wrote: Currently some of us let hg choose the 'username' automatically in our commits. Some of us use 'username=value' in

[petsc-dev] proceedings on preliminary work on PETSc with GPUs

2010-08-19 Thread Matthew Knepley
I spoke to Scott MacLachlan today (Victor's advisor at Tufts). He would like a little blurb saying how good Victor was this summer. He wants to give it to the Dean and get him to give Victor money for this semester instead of grading. I can write something to start, or feel free if you have time.

[petsc-dev] no module named cmakegen

2010-08-23 Thread Matthew Knepley
Can you explain the logPrint() problem? It is not quite clear. I tried to be careful that Framework was picklable. Matt On Mon, Aug 23, 2010 at 9:37 PM, Jed Brown jed at 59a2.org wrote: On Mon, 23 Aug 2010 16:18:27 -0500 (CDT), Satish Balay balay at mcs.anl.gov wrote: Well currently

[petsc-dev] no module named cmakegen

2010-08-23 Thread Matthew Knepley
On Mon, Aug 23, 2010 at 10:21 PM, Satish Balay balay at mcs.anl.gov wrote: On Mon, 23 Aug 2010, Satish Balay wrote: Now that there are more and more standalone scripts - perhaps we need a better consistant way to handle this - but I don't know what that is.. mercurial way is to have

[petsc-dev] no module named cmakegen

2010-08-24 Thread Matthew Knepley
On Mon, Aug 23, 2010 at 10:38 PM, Jed Brown jed at 59a2.org wrote: On Mon, 23 Aug 2010 17:21:52 -0500 (CDT), Satish Balay balay at mcs.anl.gov wrote: On Mon, 23 Aug 2010, Satish Balay wrote: Now that there are more and more standalone scripts - perhaps we need a better consistant way

[petsc-dev] What's the point of D(A/M)GetGlobalVector?

2010-08-27 Thread Matthew Knepley
Simply, in PETSc, getFoo() and restoreFoo() operate an object pool. Matt On Fri, Aug 27, 2010 at 11:23 AM, Jed Brown jed at 59a2.org wrote: On Fri, 27 Aug 2010 14:13:01 +0300, Aron Ahmadia aron.ahmadia at kaust.edu.sa wrote: What exactly is the purpose of these routines then? Is there

[petsc-dev] What's the point of D(A/M)GetGlobalVector?

2010-08-27 Thread Matthew Knepley
I agree with DACreateMatrix(). Matt On Fri, Aug 27, 2010 at 1:24 PM, Jed Brown jed at 59a2.org wrote: On Fri, 27 Aug 2010 16:13:14 +0300, Aron Ahmadia aron.ahmadia at kaust.edu.sa wrote: Not to mention the various Get routines that are actually used to create things, such as

[petsc-dev] What's the point of D(A/M)GetGlobalVector?

2010-08-28 Thread Matthew Knepley
I would support a name change to Create(). However, we should be really sure before we do stuff like that since nothing pisses people off (Wolfgang) like name changes. Matt On Fri, Aug 27, 2010 at 5:49 PM, Barry Smith bsmith at mcs.anl.gov wrote: Hmmm, petscda.h:EXTERN PetscErrorCode

[petsc-dev] [GPU] Performance on Fermi

2010-08-28 Thread Matthew Knepley
On Fri, Aug 27, 2010 at 7:19 PM, Keita Teranishi keita at cray.com wrote: Barry, Yes. It improves the performance dramatically, but the execution time for KSPSolve stays the same. MatMult 5.2 Gflops I will note that to put the matvec on the GPU you will also need -mat_type aijcuda.

[petsc-dev] [GPU] Performance of ex19

2010-08-31 Thread Matthew Knepley
On Tue, Aug 31, 2010 at 7:17 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Aug 31, 2010, at 3:14 PM, Keita Teranishi wrote: Does this PETSc use timers from CUDA? No, didn't know there were timers in CUDA. Yes, I use them when I really want to know how well I an utilizing the board,

Should BuildSystem download updated external packages when they are out of date?

2009-10-02 Thread Matthew Knepley
On Fri, Oct 2, 2009 at 9:11 AM, Richard Tran Mills rmills at climate.ornl.govwrote: BuildSystem Folks (Matt, mainly), This is not really a problem for me since I know about this behavior already, but I note that it can cause significant confusion when configure.py has been asked to download

Where are the gcov HTML pages?

2009-10-12 Thread Matthew Knepley
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL:

Subsecond sleep

2009-10-16 Thread Matthew Knepley
No. Matt On Fri, Oct 16, 2009 at 5:02 AM, Jed Brown jed at 59a2.org wrote: I sometimes want to slow down animation with a small problem size and -draw_pause is generally does what I want, but it can only take integer arguments and pausing for a whole second every time is often too much.

Fwd: People who need to solve large systems?

2009-10-21 Thread Matthew Knepley
This stuff is very cool, and based on very interesting mathematics. Anyone have matrices of this type? Thanks, Matt -- Forwarded message -- From: Gary Miller glmil...@cs.cmu.edu Date: Tue, Oct 20, 2009 at 2:46 PM Subject: People who need to solve large systems? To: Matthew

person who wrote PetscRandomGetValueReal() PetscRandomGetValueImaginary()

2009-10-23 Thread Matthew Knepley
Done. However, Igot an unholy merge on builder.py. It is definitely screwed up. I can fix it when I get home on Sunday. Matt On Fri, Oct 23, 2009 at 2:07 PM, Barry Smith bsmith at mcs.anl.gov wrote: PetscRandomGetValueReal() and PetscRandomGetValueImaginary() have no mention in their

person who wrote PetscRandomGetValueReal() PetscRandomGetValueImaginary()

2009-10-23 Thread Matthew Knepley
you ask me to put that stuff in? Matt Barry On Oct 23, 2009, at 2:26 PM, Matthew Knepley wrote: Done. However, Igot an unholy merge on builder.py. It is definitely screwed up. I can fix it when I get home on Sunday. Matt On Fri, Oct 23, 2009 at 2:07 PM, Barry Smith bsmith

PetscRandomGetValue(), PetscRandomGetValueReal(), PetscRandomGetValueImaginary()

2009-10-26 Thread Matthew Knepley
I am fine with that as long as the doc is clear. Matt On Mon, Oct 26, 2009 at 3:19 PM, Barry Smith bsmith at mcs.anl.gov wrote: There is a lack of symmetry in the meaning of PetscRandomGetValueReal() and PetscRandomGetValueImaginary(). The real returns a PetscReal with a random value in

changes for PetscMap

2009-10-27 Thread Matthew Knepley
I wanted that 10 years ago. Matt On Tue, Oct 27, 2009 at 4:37 PM, Barry Smith bsmith at mcs.anl.gov wrote: PetscMap is a strange beast in PETSc. 1) it is not a PetscObject (that is it has no PETSc header) 2) it is used as a PetscMap* instead of PetscMap like most objects 3) it does

Move src/sys/viewerimpl.h to include/private ?

2009-10-28 Thread Matthew Knepley
Nope, get rid of it. Matt On Wed, Oct 28, 2009 at 12:32 PM, Jed Brown jed at 59a2.org wrote: I don't think external projects should have to include ../src/whatever because this doesn't work if PETSc is installed and it's a screwy API (even if most users don't see it). Any objections to

ML symbols in petscksp.so

2009-10-28 Thread Matthew Knepley
You always have to link the ML library if you configure for ML, right? Matt On Wed, Oct 28, 2009 at 3:34 PM, Jed Brown jed at 59a2.org wrote: Something recently caused these symbols to get lost, meaning that users have to explicitly link the static libml.a. Jed -- What most

ML symbols in petscksp.so

2009-10-28 Thread Matthew Knepley
On Wed, Oct 28, 2009 at 4:14 PM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: You always have to link the ML library if you configure for ML, right? No, all the third party libs used to be put into the PETSc shared libs (but not the static libs). The easiest way to do

most efficient way to use local CSR matrix?

2009-10-29 Thread Matthew Knepley
On Thu, Oct 29, 2009 at 1:46 PM, Chris Kees christopher.e.kees at usace.army.mil wrote: Mark, Barry, Thanks for the help. For now I'm sticking with the approach of copying the csr matrix and using the csr data structures to do the preallocations. I'll eventually get around to writing the

missing CHKERRQ in MatNullSpaceRemove()

2009-10-29 Thread Matthew Knepley
Pushed. Matt On Thu, Oct 29, 2009 at 1:59 PM, Lisandro Dalcin dalcinl at gmail.com wrote: While trying to update petsc4py to latest changes in petsc-dev repo, I've noticed that at the very end of MatNullSpaceRemove() we have: if (sp-remove){ ierr = (*sp-remove)(sp,vec,sp-rmctx); }

ML symbols in petscksp.so

2009-10-29 Thread Matthew Knepley
On Thu, Oct 29, 2009 at 4:01 PM, Satish Balay balay at mcs.anl.gov wrote: On Thu, 29 Oct 2009, Barry Smith wrote: On Oct 29, 2009, at 3:52 PM, Lisandro Dalcin wrote: On Thu, Oct 29, 2009 at 4:45 PM, Jed Brown jed at 59a2.org wrote: I realize that the real problem here was

PETSc in double-double and quad-double...

2009-11-04 Thread Matthew Knepley
long? Matt Barry On Nov 4, 2009, at 1:33 PM, Matthew Knepley wrote: If you are alright with EVERYTHING running with that, it already works. Just set PetscScalar to that type. I would help you setup a configure switch for it if needed. If you only want it for some operations, we

PETSc in double-double and quad-double...

2009-11-04 Thread Matthew Knepley
On Wed, Nov 4, 2009 at 2:21 PM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: On Wed, Nov 4, 2009 at 2:06 PM, Barry Smith bsmith at mcs.anl.gov mailto:bsmith at mcs.anl.gov wrote: This is tricky because some C/C++ compilers long double uses 128 bits of space

who broke PCFieldSplit code in petsc-dev?????????

2009-11-11 Thread Matthew Knepley
Has this been broken for a while? It bombs in MatGetSubMatrix(), and this had an API change. The last change in fieldsplit.c is in July, by Jed. Matt On Wed, Nov 11, 2009 at 5:16 PM, Barry Smith bsmith at mcs.anl.gov wrote: See for example

who broke PCFieldSplit code in petsc-dev?????????

2009-11-11 Thread Matthew Knepley
I have 5 args for fieldsplit.c:290. Matt On Wed, Nov 11, 2009 at 5:27 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Wed, Nov 11, 2009 at 8:23 PM, Matthew Knepley knepley at gmail.com wrote: Has this been broken for a while? It bombs in MatGetSubMatrix(), and this had an API change

who broke PCFieldSplit code in petsc-dev?????????

2009-11-11 Thread Matthew Knepley
Fixed. Matt On Wed, Nov 11, 2009 at 5:46 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Wed, Nov 11, 2009 at 8:41 PM, Matthew Knepley knepley at gmail.com wrote: I have 5 args for fieldsplit.c:290. But the traceback Barry posted says arg #6 is wrong... In short, MatGetSubMatrix

who broke PCFieldSplit code in petsc-dev?????????

2009-11-11 Thread Matthew Knepley
could easily check what DAY things broke. Satish please have it save the old days stuff for each day. Barry On Nov 11, 2009, at 5:56 PM, Matthew Knepley wrote: Fixed. Matt On Wed, Nov 11, 2009 at 5:46 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Wed, Nov 11, 2009 at 8:41 PM

Fwd: [computational.science] CFP: Workshop on Large-Scale Parallel Processing (LSPP'10)

2009-11-13 Thread Matthew Knepley
-- Forwarded message -- From: Robert A. van de Geijn r...@cs.utexas.edu Date: Fri, Nov 13, 2009 at 3:21 PM Subject: Fwd: [computational.science] CFP: Workshop on Large-Scale Parallel Processing (LSPP'10) To: flame_cs at cs.utexas.edu Begin forwarded message: *From: *Darren

Why did you removed PETSC_ARCH_NAME variable in makefiles?

2009-11-27 Thread Matthew Knepley
Why not just hash the RDict.db? Matt On Fri, Nov 27, 2009 at 1:24 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Fri, Nov 27, 2009 at 2:21 AM, Satish Balay balay at mcs.anl.gov wrote: The way I look at it - the primary complexity here is from PETSc buildsystem. Because PETSc

PETSC_ARCH variable -- PETSC_CONF

2009-11-30 Thread Matthew Knepley
I see the point, but this is one of the oldest parts of PETSc, and I am hesitant to change one arbitrary name to another without a more convincing reason. Matt On Mon, Nov 30, 2009 at 2:17 PM, Barry Smith bsmith at mcs.anl.gov wrote: Lisandro has pointed out to me several times that the

PETSC_ARCH variable -- PETSC_CONF

2009-11-30 Thread Matthew Knepley
saving people's time and energy. This will be far outweighed by the number of people complaining about such a change. I agree with them. Matt Barry On Nov 30, 2009, at 3:25 PM, Matthew Knepley wrote: I see the point, but this is one of the oldest parts of PETSc, and I am hesitant

TS: high-order with Galerkin spatial discretization, semi-explicit index-1 DAE

2009-03-18 Thread Matthew Knepley
I could not understand the post. We will have to go more slowly. To begin, I cannot understand why you cannot use TSSetMatrices() http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/TS/TSSetMatrices.html#TSSetMatrices Thanks, Matt On Tue, Mar 17, 2009 at 9:15 PM,

TS: high-order with Galerkin spatial discretization, semi-explicit index-1 DAE

2009-03-18 Thread Matthew Knepley
On Wed, Mar 18, 2009 at 5:27 PM, Jed Brown jed at 59a2.org wrote: On Wed 2009-03-18 13:59, Matthew Knepley wrote: I cannot understand why you cannot use TSSetMatrices() TSSetMatrices is only for linear problems and doesn't generalize nicely to nonlinear problems, especially with matrix-free

TS: high-order with Galerkin spatial discretization, semi-explicit index-1 DAE

2009-03-19 Thread Matthew Knepley
with Tim Gowers' polymath experiment). Satish, do we have a Wiki, or could we setup one in short order? Thanks, Matt On Thu, Mar 19, 2009 at 5:28 AM, Jed Brown jed at 59a2.org wrote: On Wed 2009-03-18 18:34, Matthew Knepley wrote: Not a whole lot. Needs to go slower, since I think

building PETSc on Windows

2009-03-23 Thread Matthew Knepley
On Mon, Mar 23, 2009 at 8:23 AM, Lisandro Dalcin dalcinl at gmail.com wrote: The Cygwin way works, but let me say that it is a bit of pain... as a novice in this crapy environment, I've spent a lot of time just to get all the Cygwin stuff properly installed... I'm exploring the possibility

building PETSc on Windows

2009-03-23 Thread Matthew Knepley
On Mon, Mar 23, 2009 at 11:00 AM, Lisandro Dalcin dalcinl at gmail.com wrote: On Mon, Mar 23, 2009 at 10:57 AM, Satish Balay balay at mcs.anl.gov wrote: The main issue is - win32fe works with cygwin PATHs - not windows paths.. Sorry, I'm do not know what are you talking about ... Is

building PETSc on Windows

2009-03-23 Thread Matthew Knepley
On Mon, Mar 23, 2009 at 11:09 AM, Lisandro Dalcin dalcinl at gmail.com wrote: On Mon, Mar 23, 2009 at 11:48 AM, Matthew Knepley knepley at gmail.com wrote: 2) BuildSystem assumes that the shell is '/bin/sh'... I would need it to find sh in the $PATH instead... Where

MatMult_Scatter woes

2009-03-30 Thread Matthew Knepley
On Mon, Mar 30, 2009 at 5:52 AM, Jed Brown jed at 59a2.org wrote: The current implementation of MatMult_Scatter and MatMultTranspose_Scatter do not actually have matrix semantics unless the index sets are permutations. For example, it's generally desirable that MatMult(A,x,y); agrees

MatMult_Scatter woes

2009-03-30 Thread Matthew Knepley
, they should be added. You are right, that square matrices can still have this property. However, I do not see how the above sequence fixes this. I can't see how you would get around an addition. Matt Barry On Mar 30, 2009, at 7:04 PM, Matthew Knepley wrote: On Mon, Mar 30, 2009 at 5:52 AM

GPU related stuff

2009-07-09 Thread Matthew Knepley
On Thu, Jul 9, 2009 at 6:15 AM, Farshid Mossaiby mossaiby at yahoo.com wrote: Hi all, Some time ago on this list, there was some discussion about GPU and a GPU version of PETSc. I would like to know if there has been any progress. Also, I need some advice on preconditioners suitable for GPU

GPU related stuff

2009-07-09 Thread Matthew Knepley
On Thu, Jul 9, 2009 at 7:31 AM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: PCs which have high flop to memory access ratios look good. No surprise there. My concern here is that almost all good preconditioners are multiplicative in the fine-grained kernels or do significant

Python 2.6 and deprecated modules

2009-07-13 Thread Matthew Knepley
On Mon, Jul 13, 2009 at 4:56 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Mon, Jul 13, 2009 at 6:47 PM, Matt Knepleyknepley at gmail.com wrote: I like 1) as well. However someone needs to test the fixes on ancient RH python. No problem, I have custom installs of old Python's

Chat for coding collaboration

2009-07-20 Thread Matthew Knepley
I have found chat to be the best way to resolve simple problems when coding, especially questions about parts I did not write. It would be nice if we all had access to the same chat program to facilitate this kind of interaction in our small group. I tend to prefer GChat. Does anyone else have a

Chat for coding collaboration

2009-07-20 Thread Matthew Knepley
I use my mail reader for GChat. I think Apple can do that. Matt On Mon, Jul 20, 2009 at 2:53 PM, Barry Smith bsmith at mcs.anl.gov wrote: I want to use may mail reader as my front end Barry It happens to be Mail.app but that is not important On Jul 20, 2009, at 2:41 PM, Andy Ray

Chat for coding collaboration

2009-07-20 Thread Matthew Knepley
The steps are simple: 1) Invite partychat2 at gmail.com to chat using the contacts on the left 2) Open a chat with partychat2 3) In the window type '/join petsc' 4) Then we can all talk at once, and it will be saved I think Matt On Mon, Jul 20, 2009 at 4:02 PM, Andy Ray Terrel

Setting options on matrices obtained from a DM

2009-07-21 Thread Matthew Knepley
Actually, this seems like the same problem that Lisandro is having, just with different functions. I propose making data structures do the work for us rather than complicated organization in an imperative program. We could use the same mechanism we use in configure to handle issues of object

Setting options on matrices obtained from a DM

2009-07-22 Thread Matthew Knepley
No comments? Not even This is complete shit!? Matt On Tue, Jul 21, 2009 at 8:56 PM, Matthew Knepley knepley at gmail.com wrote: Actually, this seems like the same problem that Lisandro is having, just with different functions. I propose making data structures do the work for us rather

Setting options on matrices obtained from a DM

2009-07-23 Thread Matthew Knepley
On Thu, Jul 23, 2009 at 8:15 AM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: On Wed, Jul 22, 2009 at 3:57 PM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: No comments? Not even This is complete shit!? Heh, it might be overkill and, if I understand you correctly

Setting options on matrices obtained from a DM

2009-07-23 Thread Matthew Knepley
On Thu, Jul 23, 2009 at 9:45 AM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: MatExecute()? What about MatGetXXX() before MatExecute()? It will reregister that function, so you would need another SetUp() call, exactly the same interface we have now. But the user wants

Setting options on matrices obtained from a DM

2009-07-23 Thread Matthew Knepley
On Thu, Jul 23, 2009 at 10:21 AM, Jed Brown jed at 59a2.org wrote: Matthew Knepley wrote: It would be impossible for me to disagree more strongly. A user cannot get information that is not setup in the current model. I really do not understand your objection, and the Haskell discussion

Prefix for DACreate

2009-02-21 Thread Matthew Knepley
I will do it, but it will have to wait for weekend after next since I have to go the SIAM CSE. Matt On Sat, Feb 21, 2009 at 2:57 PM, Jed Brown jed at 59a2.org wrote: We solve equations in a high-aspect ratio domain and have both a 2- and 3-dimensional DA which need to be compatible. Due to

USFFT Mat class

2008-10-08 Thread Matthew Knepley
This means that you cannot use -mat_type fftw. Probably because no combination of options can currently initialize it correctly. It is fine to leave out USFFT until this is possible. Matt On Wed, Oct 8, 2008 at 2:20 PM, Dmitry Karpeev karpeev at mcs.anl.gov wrote: I'm writing a new Seq

request for comments/advices about mapping PetscCookie - Python Type

2008-10-10 Thread Matthew Knepley
That approach sounds fine. I think that initializing the COOKIEs was not liked by some compilers. Ugh. Matt On Fri, Oct 10, 2008 at 1:31 PM, Lisandro Dalcin dalcinl at gmail.com wrote: In petsc4py, I've just implemented support for compose()/query() in my base Object class. However, I

FW: [PETSC #18391] PETSc crash with memory allocation in ILU preconditioning

2008-10-10 Thread Matthew Knepley
to it. I think I'll send a test case that reproduces the problem. -Original Message- From: knepley at gmail.com [mailto:knepley at gmail.com] On Behalf Of Matthew Knepley Sent: Tuesday, October 07, 2008 2:21 PM To: Rhew, Jung-hoon Cc: PETSC Maintenance Subject: Re: [PETSC #18391] PETSc

Independently configuring down up smoothers in PCMG

2008-10-20 Thread Matthew Knepley
On Mon, Oct 20, 2008 at 4:18 PM, Barry Smith bsmith at mcs.anl.gov wrote: Sorry I haven't responded to this fully. I hate the idea of a hack that is just for this case and a special flag you got to set or pass in. For the _1_, _2_ etc I did the following. If the string in the code has

Independently configuring down up smoothers in PCMG

2008-10-20 Thread Matthew Knepley
library to do just this, and I think it is much more understandable than ad hoc rules from us. Matt Barry On Oct 20, 2008, at 4:56 PM, Matthew Knepley wrote: On Mon, Oct 20, 2008 at 4:18 PM, Barry Smith bsmith at mcs.anl.gov wrote: Sorry I haven't responded to this fully. I hate the idea

error in handling download external packages

2008-10-22 Thread Matthew Knepley
I have a problem with hiding the package includes. If I tell PETSc to install ParMetis, I think I should be able to use ParMetis however I want. I do not just want to use it through the PETSc interface. I have done this numerous times. I think most users think this way as well. What is the

Check it out

2008-03-09 Thread Matthew Knepley
http://en.wikipedia.org/wiki/Portable%2C_Extensible_Toolkit_for_Scientific_Computation Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

changes for next PETSc release

2008-03-17 Thread Matthew Knepley
On Sun, Mar 16, 2008 at 7:38 PM, Barry Smith bsmith at mcs.anl.gov wrote: There are two significant changes I'd like to see before the next PETSc release: 1) remove the overly complicated (from a user perspective) matrix subclassing for the various external matrix solver

changes for next PETSc release

2008-03-17 Thread Matthew Knepley
On Mon, Mar 17, 2008 at 12:27 PM, Barry Smith bsmith at mcs.anl.gov wrote: 2) properly name-space PETSc by putting a Petsc in front of all PETSc objects, function names etc (this will require changing a few names also to keep them below the 32 character limit). This will

changes for next PETSc release

2008-03-17 Thread Matthew Knepley
together. I jsut really hate the idea of putting PETSc in front of every word in the package. It is really the ugliest thing I can imagine and will make programming that much more of a slog. Matt On Mar 17, 2008, at 12:59 PM, Matthew Knepley wrote: On Mon, Mar 17, 2008 at 12:27 PM, Barry

petscvariables: hardwired build dir instead of install dir

2008-03-24 Thread Matthew Knepley
If someone explains to me in words what that sed crap is doing, I will rewrite it in python. Thanks, Matt On Mon, Mar 24, 2008 at 10:57 AM, Lisandro Dalcin dalcinl at gmail.com wrote: Barry, things are still broken. I think that at some point we have to review the 'install:' target

petscvariables: hardwired build dir instead of install dir

2008-03-24 Thread Matthew Knepley
On Mon, Mar 24, 2008 at 9:50 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Mar 24, 2008, at 10:57 AM, Lisandro Dalcin wrote: Barry, things are still broken. I think that at some point we have to review the 'install:' target more carefully. First, the 'sed' command i being called

petscvariables: hardwired build dir instead of install dir

2008-03-24 Thread Matthew Knepley
differently on different systems. Barry On Mar 24, 2008, at 10:07 PM, Matthew Knepley wrote: On Mon, Mar 24, 2008 at 9:50 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Mar 24, 2008, at 10:57 AM, Lisandro Dalcin wrote: Barry, things are still broken. I think that at some point

petscvariables: hardwired build dir instead of install dir

2008-03-24 Thread Matthew Knepley
On Mon, Mar 24, 2008 at 10:38 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Mar 24, 2008, at 10:25 PM, Matthew Knepley wrote: On Mon, Mar 24, 2008 at 10:14 PM, Barry Smith bsmith at mcs.anl.gov wrote: Matt, The sed is so trivial it is silly to even think about

petscvariables: hardwired build dir instead of install dir

2008-03-25 Thread Matthew Knepley
On Mon, Mar 24, 2008 at 10:53 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Mar 24, 2008, at 10:42 PM, Matthew Knepley wrote: On Mon, Mar 24, 2008 at 10:38 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Mar 24, 2008, at 10:25 PM, Matthew Knepley wrote: On Mon, Mar 24, 2008

This build is so broken

2008-03-31 Thread Matthew Knepley
I have spent the last hour trying to get my PETSc to work. This sucks. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

This build is so broken

2008-03-31 Thread Matthew Knepley
I pushed fixes for all this crap. Matt On Mon, Mar 31, 2008 at 9:46 AM, Matthew Knepley knepley at gmail.com wrote: I have spent the last hour trying to get my PETSc to work. This sucks. Matt -- What most experimenters take for granted before they begin their experiments

This build is so broken

2008-03-31 Thread Matthew Knepley
On Mon, Mar 31, 2008 at 10:12 AM, Lisandro Dalcin dalcinl at gmail.com wrote: Which part? The 'configure' one, or the 'make' one? Anyone who builds needs to pull BuildSystem and reconfigure, and then a bunch of the make was broken. It is fixed now. Matt On 3/31/08, Matthew Knepley knepley

config/configure.py - configure?

2008-07-17 Thread Matthew Knepley
On Thu, Jul 17, 2008 at 4:14 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Jul 17, 2008, at 3:58 PM, Boyana Norris wrote: Speaking of configure and user experiences, I just noticed that PETSc defines top_builddir (in petscconf), which conflicts with automake-generated makefiles, making it

Repository

2008-07-18 Thread Matthew Knepley
I cannot ping petsc.cs.iit.edu. I have a tutorial scheduled for 8:30am and was planning for everyone to pull and build. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. --

problem with PetscSynchronizedFPrintf (fwd)

2008-07-21 Thread Matthew Knepley
I made this change because synchronized printing was cutting off my stuff. Now the queue is dynamically allocated so it can accomodate any size string. I tried ex7 and it runs fine for me, even under valgrind. Shi, can you reproduce this with a Sieve example? Matt On Mon, Jul 21, 2008 at

problem with PetscSynchronizedFPrintf (fwd)

2008-07-21 Thread Matthew Knepley
not saying to change back the code, just to be careful and check it. Barry On Jul 21, 2008, at 2:47 PM, Matthew Knepley wrote: I made this change because synchronized printing was cutting off my stuff. Now the queue is dynamically allocated so it can accomodate any size string. I tried

problem with PetscSynchronizedFPrintf (fwd)

2008-07-21 Thread Matthew Knepley
-icc-ifort/include -D__SDIR__='src/sys/fileio/' mprint.c mprint.c(350): warning #592: variable fullLength is used before its value is set while(fullLength = next-size) { ^ On Jul 21, 2008, at 3:30 PM, Matthew Knepley wrote: Pushed a fix. I missed the Viewer when I

problem with PetscSynchronizedFPrintf (fwd)

2008-07-21 Thread Matthew Knepley
at all the code. There are no calls in PetscInitialize() and I cannot find anything else that might be called before PetscMalloc is set. Matt Barry On Jul 21, 2008, at 2:47 PM, Matthew Knepley wrote: I made this change because synchronized printing was cutting off my stuff. Now the queue

Check it out

2008-07-29 Thread Matthew Knepley
http://www.selenic.com/mercurial/wiki/index.cgi/ForestExtension Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

docs for MINRES, only for positive definite?

2008-02-28 Thread Matthew Knepley
Docs are wrong. Matt 2008/2/28 Lisandro Dalcin dalcinl at gmail.com: I've noticed that the docs for MINRES say that the operator and the preconditioner must be POSITIVE DEFINITE. But I understand MINRES is tailored for the symmetric/hermitian-indefinite case. Are the docs wrong? Or the

docs for MINRES, only for positive definite?

2008-02-28 Thread Matthew Knepley
, Matthew Knepley wrote: Docs are wrong. Matt 2008/2/28 Lisandro Dalcin dalcinl at gmail.com: I've noticed that the docs for MINRES say that the operator and the preconditioner must be POSITIVE DEFINITE. But I understand MINRES is tailored for the symmetric/hermitian-indefinite

docs for MINRES, only for positive definite?

2008-02-28 Thread Matthew Knepley
' ksp.setOperators(A) ksp.setFromOptions() ksp.solve(b,x) On 2/28/08, Barry Smith bsmith at mcs.anl.gov wrote: But does it require a positive definite preconditioner? Barry On Feb 28, 2008, at 9:32 AM, Matthew Knepley wrote: Docs are wrong. Matt

docs for MINRES, only for positive definite?

2008-02-28 Thread Matthew Knepley
. Matt Barry On Feb 28, 2008, at 11:28 AM, Matthew Knepley wrote: On Thu, Feb 28, 2008 at 11:25 AM, Lisandro Dalcin dalcinl at gmail.com wrote: Good point, the current code seems to require that... Its not the code, its the algorithm. It requires symmetry. Matt

Inner KSP never created on SNES

2008-12-02 Thread Matthew Knepley
I am not against requiring SetFromOptions() to create default objects. We already require it to set a default type. Matt On Tue, Dec 2, 2008 at 8:08 AM, Lisandro Dalcin dalcinl at gmail.com wrote: I've just noticed that now SNES does not create the inner KSP at the time SNESCreate() is

Inner KSP never created on SNES

2008-12-02 Thread Matthew Knepley
in hg view. Matt Barry On Dec 2, 2008, at 8:29 AM, Matthew Knepley wrote: I am not against requiring SetFromOptions() to create default objects. We already require it to set a default type. Matt On Tue, Dec 2, 2008 at 8:08 AM, Lisandro Dalcin dalcinl at gmail.com wrote: I've

an small and perhaps unimportant gotcha

2008-12-02 Thread Matthew Knepley
This is strange. I just tried KSP ex2 and after MatSetType() the type_name is set. What does *.setType() call? Matt On Tue, Dec 2, 2008 at 9:12 AM, Lisandro Dalcin dalcinl at gmail.com wrote: Other (really nice!) change I've noticed is that now we can call {Vec|Mat}Create() and next

some problems with 'make install'

2008-12-02 Thread Matthew Knepley
On Tue, Dec 2, 2008 at 1:31 PM, Satish Balay balay at mcs.anl.gov wrote: On Mon, 1 Dec 2008, Lisandro Dalcin wrote: I'm testing configure passing '--prefix=', and after building, I do $ su -c 'make install' But then, some files ('configure.log', 'RDict.db') are left with

fixes for VecPlaceArray/VecResetArray

2008-12-03 Thread Matthew Knepley
On Wed, Dec 3, 2008 at 5:04 PM, Lisandro Dalcin dalcinl at gmail.com wrote: I've pushed some fixes for VecPlaceArray/VecResetArray, trying to the pair a bit more smarter. However, the I'm not sure what's the use case in VecPlaceArray for accepting NULL for the array. Just in case, I did not

Fwd: [PETSC #18705] PETSc and Cygwin License (POSIX layer)

2008-12-04 Thread Matthew Knepley
I for one think it should be possible to remove 'make' from the toolchain, leaving us with only win32fe, which we distribute. Thus I think we could abandon cygwin once and for all. I would even be willing to write a \emph{make clone} to accomplish this, even though I am a committed enemy of make

Fwd: [PETSC #18705] PETSc and Cygwin License (POSIX layer)

2008-12-04 Thread Matthew Knepley
On Thu, Dec 4, 2008 at 12:41 PM, Satish Balay balay at mcs.anl.gov wrote: On Thu, 4 Dec 2008, Matthew Knepley wrote: I for one think it should be possible to remove 'make' from the toolchain, leaving us with only win32fe, which we distribute. Thus I think we could abandon cygwin once

[PETSC #18705] PETSc and Cygwin License (POSIX layer)

2008-12-04 Thread Matthew Knepley
toolbox, or maybe toolbag (buildstuff?). Matt I hate the term toolchain Barry On Dec 4, 2008, at 12:04 PM, Matthew Knepley wrote: I for one think it should be possible to remove 'make' from the toolchain, leaving us with only win32fe, which we distribute. Thus I think we could

  1   2   3   4   5   6   7   8   9   10   >