release Monday Dec 16th?

2008-12-11 Thread Matthew Knepley
Fine with me. Everything seems to be going for me. Matt On Thu, Dec 11, 2008 at 4:36 PM, Lisandro Dalcin dalcinl at gmail.com wrote: OK, fine with me... BTW, Had you any chance to look and try examples at tutorials/python ? On Thu, Dec 11, 2008 at 7:50 PM, Barry Smith bsmith at

Request related to SLEPc

2008-12-12 Thread Matthew Knepley
I need to understand what you use this variable for. If you are just linking, then you do not need it. If you want configure information, then you need it in order to get anything, so specifying in a file will not help. Matt On Fri, Dec 12, 2008 at 9:01 AM, Jose E. Roman jroman at dsic.upv.es

Request related to SLEPc

2008-12-12 Thread Matthew Knepley
On Fri, Dec 12, 2008 at 9:32 AM, Lisandro Dalcin dalcinl at gmail.com wrote: I would like to add that, despite the new buildsystem if by far better than the old one, PETSc has lost a nice feature of being able of being installed in a central location for multiple $PETSC_ARCH's . This feature

Request related to SLEPc

2008-12-12 Thread Matthew Knepley
If PETSC_ARCH is not set, I would just use the default arch generated by configure, like we do in the PETSc configure. Matt On Fri, Dec 12, 2008 at 11:29 AM, Jose E. Roman jroman at dsic.upv.es wrote: On 12/12/2008, Jose E. Roman wrote: SLEPc's configure.py uses the value of $PETSC_ARCH

how do you access the nonzero elements in SeqAIJ matrices

2008-08-04 Thread Matthew Knepley
On Mon, Aug 4, 2008 at 3:14 AM, Ahmed El Zein ahmed at azein.com wrote: I am working on a project where I would like to copy a sparse matrix in CSR format. I have tried MatGetRow() which works OK but I would really like to get pointers to the 3 arrays directly. I also tried MatGetRowIJ()

how do you access the nonzero elements in SeqAIJ matrices

2008-08-04 Thread Matthew Knepley
On Mon, Aug 4, 2008 at 9:58 AM, Ahmed El Zein ahmed at azein.com wrote: On Mon, 2008-08-04 at 04:17 -0500, Matthew Knepley wrote: On Mon, Aug 4, 2008 at 3:14 AM, Ahmed El Zein ahmed at azein.com wrote: I am working on a project where I would like to copy a sparse matrix in CSR format. I

Is c2html a required external package?

2008-08-08 Thread Matthew Knepley
On Fri, Aug 8, 2008 at 12:09 AM, Shi Jin jinzishuai at yahoo.com wrote: Hi there, I was trying to build petsc-dev on a IBM AIX but failed at c2html. I tried not to use it with --with-c2html=0 but it seemed still be needed. Why do we need it? I don't need to generate documentations for this

petsc and x-fem

2008-08-12 Thread Matthew Knepley
On Mon, Aug 4, 2008 at 6:49 AM, Techas techas at gmail.com wrote: Hi Matthew, I'm Sergio (X-FEM), we meet some weeks ago in Davis. I'm playing a little bit with petsc to evaluate how much work will take me mount my x-fem code on petsc. Sorry I am just replying now. I returned from Norway

PCApply_Shell ghosts...

2008-08-13 Thread Matthew Knepley
PETSc knows nothing about your domain, so it can't know what you might want ghosted. I think the thing to do is use a VecScatter to map the input Vec to a ghosted Vec (called a local vector in DA language). Matt On Wed, Aug 13, 2008 at 3:41 PM, Eric Chamberland Eric.Chamberland at

[petsc4py] What should be the 'default' communicator?

2008-08-25 Thread Matthew Knepley
I would still maintain that PETSC_COMM_WORLD is the correct default. There are better paradigms for embarassingly parallel operation, like Condor. PETSc is intended for parallel, domain decomposition runs. Matt On Mon, Aug 25, 2008 at 10:54 AM, Lisandro Dalcin dalcinl at gmail.com wrote:

[petsc4py] What should be the 'default' communicator?

2008-08-25 Thread Matthew Knepley
On Mon, Aug 25, 2008 at 1:22 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On Mon, Aug 25, 2008 at 1:08 PM, Matthew Knepley knepley at gmail.com wrote: I would still maintain that PETSC_COMM_WORLD is the correct default. There are better paradigms for embarassingly parallel operation, like

[petsc4py] What should be the 'default' communicator?

2008-08-25 Thread Matthew Knepley
Dalcin wrote: On Mon, Aug 25, 2008 at 6:22 PM, Matthew Knepley knepley at gmail.com wrote: I agree that people will do this, I just don't agree that it should be the default. Would you agree with the following: At petsc4py initialization (and after calling PetscInitialize()), I define

[petsc4py] What should be the 'default' communicator?

2008-08-25 Thread Matthew Knepley
Barry On Aug 25, 2008, at 4:43 PM, Lisandro Dalcin wrote: On Mon, Aug 25, 2008 at 6:22 PM, Matthew Knepley knepley at gmail.com wrote: I agree that people will do this, I just don't agree that it should be the default. Would you agree with the following: At petsc4py initialization

blocked index sets

2008-08-27 Thread Matthew Knepley
There is no concept of global for IS. They are purely serial. AO is the only global construct with indices. Matt On Wed, Aug 27, 2008 at 10:09 AM, Lisandro Dalcin dalcinl at gmail.com wrote: I believe we have to review the interface of ISBlock. Currently, ISBlockGetSize() return the number

blocked index sets

2008-08-27 Thread Matthew Knepley
the ISBlock version as well. Matt On Wed, Aug 27, 2008 at 1:13 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Aug 27, 2008, at 10:23 AM, Matthew Knepley wrote: There is no concept of global for IS. They are purely serial. AO is the only global construct with indices. This is kind

blocked index sets

2008-08-27 Thread Matthew Knepley
On Wed, Aug 27, 2008 at 1:08 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Aug 27, 2008, at 12:13 PM, Matthew Knepley wrote: On Wed, Aug 27, 2008 at 12:06 PM, Lisandro Dalcin dalcinl at gmail.com wrote: So, Do all us agree my proposed fix should be pushed? I'll wait for Matt comments

blocked index sets

2008-08-27 Thread Matthew Knepley
On Wed, Aug 27, 2008 at 2:10 PM, Barry Smith bsmith at mcs.anl.gov wrote: Even if an object (class) has NO collective operations, if, when you use that object, you must have partners are all other processes in a MPI_Comm then I think it is a good approach to have that be a parallel object

Block matrices and Schur complement preconditioning

2008-08-27 Thread Matthew Knepley
Just so I can remove this from my todo mail pile, it appears Barry is doing this. The actually Schur stuff is a new KSP thing so we do not link Mat to KSP, but it should work as you want. Matt On Mon, Aug 11, 2008 at 4:16 PM, Jed Brown jed at 59a2.org wrote: For context, I'm thinking of LNKS

Fortran integer declaration in petscdef.h

2008-04-07 Thread Matthew Knepley
On Mon, Apr 7, 2008 at 9:48 AM, Thomas DE-SOZA thomas.de-soza at edf.fr wrote: Hi, I was wondering if in $PETSC_DIR/include/finclude/petscdef.h : 50#if defined(PETSC_HAVE_MPIUNI) 51#define MPI_Comm PetscFortranInt 52#define PetscMPIInt PetscFortranInt

Preconditioning for saddle point problems

2008-04-29 Thread Matthew Knepley
1) I believe the Wathen-Elman-Silvester stuff is the best out there of the shelf. I love the review A.C. de Niet and F.W. Wubs Two preconditioners for saddle point problems in fluid flows Int. J. Num. Meth. Fluids 2007: 54: 355-377 2) Note in there that Augmented Lagrangian

Preconditioning for saddle point problems

2008-04-29 Thread Matthew Knepley
using / \hat AB \ \ B^T 0 / as a preconditioner is great since all the eigenvalues for the constraint are preserved. Matt Thanks, -- Boyce Matthew Knepley wrote: 1) I believe the Wathen-Elman-Silvester stuff is the best out there of the shelf. I love the review

Preconditioning for saddle point problems

2008-04-29 Thread Matthew Knepley
On Tue, Apr 29, 2008 at 12:54 PM, Boyce Griffith griffith at cims.nyu.edu wrote: Matthew Knepley wrote: On Tue, Apr 29, 2008 at 12:28 PM, Boyce Griffith griffith at cims.nyu.edu wrote: Hi, Matt et al. -- Do people ever use standard projection methods as preconditioners

adding SNESSetLinearSolve()

2007-11-01 Thread Matthew Knepley
Actually, making the regularization parameter independent for each process is much more efficient. Gene Golub had a poster on this at the last SIAM CSE meeting. Matt On Nov 1, 2007 9:24 AM, Lisandro Dalcin dalcinl at gmail.com wrote: On 10/31/07, Barry Smith bsmith at mcs.anl.gov wrote:

error checking macros with optimization

2007-11-23 Thread Matthew Knepley
On Nov 23, 2007 10:25 AM, Lisandro Dalcin dalcinl at gmail.com wrote: I would to propose some changes to error checking macros for optimized builds. 1.- SETERRQXXX: define them as #define SETERRQ[1|2|..](ierr,...) return ierr I think this is fine. However for it to matter, you need 2.

move maint to bin/maint

2007-11-24 Thread Matthew Knepley
Are they all executables? Matt On 11/23/07, Barry Smith bsmith at mcs.anl.gov wrote: Can we move maint to bin/maint? Barry -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments

move maint to bin/maint

2007-11-24 Thread Matthew Knepley
I thought he meant LFS. Matt On Nov 24, 2007 3:52 PM, Barry Smith bsmith at mcs.anl.gov wrote: Seems unlike that it would have the rule no subdirectory for bin. Especially since the H stands for hierarchical. Barry On Nov 24, 2007, at 3:49 PM, Dmitry Karpeev wrote: From what I

warn message in log summary

2007-11-26 Thread Matthew Knepley
On Nov 26, 2007 4:19 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On 11/26/07, Barry Smith bsmith at mcs.anl.gov wrote: I've looked long and hard for a PETSc bug that would cause this problem. No luck. It seems to happen mostly (only?) on certain machines. Ups! I've found a possilbe

warn message in log summary

2007-11-26 Thread Matthew Knepley
That seems fine. Matt On Nov 26, 2007 4:56 PM, Lisandro Dalcin dalcinl at gmail.com wrote: On 11/26/07, Matthew Knepley knepley at gmail.com wrote: On Nov 26, 2007 4:19 PM, Lisandro Dalcin dalcinl at gmail.com After surfing the code and hacked PetscLogPrintSummary(), I added

Compiling on Solaris 10 for X86

2007-05-10 Thread Matthew Knepley
Nothing is wrong with the PETSc build or libraries. However, your make (probably not gmake) has a problem with the make rule that cleans up files in bmake/common/rules:359. You can just ignore it (or tell us what is wrong). Thanks, Matt On 5/10/07, Yi-Feng Adam Zhang Adam.Zhang at

Compiling on Solaris 10 for X86

2007-05-10 Thread Matthew Knepley
On 5/10/07, Adam Zhang Adam.Zhang at sun.com wrote: Thank you for input. I will try to use gmake later. If this compilation is ok. can I release this as a Solaris package on the internet? I guess. What is it for? Thanks, Matt Regards, Adam Matthew Knepley wrote: Nothing

DDM with PETSc

2007-03-31 Thread Matthew Knepley
I would look at some PETSc examples, for instance http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/ksp/ksp/examples/tutorials/ex32.c.html Matt On 3/31/07, Waad Subber wsubber at connect.carleton.ca wrote: Hello everyone: I am new to PETSC. Reading the tutorials,

changes to PETSc-dev bmake system and library locations

2007-06-10 Thread Matthew Knepley
This is very vague concerning the structure of exernalpackages. I cannot tell where libraries are supposed to end up, and how/when/why they might be moved. It seems that the directory information coming from PETSc/package.py has also changed. Matt On 6/8/07, Barry Smith bsmith at mcs.anl.gov

changes to PETSc-dev bmake system and library locations

2007-06-11 Thread Matthew Knepley
(630) 252-5986 fax 9700 S. Cass Ave.tmunson at mcs.anl.gov Argonne, IL 60439http://www.mcs.anl.gov/~tmunson On Sun, 10 Jun 2007, Barry Smith wrote: On Sun, 10 Jun 2007, Matthew Knepley wrote: This is very vague concerning the structure

PETSc sparsity

2007-06-15 Thread Matthew Knepley
Also, I have talked to Wolfgang many times about this. I am a firm believer in eliminating the boundary values during assembly at the element level. PETSc provides an easy mechanism for this. By default, all negative indices in calls to VecSetValues and MatSetValues are ignored. Matt On

invitation to give talk at SDSC summer institute

2007-06-20 Thread Matthew Knepley
I am not going to ICIAM anymore (since international travel has become impossible), so I am available I guess. Matt On 6/20/07, Barry Smith bsmith at mcs.anl.gov wrote: Amit, Thank you for the invitation to speak. Unfortunately I am not available at that time, and I think a couple of

test fail (because of bad makefile?)

2007-07-06 Thread Matthew Knepley
Pushed a fix. Matt On 7/6/07, Lisandro Dalcin dalcinl at gmail.com wrote: Can someone look at this? testexamples_C in: /u/dalcinl/Devel/PETSc/petsc-dev/src/dm/adda make[4]: *** No rule to make target `testexamples_C', needed by `tree'. Stop. -- Lisandro Dalc?n --- Centro

configure.py script

2007-07-17 Thread Matthew Knepley
We do generate makefiles, just not all of them. Automake is an unreadable, inpenetrable mess that should never be used by anyone. Thus, we don't use it. Thanks, Matt On 7/17/07, Sumit Vaidya sumit_vaidya at persistent.co.in wrote: Hi, I have downloaded the version 2.3.3-p3.

AOApplicationToPetscPermuteInt() help needed

2007-07-17 Thread Matthew Knepley
Done. Matt On 7/17/07, Barry Smith bsmith at mcs.anl.gov wrote: Will whoever wrote the manual pages for all the AOApplicationToPetscPermute*() routines please add more documentation explaining what they do? Maybe some algebraic formula or something: Permutes an array of blocks of reals

COVERITY static source code analysis

2007-07-19 Thread Matthew Knepley
I'm for it. Barry, do you want to mail scan-admin at coverity.com? Matt On 7/19/07, Lisandro Dalcin dalcinl at gmail.com wrote: Have any of you ever consider asking PETSc being included here, as it is an open source project? http://scan.coverity.com/index.html From many sources (mainly

Cray and PETSc

2007-08-23 Thread Matthew Knepley
On 8/23/07, Adrian Tate adrian at cray.com wrote: Hi Barry Thanks for your reply. We will expect to provide the new full releases of PETSc as soon as we possibly can. Obviously there is something of a lag because Cray need to test, package and release the library, but we'll try to do so

rename ISLocalToGlobalMapping?

2007-08-27 Thread Matthew Knepley
*petscPerm; } AO_Mapping; typedef struct { PetscInt N; PetscInt *app,*petsc; /* app[i] is the partner for the ith PETSc slot */ /* petsc[j] is the partner for the jth app slot */ } AO_Basic; Barry On Mon, 27 Aug 2007, Matthew Knepley wrote: On 8/27/07, Barry

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
Yes, definitely. Go ahead and push it. Matt On 8/28/07, Lisandro Dalcin dalcinl at gmail.com wrote: Does it make sense to change KSPSkipConverged to return KSP_CONVERGED_ITS if iternum==maxit ? KSP_DIVERGED_ITS means convergence failure, but IMHO, KSPSkipConverged should not imply

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
On 8/28/07, Lisandro Dalcin dalcinl at gmail.com wrote: On 8/28/07, Matthew Knepley knepley at gmail.com wrote: Yes, definitely. Go ahead and push it. I started to try this by implementing first on petsc4py with petsc-2.3.3-p4, by solving a trivial SPD diagonal system { A_ii = 1/(i+1

KSPSkipConverged

2007-08-28 Thread Matthew Knepley
On 8/28/07, Lisandro Dalcin dalcinl at gmail.com wrote: On 8/28/07, Matthew Knepley knepley at gmail.com wrote: Todd's response was that there are reasons it can stop early, like happy breakdown, so we really need to test this with the dev version to track down this behavior. OK, my tests

PETSc/ridgeSim 3D scaling on BG/L

2006-04-02 Thread Matthew Knepley
Thats awesome. 80% true efficiency for 1K nodes on under 1 million unknowns. Thats only 1,000 unknowns per proc. If we were cranking the problems size, I am sure you'd be up in the 90s. Hopefully this summer we can deliver you much improved serial performance as well. BTW, I want to talk sometime

BuildSystem question

2005-09-25 Thread Matthew Knepley
S V N Vishwanathan vishy at mail.rsise.anu.edu.au writes: Hi! I am tying myself up in knots in trying to figure out your BuildSystem and how to adapt it to LINEAL. PETSc is one primary component of LINEAL but going forward there will be others (viz. TAO, OOQP). I want that my users should

BuildSystem question

2005-09-26 Thread Matthew Knepley
S V N Vishwanathan vishy at mail.rsise.anu.edu.au writes: Ideally, I would like to use your BuildSystem to hack up a quick build for Numerix. Then as we add and integrate more tools in, I would work on making it fancier. Cool, we can do that. Here is the plan I came up with (after

tracking petsc

2005-11-03 Thread Matthew Knepley
Simon Burton simon at arrowtheory.com writes: On Thu, 03 Nov 2005 17:49:46 -0600 Matthew Knepley knepley at mcs.anl.gov wrote: Simon Burton simon at arrowtheory.com writes: This was patched in 2.3.0, but maybe you have an unpatched version. It is easy to check. The symbol comes

more inlining ?

2005-11-03 Thread Matthew Knepley
Barry Smith bsmith at mcs.anl.gov writes: Should we inline VecSetValues, MatSetValues, and friends? If 1) at least gcc supports it (we will have to maintain the #define for those that are not C99 I think) 2) gdb handles it I also have another proposal. Dmitry and I are going

ufuncs, iterators

2005-08-18 Thread Matthew Knepley
Simon Burton simon at arrowtheory.com writes: Hi, Along the lines of python's numarray [1], we need some way of operating pointwise (and inner/outer operations) on Mat/Vec objects. In particular, we need things like the following: (a) v=add.reduce(m) (sum along rows/cols of a Mat to

ufuncs, iterators

2005-08-19 Thread Matthew Knepley
S V N Vishwanathan vishy at mail.rsise.anu.edu.au writes: Hi! BTW: possibly related note, are you using dense matrices sometimes to represent just 2-arrays; that is, not as representations of linear operators. If so, I do not think this is the correct approach! Conceptually PETSc Mat's are

ufuncs, iterators

2005-08-19 Thread Matthew Knepley
Simon Burton simon at arrowtheory.com writes: On Thu, 18 Aug 2005 19:57:02 -0500 Matthew Knepley knepley at mcs.anl.gov wrote: (d) m=exp(m) (pointwise exp) We could add pointwise operations just like the VecPointwise*(). Matt Yes, we need MatPointwiseMult

ufuncs, iterators

2005-08-20 Thread Matthew Knepley
Simon Burton simon at arrowtheory.com writes: Oh, that should probably read: exp( -1/2\sigma^{2} ||x1_{i} - x2_{j}||_{2}^{2}) And when we vectorize this operation: ||x1_i - x2_j||^2 = ||x1_i||^2 + ||x2_i||^2 + 2*(x1_i,x2_j) and the last term is the ip matrix. It seems that this

Bitkeeper going commercial-only

2005-04-12 Thread Matthew Knepley
Aron Ahmadia aron.ahmadia at gmail.com writes: So this is really crappy news, I really liked BK. How is this going to affect PETSc and what are your future short-term and long-term plans? I have to figure out what I'm doing with COW, and I was wondering how you guys were dealing with this.

--download-prometheus working again

2005-04-20 Thread Matthew Knepley
Barry Smith bsmith at mcs.anl.gov writes: After endless futzing I've gotten --download-prometheus working with both --with-clanguage=C and C++ with petsc-dev on my Mac. Should work everywhere :-). Satish, could you try adding it to a nightly build? You also need --download-parmetis

[petsc-dev] make test error

2012-05-01 Thread Matthew Knepley
On Tue, May 1, 2012 at 3:28 PM, John Mousel john.mousel at gmail.com wrote: I just pulled petsc-dev and am getting the following error in the tests. Pushed fix. Matt John Running test examples to verify correct installation Using

[petsc-dev] passing a function pointer?

2012-05-03 Thread Matthew Knepley
On Thu, May 3, 2012 at 10:33 PM, Cinstance cinstance at gmail.com wrote: I am using PETSC as the linear solver for a CFD project. It is part of an integrated frame work, and the project is developed under MS Windows. I use PETSC as prebuild libs in Visual Studio, so I cannot debug into PETSC

[petsc-dev] Petsc VecCopy code.

2012-05-04 Thread Matthew Knepley
On Fri, May 4, 2012 at 2:27 PM, Yoo, Andy yoo2 at llnl.gov wrote: Hi, ** ** I have written a two-dimensional MatVec routine ( http://dl.acm.org/citation.cfm?id=2063469) a while ago. Any chance of getting that into PETSc proper? Now, I am noticing there is a small memory leak

[petsc-dev] Error during KSPDestroy

2012-05-06 Thread Matthew Knepley
On Sun, May 6, 2012 at 7:28 AM, Alexander Grayver agrayver at gfz-potsdam.dewrote: Hello, I use KSP and random rhs to compute largest singular value: 1) Is this the whole program? If not, this can be caused by memory corruption somewhere else. This is what I suspect. 2) You can put in

[petsc-dev] Error during KSPDestroy

2012-05-06 Thread Matthew Knepley
On Sun, May 6, 2012 at 8:42 AM, Alexander Grayver agrayver at gfz-potsdam.dewrote: ** On 06.05.2012 14:27, Matthew Knepley wrote: On Sun, May 6, 2012 at 7:28 AM, Alexander Grayver agrayver at gfz-potsdam.de wrote: Hello, I use KSP and random rhs to compute largest singular value: 1

[petsc-dev] Error during KSPDestroy

2012-05-06 Thread Matthew Knepley
is not stack'd, malloc'd or (recently) free'd On 06.05.2012 15:21, Alexander Grayver wrote: On 06.05.2012 15:07, Matthew Knepley wrote: Hello, I use KSP and random rhs to compute largest singular value: 1) Is this the whole program? If not, this can be caused by memory corruption

[petsc-dev] DM RefineLevel and CoarsenLevel

2012-05-06 Thread Matthew Knepley
That sounds fine to me. Matt On Sun, May 6, 2012 at 3:44 PM, Jed Brown jedbrown at mcs.anl.gov wrote: Should the refinement level be copied over by DMCoarsen (and the coarsen level be copied by DMRefine)? It's useful for diagnostics to be able to define a universal level. If I use PCMG

[petsc-dev] DM RefineLevel and CoarsenLevel

2012-05-06 Thread Matthew Knepley
On Sun, May 6, 2012 at 5:57 PM, Jed Brown jedbrown at mcs.anl.gov wrote: Now the next round: For semi-coarsening, we used to have stuff like -da_refine_hierarchy_x 1,1,3 -da_refine_hierarchy_y 2,2,1 -da_refine_hierarchy_z 2,2,1. Two changes make this harder now: 1. You essentially got rid

[petsc-dev] Attaching a near null space to an IS

2012-05-06 Thread Matthew Knepley
On Sun, May 6, 2012 at 7:15 PM, Jed Brown jedbrown at mcs.anl.gov wrote: Matt introduced this concept, he says the IS is a better place to attach things. http://petsc.cs.iit.edu/petsc/petsc-dev/rev/2ad289ac99e0 I don't understand why the IS is better (because it's mostly immutable?). I'm

[petsc-dev] Attaching a near null space to an IS

2012-05-06 Thread Matthew Knepley
is the right thing to do. It allows other people to stick their own stuff in place of ours easily, without junking it up with a lots of wrapper objects. Matt Barry On May 6, 2012, at 6:30 PM, Matthew Knepley wrote: On Sun, May 6, 2012 at 7:15 PM, Jed Brown jedbrown at mcs.anl.gov wrote

[petsc-dev] Error during KSPDestroy

2012-05-07 Thread Matthew Knepley
, Matthew Knepley wrote: On Sun, May 6, 2012 at 9:24 AM, Alexander Grayveragrayver at gfz-potsdam. **de agrayver at gfz-potsdam.de wrote: Hm, valgrind gives a lot of output like that (see full log in previous message): Can you run this with --download-f-blas-lapack? This sounds much more

[petsc-dev] discrepancy with bcgs, mpi, cusp, and no preconditioning

2012-05-10 Thread Matthew Knepley
On Thu, May 10, 2012 at 3:45 PM, Chetan Jhurani chetan.jhurani at gmail.comwrote: Hi wizards of -ops- indirections, I'm trying to use bcgs without preconditioning, for now, and the iterations using -vec_type cusp -mat_type mpiaijcusp don't match serial or non-GPU options. I've attached the

[petsc-dev] Spelling Chebyshev

2012-05-10 Thread Matthew Knepley
On Thu, May 10, 2012 at 4:27 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Thu, May 10, 2012 at 3:23 PM, Sean Farley sean at mcs.anl.gov wrote: It would seem to me that the best way to fix this is to just spell his name the way he did: ? ???. All we need to do is

[petsc-dev] As if we need convincing

2012-05-11 Thread Matthew Knepley
http://www.250bpm.com/blog:4 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 -- next part -- An HTML attachment was scrubbed... URL:

[petsc-dev] SNESDMComputeJacobian()

2012-05-15 Thread Matthew Knepley
On Tue, May 15, 2012 at 10:11 PM, Jed Brown jedbrown at mcs.anl.gov wrote: *- ptr - pointer to a structure that must have a DM as its first entry.* *This ptr must have been passed into SNESDMComputeFunction() as the context.* This is considered very bad form now since

[petsc-dev] patch for block size in VecLoad_HDF5

2012-05-16 Thread Matthew Knepley
On Wed, May 16, 2012 at 6:27 PM, Blaise Bourdin bourdin at lsu.edu wrote: Hi, ex19 in src/vec/vec/examples/tutorials is still broken, since changing block size of a Vec is no longer supported. the mechanism in VecLoad_HDF5 was to rely on an undocumented flag (-vecload_block_size) Shall we

[petsc-dev] A PETSc question on scicomp I can't completely answer (behavior of PETScSetxxx)

2012-05-23 Thread Matthew Knepley
On Wed, May 23, 2012 at 7:23 AM, Aron Ahmadia aron at ahmadia.net wrote: http://scicomp.stackexchange.com/questions/2303/petscs-xxxsetxxx-methods-own-pointer-or-copy-values My understanding is that Set for the most part acts like a pass-by-reference, and my answer reflects that (

[petsc-dev] Multigrid is confusing

2012-05-24 Thread Matthew Knepley
On Thu, May 24, 2012 at 4:16 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Thu, May 24, 2012 at 3:10 PM, Barry Smith bsmith at mcs.anl.gov wrote: Absolutely. And if it turns out to be too much of a pain to write and maintain such a kernel there is something wrong with our programming model

[petsc-dev] Need a CFLAGS that does is *not* included in the link

2012-05-24 Thread Matthew Knepley
On Thu, May 24, 2012 at 10:09 PM, Jed Brown jedbrown at mcs.anl.gov wrote: Building PETSc with clang++ produces a warning about compiling *.c as C++ being deprecated. To silence the warning, we would need to pass -x c++ to the compiler, but NOT to the linker. CFLAGS is currently also passed to

[petsc-dev] Multigrid is confusing

2012-05-25 Thread Matthew Knepley
On Thu, May 24, 2012 at 10:25 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Thu, May 24, 2012 at 9:17 PM, Barry Smith bsmith at mcs.anl.gov wrote: Sarcasm on -- Yes, of course. It is essentially done, write it up and I want it on my desk by 8 AM Friday. -- Sarcasm off So what if

[petsc-dev] ctetgen into its own repository

2012-05-25 Thread Matthew Knepley
On Fri, May 25, 2012 at 4:03 PM, Barry Smith bsmith at mcs.anl.gov wrote: Because of the tetgen license we cannot include ctetgen directly in the PETSc tarball. Thus we have forked it off into its own repository and it is available for users as --download-ctetgen Developers may choose

[petsc-dev] OpenMP compiler options

2012-05-29 Thread Matthew Knepley
On Tue, May 29, 2012 at 3:52 PM, Jed Brown jedbrown at mcs.anl.gov wrote: The OpenMP flags do not definitively identify that OpenMP is used. In particular, IBM XL interprets Cray's option -h omp as being equivalent to -soname omp, then silently ignores the Open MP pragmas. We can perhaps fix

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 11:09 AM, Barry Smith bsmith at mcs.anl.gov wrote: Fredrik, This question belongs on petsc-dev at mcs.anl.gov since it involves additions/extensions to PETSc so I am moving the discussion over to there. We have not done the required work to have ghosted

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 12:40 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Mon, Feb 6, 2012 at 20:35, Matthew Knepley knepley at gmail.com wrote: VecCreate() VecSetType() VecSetGhosted() Is there a compelling reason for memory to be allocated before VecSetUp()? Deferring it would make

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 12:47 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Mon, Feb 6, 2012 at 21:42, Matthew Knepley knepley at gmail.com wrote: I don't like this because it would mean calling VecSetUp() all over the place. Couldn't the ghosting flag be on the same level as the sizes

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 1:11 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 12:47 PM, Jed Brown wrote: On Mon, Feb 6, 2012 at 21:42, Matthew Knepley knepley at gmail.com wrote: I don't like this because it would mean calling VecSetUp() all over the place. Couldn't

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 1:23 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 1:14 PM, Matthew Knepley wrote: On Mon, Feb 6, 2012 at 1:11 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 12:47 PM, Jed Brown wrote: On Mon, Feb 6, 2012 at 21:42, Matthew

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 1:30 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 1:27 PM, Matthew Knepley wrote: On Mon, Feb 6, 2012 at 1:23 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 1:14 PM, Matthew Knepley wrote: On Mon, Feb 6, 2012 at 1:11 PM

[petsc-dev] [petsc-users] VECMPICUSP with ghosted vector

2012-02-06 Thread Matthew Knepley
On Mon, Feb 6, 2012 at 1:42 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 1:39 PM, Matthew Knepley wrote: On Mon, Feb 6, 2012 at 1:30 PM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 6, 2012, at 1:27 PM, Matthew Knepley wrote: On Mon, Feb 6, 2012 at 1:23 PM

[petsc-dev] petsc-dev on bitbucket

2012-02-07 Thread Matthew Knepley
On Tue, Feb 7, 2012 at 8:14 PM, Sean Farley sean at mcs.anl.gov wrote: Yoo, i've switched petsc-dev over to bitbucket since petsc.cs.itt.edu is down pull/push from https://bitbucket.org/petsc/petsc-dev Sean - Barry impersonating me I have pushed to bring it up to date with the latest

[petsc-dev] FAS with TS

2012-02-07 Thread Matthew Knepley
On Tue, Feb 7, 2012 at 11:50 PM, Jed Brown jedbrown at mcs.anl.gov wrote: On Wed, Feb 8, 2012 at 04:19, Peter Brune prbrune at gmail.com wrote: On Tue, Feb 7, 2012 at 6:40 PM, Jed Brown jedbrown at mcs.anl.gov wrote: Suppose I want to solve a finite time step problem using FAS. The time

[petsc-dev] FAS with TS

2012-02-08 Thread Matthew Knepley
On Wed, Feb 8, 2012 at 12:35 AM, Jed Brown jedbrown at mcs.anl.gov wrote: On Wed, Feb 8, 2012 at 08:57, Matthew Knepley knepley at gmail.com wrote: I just want to point out that Jed envisions that coefficients (and maybe subproblems, etc) cannot be accommodated on the same DM. I agree

[petsc-dev] petsc-dev on bitbucket

2012-02-08 Thread Matthew Knepley
On Wed, Feb 8, 2012 at 7:54 AM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 8, 2012, at 6:43 AM, Satish Balay wrote: On Tue, 7 Feb 2012, Barry Smith wrote: On Feb 7, 2012, at 9:09 PM, Sean Farley wrote: I'm sure Jed (or Matt in his prime) could have run over to IIT and

[petsc-dev] petsc-dev on bitbucket

2012-02-08 Thread Matthew Knepley
On Wed, Feb 8, 2012 at 9:39 AM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 8, 2012, at 9:06 AM, Matthew Knepley wrote: On Wed, Feb 8, 2012 at 7:54 AM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 8, 2012, at 6:43 AM, Satish Balay wrote: On Tue, 7 Feb 2012, Barry Smith wrote

[petsc-dev] petsc-dev on bitbucket

2012-02-08 Thread Matthew Knepley
On Wed, Feb 8, 2012 at 10:54 AM, Sean Farley sean at mcs.anl.gov wrote: There is only a 1-level hierarchy based on a top level account. Sean created 'petsc' for our stuff. We can create many, so that we have 'petsc-release', 'petsc-private', etc. if we want. Of course, I want traditional

[petsc-dev] petsc-dev on bitbucket

2012-02-08 Thread Matthew Knepley
On Wed, Feb 8, 2012 at 11:32 AM, Sean Farley sean at mcs.anl.gov wrote: No I mean I want petsc/releases/petsc-3.1 petsc/tools/parsing/BarrysNewHTMLMunger Sounds like you want a traditional filesystem hierarchy and like I said before, your request will fall on deaf ears. This is the

[petsc-dev] MatMatMult gives different results

2012-02-09 Thread Matthew Knepley
On Thu, Feb 9, 2012 at 2:48 AM, Alexander Grayver agrayver at gfz-potsdam.dewrote: ** lib/petsc-dev hg pull -u pulling from http://petsc.cs.iit.edu/petsc/petsc-dev searching for changes no changes found Does it take some time till it's there? You have to remake Fortran stubs: make

[petsc-dev] Someone still cannot use the annotate link

2012-02-09 Thread Matthew Knepley
Clearly, karpeev was the only use who altered those lines: http://petsc.cs.iit.edu/petsc/BuildSystem/rev/42cbb1d6192f http://petsc.cs.iit.edu/petsc/BuildSystem/annotate/d295489bd56e/config/packages/MOAB.py Matt -- What most experimenters take for granted before they begin their experiments

[petsc-dev] petsc-dev on bitbucket

2012-02-09 Thread Matthew Knepley
On Thu, Feb 9, 2012 at 8:26 PM, Satish Balay balay at mcs.anl.gov wrote: On Wed, 8 Feb 2012, Sean Farley wrote: Hell, if you *really* want to, just create the account: petsc-release(s) then the URL would be http://bitbucket.org/petsc-release/petsc-3.1 Actually, it's even

[petsc-dev] MatMatMult gives different results

2012-02-09 Thread Matthew Knepley
On Thu, Feb 9, 2012 at 9:10 PM, Satish Balay balay at mcs.anl.gov wrote: On Thu, 9 Feb 2012, Jed Brown wrote: It's all workable, but I'm not seeing a clear advantage to keeping the primary at petsc.cs.iit.edu, currently the 2 reasons that are offered agains petsc.cs.iit.edu are: - ssh

[petsc-dev] DMGetMatrix -- DMGetMatrices?

2012-02-10 Thread Matthew Knepley
On Fri, Feb 10, 2012 at 12:05 AM, Dmitry Karpeev karpeev at mcs.anl.govwrote: 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

[petsc-dev] petsc-dev on bitbucket

2012-02-10 Thread Matthew Knepley
The thread has become too deep for me to read, hence the top posting. Barry's question is the right one: What do we gain by changing? 1) Reliability and Availability Barry, you should know that this crap about petsc.cs being backed up is farcical. We would have the same situation we had

[petsc-dev] petsc-dev on bitbucket

2012-02-10 Thread Matthew Knepley
On Fri, Feb 10, 2012 at 9:23 AM, Barry Smith bsmith at mcs.anl.gov wrote: On Feb 10, 2012, at 9:13 AM, Matthew Knepley wrote: The thread has become too deep for me to read, hence the top posting. Barry's question is the right one: What do we gain by changing? 1) Reliability

[petsc-dev] petsc-dev on bitbucket

2012-02-10 Thread Matthew Knepley
On Fri, Feb 10, 2012 at 9:35 AM, Satish Balay balay at mcs.anl.gov wrote: On Fri, 10 Feb 2012, Barry Smith wrote: On Feb 10, 2012, at 9:13 AM, Matthew Knepley wrote: The thread has become too deep for me to read, hence the top posting. Barry's question is the right one: What do

[petsc-dev] petsc-dev on bitbucket

2012-02-10 Thread Matthew Knepley
On Fri, Feb 10, 2012 at 9:46 AM, Satish Balay balay at mcs.anl.gov wrote: Matt, until we move to bitbucket - can you continue to push to petsc.cs.iit - and not bitbucket? [for now - pushes to petsc.cs.iit are automatically pushed to bitbucket] I thought we were switched already :) I cannot

<    1   2   3   4   5   6   7   8   9   10   >