On Thu, Nov 24, 2011 at 6:37 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Nov 24, 2011, at 4:41 PM, Matthew Knepley wrote:
On Thu, Nov 24, 2011 at 4:09 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Jed,
Let's stop arguing about whether MPI is or is not a good base for the
next
On Wed, Nov 30, 2011 at 12:41 AM, Dave Nystrom dnystrom1 at comcast.netwrote:
I have a linear system in a code that I have interfaced to petsc that is
taking about 80 percent of the run time per timestep. This linear system
is
a symmetric block banded matrix where the blocks are 2x2. The
On Wed, Nov 30, 2011 at 4:53 PM, Sean Farley sean at mcs.anl.gov wrote:
No, but perhaps I should. My understanding is that those lines only need
to be changed when ParMetis is included as a subproject, as they seem to
have configured the build system to use relative paths.
I would still
On Sat, May 14, 2011 at 11:51 AM, Barry Smith bsmith at mcs.anl.gov wrote:
I have a proposal for a moderate/major change for how we handle
knowledge of subfields/splits in PETSc.
Currently one can call PCFieldSplitSetIS() or PCFieldSplitSetFields() to
indicate the splitting of
On Sat, May 14, 2011 at 10:14 PM, Barry Smith bsmith at mcs.anl.gov wrote:
In preparation for adding field information to PetscLayout (as discussed
below in the previous email) I have fixed up PETSc-dev so that
VecDuplicate(), MatDuplicate() and MatGetVecs() results in new objects that
On Sun, May 15, 2011 at 10:22 AM, Jed Brown jed at 59a2.org wrote:
On Sun, May 15, 2011 at 17:17, Satish Balay balay at mcs.anl.gov wrote:
can't you get cc,cflags also from python distutils?
Should there be a --with-profile=numpy to have BuildSystem interrogate
numpy for compilers and
On Mon, May 16, 2011 at 8:26 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On May 16, 2011, at 5:57 PM, Matthew Knepley wrote:
On Sat, May 14, 2011 at 10:14 PM, Barry Smith bsmith at mcs.anl.gov
wrote:
In preparation for adding field information to PetscLayout (as
discussed below
casting to (int) is not right.
Matt
Barry
On May 30, 2011, at 9:53 AM, Jed Brown wrote:
On Mon, May 30, 2011 at 16:38, Matthew Knepley knepley at gmail.com
wrote:
You wrote the code that takes PetscInt p!!! If you do not want that,
change it to PetscReal.
Right, I
Each process calls FieldSplitSetIS() with the indices for THAT FIELD that
are owned by that process. Each
process call FieldSplitSetIS() n+1 times if you have n+1 fields. Identifying
fields with processes is a mistake.
Matt
On Wed, Mar 9, 2011 at 1:50 PM, Thomas Witkowski
Thomas.Witkowski at
! Well I do.
satish
On Thu, 10 Mar 2011, Matthew Knepley wrote:
If we use all Python to build on Windows, with win32fe and the Windows
compilers, what
dependencies on Cygwin do we have left?
Thanks,
Matt
--
What most experimenters take for granted before
Yes, there was a good amount of effort spent on it. Is it broken?
Thanks,
Matt
On Fri, Mar 11, 2011 at 10:54 AM, Jose E. Roman jroman at dsic.upv.es wrote:
We are about to develop an interface to pARMS in petsc-dev, similar to that
of Hypre. But there is already a pARMS.py module. Has
On Fri, Mar 11, 2011 at 11:14 AM, Jed Brown jed at 59a2.org wrote:
On Fri, Mar 11, 2011 at 18:08, Matthew Knepley knepley at gmail.com wrote:
This does not help AT ALL. How is it bad? I use it all the time just fine,
so I was not changing it unless there were complaints.
If you can just
On Mon, Mar 14, 2011 at 11:00 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Mar 14, 2011, at 10:33 AM, Jed Brown wrote:
On Mon, Mar 14, 2011 at 16:22, Barry Smith bsmith at mcs.anl.gov wrote:
What bothers me about this one is that private has no hierarchy and will
eventually become
Does anyone have a favorite? Blaise is going to have a student add JSON
support
to PetscOptions.
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
Looks like there are at least two:
http://stackoverflow.com/questions/819452/c-json-library-in-popular-linux-distros
Matt
On Mon, Mar 14, 2011 at 1:35 PM, Lisandro Dalcin dalcinl at gmail.com wrote:
On 14 March 2011 15:31, Matthew Knepley knepley at gmail.com wrote:
Does anyone have
On Tue, Mar 15, 2011 at 1:20 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 19:07, Matthew Knepley knepley at gmail.com wrote:
not checked by my editor.
Why doesn't your editor show which branch is active for the current
PETSC_ARCH?
Does yours? My Emacs does not follow
On Tue, Mar 15, 2011 at 1:22 PM, Matthew Knepley knepley at gmail.com wrote:
On Tue, Mar 15, 2011 at 1:20 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 19:07, Matthew Knepley knepley at gmail.com wrote:
not checked by my editor.
Why doesn't your editor show which branch
On Tue, Mar 15, 2011 at 1:32 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 19:22, Matthew Knepley knepley at gmail.com wrote:
Does yours? My Emacs does not follow #if, and it would not know what
branch anyway.
Semantic can parse petscconf.h, then it knows. You are familiar
On Tue, Mar 15, 2011 at 1:49 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 19:44, Matthew Knepley knepley at gmail.com wrote:
Verification of preprocess code is exactly as intensive, since you have to
compile each version. What is your point here?
You can see all possible
Here is the scenario. I have a FieldSplit solve with 4 splits. For the 4th
split, I have
a good preconditioner matrix. In the linear case, KSP, I am just doing this:
err = KSPSetUp(_ksp);
err = KSPGetPC(_ksp, pc);
err = PCFieldSplitGetSubKSP(pc, num, ksps);
MatStructure flag;
On Tue, Mar 15, 2011 at 3:46 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 21:39, Matthew Knepley knepley at gmail.com wrote:
Here is the scenario. I have a FieldSplit solve with 4 splits. For the 4th
split, I have
a good preconditioner matrix. In the linear case, KSP, I am
On Tue, Mar 15, 2011 at 3:57 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 21:51, Matthew Knepley knepley at gmail.com wrote:
Really? That seems perverse. I have to build some global PC matrix, with
pieces I do not even have. How is that going to happen? Is
this really how we
On Tue, Mar 15, 2011 at 4:10 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 22:07, Matthew Knepley knepley at gmail.com wrote:
I guess I have no problem with MatGetSubMatrix() being the interface for
this, since all solvers get the Mat object. However, my point was
really
On Tue, Mar 15, 2011 at 5:29 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 23:28, Matthew Knepley knepley at gmail.com wrote:
What I have now is a MatShell that implements MatGetSubMatrix(),
I did this for a couple years, but now we have MatNest, why not use it?
If it does
On Tue, Mar 15, 2011 at 6:00 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 23:50, Matthew Knepley knepley at gmail.com wrote:
That is not what I have here. I want to change a portion, but outside of
that I want it to fall back to a global matrix.
You were going to put your
On Tue, Mar 15, 2011 at 6:09 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Mar 15, 2011, at 6:03 PM, Matthew Knepley wrote:
On Tue, Mar 15, 2011 at 6:00 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 15, 2011 at 23:50, Matthew Knepley knepley at gmail.com
wrote:
That is not what I
On Tue, Mar 15, 2011 at 7:06 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Stuff like this in BlasLapack.py
def getSharedFlag(self,cflags):
for flag in ['-PIC', '-fPIC', '-KPIC', '-qpic']:
if cflags.find(flag) =0: return flag
return ''
def getPrecisionFlag(self,cflags):
On Wed, Mar 16, 2011 at 8:27 AM, Thomas Witkowski
thomas.witkowski at tu-dresden.de wrote:
Jed Brown wrote:
On Wed, Mar 16, 2011 at 07:37, Thomas Witkowski
thomas.witkowski at tu-dresden.de mailto:thomas.witkowski at tu-dresden.de
wrote:
Thanks for explanations! It works fine in my
1.449869994655e-13 ||r(i)||/||b|| 9.151778504469e-13
So, are there now two or one outer iterations?
1 outer (NOT INDENTED) which took two iterates. 1 inner (INDENTED) which ran
3 times.
Matt
Thomas
Matthew Knepley wrote:
On Wed, Mar 16, 2011 at 8:27 AM, Thomas Witkowski
I fixed some bugs having to do with builds that break. I am using it for
every day building now.
I am experimenting with the Mercurial interface. I have a very preliminary
thing going in config/builder2.py:
python ./config/builder2.py stubs fortran
python ./config/buidler2.py check
python
On Mon, Mar 21, 2011 at 5:36 PM, Satish Balay petsc-maint at mcs.anl.govwrote:
Hm - why did I think Dmitry added libcharacteristic.a? I can't
figureout who added this library - so
Matt, Barry - Should it go into libpetscdm.a?
1) I added it a few years ago
2) It holds the method of
On Mon, Mar 21, 2011 at 6:11 PM, Jed Brown jed at 59a2.org wrote:
On Tue, Mar 22, 2011 at 00:06, Matthew Knepley petsc-maint at mcs.anl.gov
wrote:
3) It should go into TS
Really? It depends on a specific DM implementation. If you offered up a
generic algebraic interface, then it could
On Wed, Mar 23, 2011 at 3:21 PM, Barry Smith bsmith at mcs.anl.gov wrote:
hg diff src/makefile
How come there isn't an option to have it show the most recent set of
changes? It seems I have to hunt down some revision numbers and type them
in the request to get any useful information?
On Wed, Mar 23, 2011 at 4:47 PM, Satish Balay balay at mcs.anl.gov wrote:
On Wed, 23 Mar 2011, Matthew Knepley wrote:
On Wed, Mar 23, 2011 at 3:21 PM, Barry Smith bsmith at mcs.anl.gov wrote:
hg diff src/makefile
How come there isn't an option to have it show the most
On Wed, Mar 23, 2011 at 7:12 PM, Lisandro Dalcin dalcinl at gmail.com wrote:
On 23 March 2011 20:09, Barry Smith bsmith at mcs.anl.gov wrote:
On Mar 23, 2011, at 5:58 PM, Lisandro Dalcin wrote:
Are you going to rename MatNullSpaceDestroy - MatNullSpaceDestroy_
and add the backward
Fixed.
Matt
On Fri, Mar 25, 2011 at 10:57 PM, Satish Balay balay at mcs.anl.gov wrote:
The default 'configure --help' is is missing everything except the
header for help
satish
On Fri, 25 Mar 2011, Matthew Knepley wrote:
Now you can print specific help sections:
./configure
On Fri, Mar 25, 2011 at 11:38 PM, Satish Balay balay at mcs.anl.gov wrote:
btw: If I do 'configure --help | less' - I see some special chars:
ESC[?1034hConfigure Help
I do not see that. Does anyone else?
Matt
Satish
On Fri, 25 Mar 2011, Matthew Knepley wrote:
Fixed.
Matt
On Wed, Jun 1, 2011 at 3:48 AM, Jed Brown jed at 59a2.org wrote:
On Wed, Jun 1, 2011 at 04:56, Matthew Knepley knepley at gmail.com wrote:
I did not want to downcast because PetscPow() works for (Scalar, Int) and
(Scalar, Real) so automatically casting to (int) is not right.
1
On Wed, Jun 1, 2011 at 7:29 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On May 31, 2011, at 9:56 PM, Matthew Knepley wrote:
On Tue, May 31, 2011 at 9:45 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Man that PowInt is one ugly mo-fo. Since an int is an int is an int
why
do you
On Wed, Jun 1, 2011 at 7:36 AM, fabien delalondre delalf at
scorec.rpi.eduwrote:
Hi All,
I was wondering if you could let me know about a stable version (petsc-dev
or not) with the latest developments of the parallel block jacobi
preconditionner.
Use petsc-dev.
Matt
Thank you.
On Wed, Jun 1, 2011 at 9:06 AM, Jed Brown jed at 59a2.org wrote:
On Wed, Jun 1, 2011 at 14:51, Matthew Knepley knepley at gmail.com wrote:
Isn't it casting to PowInt which is int?
Not exactly. If I define
#define PetscPowScalar(a, b) std::pow(a, (int) b)
Barry was not suggesting
On Wed, Jun 1, 2011 at 11:23 AM, Jed Brown jed at 59a2.org wrote:
On Wed, Jun 1, 2011 at 18:20, Barry Smith bsmith at mcs.anl.gov wrote:
There is no reason for the existence of PetscViewerBinaryCreate(), I am
removing it from PETSc-dev.
Are you also removing MatCreateSeqAIJ()? Removing it
Pushed.
Matt
On Fri, Jun 3, 2011 at 2:38 PM, Matthew Knepley knepley at gmail.com wrote:
I hate fucking C. Let me push a fix.
Matt
On Fri, Jun 3, 2011 at 2:35 PM, Barry Smith bsmith at mcs.anl.gov wrote:
/Users/barrysmith/Src/petsc-dev/arch-gnu/bin/mpicc -o ex58.o -c -Wall
-Wwrite
PM, Matthew Knepley wrote:
I hate fucking C. Let me push a fix.
Matt
On Fri, Jun 3, 2011 at 2:35 PM, Barry Smith bsmith at mcs.anl.gov wrote:
/Users/barrysmith/Src/petsc-dev/arch-gnu/bin/mpicc -o ex58.o -c -Wall
-Wwrite-strings -Wno-strict-aliasing -Wno-unknown-pragmas -g3
-I/Users
Pushed.
Matt
On Fri, Jun 3, 2011 at 3:31 PM, Barry Smith bsmith at mcs.anl.gov wrote:
C++ is always used? I don't think so
Barry
=
libfast in: /Users/barrysmith/Src/petsc-dev/src
libfast in: /Users/barrysmith/Src/petsc-dev/src/inline
I was going over some pages today from a prog lang/comp opt person, trying
to figure out
if any of their stuff could apply to what we do. It was very hard, because
there were a lot
of general methods (incrementalization, loop transformation) that are known
to fail in many
cases. I found myself
On Sat, Jun 4, 2011 at 8:03 AM, Bobby Philip philipb at ornl.gov wrote:
Barry:
I heartily agree with all that you have said. Thanks for your suggestions.
I do believe what you have suggested is the best both in the shorter and
longer term. As a clarification on freezing the petsc-dev
I have made a new tarball and pushed the new name to petsc-dev. We should
know whether anything
is messed up from nightly builds.
Matt
On Thu, Jun 2, 2011 at 9:14 AM, John Fettig john.fettig at gmail.com wrote:
I think I sent this to the wrong place, I've cc'd petsc-maint on this
reply
On Tue, Jun 7, 2011 at 5:30 PM, David Andrs David.Andrs at inl.gov wrote:
Hi all!
I'm trying to hook my custom callback when using SNES using
SNESLineSearchSetPostCheck and it seems to me that this callback is not
being called. I based my test case on top of example 3 in
On Wed, Jun 8, 2011 at 10:47 AM, Blaise Bourdin bourdin at lsu.edu wrote:
Hi,
After poking a bit, I found out the bug causing a segfault in the
above-mentioned functions. It seems that gfortran will gladly nullify a
bogus pointer, whereas ifort can be a bit pickier, hence the bug coming up
On Thu, Jun 9, 2011 at 5:50 PM, Jed Brown jed at 59a2.org wrote:
On Fri, Jun 10, 2011 at 00:47, Barry Smith bsmith at mcs.anl.gov wrote:
So you are suggesting totally removing PCFieldSplitSetSchurPrecondition()
and instead having the user call MatSetApproximateSchurComplement() if they
wish
On Wed, Jun 15, 2011 at 1:20 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jun 14, 2011, at 7:59 PM, Jed Brown wrote:
On Wed, Jun 15, 2011 at 02:53, Barry Smith bsmith at mcs.anl.gov wrote:
It seems we should provide a DMSDA that is built specifically for
staggered grids, with the
On Wed, Jun 15, 2011 at 1:47 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jun 14, 2011, at 8:45 PM, Matthew Knepley wrote:
On Wed, Jun 15, 2011 at 1:42 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Jun 14, 2011, at 8:36 PM, Matthew Knepley wrote:
On Wed, Jun 15, 2011 at 1:20 AM
On Sun, Jun 19, 2011 at 1:38 PM, Jed Brown jed at 59a2.org wrote:
On Fri, Jun 17, 2011 at 23:00, Barry Smith bsmith at mcs.anl.gov wrote:
It is possible that list of elements is not correct for periodic boundary
conditions; I've never used it. I'm forwarding this to petsc-dev so Matt can
On Sun, Jun 19, 2011 at 8:12 PM, Jed Brown jed at 59a2.org wrote:
On Sun, Jun 19, 2011 at 21:44, Barry Smith bsmith at mcs.anl.gov wrote:
here is the stream benchmark results that Hongzhang Shan collected on
Hopper for Nick's COE studies. The red curve shows performance when you
run
Done.
Matt
On Wed, Jun 22, 2011 at 4:09 AM, Jed Brown jed at 59a2.org wrote:
I don't know why these functions exist, but they aren't namespaced and the
man pages are distracting and inconsistent. Also, only one of the
assembleVector() functions below is documented and neither of the
http://arxiv.org/abs/1106.4426
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...
On Fri, Jun 24, 2011 at 11:30 AM, Barry Smith bsmith at mcs.anl.gov wrote:
Moving this discussion to petsc-dev since it is relevant for much more
than this one preconditioner.
In PETSc the Mat object is responsible for the layout of the matrix
(graph) data and the Vec object is
This usually comes from mixing something debug with something optimized.
'petscstack' is only defined for debug builds.
Matt
On Tue, Jun 28, 2011 at 4:17 AM, Blaise Bourdin bourdin at lsu.edu wrote:
Hi,
On a freshly pulled / reconfigured / rebuild version of petsc-dev, I get
the
Jed has given excellent advice. However, your problems sound small. You
should try using a direct solver
like MUMPS (with --download-mumps during configure, and -ksp_type preonly
-pc_type lu -pc_factor_mat_solver_package mumps
during the solve).
Matt
On Tue, Jun 28, 2011 at 2:00 PM, Jed Brown
Is this the bad interaction Jed found between PCMG and ML?
Matt
On Wed, Jun 29, 2011 at 12:45 PM, domenico.borzacchiello at univ-st-etienne.fr
wrote:
Hi,
I use ML as a preconditioner for the pressure poisson equation for Stokes
problem. I solve for the pressure at each time step with ad
On Wed, Jun 29, 2011 at 2:48 PM, tuane tuane at lncc.br wrote:
Thank you for your help.
Well, but if the redundant preconditioner uses a direct solve what means:
1. number of iterations diversify with tolerance
2. reductions iterations with the use the KSPSetInitialGuessNonzero.
When I use
On Thu, Dec 1, 2011 at 4:21 PM, Matthew Knepley petsc-maint at
mcs.anl.govwrote:
On Thu, Dec 1, 2011 at 4:03 PM, msahinae00 msahin.ae00 at mynet.com wrote:
Dear Matt,
Element partition will be enough.
Thanks, I have started fixing up the Fortran. I will have it for you very
soon.
I
On Fri, Dec 2, 2011 at 6:20 PM, Dave Nystrom
Dave.Nystrom at tachyonlogic.comwrote:
I never received any reply to this question but would very much appreciate
one. Not sure if it fell through the cracks.
you can do whatever you want in PCSHELL. I would look at sacusp, since this
does about
On Sat, Dec 3, 2011 at 1:29 PM, Tim Gallagher tim.gallagher at
gatech.eduwrote:
Hi,
Has there been any discussion about adding a viewer for XDMF data? I'm
certainly interested in having one -- if somebody else is already working
on it or if it's never been considered before, I wouldn't mind
On Sat, Dec 3, 2011 at 2:39 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Dec 3, 2011 at 14:31, Tim Gallagher tim.gallagher at
gatech.eduwrote:
I was looking for the 'PETSc' way of doing it, something that lets you
call VecView(vec, viewer) that generates the heavy data and XML file
On Sat, Dec 3, 2011 at 2:49 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Dec 3, 2011 at 14:45, Matthew Knepley knepley at gmail.com wrote:
Yes, you can write several vectors using VecView(), AND you can do
timesteps :)
I saw that in your HDF5 functions. How do you store non-uniform
On Sat, Dec 3, 2011 at 3:24 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Dec 3, 2011 at 15:20, Matthew Knepley knepley at gmail.com wrote:
Does that matter? The same things are stored at each step.
Sure (if the mesh doesn't change, which is okay with me now), but you have
to put
On Sat, Dec 3, 2011 at 3:34 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Dec 3, 2011 at 15:25, Matthew Knepley knepley at gmail.com wrote:
Brad did it, but I believe he uses HDF5 attributes
I was speaking about the PetscViewerHDF5 interface.
There the timestep is just
On Wed, Dec 7, 2011 at 3:23 AM, Alexander Grayver
agrayver at gfz-potsdam.dewrote:
Hi dev-team,
Trying to compile petsc-dev with IFC 11.1 I got:
[ 0%] Building Fortran object CMakeFiles/petsc.dir/src/sys/**
f90-mod/petscsysmod.F.o
On Thu, Dec 8, 2011 at 10:19 AM, Alexander Grayver
agrayver at gfz-potsdam.dewrote:
Hi dev-team,
Using this code from Fortran:
call PetscViewerBinaryOpen(sub%**comm3d,'S.dat',FILE_MODE_**
WRITE,dviewer,ierr)
call PetscViewerSetFormat(dviewer,**PETSC_VIEWER_NATIVE,ierr)
call
On Fri, Dec 9, 2011 at 10:31 AM, Alexander Grayver
agrayver at gfz-potsdam.dewrote:
Hi dev-team,
I have this code:
call KSPCreate(comm3d,ksp,ierr);**CHKERRQ(ierr)
call KSPSetOperators(ksp,A,A,**DIFFERENT_NONZERO_PATTERN,**
ierr);CHKERRQ(ierr)
call
On Mon, Dec 12, 2011 at 5:35 AM, Jose E. Roman jroman at dsic.upv.es wrote:
The manpage of MatShellSetOperation states that all user-provided
functions should have the same calling sequence as the usual matrix
interface routines, but this does not apply to MATOP_DESTROY due to the
change
On Mon, Dec 12, 2011 at 5:35 AM, Alexander Grayver
agrayver at gfz-potsdam.dewrote:
**
Hello,
I use PETSs with MUMPS and looking carefully at the -ksp_view -ksp_monitor
results I see:
KSP Object:(fwd_) 64 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
On Wed, Dec 14, 2011 at 10:54 AM, Alexander Grayver agrayver at gfz-potsdam.de
wrote:
Hello dev-team,
I use the same code to output MPIAIJ and MPISBAIJ matrices, however with
latter program hangs somewhere in MatView.
The code is:
call
On Wed, Dec 14, 2011 at 12:45 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
I have a nonlinear problem and am finding the the convergence test for the
interior PETSc linear solver is apparently testing against the norm of the
first linear solve (see appended). This is clearly not what
On Wed, Dec 14, 2011 at 2:07 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
Mark
On Dec 14, 2011, at 1:57 PM, Matthew Knepley wrote:
On Wed, Dec 14, 2011 at 12:45 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
I have a nonlinear problem and am finding the the convergence test
On Wed, Dec 14, 2011 at 3:13 PM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Dec 14, 2011, at 3:44 PM, Barry Smith wrote:
On Dec 14, 2011, at 12:45 PM, Mark F. Adams wrote:
I have a nonlinear problem and am finding the the convergence test for
the interior PETSc linear solver
On Wed, Dec 14, 2011 at 3:50 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Dec 14, 2011, at 3:13 PM, Mark F. Adams wrote:
On Dec 14, 2011, at 3:44 PM, Barry Smith wrote:
On Dec 14, 2011, at 12:45 PM, Mark F. Adams wrote:
I have a nonlinear problem and am finding the the
On Wed, Dec 14, 2011 at 4:07 PM, Philip, Bobby philipb at ornl.gov wrote:
Hi:
When doing the following operation in PETSc using FD approximations to a
Jacobian vector product
J'(u)*a = [J(u+h*a) - J(u)]/h
Usually, we would put in F where you have J.
is there a way for the user to
On Thu, Dec 15, 2011 at 10:23 AM, Ravi Kannan rxk at cfdrc.com wrote:
Dear All,
** **
This is Ravi Kannan from CFD Research Corporation. Recently, we are
experimenting with the BoomerAMG preconditioner for some ?stiff? CFD
problems. In that regard, all the other standard
1:18 PM
*To:* For users of the development version of PETSc
*Subject:* Re: [petsc-dev] boomerAmg scalability
** **
** **
On Dec 15, 2011, at 1:56 PM, Matthew Knepley wrote:
On Thu, Dec 15, 2011 at 10:23 AM, Ravi Kannan rxk at cfdrc.com wrote:
Dear All
On Fri, Dec 16, 2011 at 9:37 AM, Dave Nystrom dnystrom1 at comcast.net wrote:
I'm trying to figure out whether I can do a couple of things with petsc.
1. It looks like the preconditioning matrix can actually be different from
the full problem matrix. So I'm wondering if I could provide a
On Fri, Dec 16, 2011 at 5:59 PM, Dave Nystrom Dave.Nystrom at tachyonlogic.com
wrote:
Barry Smith writes:
On Dec 16, 2011, at 9:52 AM, Matthew Knepley wrote:
On Fri, Dec 16, 2011 at 9:37 AM, Dave Nystrom dnystrom1 at comcast.net
wrote:
I'm trying to figure out whether I can do
I looked at the new checkins for fast MatMatSym. It looks to me like you
want a data
structure that has
a) O(N) space, where N is the number of entries
b) fast insert
c) fast traversal
It seems like the best data structure of this is
http://en.wikipedia.org/wiki/Y-fast_trie
It has
On Mon, Dec 19, 2011 at 9:34 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Mon, Dec 19, 2011 at 19:30, Matthew Knepley knepley at gmail.com wrote:
I looked at the new checkins for fast MatMatSym. It looks to me like you
want a data
structure that has
a) O(N) space, where N
at Akamai to drastically increase the number of streams since our
problem was so large.
Matt
Barry
On Dec 19, 2011, at 9:30 PM, Matthew Knepley wrote:
I looked at the new checkins for fast MatMatSym. It looks to me like you
want a data
structure that has
a) O(N) space, where N
On Mon, Dec 19, 2011 at 10:54 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Dec 19, 2011, at 9:30 PM, Matthew Knepley wrote:
I looked at the new checkins for fast MatMatSym. It looks to me like you
want a data
structure that has
a) O(N) space, where N is the number of entries
On Thu, Dec 22, 2011 at 9:11 AM, Chekuri S. Choudary
cchoudary at rnet-tech.com wrote:
** **
I am building an RPM for petsc-3.2-p5. I would like the users of the RPM
to be able to mention the PETSC_DIR at installation time, i.e, the RPM
should be relocatable.
For example,
rpm
We need ISGlobalToLocalMappingApply(). This could be done with an STL map :)
Matt
On Fri, Dec 23, 2011 at 1:01 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
Can we enumerate the matrix operations that are known to be non-scalable
currently?
1. MatPermute() does ISAllGather() on the row
On Fri, Dec 23, 2011 at 10:48 AM, Mark F. Adams mark.adams at
columbia.eduwrote:
On Dec 23, 2011, at 10:53 AM, Jed Brown wrote:
On Fri, Dec 23, 2011 at 09:50, Mark F. Adams mark.adams at
columbia.eduwrote:
Humm, my G-S in not in PETSc and it is perfectly scalable. It does have
more
On Fri, Dec 23, 2011 at 12:55 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Dec 23, 2011 at 12:27, Mark F. Adams mark.adams at
columbia.eduwrote:
A more interesting thing is partition down to the thread level and keep
about 100 vertices per thread (this might be to big for a GPU...)
On Wed, Dec 28, 2011 at 12:56 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Wed, Dec 28, 2011 at 12:51, Dave Nystrom dnystrom1 at comcast.net wrote:
I am experimenting now with solving my various linear systems with petsc
using a separate preconditioning matrix. These linear systems are
On Thu, Dec 29, 2011 at 9:38 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Dec 23, 2011 at 08:09, Matthew Knepley knepley at gmail.com wrote:
We need ISGlobalToLocalMappingApply(). This could be done with an STL
map :)
I'm not convinced it's really necessary, except to work
On Thu, Dec 29, 2011 at 11:01 AM, Dave Nystrom dnystrom1 at comcast.netwrote:
I have recently added the capability to have a separate preconditioning
matrix in the petsc interface for the code I am working with. I have two
types of preconditioning matrices that I have questions about. One is
On Wed, Sep 1, 2010 at 12:55 PM, Jed Brown jed at 59a2.org wrote:
Two thoughts:
1. Systematic deprecation system:
--with-shared doesn't give any errors today, but the option goes
unused (because it's now spelled --with-shared-libraries). I'd like
to see either (a) recognize the old
On Wed, Sep 1, 2010 at 2:52 PM, Jed Brown jed at 59a2.org wrote:
On Wed, 1 Sep 2010 15:39:40 +0300, Aron Ahmadia aron.ahmadia at kaust.edu.sa
wrote:
cgparse/optparse do what I need.
argparse is a nice improvement, but not in the standard lib prior to
2.7/3.2. It's not clear to me why
On Wed, Sep 1, 2010 at 6:16 PM, Jed Brown jed at 59a2.org wrote:
On Wed, 1 Sep 2010 18:02:18 +0200, Matthew Knepley knepley at gmail.com
wrote:
On Wed, Sep 1, 2010 at 2:52 PM, Jed Brown jed at 59a2.org wrote:
On Wed, 1 Sep 2010 15:39:40 +0300, Aron Ahmadia
aron.ahmadia at kaust.edu.sa
On Thu, Sep 2, 2010 at 10:32 AM, Jed Brown jed at 59a2.org wrote:
On Thu, 2 Sep 2010 09:00:47 +0200, Matthew Knepley knepley at gmail.com
wrote:
We want them to be organized in certain ways, like sections for
help.
That's what OptionGroup is for.
We want to do special processing
On Thu, Sep 2, 2010 at 11:07 AM, Jed Brown jed at 59a2.org wrote:
On Thu, Sep 2, 2010 at 10:58, Matthew Knepley knepley at gmail.com wrote:
Will it do everything we want, and be easy enough to customize if we want
something else?
I think so, but it's significant work to change, so it would
901 - 1000 of 2694 matches
Mail list logo